sam file header, how is @SQ LN calculated?
I would like to understand how the LN values found in the header are calculated. The sam specification says that it is the reference sequence length, but it doesn't seem to match. I have the following in a sam file:
@SQ SN:I LN:15072434
but that doesn't match the nucleotide count of the reference sequence:
gunzip -c Caenorhabditis_elegans.WBcel235.dna_rm.chromosome.I.fa.gz | grep -e ^[^\>] | wc -c
returns 15323642
• 3,264 views
•
link
written
by
n.roggli •
1 answer
found why, I was counting newline characters too
gunzip -c Caenorhabditis_elegans.WBcel235.dna_rm.chromosome.I.fa.gz | grep -e ^[^\>] | tr -d '\n' | wc -c
returns the expected value 15072434
• 1 views
•
link
written
by
n.roggli •
Log in to answer this question.
More posts like this
-
Sam file Header problem
written by saifulislam99121 •
-
How is contig length calculated ?
written by Alewa
-
samtools: Numerical result out of range
written by Gregor Rot
-
BWA mem giving malformed SAM headers.
written by Ismaelrymy •
-
Reformat sam header issue
written by pg_canada •
-
How to load BAM files that do not have chromosome name in the header?
written by JC •
-
A constant problem with SAM file header
written by popescuiofelia •
-
Reference genome and samtools reheader option
written by shuksi1984 •
-
Sam Reference sequence length does not match with the actual fasta sequence input
written by ashishtx •
-
Picard API: Modify sam header
written by dariober
