Book Image

Bioinformatics with Python Cookbook

By : Tiago R Antao, Tiago Antao
Book Image

Bioinformatics with Python Cookbook

By: Tiago R Antao, Tiago Antao

Overview of this book

Table of Contents (16 chapters)
Bioinformatics with Python Cookbook
Credits
About the Author
About the Reviewers
www.PacktPub.com
Preface
Index

Traversing genome annotations


Having a genome sequence is interesting, but we will want to extract features from it: genes, exons, and coding sequences. This type of annotation information is made available in GFF and GTF files. GFF stands for Generic Feature Format. In this recipe, we will see how to parse and analyze GFF files, using the annotation of the Anopheles gambiae genome as an example.

Getting ready

We will use the gffutils library to process the annotation file.

If you do not use the notebook, you need to acquire the annotation file from our datasets page at https://github.com/tiagoantao/bioinf-python/blob/master/notebooks/Datasets.ipynb (file gambiae.gff3.gz) Rename the annotation file as gambiae.gff.gz. Preferably, use the 02_Genomes/Annotations.ipynb notebook, which is provided in the code bundle of the book.

How to do it...

Let's take a look at the following steps:

  1. Let's start by creating an annotation database with gffutils based on our GFF file:

    import gffutils
    import sqlite3...