[Bioperl-l] Fwd: program bp_pairwise_kaks
Jason Stajich
jason at bioperl.org
Tue Jan 20 11:11:29 EST 2009
Markus-
You need to have also installed PAML in order to do the calculation with
PAML. There is example code on the HOWTOs for running the perl-only Ka/Ks
counts methods that are in Bio::Align::DNAStatistics
Jason Stajich
jason at bioperl.org
http://bioperl.org/wiki/User:Jason
---------- Forwarded message ----------
From: Markus Liebscher <Markus.Liebscher at gmx.de>
Date: Tue, Jan 20, 2009 at 6:22 AM
Subject: program bp_pairwise_kaks
To: jason at bioperl.org
Dear Jason,
I found the Bioperl package very useful. But I failed to run the
bp_pairwise_kaks program.
I tried to run bp_pairwise_kaks as suggested using the command line $
bp_pairwise_kaks.pl -i ../worm_fam_2785.cdna.
The program starts the routine but at the end I get this:
--------------------------------------------------------
CLUSTAL W (1.82) Multiple Sequence Alignments
Sequence format is Pearson
Sequence 1: CBG10100 363 aa
Sequence 2: F22B7.13 525 aa
Sequence 3: C38C10.4 525 aa
Start of Pairwise alignments
Aligning...
Sequences (1:2) Aligned. Score: 9
Sequences (1:3) Aligned. Score: 8
Sequences (2:3) Aligned. Score: 96
Guide tree file created: [/tmp/5M8igBfS6R/5LbxDi3Xik.dnd]
Start of Multiple Alignment
There are 2 groups
Aligning...
Group 1: Sequences: 2 Score:11215
Group 2: Delayed
Sequence:1 Score:2745
Alignment Score 3187
GCG-Alignment file created [/tmp/5M8igBfS6R/JI0anpGzWn]
cannot remove directory for /tmp/jFAAdGicRB: Directory not empty at
/usr/lib/perl5/site_perl/5.10/Bi
o/Tools/Run/WrapperBase.pm line 243
SEQ1 SEQ2 Ka Ks Ka/Ks PROT_PERCENTID CDNA_PERCENTID
--------------------------------------------------------
Do you have an idea what happened here and how I can solve this problem? I
appreciate any help with this.
I am running Cygwin under Windows XP, and have installed clustalw, t-coffee,
the bioperl package and the bioperl package run with the "Build" scripts.
Best regards,
Markus.
Dr. Markus Liebscher
Martin-Luther-University Halle/Wittenberg
Dept. Natural product chemistry
Kurt-Mothes-Str.3
06120 Halle (Saale)
Germany
More information about the Bioperl-l
mailing list