Christian Probst
unread,Feb 16, 2011, 1:42:59 PM2/16/11Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to varid-community
Hi,
I am using varid 1.0.7f and I have received the following error. I
have a very small reference file, and the number of generated reads is
very large, so I think that could be the cause.
Thanks in advance and here is the code:
varid_exec -r /solid/reference/Pool02.fasta -a
Seq01_Pool002_SHRiMP.sam -o Seq01_Pool002_SNPs.pb --threads 4
Using 4 threads
- alignments: Seq01_Pool002_SHRiMP.sam
- ref-file: /solid/reference/Pool02.fasta
Detected:
Pool02:40146 nt
Splitting reads
Total reads: 1124890
Used reads: 1124890
Done: in: 15.274 secs
------------------------------
Processing fasta entry : Pool02 - 40146bp
Processing group number : 1
Loading reads ...
Done in 14.582 secs
Aligning reads ...
There are 598 potential insert sites
Done in: 5.698 secs
Starting VARiD Algorithm
Building Transition Matrix...
Done in: 0.081 secs
Trying to allocate 52967200 bytes...
covmaps built in: 18.856 secs, Used reads in the covmaps:
1124890/1124890, Skipped for 0 index:0
Computing emissions...
SNP penalty 1.000000e-03, log -6.907755e+00
Done in: 47.723 secs,
fwd algo done in: 40.533 secs
bwd algo done in: 25.573 secs
DONE VARiD Algorithm & prediction in: 137.784 secs
Starting primary prediction... Done in: 14.769 secs
<nuc2col.c:143> assertion: (outcol_tmp.data[i] != (VA_UINT8)-1)
failed, Last errno=34
<run_wrap.c:781> varid_exec exiting...
<run_wrap.c:784> varid_exec exited