Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

duplicate sequence IDs #30

Open
hildebra opened this issue Apr 29, 2017 · 0 comments
Open

duplicate sequence IDs #30

hildebra opened this issue Apr 29, 2017 · 0 comments

Comments

@hildebra
Copy link

Dear Salvador,
I get a segmentation fold in trimal, when running the following command:

trimal -in COG0185.0.faa -out COG0185.0.fna -backtrans inMSA0.fna -ignorestopcodon -gt 0.1 -cons 60

Without the -backtrans option, the program runs fine. I was hoping you could help me.

All the best,
Falk Hildebrand

Answer from Salvador:
Dear Falk,

Thanks for using trimAl and contacting me regarding this unexpected behaviour.

I played a bit with your input file and realized that you have some repeated IDs for the nucleotide files ...
2 >394503_COG0185
2 >411474_COG0185
2 >445973_COG0185
2 >699246_COG0185
2 >718252_COG0185
... and for the protein files:
2 >394503_COG0185

I just kept the first appearance of such sequences in the attached files, and everything worked as expected.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant