SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Amos/MuMMer/Minimus2 postnuc running multikmer assembly milo0615 Bioinformatics 0 11-11-2014 01:48 PM
Mixing rapid mode and high output mode libs in one analysis merilius Illumina/Solexa 0 11-04-2014 06:01 AM
Bowtie2: mapping qualities in paired end mode different with the single mode bill_clinton Bioinformatics 1 09-25-2013 11:47 AM
BLAT or BLAST: output SNPs between reference and query? bbimber Bioinformatics 1 02-23-2010 10:59 AM

Reply
 
Thread Tools
Old 10-22-2017, 04:45 PM   #1
thakurshalabh
Junior Member
 
Location: Canada

Join Date: Oct 2017
Posts: 2
Default Running Minimus2 in reference vs query mode

Hello all,

I have recently started to use minimus2 for merging contigs from different assemblies. I followed the steps shown on its tutorial page and merged sequences from two assemblies say (S1.fasta with 28095 sequences) and (S2.fasta with 31623 sequences) together into one file S1-S2.seq. I was also able to run the toAmos command on S1-S2.seq and generate an S1-S2.afg file. However, when I run minimus2 on S1-S2.afg file using command "minimus2 S1-S2 -D REFCOUNT=28095"
the program exit at step 13 showing an error "nucmer -maxmatch -c 40 S1-S2.ref.seq S1-S2.qry.seq -p S1-S2 exited with status: 1".

On investigating the runAmos.log file I found that step 13 was to run dumpread program and create S1-S2.qry.seq file (as shown by log file below). But dumpread failed to write any sequence in S1-S2.qry.seq file which may have lead to above error.

!!! 2017-10-22 18:00:54 Doing step 13
!!! 2017-10-22 18:00:54 Running: dumpreads -m 28095 S1-S2.bnk > S1-S2.qry.seq
Objects seen: 94
Objects written: 0


I'm not sure what wrong is happening here. Because when REFCOUNT=0, the program runs perfectly fine. Can anyone suggest me how to fix this or what might be the cause of this error?
thakurshalabh is offline   Reply With Quote
Old 10-23-2017, 12:47 AM   #2
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default

Should the refcount be the sum of the count of the two assemblies (i.e. 60528, not 28905)?
gringer is offline   Reply With Quote
Old 10-23-2017, 07:53 AM   #3
thakurshalabh
Junior Member
 
Location: Canada

Join Date: Oct 2017
Posts: 2
Default

Thank you for your suggestion. But as per the instructions in the tutorial, the value of the REFCOUNT parameter should be equal to the number of sequences in the reference assembly which in this case is assembly S1.

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

Also, I tried with different values for REFCOUNT (1, 10, 20, 25, ...) and it worked fine. So does that mean 28905 is too high for REFCOUNT? Then, how should I decide what value to choose?
thakurshalabh is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:16 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO