comparison FunAnn/annotate.pl @ 1:9b8afa0b6da4 draft

Uploaded

0:1886293c2243

{
   system("wget -i https://raw.githubusercontent.com/matteo14c/CorGAT_galaxy/dev/ann.txt");
   system("gzip -d GCF_009858895.2_ASM985889v3_genomic.fna.gz");
}




$gen_code="genetic_code";
die("need genetic code file in the current folder\n") unless -e "genetic_code";
open(IN,$gen_code);
while(<IN>)
{
	($triplet,$oneL)=(split());
	$code{$triplet}=$oneL;
}

$genome="GCF_009858895.2_ASM985889v3_genomic.fna";
die("need reference genome file in the current folder\n") unless -e "GCF_009858895.2_ASM985889v3_genomic.fna";
open(IN,$genome);
while(<IN>)
{
	next if $_=~/^>/;
	chomp;

1:9b8afa0b6da4

{
   system("wget -i https://raw.githubusercontent.com/matteo14c/CorGAT_galaxy/dev/ann.txt");
   system("gzip -d GCF_009858895.2_ASM985889v3_genomic.fna.gz");
}

$gen_code="genetic_code";
die("need genetic code file in the current folder\n") unless -e "genetic_code";
open(IN,$gen_code);
while(<IN>)
{
	($triplet,$oneL)=(split());
	$code{$triplet}=$oneL;
}

$genome="GCA_009858895.3_ASM985889v3_genomic.fna";
die("need reference genome file in the current folder\n") unless -e "GCA_009858895.3_ASM985889v3_genomic.fna";
open(IN,$genome);
while(<IN>)
{
	next if $_=~/^>/;
	chomp;