how to extract all reads covered at the somatic SNP site that from a BAM file in VCF file with using python?
HI , I am doing some tumor research. I used GATK to call somatic SNV . now I have the tumor bam file and VCF file . I wanna extract all reads suggested every SNV record from that VCF file . So how to complete that goal with using python ??
• 2,006 views
•
link
written
by
winter_li •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Somatic variant calling with Mutect2 using 2 matched-normals for one sample
written by Alexandros •
-
best pipeline/practice for analyzing WES for germline and somatic variant calling
written by epigene
-
Novogene somatic results with GATK mutect2
written by chrisgr •
-
BreakDancer output file doesn't have any SV detection from the Tumor BAM file
written by Hasib •
-
Empty output of allele specific expression analysis by GATK
written by zizigolu
-
the difference between the positions called by somatic caller (MuTect2) and the position called by …
written by Zoey Zou •
-
How to merge vcf files (output of VarScan2 somatic)?
written by Raheleh
-
Is it possible to extract somatic variants from paired VCF files?
written by csardas •
-
Can I use exome sequencing data from cancer patient's blood to call germline variants?
written by ShirleyDai •
-
No Somatic SNV Calling with GATK?
written by sichan
Do you know python? Or you want someone to throw up a snippet here for you? Did you follow 'GATK Best Practices' to get your VCF? If you know python, you should look more into the format of a VCF file , and then come back with your own attempt to parse these data you've generated.