...用如下的代码实现。
#!/usr/bin/perl -w
use strict;
my $biotype_file = shift @ARGV;
my $gtf = shift @ARGV;
my $biotype = shift @ARGV;
my %biotype_list;
open my $fh1, $biotype_file or die;
while (<$fh1>) {
chomp;
my @array = split /\t/, $_;
if($array[2]eq $biotype){
...
Fungal mitochondrions also contain non-coding RNA genes. The mitochondrial rRNA genes are single copy genes for the large and small ribosomal RNAs (rrl and rrs). These are significantly shorter than their counterpart in the nuclear genome and cannot be identified by RNAmmer(Lagesen, Hallin, Rodlan...