Skip to content

agat_sp_kraken_assess_lift_coverage.pl#

DESCRIPTION#

The script takes as input gtf produced by Kraken (lift-over tool). It will analyse the kraken_mapped attributes to calculate the mapped percentage of each mRNA. According to a threshold (0 by default), gene with a mapping percentage over that value will be reported. A plot nammed geneMapped_plot.pdf is performed to visualize the result. /! The script handles chimeric files (i.e containg gene part mapped on the template genome and others on the de-novo one) /!/! If the file is complete (containing kraken_mapped="TRUE" and kraken_mapped="FALSE" attributes), the script calcul the real percentage lentgh that has been mapped. Else the calcul is only based on feature with kraken_mapped="TRUE" attributes. So in this case the result most of time will be 100%. /!/!/! We met rare cases where Kraken mapped a feature to several locations of the de-novo genome. As result we could end up with mapping over > 100%. We report them as 100% mapped in the plot and a warning is raised to allow to check thoses cases.

SYNOPSIS#

agat_sp_kraken_assess_lift_coverage --gtf infile.gtf [ -o outfile ]
agat_sp_kraken_assess_lift_coverage --help

OPTIONS#

  • -gtf

    Input gtf file produced by Kraken.

  • --threshold or -t

    Gene mapping percentage over which a gene must be reported. By default the value is 0.

  • --verbose or -v

    Verbose information.

  • -o , --output , --out or --outfile

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

  • -c or --config

    String - Input agat config file. By default AGAT takes as input agat_config.yaml file from the working directory if any, otherwise it takes the orignal agat_config.yaml shipped with AGAT. To get the agat_config.yaml locally type: "agat config --expose". The --config option gives you the possibility to use your own AGAT config file (located elsewhere or named differently).

  • -h or --help

    Display this helpful text.