Phylogenetic Analysis By Maximum Likelihood (Paml) : Upon Phylip Input Yn00.Clt Program Not Running.
2
0
Entering edit mode
8.8 years ago
Nari ▴ 900

Dear all,
I want to use PAML for Checking positive selection by dn/ds calulation.
I am using Phylip alignment file. (Which PAML suggests to be used as input).
This is my input Phylip file.

 7   1065
MPREF94h99         ---------- ---------- ---------- ---------- ---------- 
MPREF9b499         ---------- ---------- ---------- ---------- ---------- 
MPREF03499         ---------- ---------- ---------- ---------- ---------- 
MPREF949l9         AGCGAAAAAA GTTTTAAACG TTGAAAGCTC CTATCTAGAT AATATGATTC 
MPREF9c499         ---------- ---------- ---------- ---------- ---------- 
MPREF949d9         ---------- ---------- ---------- ---------- ---------- 
MPREF0t349         ---------- ---------- ---------- ---------- ---------- 
                   ---------- ---------- ---------- ---------- ---------- 
                   ---------- ---------- ---------- ---------- ---------- 
                   ---------- ---------- ---------- ---------- ---------- 
                   TTTGGGGTAA ACGCCCCTTG ATAGATTCTA TCAACGATAT TAAAGATTAC 
                   ---------- ---------- ---------- ---------- ---------- 
                   ---------- ---------- ---------- ---------- ---------- 
                   -------TTG AAAATATTGT TATATTTTGA AGGTGAGAAA ATCTTAGCTA

And this is cmd error message.

C:\Users\Lenovo\Desktop\dnds\PAML\paml4.6>yn00 YN00 in paml version 4.6, August 2012

Reading options from yn00.ctl..

ns = 7 ls = 1065 Reading sequences, sequential format..

Error in sequence data file: P at 52 seq 1. Make sure to separate the sequence from its name by 2 or more spaces.

(Note : there are more than two spaces between name and sequence in my file)

C:\Users\Lenovo\Desktop\dnds\PAML\paml4.6>

I am not getting output after running file yn00.clt
What may be the problem?

paml • 8.9k views
ADD COMMENT
0
Entering edit mode

you are sure they are spaces and not tabs?

ADD REPLY
0
Entering edit mode

yeah! as in original file, there was only one space. I inserted them with notepad.

ADD REPLY
4
Entering edit mode
8.8 years ago
Whetting ★ 1.5k

It looks like PAML is looking for a sequential file "ns = 7 ls = 1065 Reading sequences, sequential format..".
However, you are providing an interleaved file. Try converting it to sequential and see what happens

ADD COMMENT
1
Entering edit mode

Changing the header to " 7 1065 I" may be an easier test of the format. That should tell the program it's interleaved.

ADD REPLY
0
Entering edit mode

PAML is really finicky...I am not sure that will work, but it is worth the shot

ADD REPLY
3
Entering edit mode

I actually run codeml and baseml using interleaved files with the header formatted like that (with an "I"). This was pointed out to me by a colleague about 5 years who had a script what would just do this simple reformatting. There may be other issues with the file, but this formatting works for me.

ADD REPLY
0
Entering edit mode

I stand corrected...and that is good to know

ADD REPLY
0
Entering edit mode

Thanks SES: i did that now I'm getting this msg:

C:\Users\Lenovo\Desktop\dnds\PAML\paml4.6>yn00 YN00 in paml version 4.6, August 2012

Reading options from yn00.ctl.. ns = 7 ls = 1065 Reading sequences, interlaved format.. Counting site patterns.. 0:00 26 patterns at 26 / 26 sites (100.0%), 0:00 stop codon TAA in seq # 1: MPREF94h99 codons in other sequences are seq # 1 MPREF94h99 TAA seq # 2 MPREF9b499 TGA seq # 3 MPREF03499 TTC seq # 4 MPREF949l9 TGC seq # 5 MPREF9c499 TGC seq # 6 MPREF949d9 GAA seq # 7 MPREF0t349 TGA

ADD REPLY
1
Entering edit mode

PAML does not like stop codons. Are these stopcodons internal or terminal? If terminal just remove them. If they are internal, you may need to rethink your analysis...

ADD REPLY
0
Entering edit mode
7.0 years ago
zeeefa ▴ 90

This might help: http://bit.ly/1utdbFg itĀ converts FASTA alignment file to a non-interleaved/sequential Phylip format :)

ADD COMMENT

Login before adding your answer.

Traffic: 1789 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6