GMOD

GMODTools TestCase

This is an example test case for GMOD, loading a GenBank format genome (chromosome) to Chado database, then dumping it out, including GenBank submission format with GenBank output to test round-trip agreement.

This is an abbreviated example, current as of May 2008. It uses a Chado-db template in Postgres, and assumes you have installed/tested GMOD components available.

 # get AnoGam chrX and load to chado db
 curl -OR ftp://bio-mirror.net/biomirror/ncbigenomes/Anopheles_gambiae/CHR_X/NC_004818.gbk.gz
 or curl -OR ftp://ftp.ncbi.nih.gov/genomes/Anopheles_gambiae/CHR_X/NC_004818.gbk.gz

 set dbname=anogam_x
 $pg/bin/createdb -T chado_01_template $dbname

 # fix Genbank FT to SO type map if needed
 vi  lib/Bio/SeqFeature/Tools/TypeMapper.pm : add pseudogenic tRNA

 # load chromosome X to chado
 gunzip -c NC_004818.gbk.gz |\
  perl bin/bp_genbank2gff3.pl -noCDS -in stdin -out stdout |\
  perl bin/gmod_bulk_load_gff3.pl -dbname $dbname -organism fromdata

 # create GMOD Bulkfiles conf/bulkfiles/anogam.xml from template.xml : dbname, etc. edits
 # see samples to create a new organism/project configuration

 # create standard Bulkfiles outputs for anogam_x
 perl -Ilib bin/bulkfiles.pl -config=anogam  -make >& log.anogam1 &

 # new in progress output, includes regurgitation of GenBank record to compare with original
 # you should check and edit conf/bulkfiles/genbanksubmit.xml before running this part
 perl -Ilib bin/bulkfiles.pl -config=anogam -format=genbanktbl -make

Outputs should include genbanksubmit/ folder with outputs of GMOD Bulkfiles:

anogam-all-anogam1.tbl    : feature table  ** new part
anogam-all-anogam1.fsa    : genome dna     == fasta/chromosome
anogam-all-anogam1.pep    : protein aa     == fasta/translation

and outputs of NCBI tbl2asn:

anogam-all-anogam1.sqn : ASN.1 record to submit to NCBI
anogam-all-anogam1.val : errors & warnings
anogam-all-anogam1.gbf : Genbank format for review

You will find these sample outputs here http://insects.eugenes.org/genome/Anopheles_gambiae_str._PEST/anogam_20080511/

There remain caveats on reproducing the original GenBank record this way. At this writing, the GenbankSubmit function is still in development, and the Genbank to GFF to Chado loading has some losses.

Categories:

Documentation

Community

Tools