seqdbio.cpp(221) assert failed: BufferCount > 0 - interesting error!
Violetta Florova
Hello!
I've got quite interesting error in 4 markers library out of 53.
I run this K00370 markers.faa file (and 3 more) on WGS SRR6131123.fastq - and have got exit code, although other 49 markers work good.
Could you provide any ideas why does it happen?
I use this marker library before and they all worked.
The only one issue of this fast file - its the from Hi_C deconvolution study, so it's bigger than usual WGS.
Thank you so much!
Tested usearch. Appears to be working.
Treating input as a wgs file...
usearch v8.0.1623_i86linux64, 24.5Gb RAM, 8 cores
(C) Copyright 2013-15 Robert C. Edgar, all rights reserved.
Licensed to:
00:00 37Mb 0.1% Reading K00370markers.faa^M00:00 37Mb 100.0% Reading K00370markers.faa
00:00 3.3Mb 0.1% Masking ^M00:00 3.3Mb 100.0% Masking
00:00 55Mb 0.1% Word stats^M00:00 55Mb 100.0% Word stats
00:00 55Mb 100.0% Alloc rows
00:00 55Mb 0.1% Build index^M00:00 55Mb 100.0% Build index
00:00 42Mb 0.0% Rows ^M00:00 42Mb Buffers
List of files in WGS set:SRR6131123.fastq
List of files in WGS set (after unpacking tarfiles):SRR6131123.fastq
Working on file 1 of 1
usearch v8.0.1623_i86linux64, 24.5Gb RAM, 8 cores
(C) Copyright 2013-15 Robert C. Edgar, all rights reserved.
00:00 42Mb 0.0% Rows^M00:01 42Mb 16.2% Rows^M00:02 42Mb 76.7% Rows^M00:02 42Mb 100.0% Rows
00:02 42Mb Reading pointers...done.
/projects/academic/mjbuck/Buck_Lab_Members/violetta/new_project/src/usearch --usearch_local tmp303221537839651544/fasta.fna --db K00370markers.faa.udb --id 0.95 --userout tmp303221537839651544/wgs_01out_$
---Fatal error---
seqdbio.cpp(221) assert failed: BufferCount > 0
('Using this version of usearch: ', u'v8.0.1623')
Traceback (most recent call last):
File "/projects/academic/mjbuck/Buck_Lab_Members/violetta/new_project/shortbred_quantify_test.py", line 545, in <module>
iThreads=args.iThreads,dID=args.dID, iAccepts=args.iMaxHits, iRejects=args.iMaxRejects,strUSEARCH=args.strUSEARCH )
File "/projects/academic/mjbuck/Buck_Lab_Members/violetta/new_project/src/quantify_functions.py", line 256, in RunUSEARCH
"--maxrejects",str(iRejects),"--threads", str(iThreads)])
File "/util/academic/biopython/1.70/anaconda-5.0.1/lib/python2.7/subprocess.py", line 186, in check_call
raise CalledProcessError(retcode, cmd)
subprocess.CalledProcessError: Command '['/projects/academic/mjbuck/Buck_Lab_Members/violetta/new_project/src/usearch', '--usearch_local', 'tmp303221537839651544/fasta.fna', '--db', 'K00370markers.faa.ud$