This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Change Residue Output Length in trimAl?

Hello,

I am looking for a way to change the length of the residues printed in the output of trimAl.

Using clustalo, I have generated .aln data that looks like so:

chr1_38402816-38409352                     ---------------------------------------------------------------------------------agatcatgaacttttgctc--ca---agaactattgcaggaa-cataccagga-aagccgagagaatccacagaccctttgaaggaagtggatt
chr1_41176265-41182362                     -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
chr1_42059276-42065372                     -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

I was able to see 175 residues per line by using the --wrap argument in clustalo:

clustalo --in=subsub.noline.fa --out=subsub.noline.aln --force --outfmt=clustal --wrap=175 --threads=6 --verbose

Now I am trimming regions using trimAl, and the output looks like this, with the standard 60 bp length:

chr1_440936-447357                        gggactggtcagtgggt-------------gg-tgtgaagtgccccctagag-------g
chr1_675912-682333                        gggactggtcagtgggt-------------gg-tgtgaagtgccccctagag-------g
chr1_25506586-25512707                    ggaacagcttacagctccagaagacgggtgatttctgaggta-tcactagagtgcca-gg

So, I'm wondering if trimAl has something like "wrap". I have searched the documentation and the supplementary documentation and can't find an answer, although I am new to this so I might have just missed it. Surely it is coded somewhere! Any help would be greatly appreciated.

TIA!

msa clustalo trimal trimal

0 answers

No answers yet.

Log in to answer this question.