ctb/khmer — Issue #7: FastaParser current_read.name[0]
Author: shandley
Created: 2011-10-27 22:12:20+00:00
State: open

Hi,

I receive the following error when implementing the filter-inexact-all.py script

filtering 32-mers exactly;
from file /home/scott/data/carpi/illumina_13_4/s_6_1_sequence.txt to file s_6_1_filtinexact.all.txt;
looking for sequences containing k-mers, all with count >= 2
python: parsers.cc:42: FastaParser::FastaParser(const string&): Assertion `current_read.name[0] == '>'' failed.
Aborted

Do I need to rename my reads?

The script quality-trim.py works fine on the same data file.

Thanks,

SAH