[Bioperl-l] [Bioperl-guts-l] Notification: incoming/1077 (fwd)

Jason Stajich jason@cgt.mc.duke.edu
Tue, 5 Feb 2002 18:27:53 -0500 (EST)


Lana -

Any reason you can't use fasta format as input to clustal?

We have basically stopped supporting PIR format because none of our
developers use it and after sending out request of interest on the list
severl months ago no one responded that were using it. If there are indeed
users out there we would be happy to add the appropriate code if we can
get good use cases and example files that do not work with the current
code.

In any event I would certainly upgrade to bioperl 0.7.2 in the near future
and the future 1.0 release as I suspect the pir handling is not so great
in that 0.7.0 release.

-jason

-- 
Jason Stajich
Duke University
jason@cgt.mc.duke.edu

---------- Forwarded message ----------
Date: Tue, 5 Feb 2002 18:02:18 -0500
From: bioperl-bugs@bioperl.org
To: bioperl-guts-l@bioperl.org
Subject: [Bioperl-guts-l] Notification: incoming/1077

JitterBug notification

new message incoming/1077

Message summary for PR#1077
	From: Lana Schaffer <lschaffe@agouron.com>
	Subject: file formats
	Date: Tue, 05 Feb 2002 15:08:37 -0800
	0 replies 	0 followups

====> ORIGINAL MESSAGE FOLLOWS <====

>From lschaffe@agouron.com Tue Feb  5 18:02:17 2002
Received: from tbone.agouron.com (tbone.agouron.com [198.182.177.3])
	by pw600a.bioperl.org (8.12.2/8.12.2) with SMTP id g15N2FPX023149
	for <bioperl-bugs@bioperl.org>; Tue, 5 Feb 2002 18:02:16 -0500
Received: from relay.agouron.com by tbone.agouron.com
          via smtpd (for pw600a.bioperl.org [199.93.107.70]) with SMTP; 5 Feb 2002 23:08:39 UT
Received: from agouron.com (rudy.agouron.com [10.0.76.106])
	by hermes.agouron.com (8.10.2+Sun/8.9.3) with ESMTP id g15N8bQ10069
	for <bioperl-bugs@bioperl.org>; Tue, 5 Feb 2002 15:08:37 -0800 (PST)
Sender: lschaffe@agouron.com
Message-ID: <3C6065F5.93DD6611@agouron.com>
Date: Tue, 05 Feb 2002 15:08:37 -0800
From: Lana Schaffer <lschaffe@agouron.com>
Organization: Agouron
X-Mailer: Mozilla 4.77C-SGI [en] (X11; I; IRIX64 6.5 IP30)
X-Accept-Language: en
MIME-Version: 1.0
To: bioperl-bugs@bioperl.org
Subject: file formats
Content-Type: text/plain; charset=us-ascii
Content-Transfer-Encoding: 7bit

 I am starting to use Bioperl and presently have Bioperl 0.7
and unix OS IRIX64 6.5.
I have used a small script to change the protein sequence format
from fasta to PIR.  The output perl file from Bioperl is:
-------------------------------------------------------------------
>P1;CDK2

>P1;CDK2

MENFQKVEKI GEGTYGVVYK ARNKLTGEVV ALKKIRLDTE TEGVPSTAIR EISLLKELNH
PNIVKLLDVI HTENKLYLVF EFLHQDLKKF MDASALTGIP LPLIKSYLFQ LLQGLAFCHS
HRVLHRDLKP QNLLINTEGA IKLADFGLAR AFGVPVRTYT HEVVTLWYRA PEILLGCKYY
STAVDIWSLG CIFAEMVTRR ALFPGDSEID QLFRIFRTLG TPDEVVWPGV TSMPDYKPSF
PKWARQDFSK VVPPLDEDGR SLLSQMLHYD PNKRISAKAA LAHPFFQDVT KPVPHLRL
-------------------------------------------------------------------

and the output/input perl file for clustalw is:

-------------------------------------------------------------------
>P1;CDK2

MENFQKVEKIGEGTYGVVYKARNKLTGEVVALKKIRLDTETEGVPSTAIREISLLKELNH
PNIVKLLDVIHTENKLYLVFEFLHQDLKKFMDASALTGIPLPLIKSYLFQLLQGLAFCHS
HRVLHRDLKPQNLLINTEGAIKLADFGLARAFGVPVRTYTHEVVTLWYRAPEILLGCKYY
STAVDIWSLGCIFAEMVTRRALFPGDSEIDQLFRIFRTLGTPDEVVWPGVTSMPDYKPSF
PKWARQDFSKVVPPLDEDGRSLLSQMLHYDPNKRISAKAALAHPFFQDVTKPVPHLRL
*
-------------------------------------------------------------------

The clustalw program will produce an error when an import is attempted
with the Bioperl format.  The "*" is missing and the ">P1;CDK2" is
repeated twice.
I don't know if this format is corrected in more updated releases of
Bioperl.  However, this format needs to be compatible with other
programs.  Please update me on this problem.

Thanks,


--
Lana Schaffer                   Computational Chemistry
schaffer@pfizer.com            Pfizer Global R&D-LJ
phone  (858) 622-3002           10777 Science Center Drive
fax    (858) 678-8244           San Diego, CA 92122-1111




_______________________________________________
Bioperl-guts-l mailing list
Bioperl-guts-l@bioperl.org
http://bioperl.org/mailman/listinfo/bioperl-guts-l