ref and alt allele count from BAM file
Hi All,
I would like to calculate the allele ratio of ref/alt at a specific position. How can I determine the ref and alt counts from the BAM file? (or is it better to extract this info from the VCF?)
Thank you
• 3,421 views
•
link
updated
by
Vitis
•
written
by
monaallouba •
•
written
by
monaallouba •
1 answer
Both BAM and VCF can be used for this. In VCF, you would be looking for the values in the AD field for each sample. For BAM, you may want to consider using PySAM's pileup engine to query the position and reads stacked at that position. Just be careful with indels when using PySAM, you'll have to add special code to distinguish no coverage and deletions at these positions.
To tap into VCF, use PyVCF
https://pyvcf.readthedocs.io/en/latest/
To tap into BAM, use PySAM
Log in to answer this question.
More posts like this
-
How can I get an allele count table from a .vcf file?
written by mrk •
-
VEP annotated vcf is not readable in bcftools
written by Nai •
-
Extract the DP and AD from VCF file along with the chromosome postion and alteration
written by manojkumarbioinfo
-
Extract the DP from VCF file along with the chromosome postion and alteration
written by manojkumarbioinfo
-
allele frequency bcf tools
written by jh •
-
Calculate allele frequency from many VCF files in specific locus
written by John
-
How to define output for ALT tag in vcf file
written by Kash
-
Allele specific expression
written by monaallouba •
-
Using samtools to assess the number of reads that cover a variant position, as well as the ratio of…
written by claudiadast •
-
Ancestral Allele for the Genus Gallus
written by lawalakinyanju •
