举个例子:从数据库中下载基因组数据,基因组.gff文件中染色体ID较为复杂(第一列)
从全基因组序列.genome.fa文件中找到染色体ID对应简写
想将gff文档中染色体ID全部替换成第二列LC*命名
use Data::Dumper;
use Getopt::Long;
use strict;...
物种gff文件的信息格式如下:ID号是第二列的 si*G******
perl脚本信息如下:
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})){
&...
Library preparation and Illumina sequencingGenomic DNA from each sample was prepared for libraries construction. The library preparations were sequenced by using an Illumina HiSeqX platform and 150bppaired-end reads were generated.
Data filtering and alignment
The recently released genome...