Get by-strand base counts at position from cram file
0
0
Entering edit mode
4 months ago

If I use igv interactively I can load a cram file, goto a position and click on that position to view the count details, including the numbers of alt and ref alleles which occur on each strand. I can use igv in batch mode to generate a snapshot showing the reads but I can't find a way to output the count details in text format, even though igv obviously knows them.

I just want some software that will do this in batch mode. It seems simple but I can't find anything which says it does this. I just want to specify a cram file and a position and for the program to tell me the counts for each allele on each strand.

All I can find is things that will tell me e.g. the total read depth, or total allele counts or allele balance. But I'd really like to see the counts by strand because if I see that there are reads supporting both alleles on both strands I'll feel more confident that the genotype call is valid.

Please can somebody point me at a solution. Thanks.

igv cram count strand details • 215 views
ADD COMMENT
0
Entering edit mode

Thanks. So to get the strand-specific information it looks like I need a tool called bam-readcount, which is available here: https://github.com/genome/bam-readcount

ADD REPLY

Login before adding your answer.

Traffic: 1712 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6