#!/usr/local/bin/perl -w
#tRANSCAN-SE_extract
die "usage: $0 tRNASCAN-SE.list\n" unless @ARGV ==1;
while (<>) {
	my ($seq, $tRNA_no, $p1, $p2, $amino, $anticodon, $stat1, $stat2, $percent) = split;
	$beg = $p1 if $p2 > $p1; 
	$beg = $p2 if $p2 < $p1;
	$end = $p1 if $p1 > $p2; 
	$end = $p2 if $p1 < $p2;
        system ("extractseq -sequence=$seq.seq -auto -stdout -separate -reg=$beg..$end");
	}