Losing all data after filter step using SGA
1
0
Entering edit mode
7.1 years ago
jolespin ▴ 130

Raw fastq > preprocess > index > correct > index = processed_reads.ec.fa (6.5 GB)

 

I've taken the processed_reads.ec.fa and ran the following command:

 

sga filter -k 31 processed_reads.ec.fa

 

stderr:

sga: QCProcess.cpp:233: DuplicateCheckResult QCProcess::performDuplicateCheck(const SequenceWorkItem&): Assertion `fwdIntervals.interval[0].isValid() || rcIntervals.interval[0].isValid()' failed.

 

 

stdout:

...

[sga] Processed 8000000 sequences (912.782961s elapsed)

[sga] Processed 8050000 sequences (918.024303s elapsed)

[sga] Processed 8100000 sequences (925.718005s elapsed)

[sga] Processed 8150000 sequences (932.858792s elapsed)

Abort (core dumped)

 

 

The output files were:

processed_reads.ec.filter.pass.fa (1.4 MB)

processed_reads.ec.discard.fa (2.0 GB)

 

Does anyone know what's happening? What are these errors? I've tried different thread sizes and different kmer sizes and have not seen any significant improvements...

SGA assemly filter kmer • 1.5k views
ADD COMMENT
1
Entering edit mode
7.1 years ago
jolespin ▴ 130

Needed the BWT file from the reindexing step

ADD COMMENT

Login before adding your answer.

Traffic: 820 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6