getting base count (coverage) based on BAM file with strand information
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
1 answer
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'
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Coverage depth map of genome from multiple sample coverage bedgraphs?
written by Agastya •
-
Should we discard a variant if it pass many filters and fails one?
written by Sharon
-
Median coverage from BAM file
written by firestar
-
Calculating Variant Allele Frequency
written by novice
-
Heterozygous SNPs: Coverage distribution.
written by Coryza
-
True or False, SNVs using bwa+gatk
written by J.F.Jiang
-
True or False, SNVs using bwa+gatk
written by J.F.Jiang
-
Variant Calling For Ultra High Coverage Region Of Whole Exome Sequencing Data
written by lyz10302012
-
Quality Control For Exome Sequencing
written by Omid
-
Filtering Strategy In Exome Sequencing And Quality Control
written by Omid
in IGV you can color the reads by strand. For your purpose, is not this enough?