[mira_talk] FW: Ssaha-Mira question

  • From: <mark.rose@xxxxxxxxxxxx>
  • To: <mira_talk@xxxxxxxxxxxxx>
  • Date: Mon, 4 Jan 2010 14:52:31 -0500

Hi Bastien

Having used the suggestions you made below as to parameterizing ssaha2 to 
identify contaminating cloning adapters, I then launched mira (2.0.0) thusly:

nohup mira --project=CRW_GEP_454_Body_10302009 --job=denovo,est,normal,454 
 454_SETTINGS -LR:mxti=yes -CL:msvs=yes >& log_assembly3.txt &

$ tail log_assembly3.txt

Merging data from XML trace info file sample1_traceinfo_in.454.xml ...Num 
reads: 1091861
Building hash table ... done.
Localtime: Mon Jan  4 08:59:59 2010

Generated 1091861 unique template ids for 1091861 valid reads.
Done merging XML data, matched 1091861 reads.

Starting clips:  done.
Localtime: Mon Jan  4 09:00:07 2010

Merging vector screen data from SSAHA results file 
...Building hash table ... done.

Here it has remained for several hours

-bash-3.1$ top
 3660 rosema1   25   0 1958m 1.9g 1352 R  101  2.2 348:50.91 mira

-bash-3.1$ strace -p 3660
-bash-3.1$ strace -p 3660
Process 3660 attached - interrupt to quit
<after many minutes of no activity I quit>
Process 3660 detached

Is something wrong here?  I found a similar post in mira_talk but saw no 
followup to resolve it

Thanks for your help.


-----Original Message-----
From: Bastien Chevreux [mailto:bach@xxxxxxxxxxxx] 
Sent: Sunday, January 03, 2010 6:59 PM
To: Rose Mark USRE
Subject: Re: Ssaha-Mira question

On Montag 04 Januar 2010 mark.rose@xxxxxxxxxxxx wrote:
> I noticed your MIRA update that permits use of ssaha2 output to mask
>  sequence for a mira assembly.  Thanks. Would you mind being a little more
>  explicit about any parameterization of ssaha2 required by MIRA (i.e.
>  output format, etc.)? Thanks Again

Does this help?

Note 1: the output format of SSAHA2 must the native output format (-output 
ssaha2). Other formats cannot be parsed by mira. 
 Note 2: I currently use the following SSAHA2 options: -kmer 8 -skip 1 -seeds 
1 -score 12 -cmatch 9 -ckmer 6 
 Note 3: the sequence vector clippings generated from SSAHA2 data do not 
replace sequence vector clippings loaded via the EXP, CAF or XML files, they 
rather extend them.

From the main documentation file. There's also a very short example in 


This message may contain confidential information. If you are not the 
designated recipient, please notify the sender immediately, and delete the 
original and any copies. Any use of the message by you is prohibited.

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 

Other related posts: