Skip to content

Between '-backtrans' and '-nogaps' parameter, Which one to use first? #63

@kishor2019

Description

@kishor2019

Hi,
I am using trimal v1.4.rev15 to trim and back-translate into codon to make the alignment suitable fo PAML. Since running it i have found that -automated1, -backtrans, -nogaps can not be run together. So i decieded to remove the gaps first from the peptide alignment and then to the back translation. My cmd lines look like:

/trimal-1.4.1/source/trimal -in OG0008206.fa.aln.out -out OG0008206.fa.aln.out_1st.trim -phylip_paml -nogaps

/mnt/genome3/Lab_Users/Kishor/DISK_2/softwares/trimal-1.4.1/source/trimal -in OG0008206.fa.aln.out_1st.trim -out OG0008206.fa.aln.out_2nd.trim -automated1 -backtrans OG0008206.fa.codon

I also tried in oposite way by backtranleting first and then remove the gap. But i received the two different set of codons(though the nuc/pep number was same).
Can i get any suggestion?
Thanks.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions