[Bioperl-l] help about Fasta file???
issam outaleb
outaleb at web.de
Mon Oct 1 07:46:08 EDT 2007
hallo all,
i have a little problem,:
im using this programm,
i got some experiment and get some results--> IPI hits,(IPI Accnum)
what i want is how can i correlate this IPI ACC Numbers with the FASTA FILE (database fasta),so the programm
has to look where is the IPI Accnum in the db and copy this include description and Sequence to a new file;
all
#!/usr/bin/perl#use warning;#use strict;use CGI qw(:all);open (IN,"C:/Documents and Settings/XXX/Desktop/Search_file") or die "Fehler beim oeffnen";open (FASTA_db,"C:/Documents and Settings/XXX/Desktop/FASTA1.fasta") or die "FASTA nicht möglich zum öffnen!!" ;open (OUT,">C:/Documents and Settings/XXX/Desktop/reslut.txt") or die "Fehler beim Anlegen der neuen Datei";#print "\nDateien zum kopieren geoeffnet\n";while (<IN>){ $i = $_; chomp $i; if(/Hit\d">([^<\/A> ]*)/)#match string from htm datei,give me such result-->IOP123234(just IPIs) { #print OUT $1."\n"; #print this IPIs in this file.
#what i thought was to push up this IPIs in the array than look at them in the fasta_db file and copy it to new file with the description
#and Sequence also. so generate a new fasta file include just my IPIs results. how???
$j = $1; push(@array,$j); } } while (defined($var=<suchdatei>)){
$var =~ /(>IPI:)([^| .]*)([^>]*)/ ;# }} close (IN);close (FASTA_db);close (OUT);print "\nDateien geschlossen, Kopiervorgang
.
Jetzt neu! Schützen Sie Ihren PC mit McAfee und WEB.DE. 3 Monate
kostenlos testen. *http://www.pc-sicherheit.web.de/startseite/?mc=022220* [http://www.pc-sicherheit.web.de/startseite/?mc=022220]
More information about the Bioperl-l
mailing list