How to calculate VAF from BAM or SAM(+pileups) files using pysam?
I assume pysam is a tool for BAM, SAM file.
But since the official document is a bit hard to understand. I am having difficulties calculating VAF values with pysam.
If you could suggest to me some scripts or documents, it would be grateful.
Thank you.
(It would be better if the method is restricted to using mainly pysam; not pysamstats.)
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNA-seq and Chip Seq alignment
written by alphaflylizard •
-
RdRp scan - identifying/detecting viruses- metagenomic workflow- need help
written by eric--carron •
-
Is there any good chromosome recombination(crossover) analysis tools that are not old?
written by askif4 •
-
Minfi preprocessNoob vs preprocessQuantile normalization for Illumina EPIC array
written by jeni •
-
BAM File Export Coverage Values from IGV or using PySam/Samtools
written by EOX123 •
-
Difficulty filtering VCF file with vcftools
written by rc16955
-
NCBI database schema
written by Joseph Hughes
-
Identify point mutations from each read in sam files
written by DVA
-
Which script / documentation can you recommend about promoters finding?
written by juanma_lace •
-
Reproducible Research: A More General Sweave?
written by David Quigley
Don't you need to do variant calling before you can calculate VAF?
I already have the location of the mutations, which means I have VCF files from BAM files I know that there are other methods to calculate VAF(using VCF files), even if it is not efficient, I want to use pysam. Thanks~
I see. Thanks for the update. I don't know how to calculate VAF using pysam - others on the forum might be able to help you.