Question: bam-readcount output: is case significant?
0
gravatar for prograde
4.8 years ago by
prograde30
prograde30 wrote:

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?

 

bam-readcount next-gen • 1.2k views
ADD COMMENTlink modified 4.8 years ago • written 4.8 years ago by prograde30
4
gravatar for prograde
4.8 years ago by
prograde30
prograde30 wrote:

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.

ADD COMMENTlink written 4.8 years ago by prograde30
1

This is correct. It's looking up the info from your fastq.

ADD REPLYlink written 4.8 years ago by Chris Miller21k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1125 users visited in the last hour