Question: multifasta length uniformity
gravatar for Unalibun
6 days ago by
Unalibun0 wrote:

Hi everyone, I am performing the phylogenetic reconstruction of my non-model organism, I am working with the BUSCOs output (single copy sequences files) form 7 different species, I have more than 1250 small multi-fasta files generated as output by Guidance2 software.

Ej EOG092D4S9W.faa.aln




I am trying to do my Bootstrapping but the sequence generated are not accepted by raxlm program I have to uniform them to the same length but I can not find a program to do that recursively.

This is how the multifasta should looks like. However, every single files has sequences with different lengths so I can not choice an standard number to trimming the sequences. I think that the solution could be calculate the length smaller sequence and then cutting the others using the calculated length as a reference.




Can anybody help to untingling this issue me please

alignment • 72 views
ADD COMMENTlink modified 6 days ago by Mensur Dlakic4.1k • written 6 days ago by Unalibun0

You should align the sequences, not trim them.

ADD REPLYlink modified 6 days ago • written 6 days ago by Asaf7.2k

Please use the formatting bar (especially the code option) to present your post better. You can use backticks for inline code (`text` becomes text), or select a chunk of text and use the highlighted button to format it as a code block. I've done it for you this time.

P.S. length is spelled leng-T-H, not leng-H-T. I've corrected those in the post too.

ADD REPLYlink written 6 days ago by RamRS26k
gravatar for Mensur Dlakic
6 days ago by
Mensur Dlakic4.1k
Mensur Dlakic4.1k wrote:

A proper way to do this is by aligning sequences first, and trimming them afterwards so they are all the same length. Some positions in the alignment will be gaps rather than letters, which is normal. I have responded recently to a similar question here and here.

ADD COMMENTlink written 6 days ago by Mensur Dlakic4.1k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1049 users visited in the last hour