ERROR: Bad base (.) at site 6 of sequence 2

270 views
Skip to first unread message

Yuen

unread,
Apr 4, 2017, 7:30:11 AM4/4/17
to raxml
Hello. I've tried several times generating RAxML-HPC2 on XSEDE with 64 of my sequences but I could not get the RAxML_bipartitions or any output files as I named. Instead, I checked the stdout.txt file and they all gave the same error which was as follows. What should I do? 

CIPRES_THREADSPP=4
CIPRES_NP=6
running:
 raxmlHPC-HYBRID -T 4 -n result -s infile.txt -c 25 -m GTRCAT -p 12345 -f a -N 1000 -x 12345 -o HM240413.1 Paraboysidia tarutao voucher MOL 116412,HM240412.1 Paraboysidia tarutao voucher MOL 116411,HM240411.1 Paraboysidia tarutao voucher MOL 116385,HM240410.1 Gyliotrachela depressispira voucher MOL 116409,HM240409.1 Gyliotrachela frequens voucher MOL 116407,HM240408.1 Gyliotrachela salpinx voucher MOL 116372,HM240407.1 Gyliotrachela salpinx voucher MOL 116371,HM240406.1 Gyliotrachela transitans voucher MOL 116382,HM240405.1 Gyliotrachela transitans voucher MOL 116381
CIPRES_THREADSPP=4
CIPRES_NP=6
running:
 raxmlHPC-HYBRID -T 4 -n result -s infile.txt -c 25 -m GTRCAT -p 12345 -f a -N 1000 -x 12345 -o HM240413.1 Paraboysidia tarutao voucher MOL 116412,HM240412.1 Paraboysidia tarutao voucher MOL 116411,HM240411.1 Paraboysidia tarutao voucher MOL 116385,HM240410.1 Gyliotrachela depressispira voucher MOL 116409,HM240409.1 Gyliotrachela frequens voucher MOL 116407,HM240408.1 Gyliotrachela salpinx voucher MOL 116372,HM240407.1 Gyliotrachela salpinx voucher MOL 116371,HM240406.1 Gyliotrachela transitans voucher MOL 116382,HM240405.1 Gyliotrachela transitans voucher MOL 116381
CIPRES_THREADSPP=4
CIPRES_NP=6
running:
 raxmlHPC-HYBRID -T 4 -n result -s infile.txt -c 25 -m GTRCAT -p 12345 -f a -N 1000 -x 12345 -o HM240413.1 Paraboysidia tarutao voucher MOL 116412,HM240412.1 Paraboysidia tarutao voucher MOL 116411,HM240411.1 Paraboysidia tarutao voucher MOL 116385,HM240410.1 Gyliotrachela depressispira voucher MOL 116409,HM240409.1 Gyliotrachela frequens voucher MOL 116407,HM240408.1 Gyliotrachela salpinx voucher MOL 116372,HM240407.1 Gyliotrachela salpinx voucher MOL 116371,HM240406.1 Gyliotrachela transitans voucher MOL 116382,HM240405.1 Gyliotrachela transitans voucher MOL 116381
CIPRES_THREADSPP=4
CIPRES_NP=6
running:
 raxmlHPC-HYBRID -T 4 -n result -s infile.txt -c 25 -m GTRCAT -p 12345 -f a -N 1000 -x 12345 -o HM240413.1 Paraboysidia tarutao voucher MOL 116412,HM240412.1 Paraboysidia tarutao voucher MOL 116411,HM240411.1 Paraboysidia tarutao voucher MOL 116385,HM240410.1 Gyliotrachela depressispira voucher MOL 116409,HM240409.1 Gyliotrachela frequens voucher MOL 116407,HM240408.1 Gyliotrachela salpinx voucher MOL 116372,HM240407.1 Gyliotrachela salpinx voucher MOL 116371,HM240406.1 Gyliotrachela transitans voucher MOL 116382,HM240405.1 Gyliotrachela transitans voucher MOL 116381
CIPRES_THREADSPP=4
CIPRES_NP=6
running:
 raxmlHPC-HYBRID -T 4 -n result -s infile.txt -c 25 -m GTRCAT -p 12345 -f a -N 1000 -x 12345 -o HM240413.1 Paraboysidia tarutao voucher MOL 116412,HM240412.1 Paraboysidia tarutao voucher MOL 116411,HM240411.1 Paraboysidia tarutao voucher MOL 116385,HM240410.1 Gyliotrachela depressispira voucher MOL 116409,HM240409.1 Gyliotrachela frequens voucher MOL 116407,HM240408.1 Gyliotrachela salpinx voucher MOL 116372,HM240407.1 Gyliotrachela salpinx voucher MOL 116371,HM240406.1 Gyliotrachela transitans voucher MOL 116382,HM240405.1 Gyliotrachela transitans voucher MOL 116381
CIPRES_THREADSPP=4
CIPRES_NP=6
running:
 raxmlHPC-HYBRID -T 4 -n result -s infile.txt -c 25 -m GTRCAT -p 12345 -f a -N 1000 -x 12345 -o HM240413.1 Paraboysidia tarutao voucher MOL 116412,HM240412.1 Paraboysidia tarutao voucher MOL 116411,HM240411.1 Paraboysidia tarutao voucher MOL 116385,HM240410.1 Gyliotrachela depressispira voucher MOL 116409,HM240409.1 Gyliotrachela frequens voucher MOL 116407,HM240408.1 Gyliotrachela salpinx voucher MOL 116372,HM240407.1 Gyliotrachela salpinx voucher MOL 116371,HM240406.1 Gyliotrachela transitans voucher MOL 116382,HM240405.1 Gyliotrachela transitans voucher MOL 116381

This is RAxML MPI Process Number: 0

This is RAxML MPI Process Number: 1

This is RAxML MPI Process Number: 2

This is RAxML MPI Process Number: 3

This is RAxML MPI Process Number: 4

This is RAxML MPI Process Number: 5
ERROR: Bad base (.) at site 6 of sequence 2
Printing error context:

CT
142_Gyliot -----..... .......... ....

Problem reading alignment file 
ERROR: Bad base (.) at site 6 of sequence 2
Printing error context:

CT
142_Gyliot -----..... .......... ....

Problem reading alignment file 
ERROR: Bad base (.) at site 6 of sequence 2
Printing error context:

CT
142_Gyliot -----..... .......... ....

Problem reading alignment file 
ERROR: Bad base (.) at site 6 of sequence 2
Printing error context:

CT
142_Gyliot -----..... .......... ....

Problem reading alignment file 
ERROR: Bad base (.) at site 6 of sequence 2
Printing error context:

CT
142_Gyliot -----..... .......... ....

Problem reading alignment file 
ERROR: Bad base (.) at site 6 of sequence 2
Printing error context:

CT
142_Gyliot -----..... .......... ....

Problem reading alignment file 

Alexandros Stamatakis

unread,
Apr 4, 2017, 1:11:04 PM4/4/17
to ra...@googlegroups.com
replace . by N in your alignent file ...

alexis
> --
> 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
> <mailto:raxml+un...@googlegroups.com>.
> For more options, visit https://groups.google.com/d/optout.

--
Alexandros (Alexis) Stamatakis

Research Group Leader, Heidelberg Institute for Theoretical Studies
Full Professor, Dept. of Informatics, Karlsruhe Institute of Technology
Adjunct Professor, Dept. of Ecology and Evolutionary Biology, University
of Arizona at Tucson

www.exelixis-lab.org

Gianmarco Virga

unread,
Jun 5, 2017, 8:34:59 AM6/5/17
to raxml
I've got the same problem: "ERROR: Bad base <2> at site 47 of sequence 1"
printing error context:
 
 0 ? ? ? ? ? ? 0 0 2 ? 0
t1                   0 0 ? 0 ?

what should i do?
Thanks in advance

Alexandros Stamatakis

unread,
Jun 8, 2017, 12:07:49 AM6/8/17
to ra...@googlegroups.com
can you please send us the full alignment file such that we can check?

alexis
Reply all
Reply to author
Forward
0 new messages