Discussion :

[linalin]

Re bonjour à tout le monde,
j'ai un fichier multifasta et je veux récupérer le tout pour l'insérer dans ma base de données et notamment la table transcript!
Ce fichier contient l'identifiant et la description du transcrit suivi de sa séquence! (il ya à peu près 35000 transcrits dans ce fichier).
Donc je voudrais un peu d'aide pour écrire un script pour pouvoir récupérer l'identifiant dans le champs identifiant, la description dans le champs description, et la séquence dans le champs séquence! Pour vs donner une petite idée du fichier fasta, voilà un exemple de son contenu pour un trancrit :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
>gi|124249439|ref|NM_001080875.1| Gallus gallus succinate dehydrogenase Ip (SDHB), mRNA
AATTCCACTAGTTACGCGACCCCCGGAACCCAAAGACTTGGGTTTCCCGGGAGCTGCCCCGGTACGGCTC
ATGAGTGGAGGACGTCATCAGCGGTGATGATTACTCGAATGGGGGATTCTATGGGGATTACAATGCGATG
GTCGACCAAAATGGCGGCGGCCGTGGTGGGAGTCTCCTTGAGGCGCGGCGTCCCCGCGCGGTTCCTGAGG
GCTGGGCTGCGGCCGGTGAGGGGACTGGAGGCGGTCCACGGGATCTGTCGAGGAGCTCAGACGGCGGCTG
CAGCTACATCACGTATCAAGAAGTTTTCCATCTACAGATGGGATCCTGATAAGCCTGGGGACAAGCCCCG
CATGCAGACCTATGAAGTGGATTTAAATAAATGTGGGCCTATGGTACTTGATGCTCTGATTAAGATTAAA
AACGAGTTGGACTCCACTCTGACCTTCCGCAGGTCATGCAGGGAAGGCATCTGTGGCTCCTGTGCAATGA
ACATTGCTGGTGGAAACACCCTGGCTTGTACTAAAAAAATTGACCCTGATCTCAGCAAGACCACTAAAAT
CTACCCTCTCCCCCACATGTATGTGGTGAAGGATCTCGTTCCAGACTTGAGTAACTTCTACGCACAGTAC
AAATCCATCGAGCCTTACCTGAAGAAGAAGGACGAGTCGAAACAGGGCAAGGAGCAGTACTTGCAATCCA
TAGAAGACCGTCAGAAACTGGACGGACTCTATGAGTGCATCCTCTGTGCCTGCTGCAGCACCAGCTGTCC
CAGTTACTGGTGGAATGGGGACAAGTACTTGGGTCCTGCTGTACTCATGCAGGCCTATCGCTGGATGATC
GACTCCAGAGATGACTACACAGAGGAACGCCTGGCACAACTTCAAGACCCATTTTCTCTCTACCGTTGTC
ACACTATCATGAACTGCACAAGGACTTGCCCAAAGGGTTTGAACCCTGGCAAAGCAATTGCTGAGATCAA
GAAGATGATGGCAACTTACAAAGAGAAGGCAGCCGCTGCATAATGCTGTTCCCCGACTGGAAATGTAACA
CACAAACATGCTTTACCTGGAAGTAACCTGCACTTGATTACACCTGGTTCCAGCAGGAATGTTGGTGTTT
TTCCCATATGCTTGTGTGTACAATAAATGCTGTAAAGAACAAAAAAAA

Merci

[Leishmaniose]

quel est le SGBD que tu utilises et surtout via quel interface (e.g. MySQL via phpmyadmin) ?

[linalin]

Je suis sur linux ubuntu et j'utilise mysql avec le terminal donc voilà tout! Tu penses que tu pourrais m'aider??

[Leishmaniose]

Bon je ne comprends réellement pas où est le problème surtout que tu utilises MySQL !!! le truc est très simple... Avec MySQL il est possible de charger les données dans une base à partir d'un fichier sachant que les données en question doivent être agencées d'une manière que le SGBD arrive à identifier chacun des enregitrements champs par champs. Pour ceci il existe ce qu'on appel les séparateurs par exemple si un fichier (data.txt) contient les données suivants :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
1;data1;data11;data12
2;data2;data21;data22
3;data3;data31;data32
....

il est possible de l'utiliser pour charger une table constituée de 4 champs et ce en utilisant le SQL suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

LOAD DATA LOCAL INFILE '/home/user/data_files/data.txt' INTO TABLE `ma_table` FIELDS TERMINATED BY ';' ESCAPED BY '\\' LINES TERMINATED BY '\r\n'

qui va permettre d'obtenir la table suivante :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
champs1 | champs2 | champs3 | champs4
1       | data1   | data11  | data12
2       | data2   | data21  | data22
3       | data3   | data31  | data32

Partant de ce principe, il est possible et même judicieux d'utiliser Perl pour mettre en forme ton fichier source.

Bon avant d'aller plus loin je veux m'assurer que ça correspond bien à ce que tu cherches à faire !!

@+

[Jasmine80]

Coucou,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
#!/usr/local/bin/perl
 
use strict;
use warnings;
use FileHandle;
 
 
 
# Fichiers d'entrée et de sortie
#----------------------------------
 
my $File_In = "Fasta";
my $InFile = "P:/Perl/scripts/Files/$File_In.txt";
my $File_Out = "FichierPourDB";
my $OutFile = FileHandle->new (">P:/Perl/scripts/Files/$File_Out.txt");
 
 
 
# Ouverture du fichier fasta
#----------------------------------
open (FileFas,"$InFile")  or die "Can't open file\n";
 
 
 
 
# Lecture du fichier fasta
#----------------------------------
 
my $on =0;   # on ne passera pas dans le if lors du 1er passage
my $LineFas;
my $SequenceTotale;
my $gi;
my $Description;
 
while ($LineFas=<FileFas>)
{
 
        if($LineFas =~ /^>gi\|(\d+)\|ref\|[\w.]*\|\s*([\w\s\)\(]*)/)
        {
                 if ($on == 1) {print $OutFile ($gi."\t".$Description."\t".$SequenceTotale."\n")};
                 ($gi) =  $1;
                 ($Description) =  $2;
                 $on = 1;
        }
 
 
        if($LineFas =~ /^(\w*)\n/)
        {
                 my ($Sequence) =  $1;
                 $SequenceTotale = $SequenceTotale.$Sequence;
        }
 
 
}
 
# pour la dernière séquence
print $OutFile ($gi."\t".$Description."\t".$SequenceTotale."\n");
 
 
close;

La fonction LOAD DATA LOCAL INFILE prend par défaut le signe de tabulation comme séparateurs. Il te suffit donc de charger le fichier comme Leishmaniose te l'a indiqué sans aucun paramètre à la fonction.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
CREATE TABLE ma_table (Gi VARCHAR(20), Description TEXT, Sequence LONGTEXT, PRIMARY KEY (Gi));
LOAD DATA LOCAL INFILE 'P:/Perl/scripts/Files/FichierPourDB.txt' INTO TABLE `ma_table`;

Bon travail,

Jasmine,



PS1: Je travaille sous windows, je ne sais pas si cela fonctionnera sous linux pour ce qui est de l'appel des fichiers, peut être faut-il utiliser l'antislash au lieu du slash

PS2: Sais-tu que tu peux aussi te connecter directement à Genbank, y faire ta requête via BioPerl puis charger les informations renvoyées dans une base de données sans passer par aucun fichier texte?

[bigplayer]

bonjour je voudrais tester un programme en perl.quel est l'editeur de perl et comment tester ce programme
merci d'avance.