Snp Contig-Synonymous & Nonsynonymous
1
0
Entering edit mode
9.1 years ago
jackuser1979 ▴ 880

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?.

snp vcf contigs annotation mutation • 3.7k views
ADD COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

Will this pipeline accept contig lengths?

ADD REPLY
0
Entering edit mode

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 REPLY
2
Entering edit mode
9.1 years ago
Bioch'Ti ★ 1.0k

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 COMMENT
0
Entering edit mode

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 REPLY
0
Entering edit mode

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 REPLY
0
Entering edit mode

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

ADD REPLY

Login before adding your answer.

Traffic: 3260 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6