getting base count (coverage) based on BAM file with strand information
1
0
Entering edit mode
8.2 years ago
J.F.Jiang ▴ 910

Hi all,

I want to know if there is an convenient way to get the base count from a BAM file with strand information,

for example, on chr1:1-1, the depth of coverage is 1000, from IGV we could know the site is heterozygous, and the depth for each allele is A: 250(+), 250(-), C: 250(+), 250(-), from which we could know the site is definitely a heterozygous variant without any strand bias or allele drop off.

Best

bam strand • 2.3k views
ADD COMMENT
0
Entering edit mode

in IGV you can color the reads by strand. For your purpose, is not this enough?

ADD REPLY
0
Entering edit mode
8.2 years ago

you could use my tool sam2tsv : https://github.com/lindenb/jvarkit/wiki/SAM2Tsv#output , use awk+bit-operation to extract the strand from the flag column, cut columns, 'sort|uniq -c'

ADD COMMENT

Login before adding your answer.

Traffic: 1793 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