MIRA for metagenome assembly
0
0
Entering edit mode
9.1 years ago
Rahul Sharma ▴ 660

Hi all,

I am working on a de novo metagenome assembly project. In this data I am expecting several species of bacteria. The sequencing data I have received have been sequenced with Illumina 100-300 bp inserts with a read length of around 150 bps. I was giving a try to MIRA genome assembler, however there are around 150+ parameters to tune a assembly. I will be very thankful if anyone can share his/her experience with MIRA for metagenome assembly and how to optimize it for several closely related bacterial species within the sample.

Here is how my manifest file looks like:

project = MyFirstAssembly
job = genome,denovo,accurate
parameters = -GE:not=4

readgroup = SomePairedEndIlluminaReadsIGotFromTheLab

data = datape*.fastq
technology = solexa
template_size = 100 300 autorefine
segment_placement = ---> <---
segment_naming = solexa

I would like to know if anyone have tested MIRA with some specific parameter for metagenomic data, where diversity is high and chances missing out is also high. It would be very great if anyone can share an example manifest file.

Many thanks in advance!

Best regards,
Rahul

Metagenome next-gen MIRA • 3.0k views
ADD COMMENT
0
Entering edit mode

Mira relies heavily on the assumption, that the coverage along the whole genome is uniform/constant. If your DNA was harvested from a single organism or from a single bacterial colony, this assumption is true for most part of the genome. In metagenomics the genome consists of several chromosomes and plasmids. I doubt, that all of them are present in your sample with the same copy number.

ADD REPLY

Login before adding your answer.

Traffic: 1982 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