Mercurial > repos > iuc > pilon

changeset 0:fe0dc27e6327 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/pilon commit 2645abbd04dd68266f995b8259e991c31388cda8
author iuc
date Sat, 13 Aug 2016 04:50:55 -0400
parents
children 11e5408fd238
files all_fasta.loc.sample pilon.xml test-data/test1-vcf-part test-data/test1.bam test-data/test1.fasta tool_data_table_conf.xml.sample
diffstat 6 files changed, 596 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/all_fasta.loc.sample	Sat Aug 13 04:50:55 2016 -0400
@@ -0,0 +1,18 @@
+#This file lists the locations and dbkeys of all the fasta files
+#under the "genome" directory (a directory that contains a directory
+#for each build). The script extract_fasta.py will generate the file
+#all_fasta.loc. This file has the format (white space characters are
+#TAB characters):
+#
+#<unique_build_id>	<dbkey>	<display_name>	<file_path>
+#
+#So, all_fasta.loc could look something like this:
+#
+#apiMel3	apiMel3	Honeybee (Apis mellifera): apiMel3	/path/to/genome/apiMel3/apiMel3.fa
+#hg19canon	hg19	Human (Homo sapiens): hg19 Canonical	/path/to/genome/hg19/hg19canon.fa
+#hg19full	hg19	Human (Homo sapiens): hg19 Full	/path/to/genome/hg19/hg19full.fa
+#
+#Your all_fasta.loc file should contain an entry for each individual
+#fasta file. So there will be multiple fasta files for each build,
+#such as with hg19 above.
+#
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/pilon.xml	Sat Aug 13 04:50:55 2016 -0400
@@ -0,0 +1,317 @@
+<tool id="pilon" name="pilon" version="0.1">
+    <description>An automated genome assembly improvement and variant detection tool</description>
+    <requirements>
+      <requirement type="package" version="1.18">pilon</requirement>
+    </requirements>
+    <command detect_errors="exit_code"><![CDATA[
+      #if $auto_selection.auto_enabled == "yes"
+        #for $bamfile in $auto_selection.bam
+          ln -f -s "$bamfile" "\$(basename $bamfile)" &&
+          ln -f -s "$bamfile.metadata.bam_index" "\$(basename $bamfile).bai" &&
+        #end for
+      #end if
+      #if $options.selection_mode == "advanced"
+        #if $options.frags_selection.frags_enabled == "yes"
+          #for $bamfile in $options.frags_selection.frags
+            ln -f -s "$bamfile" "\$(basename $bamfile)" &&
+            ln -f -s "$bamfile.metadata.bam_index" "\$(basename $bamfile).bai" &&
+          #end for
+        #end if
+        #if $options.jumps_selection.jumps_enabled == "yes"
+         #for $bamfile in $options.jumps_selection.jumps
+           ln -f -s "$bamfile" "\$(basename $bamfile)" &&
+           ln -f -s "$bamfile.metadata.bam_index" "\$(basename $bamfile).bai" &&
+         #end for
+        #end if
+        #if $options.unpaired_selection.unpaired_enabled == "yes"
+          #for $bamfile in $options.unpaired_selection.unpaired
+            ln -f -s "$bamfile" "\$(basename $bamfile)" &&
+            ln -f -s "$bamfile.metadata.bam_index" "\$(basename $bamfile).bai" &&
+          #end for
+        #end if
+      #end if
+        ln -s -f
+          #if $reference_genome.reference_genome_source == "history"
+            "$reference_genome.history_item"
+          #else
+            "$reference_genome.builtin.fields.path"
+          #end if
+          reference.fasta &&
+        pilon
+        --genome reference.fasta
+        $variant
+        $changes
+        #if $auto_selection.auto_enabled == "yes"
+          #for $bamfile in $auto_selection.bam
+            --bam "\$(basename $bamfile)"
+          #end for
+        #end if
+        #if $options.selection_mode == "advanced"
+          #if $options.frags_selection.frags_enabled == "yes"
+            #for $bamfile in $options.frags_selection.frags
+              --frags "\$(basename $bamfile)"
+            #end for
+          #end if
+          #if $options.jumps_selection.jumps_enabled == "yes"
+            #for $bamfile in $options.jumps_selection.jumps
+              --jumps "\$(basename $bamfile)"
+            #end for
+          #end if
+          #if $options.unpaired_selection.unpaired_enabled == "yes"
+            #for $bamfile in $options.unpaired_selection.unpaired
+              --unpaired "\$(basename $bamfile)"
+            #end for
+          #end if
+          $options.vcfqe
+          $options.vcf
+          $options.tracks
+          --chunksize $options.chunk_size
+          $options.diploid
+          $options.duplicates
+          $options.iupac
+          $options.nonpf
+          #if len($options.targetlist.strip()) > 0
+            --targetlist $options.targetlist
+          #end if
+          --fix $options.fixes
+          $options.verbose
+          --defaultqual $options.defaultqual
+          --flank $options.flank
+          --gapmargin $options.gapmargin
+          --K $options.kmersize
+          --mindepth $options.mindepth
+          --mingap $options.mingap
+          --minmq $options.minmq
+          --minqual $options.minqual
+          $options.nostrays
+        #end if
+        --threads \${GALAXY_SLOTS:-1}
+        --output pilon
+    ]]></command>
+    <inputs>
+      <conditional name="reference_genome">
+        <param label="Source for reference genome used for BAM alignments" name="reference_genome_source" type="select">
+          <option selected="True" value="history">Use a genome from history</option>
+          <option value="builtin">Use a built-in genome"</option>
+        </param>
+        <when value="history">
+          <param format="fasta" type="data" metadata_name="dbkey" name="history_item" />
+        </when>
+        <when value="builtin">
+          <param label="Select a reference genome" name="builtin" type="select">
+            <options from_data_table="all_fasta">
+              <filter column="2" type="sort_by" />
+              <validator message="No genomes are available for the selected input dataset" type="no_options" />
+            </options>
+          </param>
+        </when>
+      </conditional>
+      <conditional name="auto_selection">
+        <param name="auto_enabled" label="Type automatically determined by pilon"
+          type="boolean" checked="true" truevalue="yes" falsevalue="no" />
+        <when value="yes">
+          <param argument="bam" label="Input BAM file" multiple="true" type="data" format="bam"/>
+        </when>
+        <when value="no"></when>
+      </conditional>
+      <param argument="variant" type="boolean" label="Variant calling mode" checked="true" truevalue="--variant" falsevalue=""
+        help="Sets up heuristics for variant calling, as opposed to assembly improvement; equivalent to '--vcf --fix all,breaks'."/>
+      <param argument="changes" type="boolean" label="Create changes file" truevalue="--changes" falsevalue=""
+        help="If specified, a file listing changes in the &lt;output&gt;.fasta will be generated."/>
+      <conditional name="options">
+        <param label="Use advanced options" name="selection_mode" type="select">
+          <option selected="True" value="default">Use default options</option>
+          <option value="advanced">Use advanced options</option>
+        </param>
+        <when value="default"> </when>
+        <when value="advanced">
+          <conditional name="frags_selection">
+            <param name="frags_enabled" label="Paired end fragments" type="boolean" truevalue="yes" falsevalue="no" />
+            <when value="yes">
+              <param argument="frags" label="Input BAM file (paired end fragments)" multiple="true" type="data" format="bam"
+                help="BAM file consisting of fragment paired-end alignments." />
+            </when>
+            <when value="no"></when>
+          </conditional>
+          <conditional name="jumps_selection">
+            <param name="jumps_enabled" label="Mate pairs" type="boolean" truevalue="yes" falsevalue="no" />
+            <when value="yes">
+              <param argument="jumps" label="Input BAM file (mate pairs)" multiple="true" type="data" format="bam"
+                help="BAM file consisting of jump (mate pair) paired-end alignments." />
+            </when>
+            <when value="no"></when>
+          </conditional>
+          <conditional name="unpaired_selection">
+            <param name="unpaired_enabled" label="Unpaired reads"  type="boolean" truevalue="yes" falsevalue="no" />
+            <when value="yes">
+              <param argument="unpaired" label="Input BAM file (unpaired)" multiple="true" type="data" format="bam"
+                help="BAM file consisting of unpaired alignments." />
+            </when>
+            <when value="no"></when>
+          </conditional>
+          <param argument="vcf" type="boolean" checked="false" label="VCF output"
+            truevalue="--vcf" falsevalue="" help="If specified, a vcf file will be generated (even if 'Variant calling mode' is off)"/>
+          <param argument="tracks" type="boolean" checked="false" label="Output annotation tracks"
+            help="Write many track files (*.bed, *.wig) suitable for viewing in a genome browser."
+            truevalue="--tracks" falsevalue="" />
+          <param argument="chunk_size" type="integer" min="1" value="10000000" label="Chunk size"
+            help="Input FASTA elements larger than this will be processed in smaller pieces not to
+            exceed this size." />
+          <param argument="vcfqe" type="boolean" checked="false" label="QE (not QP) in VCF" help="If specified the VCF will contain a QE (quality-weighted evidence) field rather
+          than the default QP (quality-weighted percentage of evidence) field." truevalue="--vcfqe" falsevalue="" />
+          <param argument="fixes" label="Issues that pilon should try and fix" type="select" multiple="true">
+            <option value="all" selected="true">All non-experimental fixes</option>
+            <option value="bases">Individual bases and small indels</option>
+            <option value="gaps">Fill gaps</option>
+            <option value="local">Detect and fix local misassemblies</option>
+            <option value="none">Do none of these fixes (no FASTA will be written)</option>
+            <option value="amb">Fix ambigious bases in FASTA output (experimental)</option>
+            <option value="breaks">Allow local reassembly to open new gaps (experimental, requires local assembly fixing to be selected)</option>
+            <option value="novel">Assemble novel sequence from unaligned non-jump reads (experimental)</option>
+          </param>
+          <param argument="diploid" label="Organism is diploid" type="boolean" checked="false"
+            help="Sample is from diploid organism; will eventually affect calling of heterozygous SNPs"
+            truevalue="--diploid" falsevalue="" />
+          <param argument="duplicates" label="Use duplicates" type="boolean" checked="false"
+            help="Use reads marked as duplicates in the input BAMs"
+            truevalue="--duplicates" falsevalue="" />
+          <param argument="iupac" label="Use IUPAC codes in FASTA output" type="boolean" checked="false"
+            help="Output IUPAC ambiguous base codes in the output FASTA file when appropriate"
+            truevalue="--iupac" falsevalue="" />
+          <param argument="nonpf" label="Use low quality reads" type="boolean" checked="false"
+            help="Use reads which failed sequencer quality filtering"
+            truevalue="--nonpf" falsevalue="" />
+          <param argument="targetlist" label="List of targets to process (leave blank for all)" default="" type="text" length="40"
+            help="Only process the specified target(s).  Targets are comma-separated, and each target is a fasta element name optionally followed by a base range." />
+          <param argument="verbose" label="Verbose output (in tool log)" type="boolean" check="false"
+            truevalue="--verbose" falsevalue="" />
+          <param argument="defaultqual" label="Default base quality" type="integer" min="1" value="15"
+            help="Assumes bases are of this quality if quals are no present in input BAMs" />
+          <param argument="flank" label="Flanking bases to ignore" type="integer" min="1" value="10"
+            help="This many bases at each end of the good reads will be ignored." />
+          <param argument="gapmargin" label="Allowable gap margin" type="integer" min="1" value="100000"
+            help="Closed gaps must be within this number of bases of true size to be closed" />
+          <param argument="kmersize" label="Kmer size" type="integer" min="1" value="47"
+            help="Kmer size used by internal assembler" />
+          <param argument="mindepth" label="Minimum depth" type="float" value="0.1"
+            help="Minimum depth of coverage required for variants to be called. See complete documentation below." />
+          <param argument="mingap" label="Mininum gap size" type="integer" value="10"
+            help="Minimum size for unclosed gaps" />
+          <param argument="minmq" label="Minimum mapping quality" type="integer" value="0"
+            help="Minimum alignment mapping quality for a read to count in pileups" />
+          <param argument="minqual" label="Minimum base quality" type="integer" value="0"
+            help="Minimum base quality to consider for pileups" />
+          <param argument="nostrays" label="Disable 'stray read filtering'" type="boolean" checked="false"
+            help="See documentation below"
+            truevalue="--nostrays" falsevalue="" />
+        </when>
+      </conditional>
+    </inputs>
+    <outputs>
+      <data format="vcf" from_work_dir="pilon.vcf" label="VCF from ${tool.name} on ${on_string}" name="output_vcf">
+        <filter>variant or (options['selection_mode'] == 'advanced' and options['vcf'])</filter>
+      </data>
+      <data format="txt" from_work_dir="pilon.changes" label="Changes in FASTA from ${tool.name} on ${on_string}" name="output_changes">
+        <filter>changes</filter>
+      </data>
+      <data format="fasta" from_work_dir="pilon.fasta" label="FASTA from ${tool.name} on ${on_string}" name="output_fasta">
+        <filter>options['selection_mode'] == 'default' or (options['selection_mode'] == 'advanced' and 'none' not in options['fixes'])</filter>
+      </data>
+      <collection name="tracks" type="list" label="Annotation tracks from ${tool.name} on ${on_string}">
+        <filter>options['selection_mode'] == 'advanced' and options['tracks']</filter>
+        <data format="bed" from_work_dir="pilonPilon.bed" label="Features from ${tool.name} on ${on_string} (BED format)" name="output_pilon_bed" />
+        <data format="wig" from_work_dir="pilonChanges.wig" label="${tool.name} changes track on ${on_string} (WIG format)" name="output_changes_wig" />
+        <data format="wig" from_work_dir="pilonUnconfirmed.wig" label="${tool.name} unconfirmed track on ${on_string}" name="output_unconfirmed_wig" />
+        <data format="wig" from_work_dir="pilonCopyNumber.wig" label="${tool.name} copy number track on ${on_string}" name="output_copynumber_wig" />
+        <data format="wig" from_work_dir="pilonCoverage.wig" label="${tool.name} coverage track on ${on_string}" name="output_coverage_wig" />
+        <data format="wig" from_work_dir="pilonBadCoverage.wig" label="${tool.name} bad coverage track on ${on_string}" name="output_badcoverage_wig" />
+        <data format="wig" from_work_dir="pilonPctBad.wig" label="${tool.name} pct bad track on ${on_string}" name="output_pctbad_wig" />
+        <data format="wig" from_work_dir="pilonDeltaCoverage.wig" label="${tool.name} delta coverage track on ${on_string}" name="output_deltacoverage_wig" />
+        <data format="wig" from_work_dir="pilonDipCoverage.wig" label="${tool.name} dip coverage track on ${on_string}" name="output_dipcoverage_wig" />
+        <data format="wig" from_work_dir="pilonPhysicalCoverage.wig" label="${tool.name} physical coverage track on ${on_string}" name="output_physicalcoverage_wig" />
+        <data format="wig" from_work_dir="pilonClippedAlignments.wig" label="${tool.name} clipped alignments track on ${on_string}" name="output_clippedalignments_wig" />
+        <data format="wig" from_work_dir="pilonWeightedQual.wig" label="${tool.name} weighted quality track on ${on_string}" name="output_weightedqual_wig" />
+        <data format="wig" from_work_dir="pilonWeightedMq.wig" label="${tool.name} weighted MQ track on ${on_string}" name="output_weightedmq_wig" />
+        <data format="wig" from_work_dir="pilonGC.wig" label="${tool.name} GC track on ${on_string}" name="output_gc_wig" />
+      </collection>
+     </outputs>
+    <tests>
+      <test>
+        <param name="reference_genome_source" value="history" />
+        <param ftype="fasta" name="history_item" value="test1.fasta" />
+        <param name="bam" value="test1.bam" />
+        <param name="variant" value="true" />
+        <output file="test1-vcf-part" ftype="vcf" name="output_vcf" compare="contains" />
+        <output md5="352907b0d965bc926289b1b2eb9cbecb" ftype="fasta" name="output_fasta" />
+      </test>
+      <test>
+        <param name="reference_genome_source" value="history" />
+        <param ftype="fasta" name="history_item" value="test1.fasta" />
+        <param name="bam" value="test1.bam" />
+        <param name="variant" value="true" />
+        <param name="selection_mode" value="advanced" />
+        <param name="tracks" value="true" />
+        <output file="test1-vcf-part" ftype="vcf" name="output_vcf" compare="contains" />
+        <output md5="352907b0d965bc926289b1b2eb9cbecb" ftype="fasta" name="output_fasta" />
+        <output_collection name="tracks">
+          <element name="output_pilon_bed" md5="ae0518a6d641efecdcf41c808c014226" />
+          <element name="output_changes_wig" md5="26e0a9a3793e6604673e9911c4d334ab" />
+          <element name="output_unconfirmed_wig" md5="f8bb3def4547e854bd90e1a6b4f8cc66" />
+          <element name="output_copynumber_wig" md5="c1e94296fbc24d00fa7d057f30e4de8f" />
+          <element name="output_coverage_wig" md5="48f7da4bad60fcbb49aaaa992aef8a02" />
+          <element name="output_badcoverage_wig" md5="31e3ad4a5d0f8cd1e296e4b1424d6a95" />
+          <element name="output_pctbad_wig" md5="8353134113c87dbc7d9497f983d58b4b" />
+          <element name="output_deltacoverage_wig" md5="e0fb6474b851e1890f91c57c2ba1fd76" />
+          <element name="output_dipcoverage_wig" md5="8885c63df1dc7309a7bab371b9eb449b" />
+          <element name="output_physicalcoverage_wig" md5="db485f84b51499c3a3c72da70a8ef7af" />
+          <element name="output_clippedalignments_wig" md5="b6eab0827a9b6f2e925d9ece5ee4f87f" />
+          <element name="output_weightedqual_wig" md5="9987289bd0ec8cd05bd6d330bdd1d01d" />
+          <element name="output_weightedmq_wig" md5="7ec1b6420f52fc44bf3c04aa593fcdeb" />
+          <element name="output_gc_wig" md5="815af1378d016b85ed6f52667dde10c1" />
+        </output_collection>
+      </test>
+
+    </tests>
+    <help><![CDATA[
+  Pilon is a software tool which can be used to:
+
+  * Automatically improve draft assemblies
+
+  * Find variation among strains, including large event detection
+
+  Pilon requires as input a FASTA file of the genome along with one or more BAM files of reads aligned to the input FASTA file. Pilon uses read alignment analysis to identify inconsistencies between the input genome and the evidence in the reads. It then attempts to make improvements to the input genome, including:
+
+  * Single base differences
+
+  * Small indels
+
+  * Larger indel or block substitution events
+
+  * Gap filling
+
+  * Identification of local misassemblies, including optional opening of new gaps
+
+  Pilon then outputs a FASTA file containing an improved representation of the genome from the read data and an optional VCF file detailing variation seen between the read data and the input genome.
+
+  To aid manual inspection and improvement by an analyst, Pilon can optionally produce tracks that can be displayed in genome viewers such as IGV and GenomeView, and it reports other events (such as possible large collapsed repeat regions) in its standard output.
+
+  Note on **mindepth**:
+
+  Variants (snps and indels) will only be called if there is coverage of good pairs
+  at the value set for *mindepth* depth or more; if this value is >= 1, it is an absolute depth, if it is a
+  fraction < 1, then minimum depth is computed by multiplying this value by the mean
+  coverage for the region, with a minumum value of 5 (default 0.1: min depth to call
+  is 10% of mean coverage or 5, whichever is greater).
+
+  Note on **stray read filtering**
+
+  By default a pass is made through the input BAM files to identify stray pairs, that is,
+  those pairs in which both reads are aligned but not marked valid because they have
+  inconsistent orientation or separation. Identifying stray pairs can help fill gaps
+  and assemble larger insertions, especially of repeat content.  However, doing so
+  sometimes consumes considerable memory.
+    ]]></help>
+    <citations>
+        <citation type="doi">10.1371/journal.pone.0112963</citation>
+    </citations>
+</tool>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test1-vcf-part	Sat Aug 13 04:50:55 2016 -0400
@@ -0,0 +1,85 @@
+#CHROM	POS	ID	REF	ALT	QUAL	FILTER	INFO	FORMAT	SAMPLE
+gi|448814763|ref|NC_000962.3|	1	.	T	.	1276	PASS	DP=32;TD=65;BQ=40;MQ=46;QD=39;BC=0,0,0,32;QP=0,0,0,100;PC=55;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	2	.	T	.	1408	PASS	DP=36;TD=73;BQ=39;MQ=47;QD=39;BC=0,0,0,36;QP=0,0,0,100;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	3	.	G	.	1450	PASS	DP=37;TD=74;BQ=39;MQ=47;QD=39;BC=0,0,37,0;QP=0,0,100,0;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	4	.	A	.	1456	PASS	DP=37;TD=77;BQ=39;MQ=47;QD=39;BC=37,0,0,0;QP=100,0,0,0;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	5	.	C	.	1567	PASS	DP=39;TD=79;BQ=40;MQ=48;QD=40;BC=0,39,0,0;QP=0,100,0,0;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	6	.	C	.	1606	PASS	DP=40;TD=80;BQ=40;MQ=48;QD=40;BC=0,40,0,0;QP=0,100,0,0;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	7	.	G	.	1586	PASS	DP=40;TD=80;BQ=40;MQ=48;QD=39;BC=0,0,40,0;QP=0,0,100,0;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	8	.	A	.	1557	PASS	DP=40;TD=81;BQ=39;MQ=48;QD=38;BC=40,0,0,0;QP=100,0,0,0;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	9	.	T	.	1641	PASS	DP=42;TD=84;BQ=39;MQ=48;QD=39;BC=0,0,0,42;QP=0,0,0,100;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	10	.	G	.	1748	PASS	DP=45;TD=87;BQ=39;MQ=49;QD=38;BC=0,0,45,0;QP=0,0,100,0;PC=56;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	11	.	A	.	2104	PASS	DP=55;TD=104;BQ=38;MQ=51;QD=38;BC=55,0,0,0;QP=100,0,0,0;PC=58;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	12	.	C	.	2222	PASS	DP=56;TD=106;BQ=40;MQ=51;QD=39;BC=0,56,0,0;QP=0,100,0,0;PC=59;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	13	.	C	.	2231	PASS	DP=56;TD=106;BQ=40;MQ=51;QD=39;BC=0,56,0,0;QP=0,100,0,0;PC=59;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	14	.	C	.	2237	PASS	DP=56;TD=107;BQ=40;MQ=51;QD=39;BC=0,56,0,0;QP=0,100,0,0;PC=59;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	15	.	C	.	2230	PASS	DP=56;TD=107;BQ=40;MQ=51;QD=39;BC=0,56,0,0;QP=0,100,0,0;PC=60;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	16	.	G	.	2236	PASS	DP=56;TD=107;BQ=40;MQ=51;QD=39;BC=0,0,56,0;QP=0,0,100,0;PC=61;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	17	.	G	.	2213	PASS	DP=56;TD=108;BQ=40;MQ=51;QD=39;BC=0,0,56,0;QP=0,0,100,0;PC=64;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	18	.	T	.	2178	PASS	DP=56;TD=108;BQ=39;MQ=51;QD=38;BC=0,0,0,56;QP=0,0,0,100;PC=65;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	19	.	T	.	2161	PASS	DP=56;TD=110;BQ=39;MQ=51;QD=38;BC=0,0,0,56;QP=0,0,0,100;PC=66;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	20	.	C	.	2164	PASS	DP=56;TD=110;BQ=39;MQ=51;QD=38;BC=0,56,0,0;QP=0,100,0,0;PC=70;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	21	.	A	.	2188	PASS	DP=57;TD=112;BQ=38;MQ=51;QD=38;BC=57,0,0,0;QP=100,0,0,0;PC=72;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	22	.	G	.	2127	PASS	DP=56;TD=112;BQ=38;MQ=51;QD=37;BC=0,0,56,0;QP=0,0,100,0;PC=72;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	23	.	G	.	2102	PASS	DP=55;TD=112;BQ=38;MQ=51;QD=38;BC=0,0,55,0;QP=0,0,100,0;PC=73;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	24	.	C	.	2082	PASS	DP=55;TD=113;BQ=38;MQ=51;QD=37;BC=0,55,0,0;QP=0,100,0,0;PC=74;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	25	.	T	.	2138	PASS	DP=56;TD=111;BQ=38;MQ=51;QD=38;BC=0,0,0,56;QP=0,0,0,100;PC=74;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	26	.	T	.	2126	PASS	DP=55;TD=112;BQ=39;MQ=51;QD=38;BC=0,0,0,55;QP=0,0,0,100;PC=75;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	27	.	C	.	2193	PASS	DP=58;TD=115;BQ=38;MQ=51;QD=37;BC=0,58,0,0;QP=0,100,0,0;PC=75;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	28	.	A	.	2073	PASS	DP=58;TD=116;BQ=36;MQ=51;QD=35;BC=57,1,0,0;QP=100,0,0,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	29	.	C	.	2208	PASS	DP=58;TD=117;BQ=38;MQ=52;QD=38;BC=0,58,0,0;QP=0,100,0,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	30	.	C	.	2314	PASS	DP=61;TD=120;BQ=38;MQ=52;QD=37;BC=0,61,0,0;QP=0,100,0,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	31	.	A	.	2362	PASS	DP=63;TD=120;BQ=38;MQ=52;QD=37;BC=63,0,0,0;QP=100,0,0,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	32	.	C	.	2380	PASS	DP=63;TD=120;BQ=38;MQ=52;QD=37;BC=0,63,0,0;QP=0,100,0,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	33	.	A	.	2332	PASS	DP=62;TD=120;BQ=38;MQ=53;QD=37;BC=62,0,0,0;QP=100,0,0,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	34	.	G	.	2382	PASS	DP=63;TD=120;BQ=38;MQ=53;QD=37;BC=0,0,63,0;QP=0,0,100,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	35	.	T	.	2298	PASS	DP=62;TD=120;BQ=37;MQ=53;QD=37;BC=0,0,0,62;QP=0,0,0,100;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	36	.	G	.	2248	PASS	DP=60;TD=119;BQ=37;MQ=54;QD=37;BC=0,0,60,0;QP=0,0,100,0;PC=77;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	37	.	T	.	2064	PASS	DP=59;TD=118;BQ=36;MQ=54;QD=34;BC=0,0,1,58;QP=0,0,0,100;PC=78;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	38	.	G	.	2181	PASS	DP=61;TD=119;BQ=36;MQ=54;QD=35;BC=1,0,60,0;QP=0,0,100,0;PC=79;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	39	.	G	.	2217	PASS	DP=60;TD=116;BQ=37;MQ=54;QD=36;BC=0,0,60,0;QP=0,0,100,0;PC=80;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	40	.	A	.	2100	PASS	DP=56;TD=112;BQ=38;MQ=55;QD=37;BC=56,0,0,0;QP=100,0,0,0;PC=82;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	41	.	A	.	2057	PASS	DP=56;TD=110;BQ=37;MQ=55;QD=36;BC=56,0,0,0;QP=100,0,0,0;PC=83;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	42	.	C	.	2069	PASS	DP=56;TD=111;BQ=37;MQ=55;QD=36;BC=0,56,0,0;QP=0,100,0,0;PC=83;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	43	.	G	.	1987	PASS	DP=56;TD=111;BQ=35;MQ=55;QD=35;BC=0,0,56,0;QP=0,0,100,0;PC=84;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	44	.	C	.	1985	PASS	DP=56;TD=111;BQ=35;MQ=55;QD=35;BC=0,56,0,0;QP=0,100,0,0;PC=84;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	45	.	G	.	2018	PASS	DP=56;TD=104;BQ=36;MQ=55;QD=36;BC=0,0,56,0;QP=0,0,100,0;PC=85;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	46	.	G	.	2029	PASS	DP=56;TD=106;BQ=36;MQ=55;QD=36;BC=0,0,56,0;QP=0,0,100,0;PC=85;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	47	.	T	.	1935	PASS	DP=57;TD=109;BQ=34;MQ=55;QD=33;BC=0,0,0,57;QP=0,0,0,100;PC=85;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	48	.	C	.	1871	PASS	DP=57;TD=109;BQ=34;MQ=55;QD=32;BC=0,56,0,1;QP=0,99,0,1;PC=85;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	49	.	G	.	1974	PASS	DP=56;TD=109;BQ=35;MQ=56;QD=35;BC=0,0,56,0;QP=0,0,100,0;PC=86;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	50	.	T	.	2024	PASS	DP=58;TD=110;BQ=35;MQ=56;QD=34;BC=0,0,0,58;QP=0,0,0,100;PC=86;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	51	.	C	.	2041	PASS	DP=58;TD=111;BQ=35;MQ=56;QD=35;BC=0,58,0,0;QP=0,100,0,0;PC=86;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	52	.	T	.	1935	PASS	DP=54;TD=108;BQ=36;MQ=56;QD=35;BC=0,0,0,54;QP=0,0,0,100;PC=86;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	53	.	C	.	2018	PASS	DP=55;TD=111;BQ=37;MQ=56;QD=36;BC=0,55,0,0;QP=0,100,0,0;PC=88;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	54	.	C	.	1945	PASS	DP=54;TD=109;BQ=36;MQ=56;QD=36;BC=0,54,0,0;QP=0,100,0,0;PC=88;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	55	.	G	.	1957	PASS	DP=54;TD=108;BQ=36;MQ=56;QD=36;BC=0,0,54,0;QP=0,0,100,0;PC=88;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	56	.	A	.	1758	PASS	DP=51;TD=105;BQ=34;MQ=57;QD=34;BC=51,0,0,0;QP=100,0,0,0;PC=88;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	57	.	A	.	1686	PASS	DP=50;TD=104;BQ=34;MQ=57;QD=33;BC=50,0,0,0;QP=100,0,0,0;PC=89;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	58	.	C	.	1542	PASS	DP=47;TD=101;BQ=33;MQ=57;QD=32;BC=0,47,0,0;QP=0,100,0,0;PC=89;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	59	.	T	.	1656	PASS	DP=48;TD=102;BQ=35;MQ=57;QD=34;BC=0,0,0,48;QP=0,0,0,100;PC=90;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	60	.	T	.	1610	PASS	DP=45;TD=99;BQ=36;MQ=57;QD=35;BC=0,0,0,45;QP=0,0,0,100;PC=90;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	61	.	A	.	1494	PASS	DP=43;TD=96;BQ=35;MQ=57;QD=34;BC=43,0,0,0;QP=100,0,0,0;PC=91;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	62	.	A	.	1479	PASS	DP=43;TD=95;BQ=34;MQ=57;QD=34;BC=43,0,0,0;QP=100,0,0,0;PC=91;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	63	.	C	.	1569	PASS	DP=45;TD=96;BQ=35;MQ=58;QD=34;BC=0,45,0,0;QP=0,100,0,0;PC=91;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	64	.	G	.	1519	PASS	DP=44;TD=95;BQ=35;MQ=58;QD=34;BC=0,0,44,0;QP=0,0,100,0;PC=91;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	65	.	G	.	1472	PASS	DP=43;TD=92;BQ=34;MQ=58;QD=34;BC=0,0,43,0;QP=0,0,100,0;PC=91;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	66	.	C	.	1336	PASS	DP=40;TD=91;BQ=33;MQ=59;QD=33;BC=0,40,0,0;QP=0,100,0,0;PC=91;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	67	.	G	.	1338	PASS	DP=40;TD=91;BQ=33;MQ=59;QD=33;BC=0,0,40,0;QP=0,0,100,0;PC=91;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	68	.	A	.	1232	PASS	DP=38;TD=87;BQ=32;MQ=59;QD=32;BC=38,0,0,0;QP=100,0,0,0;PC=92;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	69	.	C	.	1192	PASS	DP=35;TD=86;BQ=34;MQ=59;QD=34;BC=0,35,0,0;QP=0,100,0,0;PC=93;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	70	.	C	.	1115	PASS	DP=32;TD=81;BQ=35;MQ=59;QD=34;BC=0,32,0,0;QP=0,100,0,0;PC=93;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	71	.	C	.	1093	PASS	DP=31;TD=84;BQ=35;MQ=60;QD=35;BC=0,31,0,0;QP=0,100,0,0;PC=93;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	72	.	T	.	1055	PASS	DP=30;TD=80;BQ=35;MQ=60;QD=35;BC=0,0,0,30;QP=0,0,0,100;PC=94;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	73	.	A	.	1016	PASS	DP=29;TD=80;BQ=35;MQ=60;QD=35;BC=29,0,0,0;QP=100,0,0,0;PC=99;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	74	.	A	.	994	PASS	DP=29;TD=80;BQ=34;MQ=60;QD=34;BC=29,0,0,0;QP=100,0,0,0;PC=100;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	75	.	G	.	1023	PASS	DP=29;TD=81;BQ=35;MQ=60;QD=35;BC=0,0,29,0;QP=0,0,100,0;PC=102;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	76	.	G	.	996	PASS	DP=28;TD=80;BQ=36;MQ=60;QD=35;BC=0,0,28,0;QP=0,0,100,0;PC=104;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	77	.	T	.	919	PASS	DP=26;TD=79;BQ=35;MQ=60;QD=35;BC=0,0,0,26;QP=0,0,0,100;PC=104;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	78	.	T	.	912	PASS	DP=26;TD=79;BQ=35;MQ=60;QD=35;BC=0,0,0,26;QP=0,0,0,100;PC=104;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	79	.	G	.	832	PASS	DP=25;TD=78;BQ=33;MQ=60;QD=33;BC=0,0,25,0;QP=0,0,100,0;PC=104;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	80	.	A	.	788	PASS	DP=23;TD=77;BQ=34;MQ=60;QD=34;BC=23,0,0,0;QP=100,0,0,0;PC=104;IC=0;DC=0;XC=1;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	81	.	C	.	803	PASS	DP=23;TD=75;BQ=35;MQ=60;QD=34;BC=0,23,0,0;QP=0,100,0,0;PC=104;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	82	.	G	.	807	PASS	DP=24;TD=75;BQ=34;MQ=60;QD=33;BC=0,0,24,0;QP=0,0,100,0;PC=104;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	83	.	A	.	1025	PASS	DP=29;TD=82;BQ=35;MQ=60;QD=35;BC=29,0,0,0;QP=100,0,0,0;PC=106;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
+gi|448814763|ref|NC_000962.3|	84	.	C	.	1080	PASS	DP=30;TD=84;BQ=36;MQ=60;QD=36;BC=0,30,0,0;QP=0,100,0,0;PC=106;IC=0;DC=0;XC=0;AC=0;AF=0.00	GT	0/0
Binary file test-data/test1.bam has changed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test1.fasta	Sat Aug 13 04:50:55 2016 -0400
@@ -0,0 +1,168 @@
+>gi|448814763|ref|NC_000962.3| NC_000962.3 Mycobacterium tuberculosis H37Rv, complete genome
+TTGACCGATGACCCCGGTTCAGGCTTCACCACAGTGTGGAACGCGGTCGTCTCCGAACTT
+AACGGCGACCCTAAGGTTGACGACGGACCCAGCAGTGATGCTAATCTCAGCGCTCCGCTG
+ACCCCTCAGCAAAGGGCTTGGCTCAATCTCGTCCAGCCATTGACCATCGTCGAGGGGTTT
+GCTCTGTTATCCGTGCCGAGCAGCTTTGTCCAAAACGAAATCGAGCGCCATCTGCGGGCC
+CCGATTACCGACGCTCTCAGCCGCCGACTCGGACATCAGATCCAACTCGGGGTCCGCATC
+GCTCCGCCGGCGACCGACGAAGCCGACGACACTACCGTGCCGCCTTCCGAAAATCCTGCT
+ACCACATCGCCAGACACCACAACCGACAACGACGAGATTGATGACAGCGCTGCGGCACGG
+GGCGATAACCAGCACAGTTGGCCAAGTTACTTCACCGAGCGCCCGCACAATACCGATTCC
+GCTACCGCTGGCGTAACCAGCCTTAACCGTCGCTACACCTTTGATACGTTCGTTATCGGC
+GCCTCCAACCGGTTCGCGCACGCCGCCGCCTTGGCGATCGCAGAAGCACCCGCCCGCGCT
+TACAACCCCCTGTTCATCTGGGGCGAGTCCGGTCTCGGCAAGACACACCTGCTACACGCG
+GCAGGCAACTATGCCCAACGGTTGTTCCCGGGAATGCGGGTCAAATATGTCTCCACCGAG
+GAATTCACCAACGACTTCATTAACTCGCTCCGCGATGACCGCAAGGTCGCATTCAAACGC
+AGCTACCGCGACGTAGACGTGCTGTTGGTCGACGACATCCAATTCATTGAAGGCAAAGAG
+GGTATTCAAGAGGAGTTCTTCCACACCTTCAACACCTTGCACAATGCCAACAAGCAAATC
+GTCATCTCATCTGACCGCCCACCCAAGCAGCTCGCCACCCTCGAGGACCGGCTGAGAACC
+CGCTTTGAGTGGGGGCTGATCACTGACGTACAACCACCCGAGCTGGAGACCCGCATCGCC
+ATCTTGCGCAAGAAAGCACAGATGGAACGGCTCGCGGTCCCCGACGATGTCCTCGAACTC
+ATCGCCAGCAGTATCGAACGCAATATCCGTGAACTCGAGGGCGCGCTGATCCGGGTCACC
+GCGTTCGCCTCATTGAACAAAACACCAATCGACAAAGCGCTGGCCGAGATTGTGCTTCGC
+GATCTGATCGCCGACGCCAACACCATGCAAATCAGCGCGGCGACGATCATGGCTGCCACC
+GCCGAATACTTCGACACTACCGTCGAAGAGCTTCGCGGGCCCGGCAAGACCCGAGCACTG
+GCCCAGTCACGACAGATTGCGATGTACCTGTGTCGTGAGCTCACCGATCTTTCGTTGCCC
+AAAATCGGCCAAGCGTTCGGCCGTGATCACACAACCGTCATGTACGCCCAACGCAAGATC
+CTGTCCGAGATGGCCGAGCGCCGTGAGGTCTTTGATCACGTCAAAGAACTCACCACTCGC
+ATCCGTCAGCGCTCCAAGCGCTAGCACGGCGTGTTCTTCCGACAACGTTCTTAAAAAAAC
+TTCTCTCTCCCAGGTCACACCAGTCACAGAGATTGGCTGTGAGTGTCGCTGTGCACAAAC
+CGCGCACAGACTCATACAGTCCCGGCGGTTCCGTTCACAACCCACGCCTCATCCCCACCG
+ACCCAACACACACCCCACAGTCATCGCCACCGTCATCCACAACTCCGACCGACGTCGACC
+TGCACCAAGACCAGACTGTCCCCAAACTGCACACCCTCTAATACTGTTACCGAGATTTCT
+TCGTCGTTTGTTCTTGGAAAGACAGCGCTGGGGATCGTTCGCTGGATACCACCCGCATAA
+CTGGCTCGTCGCGGTGGGTCAGAGGTCAATGATGAACTTTCAAGTTGACGTGAGAAGCTC
+TACGGTTGTTGTTCGACTGCTGTTGCGGCCGTCGTGGCGGGTCACGCGTCATGGGCATTC
+GTCGTTGGCAGTCCCCACGCTAGCGGGGCGCTAGCCACGGGATCGAACTCATCGTGAGGT
+GAAAGGGCGCAATGGACGCGGCTACGACAAGAGTTGGCCTCACCGACTTGACGTTTCGTT
+TGCTACGAGAGTCTTTCGCCGATGCGGTGTCGTGGGTGGCTAAAAATCTGCCAGCCAGGC
+CCGCGGTGCCGGTGCTCTCCGGCGTGTTGTTGACCGGCTCGGACAACGGTCTGACGATTT
+CCGGATTCGACTACGAGGTTTCCGCCGAGGCCCAGGTTGGCGCTGAAATTGTTTCTCCTG
+GAAGCGTTTTAGTTTCTGGCCGATTGTTGTCCGATATTACCCGGGCGTTGCCTAACAAGC
+CCGTAGACGTTCATGTCGAAGGTAACCGGGTCGCATTGACCTGCGGTAACGCCAGGTTTT
+CGCTACCGACGATGCCAGTCGAGGATTATCCGACGCTGCCGACGCTGCCGGAAGAGACCG
+GATTGTTGCCTGCGGAATTATTCGCCGAGGCAATCAGTCAGGTCGCTATCGCCGCCGGCC
+GGGACGACACGTTGCCTATGTTGACCGGCATCCGGGTCGAAATCCTCGGTGAGACGGTGG
+TTTTGGCCGCTACCGACAGGTTTCGCCTGGCTGTTCGAGAACTGAAGTGGTCGGCGTCGT
+CGCCAGATATCGAAGCGGCTGTGCTGGTCCCGGCCAAGACGCTGGCCGAGGCCGCCAAAG
+CGGGCATCGGCGGCTCTGACGTTCGTTTGTCGTTGGGTACTGGGCCGGGGGTGGGCAAGG
+ATGGCCTGCTCGGTATCAGTGGGAACGGCAAGCGCAGCACCACGCGACTTCTTGATGCCG
+AGTTCCCGAAGTTTCGGCAGTTGCTACCAACCGAACACACCGCGGTGGCCACCATGGACG
+TGGCCGAGTTGATCGAAGCGATCAAGCTGGTTGCGTTGGTAGCTGATCGGGGCGCGCAGG
+TGCGCATGGAGTTCGCTGATGGCAGCGTGCGGCTTTCTGCGGGTGCCGATGATGTTGGAC
+GAGCCGAGGAAGATCTTGTTGTTGACTATGCCGGTGAACCATTGACGATTGCGTTTAACC
+CAACCTATCTAACGGACGGTTTGAGTTCGTTGCGCTCGGAGCGAGTGTCTTTCGGGTTTA
+CGACTGCGGGTAAGCCTGCCTTGCTACGTCCGGTGTCCGGGGACGATCGCCCTGTGGCGG
+GTCTGAATGGCAACGGTCCGTTCCCGGCGGTGTCGACGGACTATGTCTATCTGTTGATGC
+CGGTTCGGTTGCCGGGCTGAGCACTTGGCGCCCGGGTAGGTGTACGTCCGTCATTTGGGG
+CTGCGTGACTTCCGGTCCTGGGCATGTGTAGATCTGGAATTGCATCCAGGGCGGACGGTT
+TTTGTTGGGCCTAACGGTTATGGTAAGACGAATCTTATTGAGGCACTGTGGTATTCGACG
+ACGTTAGGTTCGCACCGCGTTAGCGCCGATTTGCCGTTGATCCGGGTAGGTACCGATCGT
+GCGGTGATCTCCACGATCGTGGTGAACGACGGTAGAGAATGTGCCGTCGACCTCGAGATC
+GCCACGGGGCGAGTCAACAAAGCGCGATTGAATCGATCATCGGTCCGAAGTACACGTGAT
+GTGGTCGGAGTGCTTCGAGCTGTGTTGTTTGCCCCTGAGGATCTGGGGTTGGTTCGTGGG
+GATCCCGCTGACCGGCGGCGCTATCTGGATGATCTGGCGATCGTGCGTAGGCCTGCGATC
+GCTGCGGTACGAGCCGAATATGAGAGGGTGTTGCGCCAGCGGACGGCGTTATTGAAGTCC
+GTACCTGGAGCACGGTATCGGGGTGACCGGGGTGTGTTTGACACTCTTGAGGTATGGGAC
+AGTCGTTTGGCGGAGCACGGGGCTGAACTGGTGGCCGCCCGCATCGATTTGGTCAACCAG
+TTGGCACCGGAAGTGAAGAAGGCATACCAGCTGTTGGCGCCGGAATCGCGATCGGCGTCT
+ATCGGTTATCGGGCCAGCATGGATGTAACCGGTCCCAGCGAGCAGTCAGATATCGATCGG
+CAATTGTTAGCAGCTCGGCTGTTGGCGGCGCTGGCGGCCCGTCGGGATGCCGAACTCGAG
+CGTGGGGTTTGTCTAGTTGGTCCGCACCGTGACGACCTAATACTGCGACTAGGCGATCAA
+CCCGCGAAAGGATTTGCTAGCCATGGGGAGGCGTGGTCGTTGGCGGTGGCACTGCGGTTG
+GCGGCCTATCAACTGTTACGCGTTGATGGTGGTGAGCCGGTGTTGTTGCTCGACGACGTG
+TTCGCCGAACTGGATGTCATGCGCCGTCGAGCGTTGGCGACGGCGGCCGAGTCCGCCGAA
+CAGGTGTTGGTGACTGCCGCGGTGCTCGAGGATATTCCCGCCGGCTGGGACGCCAGGCGG
+GTGCACATCGATGTGCGTGCCGATGACACCGGATCGATGTCGGTGGTTCTGCCATGACGG
+GTTCTGTTGACCGGCCCGACCAGAATCGCGGTGAGCGATCAATGAAGTCACCAGGGTTGG
+ATTTGGTCAGGCGCACCCTGGACGAAGCTCGTGCTGCTGCCCGCGCGCGCGGACAAGACG
+CCGGTCGAGGGCGGGTCGCTTCCGTTGCGTCGGGTCGGGTGGCCGGACGGCGACGAAGCT
+GGTCGGGTCCGGGGCCCGACATTCGTGATCCACAACCGCTGGGTAAGGCCGCTCGTGAGC
+TGGCAAAGAAACGCGGCTGGTCGGTGCGGGTCGCCGAGGGTATGGTGCTCGGCCAGTGGT
+CTGCGGTGGTCGGCCACCAGATCGCCGAACATGCACGCCCGACTGCGCTAAACGACGGGG
+TGTTGAGCGTGATTGCGGAGTCGACGGCGTGGGCGACGCAGTTGAGGATCATGCAGGCCC
+AGCTTCTGGCCAAGATCGCCGCAGCGGTTGGCAACGATGTGGTGCGATCGCTAAAGATCA
+CCGGGCCGGCGGCACCATCGTGGCGCAAGGGGCCTCGCCATATTGCCGGTAGGGGTCCGC
+GCGACACCTACGGATAACACGTCGATCGGCCCAGAACAAGGCGCTCCGGTCCCGGCCTGA
+GAGCCTCGAGGACGAAGCGGATCCGTATGCCGGACGTCGGGACGCACCAGGAAGAAAGAT
+GTCCGACGCACGGCGCGGTTAGATGGGTAAAAACGAGGCCAGAAGATCGGCCCTGGCGCC
+CGATCACGGTACAGTGGTGTGCGACCCCCTGCGGCGACTCAACCGCATGCACGCAACCCC
+TGAGGAGAGTATTCGGATCGTGGCTGCCCAGAAAAAGAAGGCCCAAGACGAATACGGCGC
+TGCGTCTATCACCATTCTCGAAGGGCTGGAGGCCGTCCGCAAACGTCCCGGCATGTACAT
+TGGCTCGACCGGTGAGCGCGGTTTACACCATCTCATTTGGGAGGTGGTCGACAACGCGGT
+CGACGAGGCGATGGCCGGTTATGCAACCACAGTGAACGTAGTGCTGCTTGAGGATGGCGG
+TGTCGAGGTCGCCGACGACGGCCGCGGCATTCCGGTCGCCACCCACGCCTCCGGCATACC
+GACCGTCGACGTGGTGATGACACAACTACATGCCGGCGGCAAGTTCGACTCGGACGCGTA
+TGCGATATCTGGTGGTCTGCACGGCGTCGGCGTGTCGGTGGTTAACGCGCTATCCACCCG
+GCTCGAAGTCGAGATCAAGCGCGACGGGTACGAGTGGTCTCAGGTTTATGAGAAGTCGGA
+ACCCCTGGGCCTCAAGCAAGGGGCGCCGACCAAGAAGACGGGGTCAACGGTGCGGTTCTG
+GGCCGACCCCGCTGTTTTCGAAACCACGGAATACGACTTCGAAACCGTCGCCCGCCGGCT
+GCAAGAGATGGCGTTCCTCAACAAGGGGCTGACCATCAACCTGACCGACGAGAGGGTGAC
+CCAAGACGAGGTCGTCGACGAAGTGGTCAGCGACGTCGCCGAGGCGCCGAAGTCGGCAAG
+TGAACGCGCAGCCGAATCCACTGCACCGCACAAAGTTAAGAGCCGCACCTTTCACTATCC
+GGGTGGCCTGGTGGACTTCGTGAAACACATCAACCGCACCAAGAACGCGATTCATAGCAG
+CATCGTGGACTTTTCCGGCAAGGGCACCGGGCACGAGGTGGAGATCGCGATGCAATGGAA
+CGCCGGGTATTCGGAGTCGGTGCACACCTTCGCCAACACCATCAACACCCACGAGGGCGG
+CACCCACGAAGAGGGCTTCCGCAGCGCGCTGACGTCGGTGGTGAACAAGTACGCCAAGGA
+CCGCAAGCTACTGAAGGACAAGGACCCCAACCTCACCGGTGACGATATCCGGGAAGGCCT
+GGCCGCTGTGATCTCGGTGAAGGTCAGCGAACCGCAGTTCGAGGGCCAGACCAAGACCAA
+GTTGGGCAACACCGAGGTCAAATCGTTTGTGCAGAAGGTCTGTAACGAACAGCTGACCCA
+CTGGTTTGAAGCCAACCCCACCGACGCGAAAGTCGTTGTGAACAAGGCTGTGTCCTCGGC
+GCAAGCCCGTATCGCGGCACGTAAGGCACGAGAGTTGGTGCGGCGTAAGAGCGCCACCGA
+CATCGGTGGATTGCCCGGCAAGCTGGCCGATTGCCGTTCCACGGATCCGCGCAAGTCCGA
+ACTGTATGTCGTAGAAGGTGACTCGGCCGGCGGTTCTGCAAAAAGCGGTCGCGATTCGAT
+GTTCCAGGCGATACTTCCGCTGCGCGGCAAGATCATCAATGTGGAGAAAGCGCGCATCGA
+CCGGGTGCTAAAGAACACCGAAGTTCAGGCGATCATCACGGCGCTGGGCACCGGGATCCA
+CGACGAGTTCGATATCGGCAAGCTGCGCTACCACAAGATCGTGCTGATGGCCGACGCCGA
+TGTTGACGGCCAACATATTTCCACGCTGTTGTTGACGTTGTTGTTCCGGTTCATGCGGCC
+GCTCATCGAGAACGGGCATGTGTTTTTGGCACAACCGCCGCTGTACAAACTCAAGTGGCA
+GCGCAGTGACCCGGAATTCGCATACTCCGACCGCGAGCGCGACGGTCTGCTGGAGGCGGG
+GCTGAAGGCCGGGAAGAAGATCAACAAGGAAGACGGCATTCAGCGGTACAAGGGTCTAGG
+TGAAATGGACGCTAAGGAGTTGTGGGAGACCACCATGGATCCCTCGGTTCGTGTGTTGCG
+TCAAGTGACGCTGGACGACGCCGCCGCCGCCGACGAGTTGTTCTCCATCCTGATGGGCGA
+GGACGTCGACGCGCGGCGCAGCTTTATCACCCGCAACGCCAAGGATGTTCGGTTCCTGGA
+TGTCTAACGCAACCCTGCGTTCGATTGCAAACGAGGAATAGATGACAGACACGACGTTGC
+CGCCTGACGACTCGCTCGACCGGATCGAACCGGTTGACATCGAGCAGGAGATGCAGCGCA
+GCTACATCGACTATGCGATGAGCGTGATCGTCGGCCGCGCGCTGCCGGAGGTGCGCGACG
+GGCTCAAGCCCGTGCATCGCCGGGTGCTCTATGCAATGTTCGATTCCGGCTTCCGCCCGG
+ACCGCAGCCACGCCAAGTCGGCCCGGTCGGTTGCCGAGACCATGGGCAACTACCACCCGC
+ACGGCGACGCGTCGATCTACGACAGCCTGGTGCGCATGGCCCAGCCCTGGTCGCTGCGCT
+ACCCGCTGGTGGACGGCCAGGGCAACTTCGGCTCGCCAGGCAATGACCCACCGGCGGCGA
+TGAGGTACACCGAAGCCCGGCTGACCCCGTTGGCGATGGAGATGCTGAGGGAAATCGACG
+AGGAGACAGTCGATTTCATCCCTAACTACGACGGCCGGGTGCAAGAGCCGACGGTGCTAC
+CCAGCCGGTTCCCCAACCTGCTGGCCAACGGGTCAGGCGGCATCGCGGTCGGCATGGCAA
+CCAATATCCCGCCGCACAACCTGCGTGAGCTGGCCGACGCGGTGTTCTGGGCGCTGGAGA
+ATCACGACGCCGACGAAGAGGAGACCCTGGCCGCGGTCATGGGGCGGGTTAAAGGCCCGG
+ACTTCCCGACCGCCGGACTGATCGTCGGATCCCAGGGCACCGCTGATGCCTACAAAACTG
+GCCGCGGCTCCATTCGAATGCGCGGAGTTGTTGAGGTAGAAGAGGATTCCCGCGGTCGTA
+CCTCGCTGGTGATCACCGAGTTGCCGTATCAGGTCAACCACGACAACTTCATCACTTCGA
+TCGCCGAACAGGTCCGAGACGGCAAGCTGGCCGGCATTTCCAACATTGAGGACCAGTCTA
+GCGATCGGGTCGGTTTACGCATCGTCATCGAGATCAAGCGCGATGCGGTGGCCAAGGTGG
+TGATCAATAACCTTTACAAGCACACCCAGCTGCAGACCAGCTTTGGCGCCAACATGCTAG
+CGATCGTCGACGGGGTGCCGCGCACGCTGCGGCTGGACCAGCTGATCCGCTATTACGTTG
+ACCACCAACTCGACGTCATTGTGCGGCGCACCACCTACCGGCTGCGCAAGGCAAACGAGC
+GAGCCCACATTCTGCGCGGCCTGGTTAAAGCGCTCGACGCGCTGGACGAGGTCATTGCAC
+TGATCCGGGCGTCGGAGACCGTCGATATCGCCCGGGCCGGACTGATCGAGCTGCTCGACA
+TCGACGAGATCCAGGCCCAGGCAATCCTGGACATGCAGTTGCGGCGCCTGGCCGCACTGG
+AACGCCAGCGCATCATCGACGACCTGGCCAAAATCGAGGCCGAGATCGCCGATCTGGAAG
+ACATCCTGGCAAAACCCGAGCGGCAGCGTGGGATCGTGCGCGACGAACTCGCCGAAATCG
+TGGACAGGCACGGCGACGACCGGCGTACCCGGATCATCGCGGCCGACGGAGACGTCAGCG
+ACGAGGATTTGATCGCCCGCGAGGACGTCGTTGTCACTATCACCGAAACGGGATACGCCA
+AGCGCACCAAGACCGATCTGTATCGCAGCCAGAAACGCGGCGGCAAGGGCGTGCAGGGTG
+CGGGGTTGAAGCAGGACGACATCGTCGCGCACTTCTTCGTGTGCTCCACCCACGATTTGA
+TCCTGTTCTTCACCACCCAGGGACGGGTTTATCGGGCCAAGGCCTACGACTTGCCCGAGG
+CCTCCCGGACGGCGCGCGGGCAGCACGTGGCCAACCTGTTAGCCTTCCAGCCCGAGGAAC
+GCATCGCCCAGGTCATCCAGATTCGCGGCTACACCGACGCCCCGTACCTGGTGCTGGCCA
+CTCGCAACGGGCTGGTGAAAAAGTCCAAGCTGACCGACTTCGACTCCAATCGCTCGGGCG
+GAATCGTGGCGGTCAACCTGCGCGACAACGACGAGCTGGTCGGTGCGGTGCTGTGTTCGG
+CCGGCGACGACCTGCTGCTGGTCTCGGCCAACGGGCAGTCCATCAGGTTCTCGGCGACCG
+ACGAGGCGCTGCGGCCAATGGGTCGTGCCACCTCGGGTGTGCAGGGCATGCGGTTCAATA
+TCGACGACCGGCTGCTGTCGCTGAACGTCGTGCGTGAAGGCACCTATCTGCTGGTGGCGA
+CGTCAGGGGGCTATGCGAAACGTACCGCGATCGAGGAATACCCGGTACAGGGCCGCGGCG
+GTAAAGGTGTGCTGACGGTCATGTACGACCGCCGGCGCGGCAGGTTGGTTGGGGCGTTGA
+TTGTCGACGACGACAGCGAGCTGTATGCCGTCACTTCCGGCGGTGGCGTGATCCGCACCG
+CGGCACGCCAGGTTCGCAAGGCGGGACGGCAGACCAAGGGTGTTCGGTTGATGAATCTGG
+GCGAGGGCGACACACTGTTGGCCATCGCGCGCAACGCCGAAGAAAGTGGCGACGATAATG
+CCGTGGACGCCAACGGCGCAGACCAGACGGGCAATTAATCAGGCTCGCCCGACGACGATG
+CGGATCGCGTAGCGATCTGAGGAGGAATCGGGCAGCTAGGCTCGGCAGCCGGGTACGAGT
+GTTAGGAGTCGGGGTGACTGCACCGAACGAGCCGGGGGCGCTCAGCAAGGGCGACGGCCC
+GAATGCGGATGGCTTGGTCGACCGTGGGGGCGCACATCGG
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_data_table_conf.xml.sample	Sat Aug 13 04:50:55 2016 -0400
@@ -0,0 +1,8 @@
+<tables>
+    <!-- Locations of all fasta files under genome directory -->
+    <table name="all_fasta" comment_char="#">
+        <columns>value, dbkey, name, path</columns>
+        <file path="tool-data/all_fasta.loc" />
+    </table>
+</tables>
+