Skip to main content

A Python3 annotation program to select the best gene model in each locus

Project description

Mikado is a lightweight Python3 pipeline whose purpose is to facilitate the identification of expressed loci from RNA-Seq data * and to select the best models in each locus.

The logic of the pipeline is as follows:

  1. In a first step, the annotation (provided in GTF/GFF3 format) is parsed to locate superloci of overlapping features on the same strand.

  2. The superloci are divided into different subloci, each of which is defined as follows:

    • For multiexonic transcripts, to belong to the same sublocus they must share at least a splicing junction (i.e. an intron)
    • For monoexonic transcripts, they must overlap for at least one base pair
    • All subloci must contain either only multiexonic or only monoexonic transcripts
  3. In each sublocus, the pipeline selects the best transcript according to a user-defined prioritization scheme.

  4. The resulting monosubloci are merged together, if applicable, into monosubloci_holders

  5. The best non-overlapping transcripts are selected, in order to define the loci contained inside the superlocus.

    • At this stage, monoexonic and multiexonic transcript are checked for overlaps
    • Moreover, two multiexonic transcripts are considered to belong to the same locus if they share a splice site (not junction)
  6. Once the loci have been defined, the program backtracks and looks for transcripts which can be assigned unambiguously to a single locus and constitute valid alternative splicing isoforms of the main transcripts.

The criteria used to select the “best” transcript are left to the user’s discretion, using specific configuration files.

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Filename, size & hash SHA256 hash help File type Python version Upload date
Mikado-1.2.4-cp36-cp36m-macosx_10_7_x86_64.whl (16.7 MB) Copy SHA256 hash SHA256 Wheel 3.6 Aug 8, 2018
Mikado-1.2.4.tar.gz (19.9 MB) Copy SHA256 hash SHA256 Source None Aug 8, 2018

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN SignalFx SignalFx Supporter DigiCert DigiCert EV certificate StatusPage StatusPage Status page