bonjour
voici mon problème :
je récupère des identifiants et les séquences associées dans une base de données.
J'ai besoin de formater ces informations sous la forme d'un fichier fasta.
En utilisant Bio::SeqIO de cette manière là :
le problème c'est que j'ai ceci et non pas d'objet séquence :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6 my $in = Bio::SeqIO-> new ( -file =>">$file", -format =>'fasta' ); $in-> write_seq($seq);
ma question serait comment créer un objet séquence permettant l'utilisation des fonctions de Bio::SeqIO.
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7 foreach my $sequence (@tot_seq) { my $seq_id = $$sequence[1]; my $seq = $$sequence[2]; }
Ou bien comment à partir directement de l'identifiant et de la séquence créer un fichier multifasta.
Peut-être est-il possible également de faire une conversion à partir de ce format de fichier qui est facile d'obtenir ? ( j'ignore le nom de ce format, format brute peut-être ?)
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2>Glyma2227s00200.1 EWIYILSVSGHFVLTNLLVKKMVETAKETGVQGRIVNVSSSIHGWFSGDAISYLALISRNKRHYDATRAYALSKLANVFHTKELSRRLQQMGANVTVNCVHPGIVRTRLTREREGLLTDLVFFLASKLLKTIPQVIKYNMK*
merci par avance, de vos remarques et suggestions
Partager