ABySS error at rescaffolding with long sequences
0
0
Entering edit mode
8.7 years ago
wstfljs ▴ 100

Hi! I'm running the abyss-pe script with the 'long' parameter to enable rescaffolding of my assembly with longer contigs. Unfortunately, the process stops at this very step. The problem is most likely related to the 'bwa mem' command. When comparing how abyss-pe executed 'bwa mem' with what the command is supposed to look like, I see that one of the input files is missing.

bwa mem -a -t24 -S -P -k63 myassembly_k63-8.fa  \
        |gzip >../long_contigs.fasta-8.sam.gz
bwa mem -a -t$j -S -P -k$l $(name)-8.fa $(strip $($*)) \
        |gzip >$@

As you can see, the following fragment of the code $(strip $($*)) returns nothing and due to insufficient number of input parameters, the 'bwa mem' exits:

Usage: bwa mem [options] <idxbase> <in1.fq> [in2.fq]
abyss assembly • 2.1k views
ADD COMMENT
0
Entering edit mode

Hi Jakalssj3,

I've experienced the same thing (I just posted my question here ABySS fails with -long reads at BWA step). Did you ever get this issue resolved?

Thanks!

ADD REPLY

Login before adding your answer.

Traffic: 1490 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