Question: A5 Pipeline : Idba To Crash
0
gravatar for HG
6.7 years ago by
HG1.1k
Germany
HG1.1k wrote:

Hi everyone, i am using A5 pipeline for Miseq 250bp pairend reads. Getting some error message like this :

1876358    Sequences
9744    Tags rejected (0.5%) at 29.0% coverage cutoff (0.010000 FDR).
    Elapsed time:
        131 seconds (130.77 CPU+SYS seconds)
[a5_s2] Building contigs from mygenome.ec.fastq with IDBA
[a5_s2] Building contigs from mygenome.ec.fastq with IDBA
[a5] idba -long mygenome.s2/mygenome.ec.fasta -o mygenome.s2/mygenome --mink 25 --maxk 90
terminate called after throwing an instance of 'std::logic_error'
  what():  SequenceReader::SequenceReader() istream is invalid
Aborted (core dumped)
[a5] Error building contigs with IDBA

Please help me out

• 2.1k views
ADD COMMENTlink modified 6.7 years ago • written 6.7 years ago by HG1.1k

You should have added to your question or started a new post instead of editing the title and completely changing the question. Now, it is unclear what my answer is about (even though it solved your problem) and if someone else has this issue, it won't be as easy for them to find a solution.

ADD REPLYlink modified 6.7 years ago • written 6.7 years ago by SES8.4k
0
gravatar for SES
6.7 years ago by
SES8.4k
Vancouver, BC
SES8.4k wrote:

That message is saying the quality encoding of your data is not what is expected by the program. This could be do to not setting an argument to the program correctly, or it could be a bug. However, this message in the output:

"Check your data and re-run preprocess with the correct quality scaling flag."

seems to indicate it is the first of those two possible issues I mentioned. So, try to rerun the preprocess program with the appropriate settings. If that does not solve your problem, you may want to read the conversation about issue #5 reported to the mailing list. It appears similar issues have been seen in the past, and there are ways to solve them.

ADD COMMENTlink written 6.7 years ago by SES8.4k

Thank you for your suggestion now i have changed sub get_phred6. So far earlier issue solved but now problem with IDBA, i am using current version of IDBA V1.1.1

1876358    Sequences
9746    Tags rejected (0.5%) at 29.0% coverage cutoff (0.010000 FDR).


    Elapsed time:
        131 seconds (130.42 CPU+SYS seconds)
[a5_s2] Building contigs from mygenome.ec.fastq with IDBA
[a5_s2] Building contigs from mygenome.ec.fastq with IDBA
[a5] idba -long mygenome.s2/mygenome.ec.fasta -o mygenome.s2/mygenome --mink 10 --maxk 250
[a5] Error building contigs with IDBA
ADD REPLYlink modified 6.7 years ago • written 6.7 years ago by HG1.1k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1399 users visited in the last hour