A qusetion about denovo assembly 454 sequence using MIRA

[kentnf]

Hi,

I use mira to do assembling on my 454 sequence, but the problem is if there are two differnet nt in the same position, the contig sequence will be change.

for example:
contig1: AAAARTTTYTCCCCGGGG
member1: AAAAGTTTCTCCCCGGGG
member2: AAAAATTTTTCCCCGGGG

it means
AG =R
CT =Y
AC=M

How can change the R Y M to ATCG?

Thank you.

[coldturkey]

I suggest joining the mira_talk mailing list http://www.chevreux.org/mira_mailinglists.html
You'll get a prompt answer to your question

[kentnf]

hehe, Thank you!

[sklages]

You might want to have a look at:

force_nonIUPACconsensus_perseqtype(fnicpst)=on|yes|1, off|no|0

cheers,
Sven

[kentnf]

Thank you so much! sklages
It work when I set the -job:sanger, but not work well when set -job:454.

The new problems is how to set the ends clip in MIRA,
I try -CL:emlc and -CL:emrc, some of my sequences have different in the ends, but it not work
example is in cap3, I can set -y 15, and if the ends different lower than 15bp , it clipped.
In MIRA, it can not assembly when set -CL:emlc=no,emrc=no

[sklages]

well, it should work with -job:454 as the manual doesn't distinguish
sequencing chemistries at this point.

[kentnf]

well, it should work with -job:454 as the manual doesn't distinguish
sequencing chemistries at this point.

I use 2.9.37 version before, maybe some other parameter different when I set -job:454 or sanger, now I download 2.9.44 version, use -CLec=yes to clip the ends of sequence, but the clipping part is too long, I want set the clipping length lower than 15bp; if clip longer than 15bp,I think the two reads can not assembly. How to set it? Thank you.

[kentnf]

up up up up up