[Bioperl-l] FASTA parsing

simona bazzocchi simona_bazzocchi at yahoo.it
Tue Sep 19 22:54:53 UTC 2006


Hi, I have a problem.
I need to filter just the query sequences that match
with the first 4 bases vs library sequences.
I used FASTA to get the output (FASTA_RESULTS104.txt)
using my perl script below, but now I don't know how i
can filter that output.
I read the Bio::SearchIO manual but it's not so clear
to me.
Have you any idea?
Thank you very much
Simona 

#!/usr/bin/perl -w
use strict;
my $location="/home/kei/fasta/fasta34";
my $library="DNAUS.txt";
my $query="mature2.txt";
my $ktup=1;
my $options="-r +5/-4 -H -w 100 -m 9 -f -12 -g -4 -q";
my $command="$location $options $query $library
$ktup";
my $output = `$command`;
open FASTA,">>","FASTA_RESULTS104.txt";
print FASTA $output;
close (FASTA);


__________________________________________________
Do You Yahoo!?
Poco spazio e tanto spam? Yahoo! Mail ti protegge dallo spam e ti da tanto spazio gratuito per i tuoi file e i messaggi 
http://mail.yahoo.it 



More information about the Bioperl-l mailing list