This is a test version of Biostars. For the public version, visit https://www.biostars.org.
samtools idxstats showing accession number and genomic scaffold

Hi,

Is this an expected output for Samtools idxstats? Why is it showing accession instead of chromosome number? most of the NC_... corresponds to a chromosome but why do I have NT and NW genomic scaffold?

Samtools Manuel :

output is TAB-delimited with each line consisting of reference sequence name, sequence length, # mapped read-segments and # unmapped read-segments.

But why do I have have four columns in some?

Samtools idxstats file.bam

Screen-Shot-2021-01-02-at-7-30-07-pm
multiple image

samtools idxstats bams
• 1,673 views • link
updated by Pierre Lindenbaum Pierre Lindenbaum • written by whb • whb

1 answer

Is this an expected output for Samtools idxstats?

yes if those accession numbers were used as the name of the chromosomes in the reference genome....

you can check this using

samtools view -H your.bam | grep '^@SQ'
• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this