[mira_talk] Re: 454 mapping

Thanks a lot. Will give it a try.

On Thu, May 3, 2012 at 2:57 PM, Isabelle Lesur
<isabelle.lesur.kupin@xxxxxxxxx> wrote:
> 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
>
>



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