MIRA output for Illumina mapping giving 100% coverage!

[Kasycas]

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!

Thanks,

Kasycas

[maubp]

I'd ask on the MIRA mailing list, but Bastien also posts here.

My guess is different rules/settings for mismatch thresholds, but also MIRA isn't really doing a mapping but a guided assembly (it will extend the reference sequence(s) if it can).