SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   MIRA output for Illumina mapping giving 100% coverage! (http://seqanswers.com/forums/showthread.php?t=6798)

Kasycas 09-10-2010 02:52 AM

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!

Thanks,

Kasycas

maubp 09-10-2010 03:45 AM

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


All times are GMT -8. The time now is 08:54 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.