Home - this site is powered by TWiki(R)
TWiki webs: Main | TWiki | Sandbox   Log In or Register

Changes | Index | Search | Go
-- %TEACHINGWEB%.MarkDewitt - 10 Dec 2009

GFF processing from GenBank files

Other members of team OMGBrowse are writing scripts that can produce genome browsers directly from genbank accession numbers, and to manage them through a more user-friendly front end than JBrowse itself. My contributions to this project mostly concerned learning how to use JBrowse and to help write perl scripts that process GenBank files into GFF files that are easily readable by the JBrowse installation scripts.

The perlscripts attached below help with this automation process.

How to use it:

0. Install JBrowse, if you haven't already.

1. Get a GenBank file from GenBank, multiple files as one big file are acceptable. For example, I looked up "dengue" on NCBI's genome database, and selected GenBank output to file at the top of the window. This creates a file with four GenBank entries, one for each serotype.

2. Copy the perlscripts attached below to the bin directory in your JBrowse folder.

3. Run the gffproc.pl program:

gffproc.pl -o [outputdirectory] [inputfile]

This takes a GenBank file from the input path, creates a diectory called outputdirectory, and puts the GFFS in that directory. The script takes GFFs from the BioPerl script and splits them up by feature (the third column, which the BioPerl script writes using SOBA). So if the GenBank file contains 10 annotations (e.g. mRNA, CDS, and mature peptide annotations), each is written to a separate GFF file.

4. The GFF files are ready to be fed into bin/flatfile-to-json.pl.

I AttachmentSorted ascending Action Size Date Who Comment
Txttxt bp_genbank2gff3.pl.txt manage 36.7 K 2009-12-10 - 09:26 MarkDewitt GenBank converter, from BioPerl
Txttxt gffproc.pl.txt manage 1.8 K 2009-12-10 - 09:25 MarkDewitt GFF processing script, my me
Edit | Attach | Print version | History: r8 < r7 < r6 < r5 < r4 | Backlinks | Raw View | Raw edit | More topic actions

This site is powered by the TWiki collaboration platformCopyright © 2008-2014 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback
TWiki Appliance - Powered by TurnKey Linux