[mira_talk] Assembling with no quality file. System requesting quality file even when indicated the contrary.

  • From: Patricia Carvajal <patriciacarvajal.mxli@xxxxxxxxx>
  • To: mira_talk@xxxxxxxxxxxxx
  • Date: Mon, 5 Oct 2015 22:40:13 -0600

Hello!!

I'm trying to de novo assembly a set of mRNA sequences of a fasta (.fa)
file, there's no quality file on these sequences.

Per the assembly manaual the input file is a .fa file, the --noqualities
and -CL:pec=false parameters are stated in the manifest file as well as a
default_qual=30.

When assemblying the system is stil requesting a quality file for the fasta
file.

I'm using Mira 4.9.5_2 version (the defaul version of BioLinux 8)

############################
My manifest file is as follows...
################################

# First part: defining some basic things
project = Cqd_assembly
job = est,denovo,accurate
parameters = COMMON_SETTINGS -GENERAL:number_of_threads=2 \
-OUTPUT:output_result_caf=no
parameters= TEXT_SETTINGS -ALIGN:min_relative_score=40 \
-ASSEMBLY:minimum_read_length=50\
--noqualities\
--noclipping=all\
-CL:pec=false

# The second part
readgroup = RNAs_from_Cqd
data = /cqd-dir/Cqd-mRNAs.fa
default_qual=30
technology = text


#########################################

Here's the error I'm getting...
##########################################

********************************************************************************
* FASTA quality file expected to exist, but no quality filename given???
Refer *
* to
*
*
http://mira-assembler.sourceforge.net/docs/DefinitiveGuideToMIRA.html#sect_ref_manifest_readgroups*
* on how to load FASTA data without corresponding quality files.

####################################3

Please advise.


THANKS =)

Other related posts: