view scripts/get_sequence_based_on_id.pl @ 30:57b2221536fc

remove custom samtools package to its own package on TTS

use Bio::DB::Fasta;

my $fastaFile = shift;
my $queryFile = shift;

my $db = Bio::DB::Fasta->new( $fastaFile );
    $seq = $queryFile;
    chomp($queryFile);

    my $sequence = $db->seq($seq);
    if  (!defined( $sequence )) {
            die "Sequence $seq not found. \n" 
    }   
    print ">$seq\n", "$sequence\n";