Merging of contigs files by minimus2
1
1
Entering edit mode
5.4 years ago
waqasnayab ▴ 200

Hi,

I am following the link:

http://amos.sourceforge.net/wiki/index.php/Minimus2

for the merging of two contig files from velvet. It is written on the link: "......for merging one or two sequence sets (S1,S2)." here, seqeucnes sets means contig files in fasta format.???

I am facing the following error during merging:

command is: ./minimus2 all_velvet -D REFCOUNT=213048

error: The log file is: all_velvet.runAmos.log
Doing step 10:  Building AMOS bank & Dumping reads
Doing step 11
Doing step 12
Doing step 13
Doing step 20:  Getting overlaps
Doing step 21
Command: /usr/bin/show-coords -H -c -l -o -r -I 94 all_velvet.delta | /home/kamran/Desktop/amos-3.1.0/bin/nucmerAnnotate | egrep 'BEGIN|END|CONTAIN|IDENTITY' > all_velvet.coords  exited with status: 1

Help required,

Regards,

Waqas.

Assembly software error • 3.0k views
ADD COMMENT
0
Entering edit mode

I think there may be no over lapping sequences to connect

ADD REPLY
5
Entering edit mode
4.1 years ago
JV ▴ 430

I know this post is OLD. Just answering anyway, in case anybody else has problems with this.

I got the same error message when testing minimus2 recently, and it turns out this was NOT because no overlaps would be detected (would have hoped for a more informative message on that, anyway). Turns out, the standard error-output of minimus2 is not very helpful. You have to check the ".runAmos.log" file for actual info on what went wrong. (in your case that would be "all_velvet.runAmos.log".

In my case I found that minimus was expecting the tools delta-filter and show-coords at /usr/local/bin. However this may not be the case if you installed amos locally on a shared server, where you have no admin rights (as I did via bioconda).

Therefore you may have to manually edit the minimus2 script at lines 46 & 47 to point to the correct executables. Solved the problem for me!

ADD COMMENT
0
Entering edit mode

But still Thanksssss

ADD REPLY

Login before adding your answer.

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