[mira_talk] Re: 454 mapping
- From: Isabelle Lesur <isabelle.lesur.kupin@xxxxxxxxx>
- To: mira_talk@xxxxxxxxxxxxx
- Date: Thu, 3 May 2012 13:57:45 +0200
Hello!
I used this parameter with 454 and Sanger sequences:
-AS:epoq=no
You can give it a try...
MIRA manual says:
[enforce_presence_of_
qualities(epoq)=on|yes|1, off|no|0] Default is yes. When set to yes, MIRA
will stop the assem-
bly if any read has no quality values loaded.
However, I think your 454 sequences quality data will be ignored as well..
Isabelle
2012/5/3 davis gimode <dgimode@xxxxxxxxx>
> Hi,
> I am trying to do a mapping assembly of 454 data on to a reference
> genome. My call parameters are:
>
> mira --project=indica --job=mapping,genome,accurate,454 -AS:nop=4
>
> -SB:load_backbone=yes:startbackboneusage_inpass=3:backbone_strainname=indica_backbone_in.fasta:backbone_filetype=fasta:bbq=-1
>
> However, it requires quality scores for the backbone i.e it wants
> indica_backbone_in.fasta.qual which i dont have. How do i disable mira
> from requiring the quality files?
>
> Thanks in advance
>
> --
> Davis M. Gimode
>
> --
> You have received this mail because you are subscribed to the mira_talk
> mailing list. For information on how to subscribe or unsubscribe, please
> visit http://www.chevreux.org/mira_mailinglists.html
>
Other related posts:
