Hola.
Necesito hacer un script que extraiga regiones al azar de un archivo genbank. La región debe tener 500 nt de longitud.
Mi código es el siguiente:
use Bio::SeqIO;
sub extract_intergenic_from_genbank {
my ( $infile, $out_intergenic_file, $min_intergenic_size, $max_intergenic_size ) = @_;
my ( $n_of_intergenic, $gi, $start, $end, $length, $strand, $taxon ) = (0);
my $in = new Bio::SeqIO( -file => "NC_002755.gb", -format => 'genbank' );
open( FNA, ">$salida.FASTA" )
|| die "# extract_intergenic_from_genbank : ...
News of Foro - Perl en Español