使用gsds绘制基因结构图_从gff文件当中提取对应转录本ID的基因结构信息用于GSDS绘制结构图脚本更新...

脚本名称：get_gene_exon_from_gff.pl

用法：perl ../script/get_gene_exon_from_gff.pl -in1 WRKY_domain_new_out_emoved_redundant.txt -in2 Arabidopsis_thaliana.TAIR10.31.gff3 -out gene_exon_info.gff

输出结果(部分)：

脚本源代码：

use Getopt::Long;

my %opts;

use Data::Dumper;

GetOptions( \%opts, "in1=s", "in2=s", "out=s", "h" );

if ( !defined( $opts{in1} )

|| !defined( $opts{in2} )

|| !defined( $opts{out} )

|| defined( $opts{h} ) )

{

&USAGE;

}

open( IN1, "$opts{in1}" ) || die "open $opts{in1} failed\n";

open( IN2, "$opts{in2}" ) || die "open $opts{in2} failed\n";

open( OUT, ">$opts{out}" ) || die "open $opts{out} failed\n";

my %gffs;