<opt
  name="drosmelgb"
  relid="1"
  make_current="1"
  >

  <no_csomesplit>1</no_csomesplit>
  <valid>0</valid> <!-- set 1 to skip db checks of above variables -->
  
  <about id="drosmelgb"><![CDATA[
   See bulkfiles_template.xml for further documentation. 
   See site_defaults for included configurations.
  ]]></about>
  
  <include>site_defaults</include>
  
  <title>DrosMelGb</title>
  
  <doc id="Release.txt"><![CDATA[
  Bulk output for Chado genome databases
  Database: ${title} 
  Species:  ${species}
  Release:  ${rel}, dated ${date} as ${relfull}
  
  This is a test case for <a href="http://gmod.org/GMODTools">GMOD Bulkfiles</a>; 
  loaded from a NCBI GenBank chromosome record
  ]]></doc>
    
  <release id="1" rel="drosmelgb4"  dbname="drosmel4" date="20080512" />
  <release id="dummy" title="placeholder for xmlsimple"/>

  <relfull>${org}_${release_date}</relfull> 
  <release_url>/genome/${species}/release-${release_id}.html</release_url>

  <org>dromel</org>
  <species>Drosophila_melanogaster</species>

  <GFF_source>Genbank</GFF_source>
  <seq_ontology>sequence</seq_ontology> 
  <golden_path>'chromosome','chromosome_arm'</golden_path>

  <!-- some useful Genbank feature properties
    see tables/properties-overview; 
    not keeping: codon_start, date, old_locus_tag, organism, strain?
      translation: should be separate .aa file from fasta/ 
    -->
  <featureprops>
      'Note','comment1','eC_number','exception',
      'gene','locus_tag','inference','map',
      'product','protein_id', 'transcript_id',
      'standard_name', 'strain', 'transposon'
  </featureprops>
  
  <featdump
    path="tmp/featdump/\w+.tsv"
    config="chadofeatsql"
    tag="feature_sql"
    type="feature_table"
    splitname="chadofeat"
    >
    <target>chromosomes</target>
    <target>features</target>
    <target>matches</target>
    <target>analysis</target>
  </featdump>
  
  <!-- feature sets to make fasta bulk files -->
  <featset>chromosome</featset>
  <featset>gene</featset>
  <featset>transcript</featset>
  <featset>CDS</featset>  <!-- dna of protein coding sequence -->
  <featset>translation</featset> <!-- protein aminos -->
  <featset>CDS_translation</featset> <!-- compare to protein from db -->
  <featset>ncRNA</featset>
  <featset>pseudogene</featset>
  <featset>transposon</featset>
  <!-- featset>gene_extended2000</featset -->
  <!-- featset>intergenic</featset -->


  <doc id="TestCase.txt"><![CDATA[

  # get Drosophila_melanogaster CHR_4 and load to chado db
  curl -OR ftp://bio-mirror.net/biomirror/ncbigenomes/Drosophila_melanogaster/CHR_4/NC_004353.gbk.gz

  set dbname=drosmel4
  $pg/bin/createdb -T chado_01_template $dbname
  
  # fix Genbank FT to SO type map
  vi  lib/Bio/SeqFeature/Tools/TypeMapper.pm : add pseudogenic_tRNA
 
  # load chromosome to chado
  gunzip -c NC_004353.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/drosmelgb.xml from template.xml : dbname, etc. edits
  
  # create Bulkfiles outputs
  perl -Ilib bin/bulkfiles.pl -config=drosmelgb -debug -make >& log.drosmel1 &
  # and genbank submit table
  perl -Ilib bin/bulkfiles.pl -config=drosmelgb -format=genbanktbl -debug -make  

 problems:
  transcript : no db residues; reconfig to get froms eq
  
  ]]></doc>

</opt>
