[mira_talk] De novo assembly using Illumina and PacBio data aborts after internal error

  • From: Bibi Lujan <bibilujan@xxxxxxxxx>
  • To: mira_talk@xxxxxxxxxxxxx
  • Date: Thu, 23 Jun 2016 14:08:19 -0400

Hello,

I am using MIRA (4.9.6) to assemble a plant genome (260MB) with illumina
and corrected PacBio data. However my assembly has aborted after the
following error:

Internal logic/programming/debugging error (*sigh* this should not have
happened)
------------------------------

   - Alignment of two sequences: 81726 > 65530 bases, cannot handle,
   shouldn't *
   - have happened! *

------------------------------

->Thrown: void AlignedDualSeq::acquireSequences(const char * seq1, const
char * seq2, readid_t id1, readid_t id2, int8 id1dir, int8 id2dir, bool
enforce_clean_ends, bool dontpenalisengaps, int32 bandwidthused, int32
minbanddistance)
->Caught: void Align::fullAlign(std::list<AlignedDualSeq> * adslist, bool
enforce_clean_ends, bool dontpenalisengaps)

I have put the full log here :

https://www.dropbox.com/s/47stfpkn14vgqm5/Assembly_MIRA_June13r.txt.zip?dl=0

I have re-started the assembly only to get a similar error (Alignment of
two sequences: 66164 > 65530 bases, cannot handle, shouldn't  have
happened!)

Is there any way to continue the assembly?

Thanks,

Beatriz

Other related posts: