How do I display a single nucleotide position from reads in a BAM file? I have been looking at variation using samtools mpileup, but I want to actually just display the nucleotide at the position I am interested in. This seems like you should be able to do it but I can't figure out how.
To be clear I have a BAM file of a bunch of reads. I'm looking to do something like samtools magic reads.bam chr3:10000
where I get back something like:
T T T T T T T T T T T T T T T T T T T A A A A A A A A A A A A A.
I just want to sanity check the output of bcftools by actually looking at the base calls.
No comments:
Post a Comment