Mercurial > repos > biomonika > linkyx

annotate scripts/get_sequence_based_on_id.pl @ 9:695d28139f3e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
toolshed8
author biomonika <biomonika@psu.edu>
date Tue, 09 Sep 2014 14:31:02 -0400
parents 1955f03f092e
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
1 use Bio::DB::Fasta;
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
2
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
3 my $fastaFile = shift;
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
4 my $queryFile = shift;
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
5
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
6 my $db = Bio::DB::Fasta->new( $fastaFile );
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
7 $seq = $queryFile;
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
8 chomp($queryFile);
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
9
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
10 my $sequence = $db->seq($seq);
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
11 if (!defined( $sequence )) {
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
12 die "Sequence $seq not found. \n"
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
13 }
1955f03f092e initial commit
biomonika <biomonika@psu.edu>
parents:
diff changeset
14 print ">$seq\n", "$sequence\n";