raxmlHPC ERROR: Bad base (A) at site 1 of sequence 1

797 views
Skip to first unread message

matthe...@gmail.com

unread,
Dec 28, 2018, 2:33:13 AM12/28/18
to raxml
Hello, my raxml raised error, I google a lot and no solutions can fix this, I hope someone know how to fix. Thanks.

My command:
raxmlHPC -s COX1_alignment.phy -n COX1 -m ASC_BINCATX -b 12345 -p 12345 -# 1000

Error msg:
ERROR: Bad base (A) at site 1 of sequence 1
Printing error context:

 1614
57          ATTATTAGAC GGTTTGAAAG

Problem reading alignment file


My alignment file looks:

18   1614
57          ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
27          ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
62          ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
58          ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
82          ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
64          ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
100         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
109         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
115         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
112         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
118         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
119         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
126         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGATGTTTGC GGTGATTATA TAGGACAAAT
122         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
121         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
130         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
129         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT
128         ATTATTAGAC GGTTTGAAAG GTTTCGTTTG TGACGTTTGC GGTGATTATA TAGGACAAAT

            CATAAGGATA TTGGGACTCT ATATTTTCTT TTTTCTATTT GATCTGGTTT AATAGGAACT
            CATAAGGATA TTGGGACTTT ATACTTTCTT TTTTCTATTT GATCTGGTTT AATAGGAACT
            CATAAGGATA TTGGGACTCT ATATTTTCTT TTTTCTATTT GATCTGGTTT AATAGGAACT
            CATAAGGATA TTGGGACTCT ATACTTTCTT TTTTCTATTT GATCTGGTTT AATAGGAACT
            CATAAGGATA TTGGGACTTT ATACTTTCTT TTTTCTATTT GATCTGGTTT AATAGGAACT
            CATAAGGATA TTGGGACTCT ATACTTTCTT TTTTCTATTT GATCTGGTTT AATAGGAACT
            CATAAGGATA TTGGGACTCT ATACTTTCTT TTTTCTATTT GATCTGGTTT AATAGGAACT


Alexey Kozlov

unread,
Dec 28, 2018, 6:35:10 AM12/28/18
to ra...@googlegroups.com

You get this error because you specify binary model (-m ASC_BINCATX), but your alignment contains DNA data. You should use nucleotide model instead, e.g. ASC_GTRCATX.

Hope this helps,
Alexey

--
You received this message because you are subscribed to the Google Groups "raxml" group.
To unsubscribe from this group and stop receiving emails from it, send an email to raxml+un...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

matthe...@gmail.com

unread,
Dec 29, 2018, 5:02:17 AM12/29/18
to raxml
Thanks.

在 2018年12月28日星期五 UTC+8下午7:35:10,Alexey Kozlov写道:
Reply all
Reply to author
Forward
0 new messages