[mira_talk] Parameters not parsed

Hi mira-friends

I tried to tweak the parameters a little bit, because I'm getting very
strange results with one of my assemblies. It seems that mira is
unwilling to parse my parameters (darn). 

My project:
- sequencing a complete bacterial genome (~2Mb)
- half a 454 Titanium run (~100Mb of data, not quite as much as
advertised), plus 2 lanes from a titration run.
- about 50k paired-end mates (experimental 20kb library)

The bacteria that I'm sequencing experiences some kind of weird
alternate replication in late growth-phase, leading to an (at least) 2
fold amplification of about 10-20% of the genome, which may impair the
uniform read distribution (I tried with urd:yes, and yielded strange
results).

Here is some parts of the assembly_log:

<code>
-SB:sbuip is 3, but must be no more than 2. Setting to 2
</code>

Didn't say anything about strain or backbone... strange, but shouldn't
cause any problems...

<code>
Parsing parameters: --project=BAnh1NrPeM01 -fasta -OUT:ora=yes -GE:not=4
-AS:urd=no --job=denovo,genome,454,normal
</code>

Fine, that's what I want. But then:

<code>
Parameter settings seen for:
Sanger data (also common parameters), 454 data

Used parameter settings:
  General (-GE):
        Project name in (proin)                  : BAnh1NrPeM01
        Project name out (proout)                : BAnh1NrPeM01
        Number of threads (not)                  : 2
</code>

Uh? I said 4... But it's still OK...

<code>
  Assembly options (-AS):
        [...]
        Automatic repeat detection (ard)            : yes
            [...]
            Use uniform read distribution (urd)     : yes
              Start in pass (urdsip)                : 5
              Cutoff multiplier (urdcm)             :  [san]  1.5
                                                       [454]  1.5
</code>

Uh??? I asked not to use urd...

Any suggestion? Am I doing something wrong?

Thanks for your help... I might come here again to discuss my weird
assembly.

Lionel 

Uppsala University




-- 
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

Other related posts: