This is correct. It's looking up the info from your fastq.
• 0 views
•
link
written
by
Chris Miller
bam-readcount output: is case significant?
I've noticed that the output for bam-readcount sometimes has bases in uppercase, sometimes lowercase. This is in the reference column, and also for additional columns (after A, C, G, T, N) for indels. Sometimes an indel is even given in mixed case (like "aaaTT").
What is the significance of case in the bam-readcount output?
• 2,456 views
•
link
updated
by
Ram
•
written
by
prograde •
•
written
by
prograde •
1 answer
Answering my own question (I think): the lower case bases originate in my reference FASTA file (soft masked repeats). I re-wrote my reference as all uppercase and the bam-readcount output also became all uppercase (but otherwise the same).
Can anyone confirm or deny that bam-readcount maintains case based on the content of "-f ref.fa"? ...just for the sake of my own sanity.
• 0 views
•
link
written
by
prograde •
Log in to answer this question.
More posts like this
-
Bcftools consensus generates mismatched consensus sequence
written by Duy •
-
bam-readcount empty output
written by D_Jens •
-
Samtools mplileup detects wrong INDEL positions
written by Gama313
-
obisplit (from Obitools) --uppercase doesn't work
written by hexphe •
-
What are lowercase a,t,c,g refer to in a complete genome? Should lower case a,t,c,g to be converted…
written by saranpons3
-
bam-readcount fails to report all positions
written by mark.rose •
-
bam-readcount original output changes with larger site list file
written by twtoal •
-
Using sequence case to encode intron/exon information
written by anderspitman •
-
How To Interpret Output From Bam-Readcount
written by merbergd •
-
Can I Convert Fasta Lowercase Bases To 'N'?
written by User 8087 •
