This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Finding transitions and transversions in VCF file

Hello,

I have a VCF file in which I have to find all transitions and transversions. How do I detect if something is, for example, transition?

chr1    924533  .   A   G   50  PASS    .   GT:DP:ADALL:AD:GQ:IGT:IPS:PS    1|1:456:0,226:0,0:259:1/1:.:PATMAT

By looking at this line I know that ref allele is A and alt is G (which seems like transition), but then I also see that this is homozygous with GG genotype, so I really don't know how to tell.

vcf pysam python snp

0 answers

No answers yet.

Log in to answer this question.