Précédent   Forum du club des développeurs et IT Pro > Autres langages > Perl > Bioinformatique
Bioinformatique Toutes vos questions sur les scripts Perl associés à la bioinformatique, modules bioperl, projets bioinformatiques, etc ... Avant de poster, veuillez consulter les cours Perl et les critiques de livres.
Partagez cette discussion sur d'autres réseaux sociaux : Viadeo Twitter Google Facebook Digg Delicious MySpace Yahoo
Réponse
 
Outils de la discussion
Publicité
'
Vieux 23/02/2012, 18h40   #1
hajarita
Candidat au titre de Membre du Club
 
Femme
Étudiant
Inscription : janvier 2012
Messages : 49
Détails du profil
Informations personnelles :
Sexe : Femme
Localisation : France, Haute Garonne (Midi Pyrénées)

Informations professionnelles :
Activité : Étudiant
Secteur : Services à domicile

Informations forums :
Inscription : janvier 2012
Messages : 49
Points : 12
Points : 12
Par défaut Parser un fichier swissProt

Bonjour,
Je suis débutante en perl et je veux recupérer des informations dans le fichier swissProt:
#############################################"
Citation:
ID PYRB_LACLL
AC 3RQV1731;
DE RecName: Full=(trusted)Aspartate carbamoyltransferase;
DE EC=2.1.3.2;
DE AltName: Full=Aspartate transcarbamylase;
DE Short=ATCase;
GN Name=pyrB;
OS Lactococcus lactis subsp. lactis.
OC Bacteria; Firmicutes; Lactobacillales; Streptococcaceae; Lactococcus.
OX NCBI_TaxID=1360;
CC -!- CATALYTIC ACTIVITY: Carbamoyl phosphate + L-aspartate = phosphate
CC + N-carbamoyl-L-aspartate.
CC -!- PATHWAY: Pyrimidine metabolism; UMP biosynthesis via de novo
CC pathway; (S)-dihydroorotate from bicarbonate: step 2/3.
CC -!- SIMILARITY: Belongs to the ATCase/OTCase family.
DR GO; GO:0004070; F:aspartate carbamoyltransferase activity; IEA:HAMAP.
DR GO; GO:0006221; P:pyrimidine nucleotide biosynthetic process; IEA:HAMAP.
DR HAMAP; MF_00001; Asp_carb_tr; 1; -.
DR PROSITE; PS00097; CARBAMOYLTRANSFERASE; 1.
KW Pyrimidine biosynthesis; Transferase.
**
** ################# INTERNAL SECTION ##################
**HA SAM; Annotated by SAM 1.216; MF_00001.29; 14-FEB-2012.
**HA FAM; Method MF_00001; PYRB; Trusted match; 37.146 (+3.9).
**HA FAM; Method MF_01109; OTC; Weak match; 18.126 (-50.3).
**HR DE 1742307__1743239 (reverse), 311 amino acids.
SQ SEQUENCE 310 AA; 34633 MW; 53DDA4B1C6F32776 CRC64;
MSVKNGLVQL ENLTSMEKLS VDEVMGLIKR ASAFKAGTAE FDLEKQTFAS NLFFENSTRT
HHSFHIAERK LGLDVLEFDA QASSISKGET LYDTVLTLDA LGVDICVIRS GVEHYYEELV
NSDNIHCAIV NGGDGSGQHP SQCLLDLMTI YEEFGKFEGL KIAISGDLTH SRVAKSNMMM
LQKLGARLYF TGPAAWYSEE FDDYGHYANL DRILPELDVH MLLRVQHERH DSGESFSKEG
YHNHFGLTEE RAKMLKPTAI IMHPAPVNRD VEIADSLVEN SQSRIVQQMS NGVYTRMAIL
EAILAGKKAK
/###################################################"

J'ai trouver le package swissknif qui permt de le faire.
J'ai pu récupérer des champs comme; AC, ID, CC. MAIS, j'ai pas pu recuperer les champs DR,FT et **HA, **HR
Quelq'un peut m'aider s'il vous plait
hajarita est déconnecté   Envoyer un message privé Réponse avec citation 00
Vieux 24/02/2012, 13h39   #2
Beniou
Membre expérimenté
 
Avatar de Beniou
 
Homme
Inscription : novembre 2009
Messages : 357
Détails du profil
Informations personnelles :
Sexe : Homme
Âge : 33
Localisation : France, Nord (Nord Pas de Calais)

Informations forums :
Inscription : novembre 2009
Messages : 357
Points : 509
Points : 509
Bonjour,

As-tu regardé du coté de Bioperl et du package Bio::Seq.

Tous les formats les plus courants (fasta, blast, swissprot ...) sont pris en charge.
Beniou est déconnecté   Envoyer un message privé Réponse avec citation 00
Vieux 09/03/2012, 14h49   #3
alaninho
Membre régulier
 
Inscription : mars 2011
Messages : 175
Détails du profil
Informations forums :
Inscription : mars 2011
Messages : 175
Points : 74
Points : 74
bonjour,

Regarde bioPerl, c'est le meilleur moyen, j'ai déjà travailler avec et c'est efficace.

Alaninho
alaninho est déconnecté   Envoyer un message privé Réponse avec citation 00
Réponse
Outils de la discussion

Navigation rapide


Fuseau horaire GMT +2. Il est actuellement 08h21.


 
 
 
 
Partenaires

Hébergement Web