agat_sp_filter_gene_by_length.pl

DESCRIPTION

The script aims to filter level1 feature (e.g. gene, match, etc) by length. It will create two files. one with the feature passing the length filter, the other one with the remaining features. If the level1 feature has exon features, the size is computed by concatenating the exon together. If the level1 feature has several level2 features (e.g. mRNA) we apply the test on the longest one (the longest concatenated exon set).

Some examples: Select L1 feature shorter than 1000bp: agat_sp_filter_gene_by_length.pl --gff infile.gff --size 1000 --test "<" -o result.gff Select genes longer than 200bp: agat_sp_filter_gene_by_length.pl --gff infile.gff --size 200 --test ">" -o result.gff

SYNOPSIS

agat_sp_filter_gene_by_length.pl --gff infile.gff --test ">=" --nb 10 [ --output outfile ]
agat_sp_filter_gene_by_length.pl --help

OPTIONS

  • -f, --reffile, --gff or -ref

    Input GFF3 file that will be read

  • -s or --size

    Integer - Gene size in pb [Default 100]

  • -t or --test

    Test to apply (>, <, =, >= or <=). If you use one of these two characters >, <, please do not forget to quote your parameter like that "<=". Else your terminal will complain. [Default "="]

  • -o or --output

    Output GFF file. If no output file is specified, the output will be written to STDOUT.

  • -v

    Verbose option for debugging purpose.

  • -h or --help

    Display this helpful text.