...s (http://www.uniprot.org), using BLASTx (Altschul et al., 1997) with a cutoff E-value 1E-5. Blast results were imported into Blast2Go (Conesa et al., 2005) for GO term mapping. Mapping results were submitted to WEGO (Ye et al., 2006) for GO classification.
他的结果如下：
814 and 745 CDS...
Using standard /share/work/biosoft/perl/perl-5.22.1/lib/site_perl/5.22.1/auto/share/dist/AGAT/config.yaml file
File .longest_isoform.gff3 already exist.
...neRegions (e.g., exonic, intronic, non-coding RNA)) that one variant hits; please click here for details.Gene.refGeneGene name associated with one variantExonicFunc.refGeneExonic variant function, e.g., nonsynonymous, synonymous, frameshift insertion.please click here for details.AAChange.refGen...
...isect = intersect(@a, @b); # unique union
my @unique = unique(@a, @b); # check if arrays contain same members
if ( !array_diff(@a, @b) ) {
# do something
} # get items from array @a that are not in array @b
my @minus = array_minus( @a, @b );包的详细帮助：https://metacpan.or...
...lt;-read.table("All_DEG_venn.genes",header=TRUE,comment.char="",sep = "\t",check.names=FALSE,fill=T,na.strings="")colnames(All_DEG_venn)<-c("DAF2 vs GDAF2","DAF5 vs GDAF5","DAF11 vs GDAF11","DAF16 vs GDAF16")
#去除数据中的NA，以免多计数
All_DEP_ID=apply(All_DEG_venn,2,na.omit)datav<...