Hi all,
I called variants using GATK, and annotated the results using snpeff. Since the organism is chinese hamster, there is not much information available. I want to use provean to predict the effects of variants. The problem is that the input for provean for non human/mouse model should be amino acids sequence and amino acids variants. What I have now are variants in genomic level, does anyone know any tool that can transfer the genomic vcf files to input files for provean? Thanks.
1 answer
Unfortunately, it's not easy, because this requires mutation effect predictors to work in terminal, and most of those are old and work poorly. I've spent some time trying to get Provean to work, for example, but without any luck.
One tool that works brilliantly (and is well supported) is SIFT4G, a re-implementation of older SIFT. You can get all the instructions here: https://github.com/rvaser/sift4g
You'll need to work out the calibration (cutoffs) since databases change all the time.
Log in to answer this question.
-
Selective sweep for scaffold level reference genome
written by Mel Kile •
-
Ways to improve SnpEff prediction
written by Mohammed •
-
variant annotation issue
written by prasundutta87
-
annotate DNA missense variants with protein sequence context
written by Irsan
-
(SIFT, Polyphen, PROVEAN) vs. (ANNOVAR) vs. (LOFTEE)
written by genetic •
-
What are SNP effects? in terms of SNPeff output file
written by JstRoRR
-
Predict Impacts of variants (from VCF file) using Dog Models
written by sruthimouli88 •
-
Bioinformatic tools for identifying protein-interaction network for non model organism (plants)
written by rajaragupathy •
-
how to fetch GO categories from Chinese hamster genome
written by upendrakumar.devisetty
-
Tool: Snpeff - Annotating And Predicting The Effects Of Single Nucleotide Polymorphisms
written by Malachi Griffith
Did snpEff list the coding variants? You can use that to filter the VCF and get only the relevant variants (as the first step)
Thanks. Yes, there is coding variants, but provean also needs the whole protein sequence, After building the database using snpeff, there is only one .bin file, do you know how to get the whole protein sequence? The only way I can think of is using genome annotation file, but it would be a little bit tricky.
Standalone PROVEAN? I've only used Web based PROVEAN. If you know the start codon location, translation shouldn't be a huge problem. I can't recollect any tool that gives you protein mutations from nucleotide changes, sorry :-(
Hi @shl198 I have the same issue did you find a solution
Do not add answers if you're not answering the top level question. Use Comments instead.
Hi @shl198 I have the same issue did you find a solution