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
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
Comment