...用如下的代码实现。
#!/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...
...ER2 (for
details, see section 1.4.2 about repeat masking).
It is important to note that even though MAKER2 and BRAKER2 used the same repeat libraries, genomes are
masked differently by MAKER2 when comparing to our approach of using RepeatMasker/RepeatModeler for
BRAKER2.
MAKER2 runs RepeatMasker int...