[mira_talk] Problem with convert_project.

  • From: Vésteinn Snæbjarnarson <ves4@xxxxx>
  • To: mira_talk@xxxxxxxxxxxxx
  • Date: Thu, 21 Jul 2011 13:28:46 -0000 (GMT)

Hi

I'm running MIRA on a hybrid set of 454-reads and paired 2x36 SOLEXA
reads. I'm trying to convert the resulting caf file with convert_project
but the following error is reported for the first read in the file:

-----------------
Now loading and processing data:
 [0%] Sequence size: 38 Qualities size: 39

Internal logic/programming/debugging error (*sigh* this should not have
happened).
Please file a bug report on http://sourceforge.net/apps/trac/mira-assembler/

"qualities.size() != sequence.size() : 8:92:14157:14935:Y"

->Thrown: void Read::initialiseRead(bool preserve_originals, ...,
...<base_quality_t>, ...)
->Caught: Read & CAF::createCafRead()
----------------

This seems to be the same kind of error as noted here
//www.freelists.org/post/mira_talk/convert-project-internal-logicprogrammingdebugging-error.

I'm using the latest stable mira release, 3.2.1, with parameters
--job=denovo,genome,normal,454,solexa SOLEXA_SETTINGS
-GE:tismin=1500:tismax:4500 -SB:ads=yes COMMON_SETTINGS
-SB:load_straindata=1

Any thoughts?

Cheers,
Vesteinn


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