Tag: HTSlib

Ubuntu Manpage: samtools reheader – replaces the header in the input file

Provided by: samtools_1.13-2_amd64 NAME samtools reheader – replaces the header in the input file SYNOPSIS samtools reheader [-iP] [-c CMD | in.header.sam ] in.bam DESCRIPTION Replace the header in in.bam with the header in in.header.sam. This command is much faster than replacing the header with a BAM→SAM→BAM conversion. By default…

Continue Reading Ubuntu Manpage: samtools reheader – replaces the header in the input file

Fail to install Ensembl VEP

Fail to install Ensembl VEP 2 Hello I have been encounter some issues on installing the Ensembl VEP. I am using Ubuntu 20.04 system, with the newest perl v5.34. perl -MBio::Root::Version -e ‘print $Bio::Root::Version::VERSION,”n”‘ 1.006924 However, when I install the VEP, it failed. Hello! This installer is configured to install…

Continue Reading Fail to install Ensembl VEP

Running samtools view on bam affects the number of variants called by both haplotypecaller and deepvariant – C samtools

Thanks for getting back to me Valeriu. As you suggested, I used the latest commit from the develop branch in my pipeline, and the results look good. I was able to replicate the numbers from samtools v1.10.2 and v1.11 for both variant callers. FYI $ docker run scilifelabram/htslib:dev_proper /opt/samtools/samtools version…

Continue Reading Running samtools view on bam affects the number of variants called by both haplotypecaller and deepvariant – C samtools

Access read group tag for BAM reads – rust-htslib

I’m sure this is probably straightforward, but I can’t figure out how to access the read group information on a per-read basis within a BAM file. I can get it from the header, but I can’t determine how to use aux or any other record method to get that info….

Continue Reading Access read group tag for BAM reads – rust-htslib

kojix2/ruby-htslib: v0.0.4 | Zenodo

Zenodo DOI Badge DOI 10.5281/zenodo.5797825 Markdown [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.5797825.svg)](https://doi.org/10.5281/zenodo.5797825) reStructedText .. image:: zenodo.org/badge/DOI/10.5281/zenodo.5797825.svg :target: doi.org/10.5281/zenodo.5797825 HTML <a href=”https://doi.org/10.5281/zenodo.5797825″><img src=”https://zenodo.org/badge/DOI/10.5281/zenodo.5797825.svg” alt=”DOI”></a> Image URL zenodo.org/badge/DOI/10.5281/zenodo.5797825.svg Target URL doi.org/10.5281/zenodo.5797825 Read more here: Source link

Continue Reading kojix2/ruby-htslib: v0.0.4 | Zenodo

Issue with installing QIIME2 2021.11 on Windows 10 – Technical Support

Hi QIIME support team, I’m attempting to install QIIME2 on my Windows 10 machine. I installed Anaconda3, then set up conda to run in Git Bash: echo “. ${PWD}/conda.sh” >> ~/.bashrc Once I restarted Git Bash and activated Conda, I installed python-wget because installation of wget kept getting the following…

Continue Reading Issue with installing QIIME2 2021.11 on Windows 10 – Technical Support

Issues installing Pindel from git

Issues installing Pindel from git 1 I’m currently attempting to install Pindel on my ubuntu machine (16.04 server), but have been having issues with running the ./INSTALL file. The ./INSTALL file takes the path to htslib as an argument, so I also installed htslib from git with git clone github.com/samtools/htslib

Continue Reading Issues installing Pindel from git

cellranger count DETECT_COUNT_CHEMISTRY (failed)

cellranger count DETECT_COUNT_CHEMISTRY (failed) 0 I am learning scRNA-seq and the tutorial I follow uses dataset (1k pbmcs from healthy donor) from 10X genomics website. I downloaded fastq and reference transcriptome files and ran following command. cellranger-6.1.1/cellranger count –id pbmc_1k_v2_example –transcriptome /home/murat/Share/single_cell/refdata-gex-GRCh38-2020-A –fastqs /home/murat/Share/single_cell/pbmc_1k_v2_fastqs I get following message. Martian Runtime…

Continue Reading cellranger count DETECT_COUNT_CHEMISTRY (failed)

plotting roh from bcftools

plotting roh from bcftools 0 Heys, I am following this small tutorial on how to calculate ROHs from a vcf file using bcftools (samtools.github.io/bcftools/howtos/roh-calling.html) and I am getting this txt file: # This file was produced by: bcftools roh(1.10.2+htslib-1.10.2-3) # The command line was: bcftools roh -G30 –AF-dflt 0.4 my_file.vcf…

Continue Reading plotting roh from bcftools

User friendly (visual&interactive) VCF/BCF mining tools (2021)

What is currently the best user friendly (visual and interactive) VCF/BCF mining tool in 2021? For VCF/BCF similar to size or even larger than the 1000 human genomes VCF? I guess most organization do not have a visual and interactive mining VCF mining tool but use either: A website front-end…

Continue Reading User friendly (visual&interactive) VCF/BCF mining tools (2021)

tabix for ID column

tabix for ID column 4 Hello, I’m looking for something similar to tabix. But instead of looking for informations within a given region, I would like to use the values in the ID column for quickly lookup. So for example I would like to take the compressed dbSNP file, index…

Continue Reading tabix for ID column

install GenomicFeatures fail

install GenomicFeatures fail 1 @5b9023e7 Last seen 19 hours ago China BiocManager::install(‘GenomicFeatures’) results show ‘getOption(“repos”)’ replaces Bioconductor standard repositories, see ‘?repositories’ for details replacement repositories: CRAN: mirrors.tuna.tsinghua.edu.cn/CRAN/ Bioconductor version 3.14 (BiocManager 1.30.16), R 4.1.0 (2021-05-18) Installing package(s) ‘GenomicFeatures’ also installing the dependencies ‘Rhtslib’, ‘Rsamtools’, ‘GenomicAlignments’, ‘rtracklayer’ Packages which are only…

Continue Reading install GenomicFeatures fail

bcftools merge; retaining sample names

bcftools merge; retaining sample names 2 When I do bcftools merge, the headers do not retain the filenames.  How can I specify filenames? This is my command  bcftools merge vcf/unfiltered/*.vcf.gz -O z > msa/pooled.vcf.gz However this is the relevant part of my header, despite the filenames I gave it.  Is…

Continue Reading bcftools merge; retaining sample names

Disappearing CB, the bam tag after samtools sort -t CB

  I’ve been trying to setup an analysis pipline for RNAvelocity in AWS EC2. I used one of the 10x dataset, 10k Peripheral blood mononuclear cells (PBMCs) from a healthy donor, Single Indexed, as a test model to setup the pipeline. For speed and cost saving, I first used samtools…

Continue Reading Disappearing CB, the bam tag after samtools sort -t CB

What does samtools mean by ‘orientation’ when marking duplicates?

What does samtools mean by ‘orientation’ when marking duplicates? 0 Hello Biostars, It is my understanding that samtools marks duplicates on the basis of the 5′ position of reads and also the orientation of reads. This is based on my reading of the following: www.htslib.org/algorithms/duplicate.html However, I am not sure…

Continue Reading What does samtools mean by ‘orientation’ when marking duplicates?