View Single Post
Old 09-10-2010, 02:52 AM   #1
Location: Dublin, Ireland

Join Date: Sep 2009
Posts: 22
Default MIRA output for Illumina mapping giving 100% coverage!

Hi all,

I've been using MIRA as I have both Illumina and 454 reads that I wanted to do mapping/assembly as well as semi-hybrid and hybrid assembly. However, when I use MIRA on the Illumina data alone for mapping with

"mira --project=0800 --job=mapping,genome,normal,solexa --fastq -SB:bft=gbf >&0800_Illumina_log_assembly.txt"

I get 1 contig of length 4,798,671 which is pretty much 100% coverage of the genome! I've used bwa and bowtie as well which give me a more realistic ~68-69% coverage. Can anyone tell me why/how MIRA is doing this? I've played around with the setting and read the documentation but I'm at a loss!


Kasycas is offline   Reply With Quote