[mira_talk] mira run died on assertion
- From: Marc Lohse <Lohse@xxxxxxxxxxxxxxxxx>
- To: "mira_talk@xxxxxxxxxxxxx" <mira_talk@xxxxxxxxxxxxx>
- Date: Mon, 24 Aug 2009 10:28:04 +0200
Hi,
i am trying to assemble about 2 million 454 reads from three different plant
genotypes and am interested in finding genotype-specific SNPs (using
mira_3rc2hf1_dev_linux-gnu_x86_64_static). As i understood from reading the
manual this is a three-step process - i have started the first step using
the following cmd-line params:
--job=denovo,normal,454,esps1 -GE:not=1 -SK:mnr=yes:nrr=10 454_SETTINGS
-CL:qc=no:cpat=no
The data contained a traceinfo file providing clipping info to make sure all
5' adaptor and barcode sequences are removed. After running merrily for 4
days only consuming 12GB of memory on a 32 GB machine the run died with the
following error message:
miraSearchESTSNPs: contig.C:3631: void
Contig::addRead_wrapped(std::vector<Align, std::allocator<Align> >&, const
AlignedDualSeqFacts*, int32, int32, int32, bool, int32,
Contig::errorstatus_t&): Assertion `(xcut) >= (0) && (xcut) <
(CON_counts.size()+1)' failed.
At that time it had already build almost 6000 sensible-looking contigs so we
were in good spirits. Is there any obvious mistake we've made? Did anybody
else maybe also see this error?
Any help would be appreciated.
thanks in advance,
Marc
--
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:
