How to Detect Amino Acid Changes from a GATK VCF File and Reference Genome Using Code?
Hello,
I have a VCF file generated from GATK and the reference genome used for the analysis. I'm interested in detecting amino acid changes resulting from the variants listed in the VCF file without going for annotation like VEP or other tools. I understand that this involves interpreting the variants within coding regions and translating the resulting nucleotide changes into amino acids.
Specifically, I would like to:
- Parse the VCF file to extract variant information.
- Use the reference genome to determine the sequence context around each variant.
- Identify the codon changes and translate them to detect amino acid changes.
• 1,056 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
DareDevil
•
written
by
DareDevil
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Finding the Sequence of a Protein with Variant
written by hkarakurt
-
Translating a Fasta File
written by Hb •
-
Annotation of vcf with snpEff
written by Peerzada •
-
Find other variants causing same protein change
written by jakub •
-
Can HMMs detect partial proteins sequences?
written by Rui •
-
Finding wobble/redundant positions in an RNA string using R.
written by fawazfebin
-
Wrong predict consequence in variants using VEP Ensembl
written by theodorobiotec •
-
ANNOVAR ability to report reference and mutant codon
written by dahliashvets •
-
Effect Of Mutation On Codon
written by Gvj
-
Tool: Snpeff - Annotating And Predicting The Effects Of Single Nucleotide Polymorphisms
written by Malachi Griffith
You would need also a GFF/GTF file with gene models annotated in that reference genome. You could then
bedtools intersectyour genomic variants with CDS features and work out the codon changes considering the phase.that is basically what VEP and SNPEFF, bcftools csq, etc... do ...