Assembly of metatranscriptomic data

Aragion · 02-03-2017, 04:02 AM

Hello.
I have a set of reads produced by metatranscriptomic cDNA sequensing of the total RNA sample. I need to assemble the RNA-virus genome, which was detected in this sample by the RT-PCR-analysis, using some reference sequence from GenBank and the reads. Can you give me some advise on how to perform this task.
Thanks.

maubp · 02-10-2017, 09:28 AM

You could try a reference guided assembly, perhaps using MIRA?

You might want to first try mapping all the reads against the reference (with relaxed mapping settings) in order to divide the reads into those which might be part of your virus (mapped to the reference) or not (unmapped).

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Accounting for species abundance changes in metatranscriptomic data	lbragg	RNA Sequencing	9	08-12-2016 08:39 AM
cuffmerge assembly vs denovo assembly of RNAseq data	skm	Bioinformatics	0	10-16-2013 10:16 PM
ORF prediction in assembled NGS metagenomic/metatranscriptomic data	someperson	Bioinformatics	5	07-03-2013 09:22 AM
Metatranscriptomic data(paired-end,Illumina)mapping?	jojohan	Illumina/Solexa	2	09-04-2012 06:59 PM

02-03-2017, 04:02 AM	#1
Aragion Junior Member Location: Russia Join Date: Dec 2016 Posts: 3	Assembly of metatranscriptomic data Hello. I have a set of reads produced by metatranscriptomic cDNA sequensing of the total RNA sample. I need to assemble the RNA-virus genome, which was detected in this sample by the RT-PCR-analysis, using some reference sequence from GenBank and the reads. Can you give me some advise on how to perform this task. Thanks.

02-10-2017, 09:28 AM	#2
maubp Peter (Biopython etc) Location: Dundee, Scotland, UK Join Date: Jul 2009 Posts: 1,543	You could try a reference guided assembly, perhaps using MIRA? You might want to first try mapping all the reads against the reference (with relaxed mapping settings) in order to divide the reads into those which might be part of your virus (mapped to the reference) or not (unmapped).