Question: Snp Contig-Synonymous & Nonsynonymous
0
gravatar for jackuser1979
7.8 years ago by
jackuser1979870
US
jackuser1979870 wrote:

I am working on non-model organism. I have called SNPs by read mapping illumina reads to reference genome(assembled contigs/scaffolds). I have VCF file with SNP predicted for each contig/scaffold. Recently Ensembl released non-chromosomal genome assembly and also have GTF file(gene prediction based on Ensembl annotation). Now I am trying to find out synonymous and nonsynonymous mutation rate. I find tools that only available for chromosomal level, Is there any script or tool available already for contig level?.

vcf annotation contigs snp mutation • 3.3k views
ADD COMMENTlink written 7.8 years ago by jackuser1979870

The VAASST pipeline can annotate variants. All you need is the annotation set in GFF3 format, the reference sequence, the variants, and the lengths of the chromosomes. Converting GTF to GFF3 should also be trivial.

ADD REPLYlink modified 7.8 years ago • written 7.8 years ago by Zev.Kronenberg11k

Will this pipeline accept contig lengths?

ADD REPLYlink written 7.8 years ago by jackuser1979870

Yes the --build flag is a tab delimited file that has 3 columns: 1. seqid/contig_name/scaffold 2. start position (1). 3. length of sequence

ADD REPLYlink written 7.8 years ago by Zev.Kronenberg11k
2
gravatar for Bioch'Ti
7.8 years ago by
Bioch'Ti1.0k
France (Avignon)
Bioch'Ti1.0k wrote:

Hi, You may have also a look to snpEff (SNP effect predictor) which determine the nature (syn/non-syn) of variants. Here is the link: http://snpeff.sourceforge.net/ Hope this helps, C.

ADD COMMENTlink written 7.8 years ago by Bioch'Ti1.0k

I have tried this software. Somehow it doesnot output anything. Did anybody tried running snpEff with VCF as input file containing with contigs instead of chromosomes?

ADD REPLYlink written 7.8 years ago by jackuser1979870

You just have to read the documentation carefully, for adding a genome. It's all there. In terms of contigs or chromosomes... there isn't really a huge difference to snpEff.

ADD REPLYlink written 7.7 years ago by Lee Katz3.0k

This software is really great, and the developer is responsive. The documentation is good too.

ADD REPLYlink written 7.7 years ago by Lee Katz3.0k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1114 users visited in the last hour