The script aims to compare two annotation of the same assembly. It provided information about split/fusion of genes between the two annotations. The most common case are: 1 => 0 ( gene uniq to file1) 0 => 1 ( gene uniq to file2) 1 => 1 ( 1 gene from file 1 overlaps only 1 gene from file2) 1 => <many> ( 1 gene from file 1 overlaps <many> genes from file2) => split case (with file 1 as reference) <many> => 1 ( <many> genes from file 1 overlap only 1 gene from file2) => fusion case (with file 1 as reference)
Then you can get more complex cases: <many> => <many> (<many> genes from file 1 overlap <many> genes from file2)
agat_sp_compare_two_annotations.pl -gff1 infile.gff [ -o outfile ] agat_sp_compare_two_annotations.pl --help
Input GTF/GFF file1.
Input GTF/GFF file2.
-o , --output , --out or --outfile
Output GFF file. If no output file is specified, the output will be written to STDOUT.
Verbose option, make it easier to follow what is going on for debugging purpose.
-h or --help
Display this helpful text.