Sorry, I try to use a script that generates the manifest file for each
metagenome project and then executes MIRA.
The script is this:
cat >mira.config <<EOConfig
parameters = -DI:trt=/tmp/ -GE:not=4
-OUT:output_result_caf=no:output_result_maf=yes:output_result_fasta=yes:output_result_ace=yes:output_result_tcs=no
-NW:cmrnl=no\
454_SETTINGS = -AL:min_overlap=30:min_relative_score=75
-OUT:savesimplesingletsinproject=yes -LR:wqf=no -AS:epoq=no
-CL:qc=no:bsqc=no\
project = Mira_454
job = genome,denovo,draft
454_SETTINGS= -AL:min_overlap=30:min_relative_score=75
OUT:savesimplesingletsinproject=yes -LR:wqf=no -AS:epoq=no
-CL:qc=no:bsqc=no\
# now the unpaired data
readgroup = unpaired
data = $IN1
technology = 454
template_size = 50 1500 autorefine
default_qual
EOConfig
# Execute mira
mira mira.config >&log_assembly.txt
Thanks
Adriana M. Froes
Laboratório de Microbiologia, Instituto de Biologia, Depto de Biologia
Marinha
Universidade Federal do Rio de Janeiro
Av. Carlos Chagas Filho 373, Sala A3-202, Bloco A (Anexo) do CCS
21941-599, Ilha do Fundão, Rio de Janeiro, RJ
On Tue, May 5, 2015 at 7:03 PM, Bastien Chevreux <bach@xxxxxxxxxxxx> wrote:
On 05 May 2015, at 22:55 , Adriana Fróes <dricamfroes@xxxxxxxxx> wrote:
Thank you very much Marcelo!
I've tried one, actually, but it keeps showing the same error, as:
[…]
So you showed the error which, indeed, looks very strange. How about
showing the manifest file for us to see what the problem could be? ;-)
B.
--
You have received this mail because you are subscribed to the mira_talk
mailing list. For information on how to subscribe or unsubscribe, please
visit http://www.chevreux.org/mira_mailinglists.html