changeset 0:26d7ec4af119 draft

"planemo upload for repository https://github.com/sunnyisgalaxy/moabs commit fca680a439f168971afc9944ccbbdd9b3b65c845"
author iuc
date Fri, 06 Sep 2019 09:54:27 -0400
parents
children 8c8cc81b34cd
files all_fasta.loc.sample macros.xml moabs.xml test-data/6_all_1.fq.gz test-data/6_all_2.fq.gz test-data/8_all_1.fq.gz test-data/8_all_2.fq.gz test-data/PE_comp.g1.vs.g2.txt test-data/PE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr test-data/SE_comp.g1.vs.g2.txt test-data/SE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr test-data/SE_g1.G.bed test-data/SE_g1_r1.bam test-data/SE_g1_r2.bam test-data/SE_g2.G.bed test-data/SE_g2_r1.bam test-data/SE_g2_r2.bam test-data/TKO2PE1.fastq.gz test-data/TKO2PE2.fastq.gz test-data/WTPE1.fastq.gz test-data/WTPE2.fastq.gz test-data/all_fasta.loc test-data/chr11.fa test-data/seg.fa tool_data_table_conf.xml.sample tool_data_table_conf.xml.test
diffstat 26 files changed, 2415 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/all_fasta.loc.sample	Fri Sep 06 09:54:27 2019 -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/macros.xml	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,29 @@
+<?xml version="1.0"?>
+<macros>
+    <xml name="requirements">
+        <requirements>
+            <requirement type="package" version="@VERSION@">moabs</requirement>
+        </requirements>
+    </xml>
+
+    <token name="@VERSION@">1.3.4.6</token>
+
+    <token name="@DMRExample@">
+
+Example::
+
+ chr11	316	1705	-	35	30	5	-1
+ chr11	3052	3674	-	17	16	1	-1
+ chr11	18599	19013	-	9	7	3	-1
+ chr11	19120	19931	-	39	39	0	-1
+ chr11	22138	22618	-	14	12	3	-1
+ chr11	28140	28356	-	4	4	0	-1
+    </token>
+
+    <xml name="citations">
+        <citations>
+            <citation type="doi">10.1186/gb-2014-15-2-r38</citation>
+        </citations>
+    </xml>
+
+</macros>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/moabs.xml	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,465 @@
+<tool id="moabs" name="MOABS" profile="16.04" version="@VERSION@">
+    <description>MOdel based Analysis of Bisulfite Sequencing data</description>
+    <macros>
+        <import>macros.xml</import>
+    </macros>
+    <expand macro="requirements"/>
+    <command detect_errors="exit_code">
+        <![CDATA[
+        #if str( $mcomp_advanced.doComp.compare_selector ) == "0":
+            cp -f '$mcomp_advanced.doComp.compFile' comp.g1.vs.g2.txt &&
+        #end if
+        moabs -v 1 --def MMAP.p="\${GALAXY_SLOTS:-4}" --def MCALL.p="\${GALAXY_SLOTS:-4}" --def MCOMP.p="\${GALAXY_SLOTS:-4}" --cf '$cfg_file' &&
+        #if "1" in $output_selector:
+            cp -f dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr '$output1' &&
+        #end if
+        #if "2" in $output_selector:
+            cp -f comp.g1.vs.g2.txt '$output2' &&
+        #end if
+        echo Done
+        ]]>
+    </command>
+    <configfiles>
+        <configfile name="cfg_file">
+            #if str( $reference_source.reference_source_selector ) == "history":
+                #set $reference_fasta_filename = $reference_source.ref_file
+            #else:
+                #set $reference_fasta_filename = $reference_source.ref_file.fields.path
+            #end if
+            [INPUT]
+            #for $i, $s in enumerate( $g1_fastq )
+                #if str( $s.fastq_input.fastq_input_selector ) == "paired":
+                    s1_r${i+1}_1='$s.fastq_input.fastq_input1'
+                    s1_r${i+1}_2='$s.fastq_input.fastq_input2'
+                #elif str( $s.fastq_input.fastq_input_selector ) == "paired_collection":
+                    s1_r${i+1}_1='$s.fastq_input.fastq_input1.forward'
+                    s1_r${i+1}_2='$s.fastq_input.fastq_input1.reverse'
+                #else:
+                    s1_r${i+1}='$s.fastq_input.fastq_input1'
+                #end if
+            #end for
+
+            #for $i, $s in enumerate( $g2_fastq )
+                #if str( $s.fastq_input.fastq_input_selector ) == "paired":
+                    s2_r${i+1}_1='$s.fastq_input.fastq_input1'
+                    s2_r${i+1}_2='$s.fastq_input.fastq_input2'
+                #elif str( $s.fastq_input.fastq_input_selector ) == "paired_collection":
+                    s2_r${i+1}_1='$s.fastq_input.fastq_input1.forward'
+                    s2_r${i+1}_2='$s.fastq_input.fastq_input1.reverse'
+                #else:
+                    s2_r${i+1}='$s.fastq_input.fastq_input1'
+                #end if
+            #end for
+
+            [TASK]
+            Program=MMAP
+            Label=g1,g2
+            Parallel=NONE
+
+            [MMAP]
+            Path=bsmap
+            d='${reference_fasta_filename}'
+            #if str( $bsmap_advanced.bsmap_mismatch.bsmap_mismatch_selector ) != "0":
+                v=$bsmap_advanced.bsmap_mismatch.v
+            #end if
+            n=$bsmap_advanced.n
+            r=$bsmap_advanced.r
+            R=''
+
+            [MCALL]
+            Path=mcall
+            r='${reference_fasta_filename}'
+
+            [MCOMP]
+            Path=mcomp
+            reference='${reference_fasta_filename}'
+            doComp=$mcomp_advanced.doComp.compare_selector
+        </configfile>
+    </configfiles>
+
+    <inputs>
+        <conditional name="reference_source">
+            <param name="reference_source_selector" type="select" label="Will you select a reference genome from your history or use a cache FASTA?" help="Cached FASTA">
+                <option value="cached">Use a cached genome FASTA</option>
+                <option value="history">Use a genome FASTA from history</option>
+            </param>
+            <when value="cached">
+                <param name="ref_file" type="select" label="Using reference genome" help="Select genome from the list">
+                    <options from_data_table="all_fasta">
+                        <filter type="sort_by" column="2" />
+                        <validator type="no_options" message="No genome FASTA are available" />
+                    </options>
+                </param>
+            </when>
+            <when value="history">
+                <param name="ref_file" type="data" format="fasta" label="Use the following dataset as the reference sequence" help="You can upload a FASTA sequence to the history and use it as reference" />
+            </when>
+        </conditional>
+        <repeat name="g1_fastq" title="Group1: fastq files" min="1">
+            <conditional name="fastq_input">
+                <param name="fastq_input_selector" type="select" label="Single or Paired-end reads" help="Select between paired and single end data">
+                    <option value="single">Single</option>
+                    <option value="paired">Paired</option>
+                    <option value="paired_collection">Paired Collection</option>
+                </param>
+                <when value="paired">
+                    <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz" label="Select first set of reads" help="Specify dataset with forward reads"/>
+                    <param name="fastq_input2" type="data" format="fastqsanger,fastqsanger.gz" label="Select second set of reads" help="Specify dataset with reverse reads"/>
+                </when>
+                <when value="single">
+                    <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz" label="Select fastq dataset" help="Specify dataset with single reads"/>
+                </when>
+                <when value="paired_collection">
+                    <param name="fastq_input1" format="fastqsanger,fastqsanger.gz" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
+                </when>
+            </conditional>
+        </repeat>
+        <repeat name="g2_fastq" title="Group2: fastq files" min="1">
+            <conditional name="fastq_input">
+                <param name="fastq_input_selector" type="select" label="Single or Paired-end reads" help="Select between paired and single end data">
+                    <option value="single">Single</option>
+                    <option value="paired">Paired</option>
+                    <option value="paired_collection">Paired Collection</option>
+                </param>
+                <when value="paired">
+                    <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz" label="Select first set of reads" help="Specify dataset with forward reads"/>
+                    <param name="fastq_input2" type="data" format="fastqsanger,fastqsanger.gz" label="Select second set of reads" help="Specify dataset with reverse reads"/>
+                </when>
+                <when value="single">
+                    <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz" label="Select fastq dataset" help="Specify dataset with single reads"/>
+                </when>
+                <when value="paired_collection">
+                    <param name="fastq_input1" format="fastqsanger,fastqsanger.gz" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
+                </when>
+            </conditional>
+        </repeat>
+        <section name="bsmap_advanced" title="Advanced options for BSMAP" expanded="False">
+            <conditional name="bsmap_mismatch">
+                <param name="bsmap_mismatch_selector" type="select" label="Set the mismatch rate or number?" help="">
+                    <option value="0">Do not set</option>
+                    <option value="1">Set the mismatch rate</option>
+                    <option value="2">Set the mismatch number</option>
+                </param>
+                <when value="1">
+                    <param argument="-v" type="float" value="0.08" min="0" max="1" label="Mismatch rate" help="The mismatch rate w.r.t to the read length"/>
+                </when>
+                <when value="2">
+                    <param argument="-v" type="integer" value="3" min="0" label="Mismatch number" help="The maximum number of mismatches allowed on a read"/>
+                </when>
+            </conditional>
+            <param argument="-n" type="boolean" truevalue="1" falsevalue="0" checked="true" label="Mapping to four strands?" help="Yes: map SE or PE reads to all 4 strands, i.e. ++, +-, -+, --; No: only map to 2 forward strands, i.e. BSW(++) and BSC(-+)"/>
+            <param argument="-r" type="select" label="How to report repeat hits" help="0=none(unique hit/pair); 1=random one; 2=all(slow)">
+                <option value="0" selected="true">0</option>
+                <option value="1">1</option>
+                <option value="2">2</option>
+            </param>
+        </section>
+        <section name="mcomp_advanced" title="Advanced options for MCOMP" expanded="False">
+            <conditional name="doComp">
+                <param name="compare_selector" type="select" label="Run the comparison or not" help="Yes: compare; No: do not compare, using the comparison result by `-c`">
+                    <option value="1">Yes</option>
+                    <option value="0">No</option>
+                </param>
+                <when value="0">
+                    <param argument="-c" name="compFile" type="data" format="txt" label="Input comparison results" help="Previously generated comparison file from history"/>
+                </when>
+            </conditional>
+        </section>
+        <param name="output_selector" type="select" multiple="true" optional="true" label="Select output files" help="">
+            <option value="1"> dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr </option>
+            <option value="2"> comp.g1.vs.g2.txt </option>
+            <option value="3"> BAM files </option>
+            <option value="4"> Methylation calling BED files </option>
+        </param>
+    </inputs>
+    <outputs>
+        <data name="output1" format="interval" label="${tool.name} on ${on_string} : dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr">
+            <filter> "1" in output_selector </filter>
+        </data>
+        <data name="output2" format="interval" label="${tool.name} on ${on_string} : comp.g1.vs.g2.txt">
+            <filter> "2" in output_selector </filter>
+        </data>
+        <collection name="output_collection_bam" type="list" label="BAM files">
+            <filter> "3" in output_selector </filter>
+            <discover_datasets pattern="(?P&lt;designation&gt;.+\.bam$)" ext='bam'/>
+        </collection>
+        <collection name="output_collection_bed" type="list" label="Methylation calling BED files">
+            <filter> "4" in output_selector </filter>
+            <discover_datasets pattern="(?P&lt;designation&gt;g[12]\.G\.bed$)" ext='interval'/>
+        </collection>
+    </outputs>
+    <tests>
+        <test>
+            <!-- test single-end reads -->
+            <param name="reference_source_selector" value="history"/>
+            <param name="ref_file" ftype="fasta" value="chr11.fa"/>
+            <repeat name="g1_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="WTPE1.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <repeat name="g1_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="WTPE2.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <repeat name="g2_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="TKO2PE1.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <repeat name="g2_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="TKO2PE2.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <conditional name="doComp">
+                <param name="compare_selector" value="0"/>
+                <param name="compFile" value="SE_comp.g1.vs.g2.txt"/>
+            </conditional>
+            <!--
+            <conditional name="doComp">
+                <param name="compare_selector" value="1"/>
+            </conditional>
+            -->
+            <param name="output_selector" value="1,2,3,4"/>
+            <output name="output1" file="SE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr" ftype="interval" lines_diff="1"/>
+            <output name="output2" file="SE_comp.g1.vs.g2.txt" ftype="interval" lines_diff="1"/>
+            <output_collection name="output_collection_bam" count="4">
+                <element name="g1_r1.bam" file="SE_g1_r1.bam" compare="sim_size"/>
+                <element name="g1_r2.bam" file="SE_g1_r2.bam" compare="sim_size"/>
+                <element name="g2_r1.bam" file="SE_g2_r1.bam" compare="sim_size"/>
+                <element name="g2_r2.bam" file="SE_g2_r2.bam" compare="sim_size"/>
+            </output_collection>
+            <output_collection name="output_collection_bed" count="2">
+                <element name="g1.G.bed" file="SE_g1.G.bed" ftype="interval" lines_diff="1"/>
+                <element name="g2.G.bed" file="SE_g2.G.bed" ftype="interval" lines_diff="1"/>
+            </output_collection>
+        </test>
+        <test>
+            <!-- test paired-end reads -->
+            <param name="reference_source_selector" value="history"/>
+            <param name="ref_file" ftype="fasta" value="seg.fa"/>
+            <repeat name="g1_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="paired"/>
+                    <param name="fastq_input1" value="6_all_1.fq.gz"/>
+                    <param name="fastq_input2" value="6_all_2.fq.gz"/>
+                </conditional>
+            </repeat>
+            <repeat name="g2_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="paired"/>
+                    <param name="fastq_input1" value="8_all_1.fq.gz"/>
+                    <param name="fastq_input2" value="8_all_2.fq.gz"/>
+                </conditional>
+            </repeat>
+            <conditional name="doComp">
+                <param name="compare_selector" value="0"/>
+                <param name="compFile" value="PE_comp.g1.vs.g2.txt"/>
+            </conditional>
+            <!--
+            <conditional name="doComp">
+                <param name="compare_selector" value="1"/>
+            </conditional>
+            -->
+            <param name="output_selector" value="1,2"/>
+            <output name="output1" file="PE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr" ftype="interval" lines_diff="1"/>
+            <output name="output2" file="PE_comp.g1.vs.g2.txt" ftype="interval" lines_diff="1"/>
+        </test>
+        <test>
+            <!-- test paired collection -->
+            <param name="reference_source_selector" value="history"/>
+            <param name="ref_file" ftype="fasta" value="seg.fa"/>
+            <repeat name="g1_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="paired_collection"/>
+                    <param name="fastq_input1">
+                        <collection type="paired">
+                            <element name="forward" value="6_all_1.fq.gz" />
+                            <element name="reverse" value="6_all_2.fq.gz" />
+                        </collection>
+                    </param>
+                </conditional>
+            </repeat>
+            <repeat name="g2_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="paired_collection"/>
+                    <param name="fastq_input1">
+                        <collection type="paired">
+                            <element name="forward" value="8_all_1.fq.gz" />
+                            <element name="reverse" value="8_all_2.fq.gz" />
+                        </collection>
+                    </param>
+                </conditional>
+            </repeat>
+            <conditional name="doComp">
+                <param name="compare_selector" value="0"/>
+                <param name="compFile" value="PE_comp.g1.vs.g2.txt"/>
+            </conditional>
+            <!--
+            <conditional name="doComp">
+                <param name="compare_selector" value="1"/>
+            </conditional>
+            -->
+            <param name="output_selector" value="1,2"/>
+            <output name="output1" file="PE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr" ftype="interval" lines_diff="1"/>
+            <output name="output2" file="PE_comp.g1.vs.g2.txt" ftype="interval" lines_diff="1"/>
+        </test>
+        <test>
+            <!-- test data table reference -->
+            <param name="reference_source_selector" value="cached"/>
+            <param name="ref_file" value="chr11"/>
+            <repeat name="g1_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="WTPE1.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <repeat name="g1_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="WTPE2.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <repeat name="g2_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="TKO2PE1.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <repeat name="g2_fastq">
+                <conditional name="fastq_input">
+                    <param name="fastq_input_selector" value="single"/>
+                    <param name="fastq_input1" value="TKO2PE2.fastq.gz"/>
+                </conditional>
+            </repeat>
+            <conditional name="doComp">
+                <param name="compare_selector" value="0"/>
+                <param name="compFile" value="SE_comp.g1.vs.g2.txt"/>
+            </conditional>
+            <!--
+            <conditional name="doComp">
+                <param name="compare_selector" value="1"/>
+            </conditional>
+            -->
+            <param name="output_selector" value="1,2"/>
+            <output name="output1" file="SE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr" ftype="interval" lines_diff="1"/>
+            <output name="output2" file="SE_comp.g1.vs.g2.txt" ftype="interval" lines_diff="1"/>
+        </test>
+    </tests>
+    <help>
+        <![CDATA[
+**MOABS: MOdel based Analysis of Bisulfite Sequencing data**
+
+MOABS is a comprehensive, accurate and efficient solution for analysis of large
+scale base-resolution DNA methylation data, bisulfite sequencing or single
+molecule direct sequencing.
+
+MOABS seamlessly integrates alignment, methylation calling, identification of
+hypomethylation for one sample and differential methylation for multiple
+samples, and other downstream analysis.
+
+For more information, check https://github.com/sunnyisgalaxy/moabs.
+
+-----
+
+**Input files**
+
+MOABS needs to input Bisulfite sequencing reads in two groups of interest, e.g.
+KO vs WT. Each group of reads may have combined sequencing library, i.e.
+single-end reads and/or paired-end reads. Multiple replicates can be specified in each group.
+
+**Outputs**
+
+Four output files can be selected to report, namely
+
+    1. **DMR region file** - the major result file
+    2. **Comparison file between two groups** - the intermediate comparion result
+    3. **BAM files** - intermediate BAM files
+    4. **Methylation BED files** - intermediate methylation BED files
+
+-----
+
+MOABS will detect differential methylated regions (DMRs) using the input BS-Seq
+reads. The output file is a tab-delimited text file (not strictly a BED
+format), representing DMRs. It has 8 columns as below.
+
+chrom<TAB>start<TAB>end<TAB>methylation_state<TAB>CpGsites<TAB>DMCcount<TAB>nonDMCcount<TAB>hidden_state
+
+    1. **chrom** - The chromosome of the region.
+    2. **start** - The start genomic locus of the region.
+    3. **end** - The end genomic locus of the region.
+    4. **methylation_state** - The methylation state of the region, "+"/"-" representing hyper- or hypo-methylation regions.
+    5. **CpGsites** - Total number of CpG sites in the region.
+    6. **DMCcount** - The number of differential methylated CpG sites (DMCs) in the region.
+    7. **nonDMCcount** - The number of non-DMCs in the region.
+    8. **hidden_state** - The hidden state prediced by Hidden Markov Model (HMM), "1"/"-1" representing hyper- or hypo-methylation states.
+
+For example, six DMRs are identified in the following format.
+
+@DMRExample@
+
+-----
+
+The intermediate comparison file summarizes methylation ratio comparison
+results on CpG sites. It has 19 columns as below.
+
+    1. **chrom** - The chromosome of the GpG site.
+    2. **start** - The start position of the site.
+    3. **end** - The end position of the site.
+    4. **single** - The next two columns are attributes for the single position.
+    5. **totalC_0** - Total number of Cs in the first group.
+    6. **nominalRatio_0** - Nominal methylation ratio in the first group.
+    7. **ratioCI_0** - The confidence interval of the methylation ratio in the first group.
+    8. **single** - The next two columns are attributes for the single position.
+    9. **totalC_1** - Total number of Cs in the second group.
+    10. **nominalRatio_1** - Nominal methylation ratio in the second group.
+    11. **ratioCI_1** - The confidence interval of the methylation ratio in the second group.
+    12. **pair** - The next three columns are attributes for pairs of groups.
+    13. **nominalDif_1-0** - Nominal difference of methylation ratio between group 1 and group 0.
+    14. **credibleDif_1-0** - Credible methylation difference between group 1 and group 0.
+    15. **difCI_1-0** - Difference of confidence intervals between group 1 and group 0.
+    16. **p_sim** - The next column is the simulation p-value.
+    17. **p_sim_1_v_0** - Simulation p-value between group 1 and group 0.
+    18. **p_fet** - The next column is the FET p-value.
+    19. **p_fet_1_v_0** - FET p-value between group 1 and group 0.
+
+The comparison result file can be reused for DMR calling.
+
+-----
+
+BAM files are intermediate mapping results of input reads to the referene
+genome. These BAM files can be reused in downstream methylation analysis.
+
+-----
+
+Methylation calling BED files are intermediate methylation calling results of
+Cs in two groups of input reads. These methyation calling results can be easily
+reused in downstream DMR calling and visualization. The BED file has 15 columns
+as below.
+
+    1. **chrom** - The chromosome of the site.
+    2. **start** - The start position of the site.
+    3. **end** - The end position of the site.
+    4. **ratio** - Methylation ratio in the site
+    5. **totalC** - Total number of reads in current Cs.
+    6. **methC** - Methylated Cs.
+    7. **strand** - The strand information for prevous three columns.
+    8. **next** - The next base.
+    9. **Plus** - Next two columns are for forward strand.
+    10. **totalC** - Total number of Cs.
+    11. **methC** - Methylated Cs.
+    12. **Minus** - Next two columns are for reverse strand.
+    13. **totalC** - Total number of Cs.
+    14. **methC** - Methylated Cs.
+    15. **localSeq** - Local sequences.
+
+]]>
+    </help>
+    <expand macro="citations"/>
+</tool>
Binary file test-data/6_all_1.fq.gz has changed
Binary file test-data/6_all_2.fq.gz has changed
Binary file test-data/8_all_1.fq.gz has changed
Binary file test-data/8_all_2.fq.gz has changed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/PE_comp.g1.vs.g2.txt	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,211 @@
+#chrom	start	end	single:	totalC_0	nominalRatio_0	ratioCI_0	single:	totalC_1	nominalRatio_1	ratioCI_1	pair:	nominalDif_1-0	credibleDif_1-0	difCI_1-0	p_sim:	p_sim_1_v_0	p_fet:	p_fet_1_v_0
+seg	87	89	single:	6	1	0.652,1	single:	10	0.2	0,0.47	pair:	-0.8	-0.327	-0.857,-0.327	p_acc:	0.000318	p_fet:	0.00699
+seg	152	154	single:	4	1	0.549,1	single:	9	0.333	0.0575,0.609	pair:	-0.667	-0.113	-0.761,-0.113	p_acc:	0.0035	p_fet:	0.0699
+seg	253	255	single:	NA	NA	NA,NA	single:	9	0	0,0.259	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	255	257	single:	NA	NA	NA,NA	single:	10	0	0,0.238	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	258	260	single:	3	1	0.473,1	single:	13	0.231	0,0.466	pair:	-0.769	-0.165	-0.809,-0.165	p_acc:	0.00236	p_fet:	0.0357
+seg	331	333	single:	6	0.667	0.341,0.992	single:	22	0.227	0.0499,0.405	pair:	-0.439	-0.058	-0.66,-0.058	p_acc:	0.00636	p_fet:	0.0638
+seg	468	470	single:	NA	NA	NA,NA	single:	7	0.143	0,0.471	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	490	492	single:	NA	NA	NA,NA	single:	3	0.333	0,0.751	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	630	632	single:	3	1	0.473,1	single:	5	0	0,0.393	pair:	-1	-0.271	-1,-0.271	p_acc:	0.000954	p_fet:	0.0179
+seg	638	640	single:	3	0.667	0.249,1	single:	5	0.6	0.264,0.936	pair:	-0.0667	-0	-0.461,0.418	p_acc:	0.0286	p_fet:	1
+seg	641	643	single:	4	1	0.549,1	single:	5	0.4	0.064,0.736	pair:	-0.6	-0.0122	-0.748,-0.0122	p_acc:	0.00715	p_fet:	0.167
+seg	645	647	single:	4	1	0.549,1	single:	8	1	0.717,1	pair:	0	0	-0.183,0.374	p_acc:	0.067	p_fet:	1
+seg	666	668	single:	4	1	0.549,1	single:	9	0.667	0.391,0.942	pair:	-0.333	-0	-0.506,0.15	p_acc:	0.021	p_fet:	0.497
+seg	685	687	single:	3	1	0.473,1	single:	8	0	0,0.283	pair:	-1	-0.342	-1,-0.342	p_acc:	0.000364	p_fet:	0.00606
+seg	689	691	single:	4	1	0.549,1	single:	9	0.667	0.391,0.942	pair:	-0.333	-0	-0.506,0.15	p_acc:	0.021	p_fet:	0.497
+seg	697	699	single:	3	1	0.473,1	single:	7	0.571	0.27,0.872	pair:	-0.429	-0	-0.595,0.156	p_acc:	0.017	p_fet:	0.475
+seg	699	701	single:	3	1	0.473,1	single:	9	1	0.741,1	pair:	0	0	-0.149,0.454	p_acc:	0.0596	p_fet:	1
+seg	706	708	single:	NA	NA	NA,NA	single:	4	0.25	0,0.657	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	708	710	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	738	740	single:	NA	NA	NA,NA	single:	5	0	0,0.393	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	749	751	single:	NA	NA	NA,NA	single:	8	0.25	0,0.55	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	753	755	single:	NA	NA	NA,NA	single:	12	0.167	0,0.41	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	820	822	single:	7	0.857	0.529,1	single:	20	0.4	0.2,0.6	pair:	-0.457	-0.0729	-0.62,-0.0729	p_acc:	0.00528	p_fet:	0.0768
+seg	830	832	single:	7	0.857	0.529,1	single:	18	0.222	0.0252,0.419	pair:	-0.635	-0.236	-0.766,-0.236	p_acc:	0.000768	p_fet:	0.0068
+seg	842	844	single:	4	0.75	0.343,1	single:	10	0.3	0.035,0.565	pair:	-0.45	-0	-0.672,0.0536	p_acc:	0.0117	p_fet:	0.245
+seg	1149	1151	single:	NA	NA	NA,NA	single:	7	0	0,0.312	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1151	1153	single:	NA	NA	NA,NA	single:	10	0.1	0,0.364	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1172	1174	single:	3	1	0.473,1	single:	12	0.333	0.0889,0.578	pair:	-0.667	-0.0658	-0.738,-0.0658	p_acc:	0.005	p_fet:	0.0769
+seg	1180	1182	single:	3	1	0.473,1	single:	11	0.0909	0,0.339	pair:	-0.909	-0.286	-0.909,-0.286	p_acc:	0.000705	p_fet:	0.011
+seg	1201	1203	single:	3	0.667	0.249,1	single:	14	0.214	0,0.44	pair:	-0.452	-0	-0.7,0.0401	p_acc:	0.0118	p_fet:	0.191
+seg	1203	1205	single:	3	0.667	0.249,1	single:	13	0.308	0.0724,0.543	pair:	-0.359	-0	-0.629,0.133	p_acc:	0.0176	p_fet:	0.518
+seg	1303	1305	single:	3	1	0.473,1	single:	4	0	0,0.451	pair:	-1	-0.229	-1,-0.229	p_acc:	0.00143	p_fet:	0.0286
+seg	1333	1335	single:	4	0.25	0,0.657	single:	6	0.167	0,0.521	pair:	-0.0833	-0	-0.47,0.291	p_acc:	0.0339	p_fet:	1
+seg	1343	1345	single:	3	0	0,0.527	single:	10	0.2	0,0.47	pair:	0.2	0	-0.318,0.357	p_acc:	0.0363	p_fet:	1
+seg	1407	1409	single:	NA	NA	NA,NA	single:	4	0	0,0.451	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1409	1411	single:	4	0.75	0.343,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1412	1414	single:	3	1	0.473,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1422	1424	single:	5	0.6	0.264,0.936	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1428	1430	single:	5	0.6	0.264,0.936	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1445	1447	single:	6	0.833	0.479,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1448	1450	single:	6	0.833	0.479,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1453	1455	single:	6	0.833	0.479,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1488	1490	single:	NA	NA	NA,NA	single:	5	0.4	0.064,0.736	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1494	1496	single:	NA	NA	NA,NA	single:	5	0.4	0.064,0.736	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1510	1512	single:	NA	NA	NA,NA	single:	5	0.4	0.064,0.736	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1532	1534	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1538	1540	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1986	1988	single:	4	1	0.549,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1993	1995	single:	4	1	0.549,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	1997	1999	single:	4	1	0.549,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2002	2004	single:	4	0.75	0.343,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2015	2017	single:	4	1	0.549,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2019	2021	single:	4	0.75	0.343,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2028	2030	single:	4	1	0.549,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2032	2034	single:	4	1	0.549,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2035	2037	single:	3	1	0.473,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2040	2042	single:	4	0.75	0.343,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2310	2312	single:	8	1	0.717,1	single:	11	0.455	0.198,0.711	pair:	-0.545	-0.162	-0.69,-0.162	p_acc:	0.00184	p_fet:	0.0181
+seg	2336	2338	single:	8	1	0.717,1	single:	5	0.6	0.264,0.936	pair:	-0.4	-0.0115	-0.656,-0.0115	p_acc:	0.0099	p_fet:	0.128
+seg	2396	2398	single:	8	1	0.717,1	single:	12	0.333	0.0889,0.578	pair:	-0.667	-0.273	-0.774,-0.273	p_acc:	0.000439	p_fet:	0.0047
+seg	2413	2415	single:	NA	NA	NA,NA	single:	4	0	0,0.451	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2425	2427	single:	NA	NA	NA,NA	single:	5	0	0,0.393	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2448	2450	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2481	2483	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2488	2490	single:	3	0.333	0,0.751	single:	7	0	0,0.312	pair:	-0.333	-0	-0.669,0.0603	p_acc:	0.0175	p_fet:	0.3
+seg	2516	2518	single:	7	0.571	0.27,0.872	single:	16	0.5	0.278,0.722	pair:	-0.0714	-0	-0.372,0.269	p_acc:	0.0377	p_fet:	1
+seg	2529	2531	single:	7	0.571	0.27,0.872	single:	25	0.04	0,0.17	pair:	-0.531	-0.203	-0.745,-0.203	p_acc:	0.000552	p_fet:	0.00445
+seg	2548	2550	single:	10	0.7	0.435,0.965	single:	17	0.0588	0,0.238	pair:	-0.641	-0.303	-0.787,-0.303	p_acc:	0.000131	p_fet:	0.000939
+seg	2550	2552	single:	10	0.8	0.53,1	single:	20	0.4	0.2,0.6	pair:	-0.4	-0.0665	-0.586,-0.0665	p_acc:	0.00582	p_fet:	0.0577
+seg	2556	2558	single:	11	0.636	0.382,0.891	single:	19	0.0526	0,0.216	pair:	-0.584	-0.271	-0.746,-0.271	p_acc:	0.000167	p_fet:	0.0011
+seg	2567	2569	single:	9	0.889	0.606,1	single:	11	0.364	0.109,0.618	pair:	-0.525	-0.135	-0.697,-0.135	p_acc:	0.0027	p_fet:	0.0281
+seg	2573	2575	single:	10	0.9	0.636,1	single:	12	0.167	0,0.41	pair:	-0.733	-0.352	-0.835,-0.352	p_acc:	0.000117	p_fet:	0.00191
+seg	2578	2580	single:	9	0.889	0.606,1	single:	10	0.5	0.234,0.766	pair:	-0.389	-0.0151	-0.6,-0.0151	p_acc:	0.0092	p_fet:	0.141
+seg	2580	2582	single:	11	0.909	0.661,1	single:	9	0.222	0,0.507	pair:	-0.687	-0.287	-0.812,-0.287	p_acc:	0.00036	p_fet:	0.00452
+seg	2603	2605	single:	8	1	0.717,1	single:	3	0.333	0,0.751	pair:	-0.667	-0.119	-0.834,-0.119	p_acc:	0.00328	p_fet:	0.0545
+seg	2625	2627	single:	8	1	0.717,1	single:	5	0.6	0.264,0.936	pair:	-0.4	-0.0115	-0.656,-0.0115	p_acc:	0.0099	p_fet:	0.128
+seg	2641	2643	single:	10	0.9	0.636,1	single:	10	0.2	0,0.47	pair:	-0.7	-0.299	-0.817,-0.299	p_acc:	0.00031	p_fet:	0.00548
+seg	2645	2647	single:	9	0.667	0.391,0.942	single:	8	0.25	0,0.55	pair:	-0.417	-0	-0.642,0.00262	p_acc:	0.00968	p_fet:	0.153
+seg	2659	2661	single:	11	0.727	0.473,0.982	single:	7	0.429	0.128,0.73	pair:	-0.299	-0	-0.567,0.0906	p_acc:	0.0183	p_fet:	0.332
+seg	2684	2686	single:	9	0.778	0.493,1	single:	10	0.5	0.234,0.766	pair:	-0.278	-0	-0.529,0.0939	p_acc:	0.0198	p_fet:	0.35
+seg	2741	2743	single:	NA	NA	NA,NA	single:	13	0.308	0.0724,0.543	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2781	2783	single:	NA	NA	NA,NA	single:	10	0.2	0,0.47	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	2976	2978	single:	92	0.88	0.811,0.95	single:	49	0.347	0.217,0.477	pair:	-0.533	-0.393	-0.639,-0.393	p_acc:	7.62e-11	p_fet:	2.14e-10
+seg	2986	2988	single:	122	0.869	0.807,0.93	single:	61	0.377	0.258,0.496	pair:	-0.492	-0.368	-0.592,-0.368	p_acc:	1.41e-11	p_fet:	1.77e-11
+seg	3042	3044	single:	5	1	0.607,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3063	3065	single:	8	0.75	0.45,1	single:	6	0	0,0.348	pair:	-0.75	-0.257	-0.833,-0.257	p_acc:	0.000785	p_fet:	0.00966
+seg	3072	3074	single:	3	0.667	0.249,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3207	3209	single:	3	1	0.473,1	single:	6	0.833	0.479,1	pair:	-0.167	-0	-0.397,0.333	p_acc:	0.0373	p_fet:	1
+seg	3359	3361	single:	NA	NA	NA,NA	single:	9	0.333	0.0575,0.609	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3405	3407	single:	4	0.75	0.343,1	single:	8	0.25	0,0.55	pair:	-0.5	-0	-0.711,0.0288	p_acc:	0.00979	p_fet:	0.222
+seg	3411	3413	single:	4	0.75	0.343,1	single:	9	0.111	0,0.394	pair:	-0.639	-0.112	-0.799,-0.112	p_acc:	0.0036	p_fet:	0.0517
+seg	3427	3429	single:	4	1	0.549,1	single:	3	0.667	0.249,1	pair:	-0.333	-0	-0.639,0.17	p_acc:	0.0214	p_fet:	0.429
+seg	3429	3431	single:	5	0.8	0.418,1	single:	4	0.5	0.147,0.853	pair:	-0.3	-0	-0.611,0.205	p_acc:	0.0214	p_fet:	0.524
+seg	3433	3435	single:	5	1	0.607,1	single:	3	0.333	0,0.751	pair:	-0.667	-0.0486	-0.815,-0.0486	p_acc:	0.00572	p_fet:	0.107
+seg	3509	3511	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3667	3669	single:	3	1	0.473,1	single:	8	0.125	0,0.429	pair:	-0.875	-0.221	-0.88,-0.221	p_acc:	0.00146	p_fet:	0.0242
+seg	3676	3678	single:	NA	NA	NA,NA	single:	7	0.429	0.128,0.73	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3724	3726	single:	NA	NA	NA,NA	single:	9	0.222	0,0.507	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3754	3756	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3799	3801	single:	NA	NA	NA,NA	single:	3	0.333	0,0.751	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3801	3803	single:	NA	NA	NA,NA	single:	3	0.333	0,0.751	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3805	3807	single:	NA	NA	NA,NA	single:	3	0.333	0,0.751	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3809	3811	single:	NA	NA	NA,NA	single:	5	0.2	0,0.582	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3811	3813	single:	NA	NA	NA,NA	single:	4	0.25	0,0.657	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3814	3816	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3824	3826	single:	NA	NA	NA,NA	single:	4	0	0,0.451	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3839	3841	single:	NA	NA	NA,NA	single:	4	0	0,0.451	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3843	3845	single:	NA	NA	NA,NA	single:	5	0.2	0,0.582	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3849	3851	single:	NA	NA	NA,NA	single:	4	0.25	0,0.657	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3851	3853	single:	NA	NA	NA,NA	single:	5	0.2	0,0.582	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3861	3863	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3865	3867	single:	NA	NA	NA,NA	single:	4	0	0,0.451	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3867	3869	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	3882	3884	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	4025	4027	single:	NA	NA	NA,NA	single:	8	0.125	0,0.429	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	4037	4039	single:	3	1	0.473,1	single:	8	0.625	0.337,0.913	pair:	-0.375	-0	-0.539,0.191	p_acc:	0.0204	p_fet:	0.491
+seg	4046	4048	single:	3	1	0.473,1	single:	7	0.143	0,0.471	pair:	-0.857	-0.189	-0.87,-0.189	p_acc:	0.00194	p_fet:	0.0333
+seg	4057	4059	single:	4	1	0.549,1	single:	7	0.429	0.128,0.73	pair:	-0.571	-0.0173	-0.71,-0.0173	p_acc:	0.00708	p_fet:	0.194
+seg	4060	4062	single:	4	1	0.549,1	single:	6	0.333	0.0079,0.659	pair:	-0.667	-0.0776	-0.778,-0.0776	p_acc:	0.00455	p_fet:	0.0762
+seg	4068	4070	single:	NA	NA	NA,NA	single:	8	1	0.717,1	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	4087	4089	single:	NA	NA	NA,NA	single:	7	0	0,0.312	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	4094	4096	single:	NA	NA	NA,NA	single:	4	0	0,0.451	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	4104	4106	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	4113	4115	single:	NA	NA	NA,NA	single:	4	0	0,0.451	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	4277	4279	single:	5	0.8	0.418,1	single:	3	0	0,0.527	pair:	-0.8	-0.0985	-0.843,-0.0985	p_acc:	0.00381	p_fet:	0.143
+seg	4280	4282	single:	5	0.6	0.264,0.936	single:	3	0.333	0,0.751	pair:	-0.267	-0	-0.594,0.282	p_acc:	0.0229	p_fet:	1
+seg	4283	4285	single:	5	0.6	0.264,0.936	single:	4	0	0,0.451	pair:	-0.6	-0.0122	-0.748,-0.0122	p_acc:	0.00715	p_fet:	0.167
+seg	4320	4322	single:	4	0.75	0.343,1	single:	4	0	0,0.451	pair:	-0.75	-0.0938	-0.834,-0.0938	p_acc:	0.00397	p_fet:	0.143
+seg	4327	4329	single:	4	1	0.549,1	single:	3	0	0,0.527	pair:	-1	-0.229	-1,-0.229	p_acc:	0.00143	p_fet:	0.0286
+seg	4333	4335	single:	4	1	0.549,1	single:	3	0	0,0.527	pair:	-1	-0.229	-1,-0.229	p_acc:	0.00143	p_fet:	0.0286
+seg	4340	4342	single:	3	1	0.473,1	single:	3	0	0,0.527	pair:	-1	-0.17	-1,-0.17	p_acc:	0.00229	p_fet:	0.1
+seg	4505	4507	single:	5	0.4	0.064,0.736	single:	16	0	0,0.162	pair:	-0.4	-0.0844	-0.683,-0.0844	p_acc:	0.00444	p_fet:	0.0476
+seg	4511	4513	single:	6	0.667	0.341,0.992	single:	18	0.167	0,0.359	pair:	-0.5	-0.107	-0.709,-0.107	p_acc:	0.00377	p_fet:	0.0381
+seg	4516	4518	single:	6	0.5	0.184,0.816	single:	21	0.0476	0,0.198	pair:	-0.452	-0.122	-0.702,-0.122	p_acc:	0.00264	p_fet:	0.0248
+seg	4550	4552	single:	9	0.667	0.391,0.942	single:	23	0.13	0,0.292	pair:	-0.536	-0.205	-0.719,-0.205	p_acc:	0.000775	p_fet:	0.00564
+seg	4605	4607	single:	8	0.75	0.45,1	single:	26	0	0,0.105	pair:	-0.75	-0.408	-0.874,-0.408	p_acc:	2.94e-06	p_fet:	2.08e-05
+seg	4648	4650	single:	15	0.667	0.445,0.889	single:	39	0.179	0.0544,0.305	pair:	-0.487	-0.232	-0.653,-0.232	p_acc:	0.00023	p_fet:	0.00108
+seg	4738	4740	single:	26	0.577	0.398,0.756	single:	53	0.151	0.0496,0.252	pair:	-0.426	-0.232	-0.576,-0.232	p_acc:	5.23e-05	p_fet:	0.000172
+seg	4797	4799	single:	15	0.6	0.374,0.826	single:	45	0.0889	0,0.188	pair:	-0.511	-0.271	-0.68,-0.271	p_acc:	3.55e-05	p_fet:	0.000153
+seg	4829	4831	single:	12	0.75	0.505,0.995	single:	38	0.105	0,0.22	pair:	-0.645	-0.364	-0.782,-0.364	p_acc:	9.24e-06	p_fet:	4.74e-05
+seg	4846	4848	single:	14	0.929	0.721,1	single:	38	0.105	0,0.22	pair:	-0.823	-0.574	-0.884,-0.574	p_acc:	1.06e-08	p_fet:	4.75e-08
+seg	4922	4924	single:	13	0.538	0.298,0.779	single:	33	0.273	0.123,0.422	pair:	-0.266	-0.00424	-0.484,-0.00424	p_acc:	0.0135	p_fet:	0.167
+seg	4925	4927	single:	13	0.615	0.377,0.854	single:	33	0.303	0.15,0.456	pair:	-0.312	-0.0403	-0.517,-0.0403	p_acc:	0.00858	p_fet:	0.0916
+seg	4932	4934	single:	15	0.8	0.583,1	single:	29	0.276	0.116,0.435	pair:	-0.524	-0.251	-0.671,-0.251	p_acc:	0.000238	p_fet:	0.00138
+seg	4949	4951	single:	16	0.625	0.407,0.843	single:	23	0.13	0,0.292	pair:	-0.495	-0.224	-0.661,-0.224	p_acc:	0.000359	p_fet:	0.00203
+seg	4982	4984	single:	15	0.6	0.374,0.826	single:	11	0.455	0.198,0.711	pair:	-0.145	-0	-0.414,0.169	p_acc:	0.034	p_fet:	0.692
+seg	4984	4986	single:	14	0.571	0.338,0.804	single:	11	0.455	0.198,0.711	pair:	-0.117	-0	-0.393,0.198	p_acc:	0.0369	p_fet:	0.695
+seg	5011	5013	single:	16	0.375	0.157,0.593	single:	12	0.25	0.00535,0.495	pair:	-0.125	-0	-0.364,0.168	p_acc:	0.0389	p_fet:	0.687
+seg	5013	5015	single:	16	0.375	0.157,0.593	single:	11	0.455	0.198,0.711	pair:	0.0795	0	-0.215,0.359	p_acc:	0.0413	p_fet:	0.71
+seg	5049	5051	single:	6	0.333	0.0079,0.659	single:	8	0.125	0,0.429	pair:	-0.208	-0	-0.509,0.156	p_acc:	0.0277	p_fet:	0.538
+seg	5057	5059	single:	4	0.75	0.343,1	single:	7	0	0,0.312	pair:	-0.75	-0.191	-0.856,-0.191	p_acc:	0.00162	p_fet:	0.0242
+seg	5068	5070	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	5070	5072	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	5128	5130	single:	3	1	0.473,1	single:	9	0	0,0.259	pair:	-1	-0.356	-1,-0.356	p_acc:	0.00028	p_fet:	0.00455
+seg	5160	5162	single:	4	1	0.549,1	single:	11	0	0,0.221	pair:	-1	-0.45	-1,-0.45	p_acc:	5.52e-05	p_fet:	0.000733
+seg	5162	5164	single:	4	1	0.549,1	single:	7	0.143	0,0.471	pair:	-0.857	-0.255	-0.88,-0.255	p_acc:	0.00101	p_fet:	0.0152
+seg	5182	5184	single:	3	1	0.473,1	single:	5	0	0,0.393	pair:	-1	-0.271	-1,-0.271	p_acc:	0.000954	p_fet:	0.0179
+seg	5192	5194	single:	3	1	0.473,1	single:	6	0.167	0,0.521	pair:	-0.833	-0.149	-0.858,-0.149	p_acc:	0.00267	p_fet:	0.0476
+seg	5201	5203	single:	NA	NA	NA,NA	single:	6	0	0,0.348	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	5333	5335	single:	14	0.857	0.637,1	single:	21	0.143	0,0.316	pair:	-0.714	-0.42	-0.819,-0.42	p_acc:	6.9e-06	p_fet:	5.49e-05
+seg	5385	5387	single:	9	0.778	0.493,1	single:	14	0.357	0.127,0.587	pair:	-0.421	-0.0503	-0.623,-0.0503	p_acc:	0.00667	p_fet:	0.0894
+seg	5537	5539	single:	5	1	0.607,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	5587	5589	single:	8	1	0.717,1	single:	18	0.778	0.581,0.975	pair:	-0.222	-0	-0.356,0.0752	p_acc:	0.026	p_fet:	0.277
+seg	5857	5859	single:	4	0.75	0.343,1	single:	5	0.2	0,0.582	pair:	-0.55	-0	-0.745,0.0401	p_acc:	0.00953	p_fet:	0.206
+seg	5869	5871	single:	4	0.75	0.343,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	5886	5888	single:	3	0.667	0.249,1	single:	4	0.25	0,0.657	pair:	-0.417	-0	-0.685,0.196	p_acc:	0.0171	p_fet:	0.486
+seg	5896	5898	single:	3	0.667	0.249,1	single:	4	1	0.549,1	pair:	0.333	0	-0.17,0.639	p_acc:	0.0214	p_fet:	0.429
+seg	6182	6184	single:	4	0.5	0.147,0.853	single:	14	0.571	0.338,0.804	pair:	0.0714	0	-0.307,0.43	p_acc:	0.0325	p_fet:	1
+seg	6293	6295	single:	4	1	0.549,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6304	6306	single:	6	0.5	0.184,0.816	single:	5	0	0,0.393	pair:	-0.5	-0.00245	-0.684,-0.00245	p_acc:	0.00849	p_fet:	0.182
+seg	6332	6334	single:	8	0.75	0.45,1	single:	6	0.167	0,0.521	pair:	-0.583	-0.096	-0.751,-0.096	p_acc:	0.00412	p_fet:	0.103
+seg	6403	6405	single:	7	0.857	0.529,1	single:	3	0	0,0.527	pair:	-0.857	-0.189	-0.87,-0.189	p_acc:	0.00194	p_fet:	0.0333
+seg	6411	6413	single:	6	1	0.652,1	single:	4	0.75	0.343,1	pair:	-0.25	-0	-0.567,0.123	p_acc:	0.0255	p_fet:	0.4
+seg	6420	6422	single:	4	1	0.549,1	single:	3	0.667	0.249,1	pair:	-0.333	-0	-0.639,0.17	p_acc:	0.0214	p_fet:	0.429
+seg	6487	6489	single:	NA	NA	NA,NA	single:	5	0	0,0.393	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6499	6501	single:	NA	NA	NA,NA	single:	5	0	0,0.393	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6533	6535	single:	3	0.667	0.249,1	single:	14	0.0714	0,0.279	pair:	-0.595	-0.101	-0.805,-0.101	p_acc:	0.00412	p_fet:	0.0632
+seg	6544	6546	single:	4	0.75	0.343,1	single:	15	0	0,0.171	pair:	-0.75	-0.272	-0.881,-0.272	p_acc:	0.000332	p_fet:	0.00413
+seg	6578	6580	single:	NA	NA	NA,NA	single:	7	0	0,0.312	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6630	6632	single:	5	1	0.607,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6633	6635	single:	5	1	0.607,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6640	6642	single:	7	1	0.688,1	single:	3	0.667	0.249,1	pair:	-0.333	-0	-0.669,0.0603	p_acc:	0.0175	p_fet:	0.3
+seg	6654	6656	single:	8	0.875	0.571,1	single:	4	0.25	0,0.657	pair:	-0.625	-0.0857	-0.789,-0.0857	p_acc:	0.00448	p_fet:	0.0667
+seg	6671	6673	single:	9	0.778	0.493,1	single:	9	0.222	0,0.507	pair:	-0.556	-0.135	-0.731,-0.135	p_acc:	0.00281	p_fet:	0.0567
+seg	6732	6734	single:	3	1	0.473,1	single:	9	0.111	0,0.394	pair:	-0.889	-0.247	-0.889,-0.247	p_acc:	0.00112	p_fet:	0.0182
+seg	6774	6776	single:	6	0.833	0.479,1	single:	7	0	0,0.312	pair:	-0.833	-0.318	-0.887,-0.318	p_acc:	0.000374	p_fet:	0.00466
+seg	6789	6791	single:	6	1	0.652,1	single:	5	0	0,0.393	pair:	-1	-0.417	-1,-0.417	p_acc:	0.000152	p_fet:	0.00216
+seg	6794	6796	single:	9	1	0.741,1	single:	7	0	0,0.312	pair:	-1	-0.549	-1,-0.549	p_acc:	8.26e-06	p_fet:	8.74e-05
+seg	6839	6841	single:	3	1	0.473,1	single:	4	0.25	0,0.657	pair:	-0.75	-0.0327	-0.822,-0.0327	p_acc:	0.00572	p_fet:	0.143
+seg	6845	6847	single:	4	1	0.549,1	single:	5	0.2	0,0.582	pair:	-0.8	-0.161	-0.854,-0.161	p_acc:	0.00238	p_fet:	0.0476
+seg	6855	6857	single:	4	1	0.549,1	single:	5	0.2	0,0.582	pair:	-0.8	-0.161	-0.854,-0.161	p_acc:	0.00238	p_fet:	0.0476
+seg	6861	6863	single:	6	1	0.652,1	single:	4	0.25	0,0.657	pair:	-0.75	-0.168	-0.85,-0.168	p_acc:	0.00212	p_fet:	0.0333
+seg	6901	6903	single:	6	0.833	0.479,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6909	6911	single:	6	1	0.652,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6927	6929	single:	6	1	0.652,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6940	6942	single:	3	1	0.473,1	single:	3	0	0,0.527	pair:	-1	-0.17	-1,-0.17	p_acc:	0.00229	p_fet:	0.1
+seg	6942	6944	single:	4	1	0.549,1	single:	4	0	0,0.451	pair:	-1	-0.291	-1,-0.291	p_acc:	0.000795	p_fet:	0.0286
+seg	6950	6952	single:	4	1	0.549,1	single:	4	0	0,0.451	pair:	-1	-0.291	-1,-0.291	p_acc:	0.000795	p_fet:	0.0286
+seg	6954	6956	single:	NA	NA	NA,NA	single:	5	0	0,0.393	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	6983	6985	single:	3	1	0.473,1	single:	5	0	0,0.393	pair:	-1	-0.271	-1,-0.271	p_acc:	0.000954	p_fet:	0.0179
+seg	7006	7008	single:	3	1	0.473,1	single:	4	0	0,0.451	pair:	-1	-0.229	-1,-0.229	p_acc:	0.00143	p_fet:	0.0286
+seg	7188	7190	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	7204	7206	single:	NA	NA	NA,NA	single:	3	0	0,0.527	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	7546	7548	single:	NA	NA	NA,NA	single:	7	0.857	0.529,1	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	7805	7807	single:	NA	NA	NA,NA	single:	8	0.125	0,0.429	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	7987	7989	single:	NA	NA	NA,NA	single:	13	0	0,0.193	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+seg	8052	8054	single:	4	1	0.549,1	single:	9	0.111	0,0.394	pair:	-0.889	-0.315	-0.896,-0.315	p_acc:	0.0005	p_fet:	0.00699
+seg	8058	8060	single:	4	1	0.549,1	single:	6	0.333	0.0079,0.659	pair:	-0.667	-0.0776	-0.778,-0.0776	p_acc:	0.00455	p_fet:	0.0762
+seg	8089	8091	single:	3	0.667	0.249,1	single:	3	0	0,0.527	pair:	-0.667	-0	-0.791,0.054	p_acc:	0.00915	p_fet:	0.4
+seg	8094	8096	single:	3	1	0.473,1	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/PE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,3 @@
+seg	152	1305	-	20	15	5	1
+seg	2310	7008	-	95	77	18	1
+seg	8052	8091	-	3	3	0	1
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/SE_comp.g1.vs.g2.txt	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,310 @@
+#chrom	start	end	single:	totalC_0	nominalRatio_0	ratioCI_0	single:	totalC_1	nominalRatio_1	ratioCI_1	pair:	nominalDif_1-0	credibleDif_1-0	difCI_1-0	p_sim:	p_sim_1_v_0	p_fet:	p_fet_1_v_0
+chr11	41	43	single:	3	0	0,0.527	single:	3	0	0,0.527	pair:	0	0	-0.391,0.391	p_acc:	0.0449	p_fet:	1
+chr11	48	50	single:	5	0	0,0.393	single:	3	0	0,0.527	pair:	0	0	-0.267,0.422	p_acc:	0.0521	p_fet:	1
+chr11	56	58	single:	6	0	0,0.348	single:	4	0	0,0.451	pair:	0	0	-0.241,0.359	p_acc:	0.0619	p_fet:	1
+chr11	59	61	single:	8	0	0,0.283	single:	5	0.4	0.064,0.736	pair:	0.4	0.0115	0.0115,0.656	p_acc:	0.0099	p_fet:	0.128
+chr11	101	103	single:	13	0.0769	0,0.297	single:	7	0.143	0,0.471	pair:	0.0659	0	-0.152,0.364	p_acc:	0.051	p_fet:	1
+chr11	140	142	single:	16	0	0,0.162	single:	9	0	0,0.259	pair:	0	0	-0.107,0.216	p_acc:	0.123	p_fet:	1
+chr11	144	146	single:	14	0	0,0.181	single:	10	0.1	0,0.364	pair:	0.1	0	-0.0722,0.317	p_acc:	0.0549	p_fet:	0.417
+chr11	210	212	single:	15	0.133	0,0.344	single:	12	0.417	0.169,0.664	pair:	0.283	0	-0.00696,0.507	p_acc:	0.0139	p_fet:	0.185
+chr11	215	217	single:	19	0.263	0.0688,0.458	single:	16	0.25	0.039,0.461	pair:	-0.0132	-0	-0.238,0.226	p_acc:	0.0566	p_fet:	1
+chr11	316	318	single:	24	0.0833	0,0.231	single:	26	0.846	0.692,1	pair:	0.763	0.538	0.538,0.845	p_acc:	9.21e-09	p_fet:	3.52e-08
+chr11	331	333	single:	27	0.037	0,0.159	single:	26	0.808	0.649,0.966	pair:	0.771	0.558	0.558,0.85	p_acc:	1.11e-09	p_fet:	4.05e-09
+chr11	400	402	single:	34	0.0588	0,0.169	single:	34	0.412	0.254,0.57	pair:	0.353	0.18	0.18,0.486	p_acc:	0.000192	p_fet:	0.00115
+chr11	407	409	single:	33	0	0,0.0843	single:	34	0.588	0.43,0.746	pair:	0.588	0.41	0.41,0.692	p_acc:	8.72e-09	p_fet:	3.39e-08
+chr11	445	447	single:	25	0.08	0,0.223	single:	33	0.606	0.447,0.766	pair:	0.526	0.316	0.316,0.648	p_acc:	9.32e-06	p_fet:	3.59e-05
+chr11	469	471	single:	26	0.0385	0,0.164	single:	34	0.588	0.43,0.746	pair:	0.55	0.352	0.352,0.661	p_acc:	1.43e-06	p_fet:	5.87e-06
+chr11	546	548	single:	20	0.1	0,0.271	single:	32	0.812	0.672,0.953	pair:	0.713	0.482	0.482,0.806	p_acc:	1.08e-07	p_fet:	4.38e-07
+chr11	666	668	single:	20	0.05	0,0.207	single:	38	0.553	0.401,0.704	pair:	0.503	0.29	0.29,0.612	p_acc:	2.88e-05	p_fet:	0.000144
+chr11	694	696	single:	22	0.318	0.132,0.504	single:	41	0.902	0.795,1	pair:	0.584	0.368	0.368,0.716	p_acc:	8.72e-07	p_fet:	2.95e-06
+chr11	712	714	single:	28	0.0714	0,0.202	single:	42	0.881	0.771,0.991	pair:	0.81	0.63	0.63,0.873	p_acc:	1.66e-12	p_fet:	4.55e-12
+chr11	721	723	single:	28	0.0357	0,0.153	single:	39	0.923	0.817,1	pair:	0.887	0.719	0.719,0.926	p_acc:	9.44e-15	p_fet:	2.57e-14
+chr11	737	739	single:	27	0.037	0,0.159	single:	37	0.622	0.471,0.772	pair:	0.585	0.393	0.393,0.687	p_acc:	2.21e-07	p_fet:	9.15e-07
+chr11	753	755	single:	23	0.087	0,0.24	single:	41	0.78	0.653,0.908	pair:	0.694	0.488	0.488,0.783	p_acc:	1.74e-08	p_fet:	5.8e-08
+chr11	762	764	single:	27	0.037	0,0.159	single:	42	0.619	0.477,0.761	pair:	0.582	0.397	0.397,0.68	p_acc:	1.42e-07	p_fet:	5.56e-07
+chr11	792	794	single:	37	0.324	0.178,0.471	single:	38	0.842	0.719,0.966	pair:	0.518	0.329	0.329,0.642	p_acc:	1.8e-06	p_fet:	6.02e-06
+chr11	798	800	single:	36	0.306	0.159,0.452	single:	40	0.875	0.761,0.989	pair:	0.569	0.384	0.384,0.685	p_acc:	1.25e-07	p_fet:	3.94e-07
+chr11	816	818	single:	33	0.0303	0,0.132	single:	34	0.912	0.793,1	pair:	0.881	0.714	0.714,0.924	p_acc:	5.66e-15	p_fet:	1.61e-14
+chr11	830	832	single:	37	0.405	0.254,0.557	single:	35	0.914	0.798,1	pair:	0.509	0.324	0.324,0.629	p_acc:	1.55e-06	p_fet:	5.05e-06
+chr11	954	956	single:	47	0	0,0.0605	single:	33	0.485	0.323,0.647	pair:	0.485	0.324	0.324,0.606	p_acc:	1.85e-08	p_fet:	4.33e-08
+chr11	977	979	single:	47	0.255	0.132,0.379	single:	37	0.703	0.559,0.847	pair:	0.447	0.264	0.264,0.579	p_acc:	1.7e-05	p_fet:	5.98e-05
+chr11	1006	1008	single:	38	0	0,0.0739	single:	38	0.526	0.374,0.678	pair:	0.526	0.364	0.364,0.633	p_acc:	1.27e-08	p_fet:	6.16e-08
+chr11	1066	1068	single:	30	0.367	0.202,0.532	single:	45	0.756	0.63,0.881	pair:	0.389	0.192	0.192,0.538	p_acc:	0.00027	p_fet:	0.00152
+chr11	1116	1118	single:	22	0.182	0.0087,0.355	single:	35	0.771	0.632,0.911	pair:	0.59	0.362	0.362,0.713	p_acc:	4.09e-06	p_fet:	2.32e-05
+chr11	1145	1147	single:	18	0.556	0.345,0.766	single:	26	0.538	0.358,0.719	pair:	-0.0171	-0	-0.248,0.222	p_acc:	0.0547	p_fet:	1
+chr11	1260	1262	single:	18	0.0556	0,0.226	single:	24	0.375	0.192,0.558	pair:	0.319	0.0943	0.0943,0.47	p_acc:	0.00355	p_fet:	0.0264
+chr11	1320	1322	single:	17	0.118	0,0.31	single:	25	0.64	0.461,0.819	pair:	0.522	0.26	0.26,0.661	p_acc:	0.000172	p_fet:	0.00122
+chr11	1559	1561	single:	36	0.194	0.0615,0.327	single:	27	0.852	0.702,1	pair:	0.657	0.451	0.451,0.762	p_acc:	5.66e-08	p_fet:	1.9e-07
+chr11	1574	1576	single:	33	0.152	0.0184,0.285	single:	28	0.607	0.435,0.779	pair:	0.456	0.245	0.245,0.601	p_acc:	7.43e-05	p_fet:	0.000388
+chr11	1592	1594	single:	35	0.343	0.191,0.495	single:	33	0.697	0.544,0.85	pair:	0.354	0.149	0.149,0.509	p_acc:	0.001	p_fet:	0.00421
+chr11	1609	1611	single:	34	0.0588	0,0.169	single:	41	0.244	0.113,0.375	pair:	0.185	0.0419	0.0419,0.305	p_acc:	0.0094	p_fet:	0.0545
+chr11	1632	1634	single:	36	0.0556	0,0.161	single:	45	0.2	0.0811,0.319	pair:	0.144	0.0144	0.0144,0.255	p_acc:	0.0193	p_fet:	0.1
+chr11	1642	1644	single:	32	0.188	0.0468,0.328	single:	41	0.268	0.134,0.402	pair:	0.0808	0	-0.0869,0.229	p_acc:	0.0606	p_fet:	0.578
+chr11	1680	1682	single:	34	0.0882	0,0.207	single:	54	0.611	0.485,0.738	pair:	0.523	0.354	0.354,0.626	p_acc:	2.25e-07	p_fet:	6.09e-07
+chr11	1689	1691	single:	35	0.2	0.0642,0.336	single:	51	0.471	0.338,0.603	pair:	0.271	0.095	0.095,0.408	p_acc:	0.00289	p_fet:	0.0124
+chr11	1703	1705	single:	42	0.262	0.13,0.393	single:	49	0.755	0.635,0.875	pair:	0.493	0.319	0.319,0.614	p_acc:	1.06e-06	p_fet:	2.93e-06
+chr11	1733	1735	single:	43	0.512	0.368,0.655	single:	51	0.745	0.626,0.864	pair:	0.233	0.0666	0.0666,0.379	p_acc:	0.00549	p_fet:	0.0305
+chr11	1766	1768	single:	46	0.13	0.0243,0.237	single:	48	0.312	0.184,0.441	pair:	0.182	0.0381	0.0381,0.31	p_acc:	0.0105	p_fet:	0.0473
+chr11	1869	1871	single:	45	0.244	0.119,0.37	single:	38	0.263	0.125,0.402	pair:	0.0187	0	-0.134,0.175	p_acc:	0.0834	p_fet:	1
+chr11	2040	2042	single:	37	0	0,0.0758	single:	34	0	0,0.082	pair:	0	0	-0.057,0.0639	p_acc:	0.31	p_fet:	1
+chr11	2104	2106	single:	37	0.027	0,0.119	single:	39	0.0513	0,0.149	pair:	0.0243	0	-0.0641,0.111	p_acc:	0.153	p_fet:	1
+chr11	2148	2150	single:	31	0.613	0.449,0.777	single:	34	0.588	0.43,0.746	pair:	-0.0247	-0	-0.214,0.17	p_acc:	0.0664	p_fet:	1
+chr11	2360	2362	single:	42	0.262	0.13,0.393	single:	26	0.538	0.358,0.719	pair:	0.277	0.0726	0.0726,0.448	p_acc:	0.00528	p_fet:	0.0373
+chr11	2384	2386	single:	34	0.382	0.226,0.539	single:	27	0.593	0.417,0.768	pair:	0.21	0	-0.00431,0.392	p_acc:	0.0176	p_fet:	0.126
+chr11	2454	2456	single:	41	0.171	0.0504,0.291	single:	28	0.107	0,0.246	pair:	-0.0636	-0	-0.189,0.0894	p_acc:	0.0753	p_fet:	0.729
+chr11	2621	2623	single:	9	0	0,0.259	single:	12	0	0,0.206	pair:	0	0	-0.206,0.147	p_acc:	0.112	p_fet:	1
+chr11	2648	2650	single:	10	0	0,0.238	single:	20	0	0,0.133	pair:	0	0	-0.202,0.0851	p_acc:	0.138	p_fet:	1
+chr11	2665	2667	single:	13	0	0,0.193	single:	19	0.105	0,0.283	pair:	0.105	0	-0.0832,0.237	p_acc:	0.0587	p_fet:	0.502
+chr11	2667	2669	single:	17	0	0,0.153	single:	20	0.1	0,0.271	pair:	0.1	0	-0.0533,0.233	p_acc:	0.057	p_fet:	0.489
+chr11	2673	2675	single:	13	0	0,0.193	single:	21	0.19	0.0119,0.369	pair:	0.19	0	-0.0203,0.323	p_acc:	0.0228	p_fet:	0.144
+chr11	2676	2678	single:	17	0	0,0.153	single:	26	0.0385	0,0.164	pair:	0.0385	0	-0.0986,0.129	p_acc:	0.132	p_fet:	1
+chr11	2685	2687	single:	16	0.0625	0,0.25	single:	28	0.143	0,0.288	pair:	0.0804	0	-0.113,0.211	p_acc:	0.0661	p_fet:	0.638
+chr11	2688	2690	single:	19	0	0,0.139	single:	28	0	0,0.0981	pair:	0	0	-0.114,0.0684	p_acc:	0.215	p_fet:	1
+chr11	2691	2693	single:	20	0	0,0.133	single:	29	0.276	0.116,0.435	pair:	0.276	0.0981	0.0981,0.397	p_acc:	0.00242	p_fet:	0.0148
+chr11	2713	2715	single:	22	0	0,0.122	single:	26	0.154	6.6e-05,0.308	pair:	0.154	0.0107	0.0107,0.277	p_acc:	0.0197	p_fet:	0.114
+chr11	2723	2725	single:	19	0	0,0.139	single:	25	0.12	0,0.272	pair:	0.12	0	-0.0272,0.238	p_acc:	0.0404	p_fet:	0.247
+chr11	2727	2729	single:	20	0	0,0.133	single:	20	0.35	0.153,0.547	pair:	0.35	0.143	0.143,0.501	p_acc:	0.000903	p_fet:	0.00832
+chr11	2731	2733	single:	17	0.0588	0,0.238	single:	19	0.211	0.02,0.401	pair:	0.152	0	-0.0525,0.322	p_acc:	0.0341	p_fet:	0.342
+chr11	2738	2740	single:	18	0	0,0.146	single:	22	0.273	0.0909,0.455	pair:	0.273	0.0783	0.0783,0.414	p_acc:	0.00417	p_fet:	0.0243
+chr11	2767	2769	single:	15	0	0,0.171	single:	21	0.333	0.142,0.525	pair:	0.333	0.107	0.107,0.473	p_acc:	0.00267	p_fet:	0.0273
+chr11	2778	2780	single:	12	0.0833	0,0.316	single:	22	0.455	0.261,0.648	pair:	0.371	0.0842	0.0842,0.526	p_acc:	0.00465	p_fet:	0.0527
+chr11	2788	2790	single:	8	0	0,0.283	single:	18	0.167	0,0.359	pair:	0.167	0	-0.118,0.297	p_acc:	0.0398	p_fet:	0.529
+chr11	2792	2794	single:	7	0.143	0,0.471	single:	22	0.318	0.132,0.504	pair:	0.175	0	-0.175,0.358	p_acc:	0.0341	p_fet:	0.635
+chr11	2810	2812	single:	3	0	0,0.527	single:	19	0.105	0,0.283	pair:	0.105	0	-0.4,0.187	p_acc:	0.0515	p_fet:	1
+chr11	2825	2827	single:	4	0	0,0.451	single:	15	0.133	0,0.344	pair:	0.133	0	-0.3,0.253	p_acc:	0.0491	p_fet:	1
+chr11	3031	3033	single:	13	0.0769	0,0.297	single:	15	0.4	0.174,0.626	pair:	0.323	0.038	0.038,0.511	p_acc:	0.00851	p_fet:	0.0836
+chr11	3034	3036	single:	13	0	0,0.193	single:	14	0.286	0.059,0.512	pair:	0.286	0.0412	0.0412,0.462	p_acc:	0.00859	p_fet:	0.0978
+chr11	3052	3054	single:	15	0.0667	0,0.264	single:	13	0.769	0.534,1	pair:	0.703	0.377	0.377,0.816	p_acc:	3.11e-05	p_fet:	0.000267
+chr11	3095	3097	single:	11	0.273	0.0181,0.527	single:	10	0.8	0.53,1	pair:	0.527	0.141	0.141,0.706	p_acc:	0.00253	p_fet:	0.03
+chr11	3117	3119	single:	12	0.0833	0,0.316	single:	14	0.571	0.338,0.804	pair:	0.488	0.161	0.161,0.655	p_acc:	0.00167	p_fet:	0.0145
+chr11	3273	3275	single:	38	0.105	0,0.22	single:	25	0.48	0.296,0.664	pair:	0.375	0.178	0.178,0.532	p_acc:	0.000337	p_fet:	0.00124
+chr11	3283	3285	single:	37	0.027	0,0.119	single:	23	0.478	0.288,0.669	pair:	0.451	0.257	0.257,0.6	p_acc:	1.1e-05	p_fet:	3.67e-05
+chr11	3323	3325	single:	35	0.0571	0,0.165	single:	24	0.625	0.442,0.808	pair:	0.568	0.358	0.358,0.698	p_acc:	8.63e-07	p_fet:	2.9e-06
+chr11	3424	3426	single:	30	0.0667	0,0.189	single:	21	0.0952	0,0.259	pair:	0.0286	0	-0.0969,0.184	p_acc:	0.0956	p_fet:	1
+chr11	3452	3454	single:	27	0	0,0.101	single:	17	0.588	0.374,0.803	pair:	0.588	0.349	0.349,0.729	p_acc:	1.79e-06	p_fet:	7.84e-06
+chr11	3484	3486	single:	25	0.04	0,0.17	single:	18	0.556	0.345,0.766	pair:	0.516	0.275	0.275,0.667	p_acc:	4.33e-05	p_fet:	0.000196
+chr11	3578	3580	single:	19	0.263	0.0688,0.458	single:	19	0.737	0.542,0.931	pair:	0.474	0.194	0.194,0.641	p_acc:	0.000789	p_fet:	0.00856
+chr11	3595	3597	single:	20	0	0,0.133	single:	21	0.238	0.0554,0.421	pair:	0.238	0.0598	0.0598,0.385	p_acc:	0.00603	p_fet:	0.0478
+chr11	3611	3613	single:	19	0.0526	0,0.216	single:	23	0.826	0.658,0.994	pair:	0.773	0.525	0.525,0.852	p_acc:	7.44e-08	p_fet:	3.4e-07
+chr11	3617	3619	single:	19	0	0,0.139	single:	21	0.762	0.579,0.945	pair:	0.762	0.514	0.514,0.842	p_acc:	7.07e-08	p_fet:	3.39e-07
+chr11	3639	3641	single:	16	0.125	0,0.326	single:	20	0.6	0.4,0.8	pair:	0.475	0.194	0.194,0.633	p_acc:	0.000772	p_fet:	0.00572
+chr11	3646	3648	single:	17	0.0588	0,0.238	single:	23	0.391	0.204,0.579	pair:	0.332	0.0973	0.0973,0.485	p_acc:	0.00347	p_fet:	0.0257
+chr11	3665	3667	single:	13	0.0769	0,0.297	single:	21	0.571	0.375,0.768	pair:	0.495	0.203	0.203,0.637	p_acc:	0.000728	p_fet:	0.00476
+chr11	3672	3674	single:	13	0	0,0.193	single:	23	0.348	0.163,0.532	pair:	0.348	0.106	0.106,0.476	p_acc:	0.00296	p_fet:	0.0319
+chr11	3677	3679	single:	14	0	0,0.181	single:	23	0.261	0.0841,0.438	pair:	0.261	0.0448	0.0448,0.391	p_acc:	0.00827	p_fet:	0.0645
+chr11	3693	3695	single:	17	0	0,0.153	single:	26	0.231	0.0678,0.394	pair:	0.231	0.0451	0.0451,0.355	p_acc:	0.0084	p_fet:	0.0662
+chr11	3704	3706	single:	20	0	0,0.133	single:	25	0.2	0.0371,0.363	pair:	0.2	0.0359	0.0359,0.329	p_acc:	0.0104	p_fet:	0.0562
+chr11	3707	3709	single:	21	0	0,0.127	single:	25	0.16	0.00176,0.318	pair:	0.16	0.0108	0.0108,0.286	p_acc:	0.019	p_fet:	0.114
+chr11	3717	3719	single:	22	0.0455	0,0.19	single:	21	0.286	0.0986,0.473	pair:	0.24	0.0455	0.0455,0.403	p_acc:	0.00857	p_fet:	0.0459
+chr11	3960	3962	single:	21	0	0,0.127	single:	18	0.167	0,0.359	pair:	0.167	0.011	0.011,0.326	p_acc:	0.0188	p_fet:	0.0893
+chr11	3991	3993	single:	20	0	0,0.133	single:	24	0.167	0.00373,0.33	pair:	0.167	0.0109	0.0109,0.296	p_acc:	0.0184	p_fet:	0.114
+chr11	4056	4058	single:	22	0	0,0.122	single:	24	0.292	0.115,0.468	pair:	0.292	0.111	0.111,0.431	p_acc:	0.0016	p_fet:	0.00965
+chr11	4163	4165	single:	26	0.0769	0,0.215	single:	30	0.0333	0,0.144	pair:	-0.0436	-0	-0.167,0.067	p_acc:	0.103	p_fet:	0.592
+chr11	4195	4197	single:	23	0.0435	0,0.183	single:	26	0.231	0.0678,0.394	pair:	0.187	0.0142	0.0142,0.331	p_acc:	0.0161	p_fet:	0.103
+chr11	4202	4204	single:	21	0	0,0.127	single:	30	0.0667	0,0.189	pair:	0.0667	0	-0.056,0.159	p_acc:	0.09	p_fet:	0.506
+chr11	4207	4209	single:	18	0.0556	0,0.226	single:	31	0.194	0.0495,0.338	pair:	0.138	0	-0.0495,0.266	p_acc:	0.0376	p_fet:	0.238
+chr11	4243	4245	single:	30	0.133	0,0.271	single:	27	0.111	0,0.254	pair:	-0.0222	-0	-0.164,0.129	p_acc:	0.0905	p_fet:	1
+chr11	4248	4250	single:	32	0.0312	0,0.136	single:	31	0.0323	0,0.14	pair:	0.00101	0	-0.0916,0.0963	p_acc:	0.165	p_fet:	1
+chr11	4284	4286	single:	34	0.0882	0,0.207	single:	31	0.226	0.0775,0.374	pair:	0.138	0	-0.014,0.281	p_acc:	0.0299	p_fet:	0.174
+chr11	4289	4291	single:	33	0	0,0.0843	single:	27	0.0741	0,0.208	pair:	0.0741	0	-0.0153,0.186	p_acc:	0.0629	p_fet:	0.198
+chr11	4302	4304	single:	34	0	0,0.082	single:	30	0	0,0.0921	pair:	0	0	-0.0611,0.0723	p_acc:	0.285	p_fet:	1
+chr11	4307	4309	single:	31	0	0,0.0894	single:	34	0.0294	0,0.129	pair:	0.0294	0	-0.0497,0.107	p_acc:	0.175	p_fet:	1
+chr11	4314	4316	single:	27	0	0,0.101	single:	33	0.0606	0,0.174	pair:	0.0606	0	-0.0384,0.149	p_acc:	0.094	p_fet:	0.497
+chr11	4320	4322	single:	25	0.04	0,0.17	single:	34	0.0882	0,0.207	pair:	0.0482	0	-0.082,0.153	p_acc:	0.0995	p_fet:	0.63
+chr11	4330	4332	single:	23	0	0,0.117	single:	36	0.0278	0,0.122	pair:	0.0278	0	-0.0768,0.0958	p_acc:	0.177	p_fet:	1
+chr11	4340	4342	single:	22	0.0455	0,0.19	single:	31	0.0323	0,0.14	pair:	-0.0132	-0	-0.143,0.0832	p_acc:	0.134	p_fet:	1
+chr11	4365	4367	single:	20	0	0,0.133	single:	33	0.0303	0,0.132	pair:	0.0303	0	-0.0885,0.103	p_acc:	0.162	p_fet:	1
+chr11	4393	4395	single:	18	0	0,0.146	single:	28	0.143	0,0.288	pair:	0.143	0	-0.0153,0.253	p_acc:	0.0296	p_fet:	0.144
+chr11	4403	4405	single:	16	0	0,0.162	single:	25	0	0,0.109	pair:	0	0	-0.133,0.0747	p_acc:	0.19	p_fet:	1
+chr11	4405	4407	single:	17	0	0,0.153	single:	24	0	0,0.113	pair:	0	0	-0.125,0.0796	p_acc:	0.193	p_fet:	1
+chr11	4420	4422	single:	14	0.0714	0,0.279	single:	22	0.136	0,0.304	pair:	0.0649	0	-0.145,0.216	p_acc:	0.0682	p_fet:	1
+chr11	4456	4458	single:	14	0	0,0.181	single:	15	0	0,0.171	pair:	0	0	-0.14,0.128	p_acc:	0.149	p_fet:	1
+chr11	4458	4460	single:	14	0	0,0.181	single:	16	0.0625	0,0.25	pair:	0.0625	0	-0.101,0.206	p_acc:	0.0874	p_fet:	1
+chr11	4469	4471	single:	14	0.0714	0,0.279	single:	13	0.154	0,0.385	pair:	0.0824	0	-0.132,0.291	p_acc:	0.0559	p_fet:	0.596
+chr11	4482	4484	single:	15	0.0667	0,0.264	single:	19	0	0,0.139	pair:	-0.0667	-0	-0.228,0.0612	p_acc:	0.0804	p_fet:	0.441
+chr11	4499	4501	single:	16	0.0625	0,0.25	single:	24	0.0417	0,0.176	pair:	-0.0208	-0	-0.189,0.102	p_acc:	0.101	p_fet:	1
+chr11	4519	4521	single:	19	0.105	0,0.283	single:	23	0.304	0.123,0.485	pair:	0.199	0	-0.019,0.37	p_acc:	0.021	p_fet:	0.149
+chr11	4523	4525	single:	20	0.05	0,0.207	single:	25	0.28	0.108,0.452	pair:	0.23	0.0331	0.0331,0.378	p_acc:	0.0106	p_fet:	0.0592
+chr11	4532	4534	single:	20	0	0,0.133	single:	23	0.174	0.00602,0.342	pair:	0.174	0.0155	0.0155,0.308	p_acc:	0.0166	p_fet:	0.111
+chr11	4549	4551	single:	21	0.0476	0,0.198	single:	26	0.231	0.0678,0.394	pair:	0.183	0.00182	0.00182,0.327	p_acc:	0.0191	p_fet:	0.112
+chr11	4558	4560	single:	21	0	0,0.127	single:	25	0.2	0.0371,0.363	pair:	0.2	0.0399	0.0399,0.33	p_acc:	0.00952	p_fet:	0.0536
+chr11	4588	4590	single:	21	0.0476	0,0.198	single:	15	0	0,0.171	pair:	-0.0476	-0	-0.159,0.105	p_acc:	0.11	p_fet:	1
+chr11	4594	4596	single:	20	0	0,0.133	single:	18	0.0556	0,0.226	pair:	0.0556	0	-0.0646,0.193	p_acc:	0.0964	p_fet:	0.474
+chr11	4633	4635	single:	15	0.0667	0,0.264	single:	19	0.263	0.0688,0.458	pair:	0.196	0	-0.0347,0.369	p_acc:	0.0237	p_fet:	0.196
+chr11	4692	4694	single:	15	0.0667	0,0.264	single:	16	0.188	0,0.396	pair:	0.121	0	-0.0941,0.308	p_acc:	0.0454	p_fet:	0.6
+chr11	4723	4725	single:	16	0	0,0.162	single:	17	0.235	0.0315,0.439	pair:	0.235	0.0328	0.0328,0.398	p_acc:	0.0105	p_fet:	0.103
+chr11	4744	4746	single:	15	0.0667	0,0.264	single:	21	0.0952	0,0.259	pair:	0.0286	0	-0.161,0.177	p_acc:	0.0837	p_fet:	1
+chr11	4756	4758	single:	19	0	0,0.139	single:	24	0.208	0.0409,0.376	pair:	0.208	0.0368	0.0368,0.34	p_acc:	0.0101	p_fet:	0.0562
+chr11	4785	4787	single:	19	0.0526	0,0.216	single:	22	0.182	0.0087,0.355	pair:	0.129	0	-0.0537,0.284	p_acc:	0.0407	p_fet:	0.35
+chr11	4789	4791	single:	21	0.0476	0,0.198	single:	24	0.0833	0,0.231	pair:	0.0357	0	-0.11,0.167	p_acc:	0.0974	p_fet:	1
+chr11	4791	4793	single:	21	0.19	0.0119,0.369	single:	23	0.261	0.0841,0.438	pair:	0.0704	0	-0.139,0.262	p_acc:	0.0571	p_fet:	0.724
+chr11	4799	4801	single:	24	0.0833	0,0.231	single:	25	0	0,0.109	pair:	-0.0833	-0	-0.203,0.0283	p_acc:	0.0616	p_fet:	0.235
+chr11	4845	4847	single:	24	0.0417	0,0.176	single:	23	0	0,0.117	pair:	-0.0417	-0	-0.147,0.063	p_acc:	0.127	p_fet:	1
+chr11	4860	4862	single:	24	0	0,0.113	single:	17	0	0,0.153	pair:	0	0	-0.0796,0.125	p_acc:	0.193	p_fet:	1
+chr11	4869	4871	single:	24	0	0,0.113	single:	21	0.238	0.0554,0.421	pair:	0.238	0.0726	0.0726,0.39	p_acc:	0.00403	p_fet:	0.0167
+chr11	4872	4874	single:	26	0	0,0.105	single:	20	0.1	0,0.271	pair:	0.1	0	-0.0148,0.243	p_acc:	0.0448	p_fet:	0.184
+chr11	4894	4896	single:	25	0	0,0.109	single:	24	0.0833	0,0.231	pair:	0.0833	0	-0.0283,0.203	p_acc:	0.0616	p_fet:	0.235
+chr11	4908	4910	single:	26	0.0385	0,0.164	single:	27	0.0741	0,0.208	pair:	0.0356	0	-0.0867,0.155	p_acc:	0.109	p_fet:	1
+chr11	4923	4925	single:	26	0	0,0.105	single:	28	0	0,0.0981	pair:	0	0	-0.0815,0.0739	p_acc:	0.249	p_fet:	1
+chr11	4942	4944	single:	25	0	0,0.109	single:	22	0.136	0,0.304	pair:	0.136	0.00635	0.00635,0.275	p_acc:	0.0239	p_fet:	0.095
+chr11	5009	5011	single:	21	0.0476	0,0.198	single:	25	0.08	0,0.223	pair:	0.0324	0	-0.113,0.16	p_acc:	0.101	p_fet:	1
+chr11	5015	5017	single:	19	0.0526	0,0.216	single:	28	0.107	0,0.246	pair:	0.0545	0	-0.109,0.179	p_acc:	0.0842	p_fet:	0.638
+chr11	5025	5027	single:	20	0.05	0,0.207	single:	26	0.192	0.0338,0.351	pair:	0.142	0	-0.0356,0.283	p_acc:	0.034	p_fet:	0.212
+chr11	5036	5038	single:	16	0.125	0,0.326	single:	25	0.04	0,0.17	pair:	-0.085	-0	-0.266,0.057	p_acc:	0.0592	p_fet:	0.55
+chr11	5047	5049	single:	17	0	0,0.153	single:	22	0	0,0.122	pair:	0	0	-0.123,0.088	p_acc:	0.188	p_fet:	1
+chr11	5051	5053	single:	14	0	0,0.181	single:	21	0.0476	0,0.198	pair:	0.0476	0	-0.115,0.157	p_acc:	0.109	p_fet:	1
+chr11	5053	5055	single:	16	0	0,0.162	single:	21	0	0,0.127	pair:	0	0	-0.13,0.0914	p_acc:	0.179	p_fet:	1
+chr11	5097	5099	single:	20	0	0,0.133	single:	24	0	0,0.113	pair:	0	0	-0.105,0.0829	p_acc:	0.209	p_fet:	1
+chr11	5099	5101	single:	18	0	0,0.146	single:	25	0.04	0,0.17	pair:	0.04	0	-0.0902,0.136	p_acc:	0.129	p_fet:	1
+chr11	5104	5106	single:	18	0	0,0.146	single:	23	0.0435	0,0.183	pair:	0.0435	0	-0.0867,0.149	p_acc:	0.121	p_fet:	1
+chr11	5106	5108	single:	17	0.0588	0,0.238	single:	22	0	0,0.122	pair:	-0.0588	-0	-0.206,0.053	p_acc:	0.0899	p_fet:	0.436
+chr11	5114	5116	single:	18	0	0,0.146	single:	19	0	0,0.139	pair:	0	0	-0.112,0.105	p_acc:	0.182	p_fet:	1
+chr11	5120	5122	single:	17	0.0588	0,0.238	single:	20	0.1	0,0.271	pair:	0.0412	0	-0.133,0.194	p_acc:	0.0825	p_fet:	1
+chr11	5124	5126	single:	17	0.0588	0,0.238	single:	18	0	0,0.146	pair:	-0.0588	-0	-0.201,0.073	p_acc:	0.0919	p_fet:	0.486
+chr11	5140	5142	single:	20	0.1	0,0.271	single:	19	0.0526	0,0.216	pair:	-0.0474	-0	-0.2,0.114	p_acc:	0.0828	p_fet:	1
+chr11	5144	5146	single:	18	0.0556	0,0.226	single:	18	0	0,0.146	pair:	-0.0556	-0	-0.19,0.0758	p_acc:	0.0971	p_fet:	1
+chr11	5149	5151	single:	17	0.0588	0,0.238	single:	16	0	0,0.162	pair:	-0.0588	-0	-0.198,0.0867	p_acc:	0.0923	p_fet:	1
+chr11	5151	5153	single:	16	0.0625	0,0.25	single:	15	0	0,0.171	pair:	-0.0625	-0	-0.208,0.0916	p_acc:	0.0874	p_fet:	1
+chr11	5157	5159	single:	16	0.0625	0,0.25	single:	15	0.133	0,0.344	pair:	0.0708	0	-0.12,0.26	p_acc:	0.0635	p_fet:	0.6
+chr11	5159	5161	single:	17	0.0588	0,0.238	single:	15	0.0667	0,0.264	pair:	0.00784	0	-0.153,0.185	p_acc:	0.0894	p_fet:	1
+chr11	5166	5168	single:	14	0.143	0,0.363	single:	16	0	0,0.162	pair:	-0.143	-0	-0.322,0.0318	p_acc:	0.0346	p_fet:	0.209
+chr11	5180	5182	single:	18	0	0,0.146	single:	13	0.0769	0,0.297	pair:	0.0769	0	-0.0593,0.258	p_acc:	0.0696	p_fet:	0.419
+chr11	5185	5187	single:	14	0.0714	0,0.279	single:	12	0	0,0.206	pair:	-0.0714	-0	-0.23,0.115	p_acc:	0.0775	p_fet:	1
+chr11	5190	5192	single:	16	0.25	0.039,0.461	single:	10	0	0,0.238	pair:	-0.25	-0	-0.404,0.017	p_acc:	0.0169	p_fet:	0.136
+chr11	5197	5199	single:	15	0.0667	0,0.264	single:	10	0	0,0.238	pair:	-0.0667	-0	-0.209,0.149	p_acc:	0.081	p_fet:	1
+chr11	5209	5211	single:	15	0.133	0,0.344	single:	9	0.111	0,0.394	pair:	-0.0222	-0	-0.22,0.252	p_acc:	0.0597	p_fet:	1
+chr11	5211	5213	single:	17	0.0588	0,0.238	single:	9	0	0,0.259	pair:	-0.0588	-0	-0.181,0.176	p_acc:	0.0868	p_fet:	1
+chr11	5227	5229	single:	18	0	0,0.146	single:	9	0.111	0,0.394	pair:	0.111	0	-0.0382,0.354	p_acc:	0.0452	p_fet:	0.333
+chr11	5231	5233	single:	17	0.0588	0,0.238	single:	11	0.0909	0,0.339	pair:	0.0321	0	-0.134,0.256	p_acc:	0.0735	p_fet:	1
+chr11	5245	5247	single:	18	0	0,0.146	single:	13	0	0,0.193	pair:	0	0	-0.103,0.155	p_acc:	0.154	p_fet:	1
+chr11	5282	5284	single:	15	0	0,0.171	single:	12	0	0,0.206	pair:	0	0	-0.124,0.163	p_acc:	0.139	p_fet:	1
+chr11	5285	5287	single:	17	0	0,0.153	single:	12	0.167	0,0.41	pair:	0.167	0	-0.012,0.369	p_acc:	0.0255	p_fet:	0.163
+chr11	5291	5293	single:	18	0.0556	0,0.226	single:	13	0.0769	0,0.297	pair:	0.0214	0	-0.135,0.219	p_acc:	0.0833	p_fet:	1
+chr11	5300	5302	single:	21	0	0,0.127	single:	13	0	0,0.193	pair:	0	0	-0.0866,0.159	p_acc:	0.161	p_fet:	1
+chr11	5318	5320	single:	22	0.136	0,0.304	single:	11	0.364	0.109,0.618	pair:	0.227	0	-0.0239,0.469	p_acc:	0.0193	p_fet:	0.186
+chr11	5342	5344	single:	21	0.0476	0,0.198	single:	12	0.167	0,0.41	pair:	0.119	0	-0.0554,0.339	p_acc:	0.0427	p_fet:	0.538
+chr11	5351	5353	single:	20	0.3	0.107,0.493	single:	14	0.429	0.196,0.662	pair:	0.129	0	-0.136,0.373	p_acc:	0.0382	p_fet:	0.487
+chr11	5379	5381	single:	16	0.125	0,0.326	single:	16	0	0,0.162	pair:	-0.125	-0	-0.285,0.043	p_acc:	0.0426	p_fet:	0.484
+chr11	5381	5383	single:	16	0.0625	0,0.25	single:	16	0	0,0.162	pair:	-0.0625	-0	-0.21,0.0836	p_acc:	0.0872	p_fet:	1
+chr11	5394	5396	single:	15	0.0667	0,0.264	single:	15	0.0667	0,0.264	pair:	0	0	-0.178,0.178	p_acc:	0.0851	p_fet:	1
+chr11	5410	5412	single:	14	0	0,0.181	single:	21	0	0,0.127	pair:	0	0	-0.148,0.0884	p_acc:	0.168	p_fet:	1
+chr11	5414	5416	single:	14	0.0714	0,0.279	single:	21	0.143	0,0.316	pair:	0.0714	0	-0.14,0.228	p_acc:	0.0651	p_fet:	0.635
+chr11	5417	5419	single:	15	0	0,0.171	single:	20	0	0,0.133	pair:	0	0	-0.137,0.0951	p_acc:	0.171	p_fet:	1
+chr11	5472	5474	single:	22	0	0,0.122	single:	23	0.174	0.00602,0.342	pair:	0.174	0.0232	0.0232,0.31	p_acc:	0.0144	p_fet:	0.109
+chr11	5498	5500	single:	23	0.0435	0,0.183	single:	21	0.143	0,0.316	pair:	0.0994	0	-0.0544,0.254	p_acc:	0.0529	p_fet:	0.335
+chr11	5532	5534	single:	30	0.0333	0,0.144	single:	26	0.115	0,0.263	pair:	0.0821	0	-0.0401,0.214	p_acc:	0.0625	p_fet:	0.328
+chr11	5535	5537	single:	31	0.0323	0,0.14	single:	26	0.231	0.0678,0.394	pair:	0.199	0.0474	0.0474,0.343	p_acc:	0.00812	p_fet:	0.0394
+chr11	5538	5540	single:	33	0.242	0.0967,0.388	single:	26	0.308	0.136,0.479	pair:	0.0653	0	-0.121,0.252	p_acc:	0.0604	p_fet:	0.769
+chr11	5551	5553	single:	29	0.138	0,0.28	single:	23	0.348	0.163,0.532	pair:	0.21	0.0113	0.0113,0.388	p_acc:	0.0154	p_fet:	0.102
+chr11	5564	5566	single:	27	0.037	0,0.159	single:	19	0	0,0.139	pair:	-0.037	-0	-0.127,0.087	p_acc:	0.138	p_fet:	1
+chr11	5576	5578	single:	28	0.0357	0,0.153	single:	21	0.0476	0,0.198	pair:	0.0119	0	-0.0935,0.146	p_acc:	0.126	p_fet:	1
+chr11	5600	5602	single:	25	0.08	0,0.223	single:	24	0.0833	0,0.231	pair:	0.00333	0	-0.136,0.146	p_acc:	0.101	p_fet:	1
+chr11	5622	5624	single:	21	0.0476	0,0.198	single:	19	0.211	0.02,0.401	pair:	0.163	0	-0.0196,0.333	p_acc:	0.0266	p_fet:	0.172
+chr11	5624	5626	single:	20	0.2	0.0156,0.384	single:	17	0.235	0.0315,0.439	pair:	0.0353	0	-0.179,0.255	p_acc:	0.0594	p_fet:	1
+chr11	5641	5643	single:	19	0	0,0.139	single:	17	0.0588	0,0.238	pair:	0.0588	0	-0.0672,0.202	p_acc:	0.0915	p_fet:	0.472
+chr11	5667	5669	single:	16	0	0,0.162	single:	24	0.0417	0,0.176	pair:	0.0417	0	-0.103,0.139	p_acc:	0.123	p_fet:	1
+chr11	5720	5722	single:	22	0	0,0.122	single:	24	0.125	0,0.282	pair:	0.125	0	-0.0104,0.251	p_acc:	0.0329	p_fet:	0.235
+chr11	5816	5818	single:	14	0	0,0.181	single:	21	0.238	0.0554,0.421	pair:	0.238	0.0251	0.0251,0.376	p_acc:	0.0115	p_fet:	0.0689
+chr11	5828	5830	single:	12	0	0,0.206	single:	21	0.238	0.0554,0.421	pair:	0.238	0.00538	0.00538,0.371	p_acc:	0.0145	p_fet:	0.133
+chr11	5830	5832	single:	13	0	0,0.193	single:	22	0.0455	0,0.19	pair:	0.0455	0	-0.128,0.148	p_acc:	0.112	p_fet:	1
+chr11	5843	5845	single:	11	0.0909	0,0.339	single:	18	0.0556	0,0.226	pair:	-0.0354	-0	-0.259,0.125	p_acc:	0.0739	p_fet:	1
+chr11	5905	5907	single:	3	0.333	0,0.751	single:	NA	NA	NA,NA	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
+chr11	5969	5971	single:	10	0.3	0.035,0.565	single:	15	0.267	0.0481,0.485	pair:	-0.0333	-0	-0.322,0.235	p_acc:	0.0461	p_fet:	1
+chr11	5972	5974	single:	11	0.0909	0,0.339	single:	15	0	0,0.171	pair:	-0.0909	-0	-0.294,0.0699	p_acc:	0.0601	p_fet:	0.423
+chr11	6023	6025	single:	11	0.0909	0,0.339	single:	18	0.333	0.128,0.538	pair:	0.242	0	-0.0451,0.422	p_acc:	0.0199	p_fet:	0.202
+chr11	6044	6046	single:	10	0	0,0.238	single:	12	0	0,0.206	pair:	0	0	-0.187,0.15	p_acc:	0.117	p_fet:	1
+chr11	6049	6051	single:	10	0	0,0.238	single:	12	0.0833	0,0.316	pair:	0.0833	0	-0.134,0.259	p_acc:	0.0677	p_fet:	1
+chr11	6076	6078	single:	11	0.0909	0,0.339	single:	7	0	0,0.312	pair:	-0.0909	-0	-0.266,0.195	p_acc:	0.0616	p_fet:	1
+chr11	6089	6091	single:	9	0	0,0.259	single:	10	0	0,0.238	pair:	0	0	-0.2,0.177	p_acc:	0.105	p_fet:	1
+chr11	6099	6101	single:	6	0	0,0.348	single:	9	0.111	0,0.394	pair:	0.111	0	-0.21,0.312	p_acc:	0.0524	p_fet:	1
+chr11	6148	6150	single:	9	0	0,0.259	single:	13	0	0,0.193	pair:	0	0	-0.209,0.135	p_acc:	0.115	p_fet:	1
+chr11	6156	6158	single:	8	0	0,0.283	single:	14	0	0,0.181	pair:	0	0	-0.235,0.121	p_acc:	0.111	p_fet:	1
+chr11	6164	6166	single:	8	0	0,0.283	single:	12	0	0,0.206	pair:	0	0	-0.23,0.143	p_acc:	0.106	p_fet:	1
+chr11	6168	6170	single:	8	0.125	0,0.429	single:	14	0	0,0.181	pair:	-0.125	-0	-0.38,0.0569	p_acc:	0.0427	p_fet:	0.364
+chr11	6208	6210	single:	10	0.1	0,0.364	single:	10	0.2	0,0.47	pair:	0.1	0	-0.177,0.347	p_acc:	0.0455	p_fet:	1
+chr11	6226	6228	single:	10	0	0,0.238	single:	10	0	0,0.238	pair:	0	0	-0.181,0.181	p_acc:	0.109	p_fet:	1
+chr11	6239	6241	single:	5	0	0,0.393	single:	12	0	0,0.206	pair:	0	0	-0.336,0.126	p_acc:	0.0829	p_fet:	1
+chr11	6247	6249	single:	11	0.0909	0,0.339	single:	11	0.0909	0,0.339	pair:	0	0	-0.226,0.226	p_acc:	0.0655	p_fet:	1
+chr11	6251	6253	single:	13	0.0769	0,0.297	single:	9	0.222	0,0.507	pair:	0.145	0	-0.105,0.403	p_acc:	0.037	p_fet:	0.544
+chr11	6253	6255	single:	13	0.0769	0,0.297	single:	9	0.111	0,0.394	pair:	0.0342	0	-0.172,0.291	p_acc:	0.0615	p_fet:	1
+chr11	6258	6260	single:	12	0	0,0.206	single:	11	0.0909	0,0.339	pair:	0.0909	0	-0.0993,0.287	p_acc:	0.0621	p_fet:	0.478
+chr11	6269	6271	single:	15	0	0,0.171	single:	12	0.0833	0,0.316	pair:	0.0833	0	-0.0753,0.272	p_acc:	0.0659	p_fet:	0.444
+chr11	6272	6274	single:	15	0.2	0,0.417	single:	11	0.0909	0,0.339	pair:	-0.109	-0	-0.307,0.152	p_acc:	0.0476	p_fet:	0.614
+chr11	6282	6284	single:	15	0.0667	0,0.264	single:	10	0.1	0,0.364	pair:	0.0333	0	-0.151,0.272	p_acc:	0.0675	p_fet:	1
+chr11	6293	6295	single:	16	0	0,0.162	single:	14	0.143	0,0.363	pair:	0.143	0	-0.0318,0.322	p_acc:	0.0346	p_fet:	0.209
+chr11	6297	6299	single:	18	0	0,0.146	single:	13	0	0,0.193	pair:	0	0	-0.103,0.155	p_acc:	0.154	p_fet:	1
+chr11	6300	6302	single:	15	0.133	0,0.344	single:	13	0.154	0,0.385	pair:	0.0205	0	-0.195,0.249	p_acc:	0.0617	p_fet:	1
+chr11	6304	6306	single:	17	0.118	0,0.31	single:	12	0.0833	0,0.316	pair:	-0.0343	-0	-0.21,0.192	p_acc:	0.0696	p_fet:	1
+chr11	6312	6314	single:	20	0.05	0,0.207	single:	12	0	0,0.206	pair:	-0.05	-0	-0.161,0.135	p_acc:	0.103	p_fet:	1
+chr11	6314	6316	single:	20	0.05	0,0.207	single:	12	0.167	0,0.41	pair:	0.117	0	-0.0623,0.337	p_acc:	0.044	p_fet:	0.54
+chr11	6318	6320	single:	20	0.05	0,0.207	single:	12	0	0,0.206	pair:	-0.05	-0	-0.161,0.135	p_acc:	0.103	p_fet:	1
+chr11	6321	6323	single:	19	0	0,0.139	single:	13	0	0,0.193	pair:	0	0	-0.0972,0.157	p_acc:	0.156	p_fet:	1
+chr11	6323	6325	single:	20	0	0,0.133	single:	12	0	0,0.206	pair:	0	0	-0.0897,0.171	p_acc:	0.153	p_fet:	1
+chr11	6350	6352	single:	18	0.0556	0,0.226	single:	12	0	0,0.206	pair:	-0.0556	-0	-0.179,0.13	p_acc:	0.0952	p_fet:	1
+chr11	6356	6358	single:	19	0.0526	0,0.216	single:	10	0.1	0,0.364	pair:	0.0474	0	-0.11,0.287	p_acc:	0.0684	p_fet:	1
+chr11	6358	6360	single:	20	0.05	0,0.207	single:	11	0	0,0.221	pair:	-0.05	-0	-0.158,0.149	p_acc:	0.101	p_fet:	1
+chr11	6365	6367	single:	19	0	0,0.139	single:	11	0	0,0.221	pair:	0	0	-0.0929,0.184	p_acc:	0.143	p_fet:	1
+chr11	6371	6373	single:	18	0.0556	0,0.226	single:	13	0.0769	0,0.297	pair:	0.0214	0	-0.135,0.219	p_acc:	0.0833	p_fet:	1
+chr11	6430	6432	single:	14	0	0,0.181	single:	19	0.105	0,0.283	pair:	0.105	0	-0.0731,0.239	p_acc:	0.0573	p_fet:	0.496
+chr11	6432	6434	single:	16	0.0625	0,0.25	single:	19	0.105	0,0.283	pair:	0.0428	0	-0.141,0.202	p_acc:	0.0788	p_fet:	1
+chr11	6442	6444	single:	20	0.1	0,0.271	single:	22	0.0909	0,0.249	pair:	-0.00909	-0	-0.174,0.147	p_acc:	0.0879	p_fet:	1
+chr11	6449	6451	single:	19	0	0,0.139	single:	24	0.0417	0,0.176	pair:	0.0417	0	-0.0824,0.143	p_acc:	0.126	p_fet:	1
+chr11	6491	6493	single:	25	0.04	0,0.17	single:	18	0.278	0.0771,0.478	pair:	0.238	0.0497	0.0497,0.416	p_acc:	0.00794	p_fet:	0.0672
+chr11	6526	6528	single:	35	0.0286	0,0.125	single:	26	0.115	0,0.263	pair:	0.0868	0	-0.0247,0.219	p_acc:	0.0549	p_fet:	0.303
+chr11	6560	6562	single:	25	0.24	0.0727,0.407	single:	30	0.233	0.0816,0.385	pair:	-0.00667	-0	-0.195,0.173	p_acc:	0.0714	p_fet:	1
+chr11	6596	6598	single:	38	0.421	0.27,0.572	single:	30	0.6	0.433,0.767	pair:	0.179	0	-0.0236,0.356	p_acc:	0.0237	p_fet:	0.222
+chr11	6730	6732	single:	41	0.366	0.223,0.509	single:	31	0.516	0.349,0.683	pair:	0.15	0	-0.0438,0.327	p_acc:	0.0316	p_fet:	0.236
+chr11	6857	6859	single:	33	0.212	0.0704,0.354	single:	46	0.239	0.116,0.362	pair:	0.027	0	-0.135,0.173	p_acc:	0.0826	p_fet:	1
+chr11	6900	6902	single:	29	0.862	0.72,1	single:	46	0.804	0.687,0.921	pair:	-0.0577	-0	-0.187,0.0999	p_acc:	0.0766	p_fet:	0.755
+chr11	6902	6904	single:	29	0.483	0.311,0.655	single:	46	0.804	0.687,0.921	pair:	0.322	0.132	0.132,0.479	p_acc:	0.00124	p_fet:	0.00524
+chr11	6979	6981	single:	25	0.28	0.108,0.452	single:	39	0.487	0.337,0.638	pair:	0.207	0	-0.00347,0.377	p_acc:	0.0177	p_fet:	0.123
+chr11	7085	7087	single:	26	0.346	0.172,0.521	single:	27	0.704	0.537,0.871	pair:	0.358	0.124	0.124,0.528	p_acc:	0.00216	p_fet:	0.0135
+chr11	7166	7168	single:	19	0.263	0.0688,0.458	single:	27	0.259	0.0961,0.422	pair:	-0.0039	-0	-0.221,0.195	p_acc:	0.0632	p_fet:	1
+chr11	7336	7338	single:	17	0.294	0.0867,0.502	single:	13	0.769	0.534,1	pair:	0.475	0.154	0.154,0.653	p_acc:	0.00186	p_fet:	0.0253
+chr11	7583	7585	single:	26	0.154	6.6e-05,0.308	single:	34	0.176	0.0419,0.311	pair:	0.0226	0	-0.145,0.172	p_acc:	0.0825	p_fet:	1
+chr11	7633	7635	single:	21	0.0952	0,0.259	single:	26	0.154	6.6e-05,0.308	pair:	0.0586	0	-0.116,0.209	p_acc:	0.0723	p_fet:	0.678
+chr11	7655	7657	single:	15	0	0,0.171	single:	21	0.143	0,0.316	pair:	0.143	0	-0.0375,0.275	p_acc:	0.0356	p_fet:	0.25
+chr11	7673	7675	single:	16	0	0,0.162	single:	17	0	0,0.153	pair:	0	0	-0.125,0.115	p_acc:	0.165	p_fet:	1
+chr11	7679	7681	single:	17	0	0,0.153	single:	21	0.0476	0,0.198	pair:	0.0476	0	-0.0895,0.162	p_acc:	0.111	p_fet:	1
+chr11	7688	7690	single:	15	0	0,0.171	single:	18	0.167	0,0.359	pair:	0.167	0	-0.0215,0.317	p_acc:	0.0269	p_fet:	0.233
+chr11	7695	7697	single:	15	0.0667	0,0.264	single:	21	0	0,0.127	pair:	-0.0667	-0	-0.23,0.0515	p_acc:	0.079	p_fet:	0.417
+chr11	7729	7731	single:	14	0	0,0.181	single:	21	0.19	0.0119,0.369	pair:	0.19	0	-0.0108,0.326	p_acc:	0.021	p_fet:	0.133
+chr11	7748	7750	single:	15	0	0,0.171	single:	23	0.0435	0,0.183	pair:	0.0435	0	-0.11,0.144	p_acc:	0.118	p_fet:	1
+chr11	7825	7827	single:	18	0	0,0.146	single:	20	0.05	0,0.207	pair:	0.05	0	-0.0806,0.171	p_acc:	0.107	p_fet:	1
+chr11	7866	7868	single:	17	0	0,0.153	single:	22	0.0455	0,0.19	pair:	0.0455	0	-0.0916,0.154	p_acc:	0.116	p_fet:	1
+chr11	7878	7880	single:	15	0.0667	0,0.264	single:	21	0.0952	0,0.259	pair:	0.0286	0	-0.161,0.177	p_acc:	0.0837	p_fet:	1
+chr11	7882	7884	single:	13	0.0769	0,0.297	single:	20	0.3	0.107,0.493	pair:	0.223	0	-0.0326,0.392	p_acc:	0.0204	p_fet:	0.202
+chr11	8027	8029	single:	26	0.0385	0,0.164	single:	18	0.167	0,0.359	pair:	0.128	0	-0.0241,0.302	p_acc:	0.0357	p_fet:	0.289
+chr11	8033	8035	single:	27	0.037	0,0.159	single:	16	0.0625	0,0.25	pair:	0.0255	0	-0.0862,0.195	p_acc:	0.103	p_fet:	1
+chr11	8303	8305	single:	36	0.0556	0,0.161	single:	23	0.0435	0,0.183	pair:	-0.0121	-0	-0.106,0.119	p_acc:	0.131	p_fet:	1
+chr11	8414	8416	single:	16	0.0625	0,0.25	single:	27	0.444	0.268,0.621	pair:	0.382	0.14	0.14,0.521	p_acc:	0.00165	p_fet:	0.0143
+chr11	8583	8585	single:	24	0.0833	0,0.231	single:	27	0.148	0,0.298	pair:	0.0648	0	-0.0948,0.209	p_acc:	0.0723	p_fet:	0.671
+chr11	8678	8680	single:	19	0	0,0.139	single:	17	0	0,0.153	pair:	0	0	-0.104,0.119	p_acc:	0.178	p_fet:	1
+chr11	8692	8694	single:	23	0.0435	0,0.183	single:	19	0.0526	0,0.216	pair:	0.00915	0	-0.116,0.155	p_acc:	0.113	p_fet:	1
+chr11	8755	8757	single:	26	0	0,0.105	single:	24	0.0417	0,0.176	pair:	0.0417	0	-0.0523,0.15	p_acc:	0.126	p_fet:	0.48
+chr11	8786	8788	single:	25	0	0,0.109	single:	19	0	0,0.139	pair:	0	0	-0.0781,0.112	p_acc:	0.207	p_fet:	1
+chr11	8812	8814	single:	22	0.0909	0,0.249	single:	21	0	0,0.127	pair:	-0.0909	-0	-0.217,0.038	p_acc:	0.0592	p_fet:	0.488
+chr11	8829	8831	single:	20	0	0,0.133	single:	22	0.136	0,0.304	pair:	0.136	0	-0.0114,0.27	p_acc:	0.0304	p_fet:	0.233
+chr11	8836	8838	single:	21	0	0,0.127	single:	21	0.0476	0,0.198	pair:	0.0476	0	-0.0664,0.167	p_acc:	0.112	p_fet:	1
+chr11	8855	8857	single:	17	0.0588	0,0.238	single:	22	0.0909	0,0.249	pair:	0.0321	0	-0.14,0.174	p_acc:	0.0887	p_fet:	1
+chr11	8902	8904	single:	17	0	0,0.153	single:	15	0.333	0.111,0.555	pair:	0.333	0.104	0.104,0.508	p_acc:	0.00262	p_fet:	0.0149
+chr11	8914	8916	single:	20	0.35	0.153,0.547	single:	13	0.308	0.0724,0.543	pair:	-0.0423	-0	-0.281,0.229	p_acc:	0.0495	p_fet:	1
+chr11	8920	8922	single:	22	0	0,0.122	single:	12	0	0,0.206	pair:	0	0	-0.0801,0.173	p_acc:	0.157	p_fet:	1
+chr11	8948	8950	single:	33	0	0,0.0843	single:	17	0.0588	0,0.238	pair:	0.0588	0	-0.0234,0.214	p_acc:	0.0814	p_fet:	0.34
+chr11	8975	8977	single:	42	0.0476	0,0.139	single:	19	0	0,0.139	pair:	-0.0476	-0	-0.11,0.0832	p_acc:	0.131	p_fet:	1
+chr11	9127	9129	single:	41	0.0244	0,0.108	single:	31	0.0968	0,0.225	pair:	0.0724	0	-0.0228,0.187	p_acc:	0.0662	p_fet:	0.308
+chr11	9136	9138	single:	42	0.0476	0,0.139	single:	33	0.0606	0,0.174	pair:	0.013	0	-0.0771,0.12	p_acc:	0.143	p_fet:	1
+chr11	9227	9229	single:	45	0.0444	0,0.131	single:	38	0.0526	0,0.153	pair:	0.00819	0	-0.0754,0.103	p_acc:	0.16	p_fet:	1
+chr11	9270	9272	single:	40	0.025	0,0.111	single:	33	0.0909	0,0.213	pair:	0.0659	0	-0.0287,0.175	p_acc:	0.0758	p_fet:	0.322
+chr11	9301	9303	single:	40	0.175	0.0523,0.298	single:	33	0.0303	0,0.132	pair:	-0.145	-0.0191	-0.253,-0.0191	p_acc:	0.0174	p_fet:	0.0652
+chr11	9371	9373	single:	25	0	0,0.109	single:	32	0	0,0.0868	pair:	0	0	-0.0872,0.0628	p_acc:	0.257	p_fet:	1
+chr11	9380	9382	single:	25	0.04	0,0.17	single:	31	0	0,0.0894	pair:	-0.04	-0	-0.147,0.0405	p_acc:	0.129	p_fet:	0.446
+chr11	9447	9449	single:	17	0	0,0.153	single:	30	0	0,0.0921	pair:	0	0	-0.129,0.0611	p_acc:	0.207	p_fet:	1
+chr11	9478	9480	single:	16	0	0,0.162	single:	21	0.0476	0,0.198	pair:	0.0476	0	-0.097,0.16	p_acc:	0.111	p_fet:	1
+chr11	9510	9512	single:	14	0	0,0.181	single:	23	0.0435	0,0.183	pair:	0.0435	0	-0.119,0.143	p_acc:	0.117	p_fet:	1
+chr11	9588	9590	single:	20	0	0,0.133	single:	18	0.111	0,0.296	pair:	0.111	0	-0.0292,0.261	p_acc:	0.0448	p_fet:	0.218
+chr11	9677	9679	single:	11	0.0909	0,0.339	single:	13	0.231	0,0.466	pair:	0.14	0	-0.131,0.354	p_acc:	0.0398	p_fet:	0.596
+chr11	9697	9699	single:	14	0	0,0.181	single:	15	0.2	0,0.417	pair:	0.2	0	-0.00842,0.371	p_acc:	0.02	p_fet:	0.224
+chr11	9709	9711	single:	11	0	0,0.221	single:	18	0.111	0,0.296	pair:	0.111	0	-0.104,0.245	p_acc:	0.0574	p_fet:	0.512
+chr11	9729	9731	single:	6	0.167	0,0.521	single:	16	0	0,0.162	pair:	-0.167	-0	-0.475,0.0224	p_acc:	0.0282	p_fet:	0.273
+chr11	9799	9801	single:	NA	NA	NA,NA	single:	7	0	0,0.312	pair:	NA	NA	NA,NA	p_acc:	NA	p_fet:	NA
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/SE_dmr_M3_g1.G.bed_vs_g2.G.bed.txt.dmr	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,2 @@
+chr11	316	1705	-	35	30	5	-1
+chr11	3052	3674	-	17	16	1	-1
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/SE_g1.G.bed	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,309 @@
+#chrom	start	end	ratio	totalC	methC	strand	next	Plus	totalC	methC	Minus	totalC	methC	localSeq
+chr11	41	43	0	3	0	B	G	+	1	0	-	2	0	AGCGG
+chr11	48	50	0	5	0	B	G	+	2	0	-	3	0	CCCGG
+chr11	56	58	0	6	0	B	G	+	2	0	-	4	0	CCCGC
+chr11	59	61	0	8	0	B	G	+	2	0	-	6	0	GCCGT
+chr11	101	103	0.0769	13	1	B	G	+	3	0	-	10	1	GTCGC
+chr11	140	142	0	16	0	B	G	+	3	0	-	13	0	GGCGA
+chr11	144	146	0	14	0	B	G	+	2	0	-	12	0	AACGG
+chr11	210	212	0.133	15	2	B	G	+	2	0	-	13	2	CCCGA
+chr11	215	217	0.263	19	5	B	G	+	3	1	-	16	4	CTCGC
+chr11	316	318	0.0833	24	2	B	G	+	4	0	-	20	2	CTCGC
+chr11	331	333	0.037	27	1	B	G	+	5	1	-	22	0	TGCGT
+chr11	400	402	0.0588	34	2	B	G	+	3	0	-	31	2	CACGC
+chr11	407	409	0	33	0	B	G	+	3	0	-	30	0	CTCGT
+chr11	445	447	0.08	25	2	B	G	+	3	0	-	22	2	CACGC
+chr11	469	471	0.0385	26	1	B	G	+	4	0	-	22	1	ACCGC
+chr11	546	548	0.1	20	2	B	G	+	9	1	-	11	1	GTCGA
+chr11	666	668	0.05	20	1	B	G	+	9	1	-	11	0	TGCGT
+chr11	694	696	0.318	22	7	B	G	+	10	4	-	12	3	TGCGT
+chr11	712	714	0.0714	28	2	B	G	+	13	0	-	15	2	ACCGA
+chr11	721	723	0.0357	28	1	B	G	+	11	0	-	17	1	CCCGG
+chr11	737	739	0.037	27	1	B	G	+	9	0	-	18	1	GGCGG
+chr11	753	755	0.087	23	2	B	G	+	8	1	-	15	1	TCCGG
+chr11	762	764	0.037	27	1	B	G	+	9	1	-	18	0	CCCGT
+chr11	792	794	0.324	37	12	B	G	+	11	1	-	26	11	GCCGT
+chr11	798	800	0.306	36	11	B	G	+	10	2	-	26	9	TGCGC
+chr11	816	818	0.0303	33	1	B	G	+	8	1	-	25	0	TTCGC
+chr11	830	832	0.405	37	15	B	G	+	7	1	-	30	14	TGCGA
+chr11	954	956	0	47	0	B	G	+	17	0	-	30	0	CCCGT
+chr11	977	979	0.255	47	12	B	G	+	20	5	-	27	7	CCCGG
+chr11	1006	1008	0	38	0	B	G	+	19	0	-	19	0	ACCGT
+chr11	1066	1068	0.367	30	11	B	G	+	17	7	-	13	4	GGCGG
+chr11	1116	1118	0.182	22	4	B	G	+	13	1	-	9	3	CACGG
+chr11	1145	1147	0.556	18	10	B	G	+	8	5	-	10	5	CCCGA
+chr11	1260	1262	0.0556	18	1	B	G	+	12	1	-	6	0	GACGA
+chr11	1320	1322	0.118	17	2	B	G	+	11	1	-	6	1	TCCGC
+chr11	1559	1561	0.194	36	7	B	G	+	25	4	-	11	3	TACGC
+chr11	1574	1576	0.152	33	5	B	G	+	20	4	-	13	1	ATCGT
+chr11	1592	1594	0.343	35	12	B	G	+	22	6	-	13	6	CTCGG
+chr11	1609	1611	0.0588	34	2	B	G	+	23	2	-	11	0	CTCGG
+chr11	1632	1634	0.0556	36	2	B	G	+	23	2	-	13	0	GACGG
+chr11	1642	1644	0.188	32	6	B	G	+	18	5	-	14	1	CCCGC
+chr11	1680	1682	0.0882	34	3	B	G	+	21	3	-	13	0	GTCGC
+chr11	1689	1691	0.2	35	7	B	G	+	21	4	-	14	3	CTCGG
+chr11	1703	1705	0.262	42	11	B	G	+	22	7	-	20	4	ACCGA
+chr11	1733	1735	0.512	43	22	B	G	+	25	14	-	18	8	CTCGT
+chr11	1766	1768	0.13	46	6	B	G	+	25	5	-	21	1	TCCGT
+chr11	1869	1871	0.244	45	11	B	G	+	21	5	-	24	6	TACGT
+chr11	2040	2042	0	37	0	B	G	+	18	0	-	19	0	AACGA
+chr11	2104	2106	0.027	37	1	B	G	+	19	0	-	18	1	AACGA
+chr11	2148	2150	0.613	31	19	B	G	+	14	11	-	17	8	GGCGT
+chr11	2360	2362	0.262	42	11	B	G	+	23	6	-	19	5	GCCGG
+chr11	2384	2386	0.382	34	13	B	G	+	18	6	-	16	7	TACGG
+chr11	2454	2456	0.171	41	7	B	G	+	29	4	-	12	3	TCCGC
+chr11	2621	2623	0	9	0	B	G	+	6	0	-	3	0	ATCGT
+chr11	2648	2650	0	10	0	B	G	+	7	0	-	3	0	GGCGT
+chr11	2665	2667	0	13	0	B	G	+	9	0	-	4	0	GCCGC
+chr11	2667	2669	0	17	0	B	G	+	13	0	-	4	0	CGCGG
+chr11	2673	2675	0	13	0	B	G	+	9	0	-	4	0	GGCGG
+chr11	2676	2678	0	17	0	B	G	+	13	0	-	4	0	GGCGA
+chr11	2685	2687	0.0625	16	1	B	G	+	12	1	-	4	0	GGCGG
+chr11	2688	2690	0	19	0	B	G	+	15	0	-	4	0	GGCGC
+chr11	2691	2693	0	20	0	B	G	+	16	0	-	4	0	GCCGG
+chr11	2713	2715	0	22	0	B	G	+	17	0	-	5	0	TTCGG
+chr11	2723	2725	0	19	0	B	G	+	16	0	-	3	0	GTCGC
+chr11	2727	2729	0	20	0	B	G	+	16	0	-	4	0	CACGG
+chr11	2731	2733	0.0588	17	1	B	G	+	13	0	-	4	1	GCCGG
+chr11	2738	2740	0	18	0	B	G	+	14	0	-	4	0	CCCGC
+chr11	2767	2769	0	15	0	B	G	+	11	0	-	4	0	GCCGG
+chr11	2778	2780	0.0833	12	1	B	G	+	8	0	-	4	1	TCCGG
+chr11	2788	2790	0	8	0	B	G	+	5	0	-	3	0	GCCGC
+chr11	2792	2794	0.143	7	1	B	G	+	4	0	-	3	1	CCCGT
+chr11	2810	2812	0	3	0	B	G	+	1	0	-	2	0	ACCGA
+chr11	2825	2827	0	4	0	B	G	+	3	0	-	1	0	TCCGC
+chr11	3031	3033	0.0769	13	1	B	G	+	6	0	-	7	1	CTCGG
+chr11	3034	3036	0	13	0	B	G	+	6	0	-	7	0	GGCGG
+chr11	3052	3054	0.0667	15	1	B	G	+	6	0	-	9	1	ATCGA
+chr11	3095	3097	0.273	11	3	B	G	+	4	2	-	7	1	ATCGG
+chr11	3117	3119	0.0833	12	1	B	G	+	3	0	-	9	1	ACCGT
+chr11	3273	3275	0.105	38	4	B	G	+	11	2	-	27	2	AGCGC
+chr11	3283	3285	0.027	37	1	B	G	+	11	1	-	26	0	ACCGG
+chr11	3323	3325	0.0571	35	2	B	G	+	13	2	-	22	0	ACCGC
+chr11	3424	3426	0.0667	30	2	B	G	+	16	0	-	14	2	GACGC
+chr11	3452	3454	0	27	0	B	G	+	17	0	-	10	0	ACCGT
+chr11	3484	3486	0.04	25	1	B	G	+	21	1	-	4	0	ATCGA
+chr11	3578	3580	0.263	19	5	B	G	+	14	5	-	5	0	TCCGC
+chr11	3595	3597	0	20	0	B	G	+	15	0	-	5	0	GACGT
+chr11	3611	3613	0.0526	19	1	B	G	+	12	0	-	7	1	CCCGG
+chr11	3617	3619	0	19	0	B	G	+	12	0	-	7	0	CTCGG
+chr11	3639	3641	0.125	16	2	B	G	+	8	1	-	8	1	TCCGG
+chr11	3646	3648	0.0588	17	1	B	G	+	9	0	-	8	1	ACCGC
+chr11	3665	3667	0.0769	13	1	B	G	+	8	0	-	5	1	GCCGG
+chr11	3672	3674	0	13	0	B	G	+	8	0	-	5	0	GTCGG
+chr11	3677	3679	0	14	0	B	G	+	10	0	-	4	0	TTCGC
+chr11	3693	3695	0	17	0	B	G	+	11	0	-	6	0	CTCGC
+chr11	3704	3706	0	20	0	B	G	+	12	0	-	8	0	TCCGC
+chr11	3707	3709	0	21	0	B	G	+	13	0	-	8	0	GCCGC
+chr11	3717	3719	0.0455	22	1	B	G	+	12	1	-	10	0	TCCGC
+chr11	3960	3962	0	21	0	B	G	+	2	0	-	19	0	CCCGA
+chr11	3991	3993	0	20	0	B	G	+	2	0	-	18	0	GACGT
+chr11	4056	4058	0	22	0	B	G	+	7	0	-	15	0	GACGC
+chr11	4163	4165	0.0769	26	2	B	G	+	17	1	-	9	1	CACGT
+chr11	4195	4197	0.0435	23	1	B	G	+	14	0	-	9	1	GGCGG
+chr11	4202	4204	0	21	0	B	G	+	12	0	-	9	0	GTCGC
+chr11	4207	4209	0.0556	18	1	B	G	+	11	0	-	7	1	TCCGG
+chr11	4243	4245	0.133	30	4	B	G	+	21	3	-	9	1	GTCGG
+chr11	4248	4250	0.0312	32	1	B	G	+	22	1	-	10	0	GGCGC
+chr11	4284	4286	0.0882	34	3	B	G	+	22	1	-	12	2	TCCGC
+chr11	4289	4291	0	33	0	B	G	+	19	0	-	14	0	TCCGG
+chr11	4302	4304	0	34	0	B	G	+	18	0	-	16	0	ATCGT
+chr11	4307	4309	0	31	0	B	G	+	16	0	-	15	0	CACGA
+chr11	4314	4316	0	27	0	B	G	+	12	0	-	15	0	AACGG
+chr11	4320	4322	0.04	25	1	B	G	+	11	0	-	14	1	ACCGT
+chr11	4330	4332	0	23	0	B	G	+	8	0	-	15	0	AACGA
+chr11	4340	4342	0.0455	22	1	B	G	+	6	0	-	16	1	TCCGT
+chr11	4365	4367	0	20	0	B	G	+	6	0	-	14	0	GACGA
+chr11	4393	4395	0	18	0	B	G	+	6	0	-	12	0	TTCGA
+chr11	4403	4405	0	16	0	B	G	+	6	0	-	10	0	GCCGC
+chr11	4405	4407	0	17	0	B	G	+	6	0	-	11	0	CGCGT
+chr11	4420	4422	0.0714	14	1	B	G	+	5	0	-	9	1	GCCGA
+chr11	4456	4458	0	14	0	B	G	+	6	0	-	8	0	GGCGC
+chr11	4458	4460	0	14	0	B	G	+	7	0	-	7	0	CGCGG
+chr11	4469	4471	0.0714	14	1	B	G	+	7	1	-	7	0	GCCGG
+chr11	4482	4484	0.0667	15	1	B	G	+	8	1	-	7	0	CTCGA
+chr11	4499	4501	0.0625	16	1	B	G	+	7	0	-	9	1	GGCGT
+chr11	4519	4521	0.105	19	2	B	G	+	8	0	-	11	2	TCCGG
+chr11	4523	4525	0.05	20	1	B	G	+	9	0	-	11	1	GCCGG
+chr11	4532	4534	0	20	0	B	G	+	12	0	-	8	0	AGCGC
+chr11	4549	4551	0.0476	21	1	B	G	+	11	0	-	10	1	TGCGG
+chr11	4558	4560	0	21	0	B	G	+	12	0	-	9	0	ACCGC
+chr11	4588	4590	0.0476	21	1	B	G	+	12	1	-	9	0	CACGT
+chr11	4594	4596	0	20	0	B	G	+	11	0	-	9	0	TCCGG
+chr11	4633	4635	0.0667	15	1	B	G	+	4	1	-	11	0	TTCGC
+chr11	4692	4694	0.0667	15	1	B	G	+	8	1	-	7	0	CCCGG
+chr11	4723	4725	0	16	0	B	G	+	8	0	-	8	0	ACCGA
+chr11	4744	4746	0.0667	15	1	B	G	+	5	1	-	10	0	TTCGG
+chr11	4756	4758	0	19	0	B	G	+	7	0	-	12	0	CACGG
+chr11	4785	4787	0.0526	19	1	B	G	+	4	1	-	15	0	GGCGG
+chr11	4789	4791	0.0476	21	1	B	G	+	5	0	-	16	1	GGCGC
+chr11	4791	4793	0.19	21	4	B	G	+	5	4	-	16	0	CGCGT
+chr11	4799	4801	0.0833	24	2	B	G	+	5	0	-	19	2	GCCGC
+chr11	4845	4847	0.0417	24	1	B	G	+	8	0	-	16	1	AACGT
+chr11	4860	4862	0	24	0	B	G	+	11	0	-	13	0	TCCGC
+chr11	4869	4871	0	24	0	B	G	+	12	0	-	12	0	ACCGC
+chr11	4872	4874	0	26	0	B	G	+	12	0	-	14	0	GCCGG
+chr11	4894	4896	0	25	0	B	G	+	13	0	-	12	0	GACGC
+chr11	4908	4910	0.0385	26	1	B	G	+	13	1	-	13	0	AGCGC
+chr11	4923	4925	0	26	0	B	G	+	14	0	-	12	0	ATCGA
+chr11	4942	4944	0	25	0	B	G	+	16	0	-	9	0	AGCGC
+chr11	5009	5011	0.0476	21	1	B	G	+	15	0	-	6	1	GGCGC
+chr11	5015	5017	0.0526	19	1	B	G	+	14	1	-	5	0	CTCGC
+chr11	5025	5027	0.05	20	1	B	G	+	14	1	-	6	0	GGCGA
+chr11	5036	5038	0.125	16	2	B	G	+	11	2	-	5	0	TGCGC
+chr11	5047	5049	0	17	0	B	G	+	12	0	-	5	0	CTCGA
+chr11	5051	5053	0	14	0	B	G	+	12	0	-	2	0	ACCGC
+chr11	5053	5055	0	16	0	B	G	+	14	0	-	2	0	CGCGA
+chr11	5097	5099	0	20	0	B	G	+	13	0	-	7	0	TCCGC
+chr11	5099	5101	0	18	0	B	G	+	11	0	-	7	0	CGCGC
+chr11	5104	5106	0	18	0	B	G	+	11	0	-	7	0	CTCGC
+chr11	5106	5108	0.0588	17	1	B	G	+	10	1	-	7	0	CGCGC
+chr11	5114	5116	0	18	0	B	G	+	11	0	-	7	0	CCCGG
+chr11	5120	5122	0.0588	17	1	B	G	+	10	1	-	7	0	AACGA
+chr11	5124	5126	0.0588	17	1	B	G	+	10	1	-	7	0	AGCGG
+chr11	5140	5142	0.1	20	2	B	G	+	12	2	-	8	0	GCCGG
+chr11	5144	5146	0.0556	18	1	B	G	+	10	1	-	8	0	GCCGG
+chr11	5149	5151	0.0588	17	1	B	G	+	9	1	-	8	0	CCCGC
+chr11	5151	5153	0.0625	16	1	B	G	+	9	1	-	7	0	CGCGA
+chr11	5157	5159	0.0625	16	1	B	G	+	7	1	-	9	0	CCCGC
+chr11	5159	5161	0.0588	17	1	B	G	+	8	0	-	9	1	CGCGC
+chr11	5166	5168	0.143	14	2	B	G	+	8	1	-	6	1	CTCGG
+chr11	5180	5182	0	18	0	B	G	+	9	0	-	9	0	GCCGG
+chr11	5185	5187	0.0714	14	1	B	G	+	6	1	-	8	0	GCCGC
+chr11	5190	5192	0.25	16	4	B	G	+	7	2	-	9	2	AGCGA
+chr11	5197	5199	0.0667	15	1	B	G	+	6	1	-	9	0	AACGC
+chr11	5209	5211	0.133	15	2	B	G	+	5	1	-	10	1	TGCGC
+chr11	5211	5213	0.0588	17	1	B	G	+	7	1	-	10	0	CGCGC
+chr11	5227	5229	0	18	0	B	G	+	7	0	-	11	0	CCCGC
+chr11	5231	5233	0.0588	17	1	B	G	+	6	0	-	11	1	CCCGA
+chr11	5245	5247	0	18	0	B	G	+	6	0	-	12	0	TGCGG
+chr11	5282	5284	0	15	0	B	G	+	6	0	-	9	0	GCCGC
+chr11	5285	5287	0	17	0	B	G	+	7	0	-	10	0	GCCGA
+chr11	5291	5293	0.0556	18	1	B	G	+	7	0	-	11	1	GGCGA
+chr11	5300	5302	0	21	0	B	G	+	8	0	-	13	0	ACCGG
+chr11	5318	5320	0.136	22	3	B	G	+	8	2	-	14	1	TCCGG
+chr11	5342	5344	0.0476	21	1	B	G	+	7	1	-	14	0	CCCGC
+chr11	5351	5353	0.3	20	6	B	G	+	5	1	-	15	5	CTCGC
+chr11	5379	5381	0.125	16	2	B	G	+	5	1	-	11	1	TACGC
+chr11	5381	5383	0.0625	16	1	B	G	+	6	1	-	10	0	CGCGC
+chr11	5394	5396	0.0667	15	1	B	G	+	6	1	-	9	0	CCCGG
+chr11	5410	5412	0	14	0	B	G	+	5	0	-	9	0	GGCGA
+chr11	5414	5416	0.0714	14	1	B	G	+	6	1	-	8	0	ACCGG
+chr11	5417	5419	0	15	0	B	G	+	6	0	-	9	0	GGCGT
+chr11	5472	5474	0	22	0	B	G	+	11	0	-	11	0	TTCGG
+chr11	5498	5500	0.0435	23	1	B	G	+	13	1	-	10	0	CCCGG
+chr11	5532	5534	0.0333	30	1	B	G	+	14	0	-	16	1	CCCGA
+chr11	5535	5537	0.0323	31	1	B	G	+	16	1	-	15	0	GACGC
+chr11	5538	5540	0.242	33	8	B	G	+	17	5	-	16	3	GCCGT
+chr11	5551	5553	0.138	29	4	B	G	+	16	1	-	13	3	ACCGT
+chr11	5564	5566	0.037	27	1	B	G	+	15	1	-	12	0	ACCGG
+chr11	5576	5578	0.0357	28	1	B	G	+	16	1	-	12	0	TTCGC
+chr11	5600	5602	0.08	25	2	B	G	+	14	2	-	11	0	GTCGG
+chr11	5622	5624	0.0476	21	1	B	G	+	15	1	-	6	0	GCCGC
+chr11	5624	5626	0.2	20	4	B	G	+	15	2	-	5	2	CGCGT
+chr11	5641	5643	0	19	0	B	G	+	13	0	-	6	0	CACGA
+chr11	5667	5669	0	16	0	B	G	+	12	0	-	4	0	CACGG
+chr11	5720	5722	0	22	0	B	G	+	9	0	-	13	0	AGCGC
+chr11	5816	5818	0	14	0	B	G	+	6	0	-	8	0	ACCGT
+chr11	5828	5830	0	12	0	B	G	+	5	0	-	7	0	CCCGC
+chr11	5830	5832	0	13	0	B	G	+	5	0	-	8	0	CGCGC
+chr11	5843	5845	0.0909	11	1	B	G	+	3	0	-	8	1	AGCGA
+chr11	5905	5907	0.333	3	1	B	G	+	1	0	-	2	1	CCCGT
+chr11	5969	5971	0.3	10	3	B	G	+	8	2	-	2	1	CCCGC
+chr11	5972	5974	0.0909	11	1	B	G	+	8	1	-	3	0	GCCGG
+chr11	6023	6025	0.0909	11	1	B	G	+	7	0	-	4	1	CTCGC
+chr11	6044	6046	0	10	0	B	G	+	6	0	-	4	0	TTCGC
+chr11	6049	6051	0	10	0	B	G	+	6	0	-	4	0	CCCGC
+chr11	6076	6078	0.0909	11	1	B	G	+	6	0	-	5	1	TCCGG
+chr11	6089	6091	0	9	0	B	G	+	5	0	-	4	0	CCCGG
+chr11	6099	6101	0	6	0	B	G	+	3	0	-	3	0	GGCGA
+chr11	6148	6150	0	9	0	B	G	+	4	0	-	5	0	ACCGC
+chr11	6156	6158	0	8	0	B	G	+	4	0	-	4	0	TTCGA
+chr11	6164	6166	0	8	0	B	G	+	3	0	-	5	0	CTCGG
+chr11	6168	6170	0.125	8	1	B	G	+	3	0	-	5	1	GGCGG
+chr11	6208	6210	0.1	10	1	B	G	+	4	1	-	6	0	ATCGA
+chr11	6226	6228	0	10	0	B	G	+	6	0	-	4	0	AACGC
+chr11	6239	6241	0	5	0	B	G	+	3	0	-	2	0	CCCGG
+chr11	6247	6249	0.0909	11	1	B	G	+	6	1	-	5	0	GCCGG
+chr11	6251	6253	0.0769	13	1	B	G	+	7	0	-	6	1	GGCGC
+chr11	6253	6255	0.0769	13	1	B	G	+	7	1	-	6	0	CGCGG
+chr11	6258	6260	0	12	0	B	G	+	7	0	-	5	0	AGCGG
+chr11	6269	6271	0	15	0	B	G	+	7	0	-	8	0	GCCGC
+chr11	6272	6274	0.2	15	3	B	G	+	7	3	-	8	0	GCCGC
+chr11	6282	6284	0.0667	15	1	B	G	+	7	0	-	8	1	CCCGG
+chr11	6293	6295	0	16	0	B	G	+	8	0	-	8	0	GTCGG
+chr11	6297	6299	0	18	0	B	G	+	10	0	-	8	0	GTCGG
+chr11	6300	6302	0.133	15	2	B	G	+	8	2	-	7	0	GGCGG
+chr11	6304	6306	0.118	17	2	B	G	+	10	2	-	7	0	GCCGG
+chr11	6312	6314	0.05	20	1	B	G	+	11	1	-	9	0	AGCGC
+chr11	6314	6316	0.05	20	1	B	G	+	10	1	-	10	0	CGCGG
+chr11	6318	6320	0.05	20	1	B	G	+	10	1	-	10	0	GGCGT
+chr11	6321	6323	0	19	0	B	G	+	9	0	-	10	0	GTCGC
+chr11	6323	6325	0	20	0	B	G	+	11	0	-	9	0	CGCGA
+chr11	6350	6352	0.0556	18	1	B	G	+	10	1	-	8	0	CTCGG
+chr11	6356	6358	0.0526	19	1	B	G	+	10	1	-	9	0	TCCGC
+chr11	6358	6360	0.05	20	1	B	G	+	11	0	-	9	1	CGCGC
+chr11	6365	6367	0	19	0	B	G	+	11	0	-	8	0	CCCGA
+chr11	6371	6373	0.0556	18	1	B	G	+	11	1	-	7	0	CCCGA
+chr11	6430	6432	0	14	0	B	G	+	5	0	-	9	0	GCCGC
+chr11	6432	6434	0.0625	16	1	B	G	+	6	1	-	10	0	CGCGG
+chr11	6442	6444	0.1	20	2	B	G	+	9	2	-	11	0	CACGG
+chr11	6449	6451	0	19	0	B	G	+	10	0	-	9	0	CTCGA
+chr11	6491	6493	0.04	25	1	B	G	+	13	0	-	12	1	CTCGG
+chr11	6526	6528	0.0286	35	1	B	G	+	15	0	-	20	1	ACCGA
+chr11	6560	6562	0.24	25	6	B	G	+	7	2	-	18	4	TCCGT
+chr11	6596	6598	0.421	38	16	B	G	+	12	3	-	26	13	GGCGA
+chr11	6730	6732	0.366	41	15	B	G	+	15	7	-	26	8	AGCGA
+chr11	6857	6859	0.212	33	7	B	G	+	16	2	-	17	5	CTCGG
+chr11	6900	6902	0.862	29	25	B	G	+	14	12	-	15	13	TGCGC
+chr11	6902	6904	0.483	29	14	B	G	+	15	6	-	14	8	CGCGC
+chr11	6979	6981	0.28	25	7	B	G	+	14	3	-	11	4	TACGT
+chr11	7085	7087	0.346	26	9	B	G	+	15	6	-	11	3	AGCGG
+chr11	7166	7168	0.263	19	5	B	G	+	13	3	-	6	2	TTCGG
+chr11	7336	7338	0.294	17	5	B	G	+	9	4	-	8	1	GCCGG
+chr11	7583	7585	0.154	26	4	B	G	+	5	0	-	21	4	GGCGT
+chr11	7633	7635	0.0952	21	2	B	G	+	2	0	-	19	2	CTCGC
+chr11	7655	7657	0	15	0	B	G	+	2	0	-	13	0	TACGA
+chr11	7673	7675	0	16	0	B	G	+	2	0	-	14	0	AGCGT
+chr11	7679	7681	0	17	0	B	G	+	4	0	-	13	0	GTCGA
+chr11	7688	7690	0	15	0	B	G	+	4	0	-	11	0	GCCGA
+chr11	7695	7697	0.0667	15	1	B	G	+	4	0	-	11	1	GACGG
+chr11	7729	7731	0	14	0	B	G	+	6	0	-	8	0	CTCGA
+chr11	7748	7750	0	15	0	B	G	+	5	0	-	10	0	ACCGA
+chr11	7825	7827	0	18	0	B	G	+	8	0	-	10	0	GACGC
+chr11	7866	7868	0	17	0	B	G	+	10	0	-	7	0	CCCGA
+chr11	7878	7880	0.0667	15	1	B	G	+	9	1	-	6	0	CCCGT
+chr11	7882	7884	0.0769	13	1	B	G	+	8	1	-	5	0	TACGA
+chr11	8027	8029	0.0385	26	1	B	G	+	22	1	-	4	0	AACGG
+chr11	8033	8035	0.037	27	1	B	G	+	22	1	-	5	0	ACCGT
+chr11	8303	8305	0.0556	36	2	B	G	+	9	1	-	27	1	ACCGA
+chr11	8414	8416	0.0625	16	1	B	G	+	5	0	-	11	1	CCCGG
+chr11	8583	8585	0.0833	24	2	B	G	+	16	2	-	8	0	CACGC
+chr11	8678	8680	0	19	0	B	G	+	9	0	-	10	0	CTCGT
+chr11	8692	8694	0.0435	23	1	B	G	+	10	0	-	13	1	TTCGA
+chr11	8755	8757	0	26	0	B	G	+	9	0	-	17	0	CTCGC
+chr11	8786	8788	0	25	0	B	G	+	12	0	-	13	0	AACGT
+chr11	8812	8814	0.0909	22	2	B	G	+	10	2	-	12	0	CCCGG
+chr11	8829	8831	0	20	0	B	G	+	11	0	-	9	0	ATCGA
+chr11	8836	8838	0	21	0	B	G	+	11	0	-	10	0	CCCGG
+chr11	8855	8857	0.0588	17	1	B	G	+	8	0	-	9	1	ACCGG
+chr11	8902	8904	0	17	0	B	G	+	11	0	-	6	0	CCCGA
+chr11	8914	8916	0.35	20	7	B	G	+	10	3	-	10	4	TACGG
+chr11	8920	8922	0	22	0	B	G	+	12	0	-	10	0	CTCGC
+chr11	8948	8950	0	33	0	B	G	+	12	0	-	21	0	AACGA
+chr11	8975	8977	0.0476	42	2	B	G	+	12	0	-	30	2	GTCGC
+chr11	9127	9129	0.0244	41	1	B	G	+	18	0	-	23	1	CTCGG
+chr11	9136	9138	0.0476	42	2	B	G	+	18	1	-	24	1	CCCGC
+chr11	9227	9229	0.0444	45	2	B	G	+	18	0	-	27	2	CACGC
+chr11	9270	9272	0.025	40	1	B	G	+	15	0	-	25	1	ACCGC
+chr11	9301	9303	0.175	40	7	B	G	+	19	3	-	21	4	CTCGC
+chr11	9371	9373	0	25	0	B	G	+	16	0	-	9	0	TTCGA
+chr11	9380	9382	0.04	25	1	B	G	+	16	1	-	9	0	AACGT
+chr11	9447	9449	0	17	0	B	G	+	13	0	-	4	0	CACGT
+chr11	9478	9480	0	16	0	B	G	+	11	0	-	5	0	ATCGC
+chr11	9510	9512	0	14	0	B	G	+	7	0	-	7	0	CTCGA
+chr11	9588	9590	0	20	0	B	G	+	14	0	-	6	0	GGCGC
+chr11	9677	9679	0.0909	11	1	B	G	+	7	0	-	4	1	ATCGA
+chr11	9697	9699	0	14	0	B	G	+	10	0	-	4	0	TCCGA
+chr11	9709	9711	0	11	0	B	G	+	7	0	-	4	0	GCCGG
+chr11	9729	9731	0.167	6	1	B	G	+	5	1	-	1	0	GACGA
Binary file test-data/SE_g1_r1.bam has changed
Binary file test-data/SE_g1_r2.bam has changed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/SE_g2.G.bed	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,310 @@
+#chrom	start	end	ratio	totalC	methC	strand	next	Plus	totalC	methC	Minus	totalC	methC	localSeq
+chr11	41	43	0	3	0	-	G	+	0	0	-	3	0	GCGGC
+chr11	48	50	0	3	0	-	G	+	0	0	-	3	0	CCGGC
+chr11	56	58	0	4	0	-	G	+	0	0	-	4	0	CCGCC
+chr11	59	61	0.4	5	2	B	G	+	1	1	-	4	1	GCCGT
+chr11	101	103	0.143	7	1	B	G	+	2	0	-	5	1	GTCGC
+chr11	140	142	0	9	0	B	G	+	4	0	-	5	0	GGCGA
+chr11	144	146	0.1	10	1	B	G	+	3	1	-	7	0	AACGG
+chr11	210	212	0.417	12	5	B	G	+	3	2	-	9	3	CCCGA
+chr11	215	217	0.25	16	4	B	G	+	3	2	-	13	2	CTCGC
+chr11	316	318	0.846	26	22	B	G	+	5	4	-	21	18	CTCGC
+chr11	331	333	0.808	26	21	B	G	+	3	2	-	23	19	TGCGT
+chr11	400	402	0.412	34	14	B	G	+	5	3	-	29	11	CACGC
+chr11	407	409	0.588	34	20	B	G	+	6	5	-	28	15	CTCGT
+chr11	445	447	0.606	33	20	B	G	+	9	6	-	24	14	CACGC
+chr11	469	471	0.588	34	20	B	G	+	9	5	-	25	15	ACCGC
+chr11	546	548	0.812	32	26	B	G	+	16	10	-	16	16	GTCGA
+chr11	666	668	0.553	38	21	B	G	+	22	13	-	16	8	TGCGT
+chr11	694	696	0.902	41	37	B	G	+	29	27	-	12	10	TGCGT
+chr11	712	714	0.881	42	37	B	G	+	23	20	-	19	17	ACCGA
+chr11	721	723	0.923	39	36	B	G	+	20	20	-	19	16	CCCGG
+chr11	737	739	0.622	37	23	B	G	+	21	13	-	16	10	GGCGG
+chr11	753	755	0.78	41	32	B	G	+	20	19	-	21	13	TCCGG
+chr11	762	764	0.619	42	26	B	G	+	21	16	-	21	10	CCCGT
+chr11	792	794	0.842	38	32	B	G	+	15	11	-	23	21	GCCGT
+chr11	798	800	0.875	40	35	B	G	+	16	13	-	24	22	TGCGC
+chr11	816	818	0.912	34	31	B	G	+	15	13	-	19	18	TTCGC
+chr11	830	832	0.914	35	32	B	G	+	17	15	-	18	17	TGCGA
+chr11	954	956	0.485	33	16	B	G	+	15	8	-	18	8	CCCGT
+chr11	977	979	0.703	37	26	B	G	+	21	16	-	16	10	CCCGG
+chr11	1006	1008	0.526	38	20	B	G	+	21	13	-	17	7	ACCGT
+chr11	1066	1068	0.756	45	34	B	G	+	29	21	-	16	13	GGCGG
+chr11	1116	1118	0.771	35	27	B	G	+	24	19	-	11	8	CACGG
+chr11	1145	1147	0.538	26	14	B	G	+	15	10	-	11	4	CCCGA
+chr11	1260	1262	0.375	24	9	B	G	+	12	5	-	12	4	GACGA
+chr11	1320	1322	0.64	25	16	B	G	+	13	11	-	12	5	TCCGC
+chr11	1559	1561	0.852	27	23	B	G	+	15	13	-	12	10	TACGC
+chr11	1574	1576	0.607	28	17	B	G	+	14	9	-	14	8	ATCGT
+chr11	1592	1594	0.697	33	23	B	G	+	14	7	-	19	16	CTCGG
+chr11	1609	1611	0.244	41	10	B	G	+	18	5	-	23	5	CTCGG
+chr11	1632	1634	0.2	45	9	B	G	+	20	5	-	25	4	GACGG
+chr11	1642	1644	0.268	41	11	B	G	+	18	6	-	23	5	CCCGC
+chr11	1680	1682	0.611	54	33	B	G	+	20	12	-	34	21	GTCGC
+chr11	1689	1691	0.471	51	24	B	G	+	20	8	-	31	16	CTCGG
+chr11	1703	1705	0.755	49	37	B	G	+	17	12	-	32	25	ACCGA
+chr11	1733	1735	0.745	51	38	B	G	+	20	15	-	31	23	CTCGT
+chr11	1766	1768	0.312	48	15	B	G	+	19	4	-	29	11	TCCGT
+chr11	1869	1871	0.263	38	10	B	G	+	21	5	-	17	5	TACGT
+chr11	2040	2042	0	34	0	B	G	+	16	0	-	18	0	AACGA
+chr11	2104	2106	0.0513	39	2	B	G	+	15	0	-	24	2	AACGA
+chr11	2148	2150	0.588	34	20	B	G	+	12	3	-	22	17	GGCGT
+chr11	2360	2362	0.538	26	14	B	G	+	22	11	-	4	3	GCCGG
+chr11	2384	2386	0.593	27	16	B	G	+	21	11	-	6	5	TACGG
+chr11	2454	2456	0.107	28	3	B	G	+	23	2	-	5	1	TCCGC
+chr11	2621	2623	0	12	0	B	G	+	11	0	-	1	0	ATCGT
+chr11	2648	2650	0	20	0	B	G	+	17	0	-	3	0	GGCGT
+chr11	2665	2667	0.105	19	2	B	G	+	16	0	-	3	2	GCCGC
+chr11	2667	2669	0.1	20	2	B	G	+	17	0	-	3	2	CGCGG
+chr11	2673	2675	0.19	21	4	B	G	+	18	2	-	3	2	GGCGG
+chr11	2676	2678	0.0385	26	1	B	G	+	23	1	-	3	0	GGCGA
+chr11	2685	2687	0.143	28	4	B	G	+	25	2	-	3	2	GGCGG
+chr11	2688	2690	0	28	0	B	G	+	24	0	-	4	0	GGCGC
+chr11	2691	2693	0.276	29	8	B	G	+	24	5	-	5	3	GCCGG
+chr11	2713	2715	0.154	26	4	B	G	+	20	3	-	6	1	TTCGG
+chr11	2723	2725	0.12	25	3	B	G	+	20	3	-	5	0	GTCGC
+chr11	2727	2729	0.35	20	7	B	G	+	16	6	-	4	1	CACGG
+chr11	2731	2733	0.211	19	4	B	G	+	15	3	-	4	1	GCCGG
+chr11	2738	2740	0.273	22	6	B	G	+	18	4	-	4	2	CCCGC
+chr11	2767	2769	0.333	21	7	B	G	+	14	5	-	7	2	GCCGG
+chr11	2778	2780	0.455	22	10	B	G	+	14	5	-	8	5	TCCGG
+chr11	2788	2790	0.167	18	3	B	G	+	14	3	-	4	0	GCCGC
+chr11	2792	2794	0.318	22	7	B	G	+	16	4	-	6	3	CCCGT
+chr11	2810	2812	0.105	19	2	B	G	+	14	1	-	5	1	ACCGA
+chr11	2825	2827	0.133	15	2	B	G	+	11	2	-	4	0	TCCGC
+chr11	3031	3033	0.4	15	6	B	G	+	11	4	-	4	2	CTCGG
+chr11	3034	3036	0.286	14	4	B	G	+	10	2	-	4	2	GGCGG
+chr11	3052	3054	0.769	13	10	B	G	+	7	5	-	6	5	ATCGA
+chr11	3095	3097	0.8	10	8	B	G	+	5	5	-	5	3	ATCGG
+chr11	3117	3119	0.571	14	8	B	G	+	4	2	-	10	6	ACCGT
+chr11	3273	3275	0.48	25	12	B	G	+	9	3	-	16	9	AGCGC
+chr11	3283	3285	0.478	23	11	B	G	+	9	2	-	14	9	ACCGG
+chr11	3323	3325	0.625	24	15	B	G	+	10	6	-	14	9	ACCGC
+chr11	3424	3426	0.0952	21	2	B	G	+	12	1	-	9	1	GACGC
+chr11	3452	3454	0.588	17	10	B	G	+	14	9	-	3	1	ACCGT
+chr11	3484	3486	0.556	18	10	B	G	+	14	7	-	4	3	ATCGA
+chr11	3578	3580	0.737	19	14	B	G	+	12	8	-	7	6	TCCGC
+chr11	3595	3597	0.238	21	5	B	G	+	14	3	-	7	2	GACGT
+chr11	3611	3613	0.826	23	19	B	G	+	16	12	-	7	7	CCCGG
+chr11	3617	3619	0.762	21	16	B	G	+	13	9	-	8	7	CTCGG
+chr11	3639	3641	0.6	20	12	B	G	+	10	6	-	10	6	TCCGG
+chr11	3646	3648	0.391	23	9	B	G	+	10	6	-	13	3	ACCGC
+chr11	3665	3667	0.571	21	12	B	G	+	10	7	-	11	5	GCCGG
+chr11	3672	3674	0.348	23	8	B	G	+	10	6	-	13	2	GTCGG
+chr11	3677	3679	0.261	23	6	B	G	+	10	3	-	13	3	TTCGC
+chr11	3693	3695	0.231	26	6	B	G	+	10	2	-	16	4	CTCGC
+chr11	3704	3706	0.2	25	5	B	G	+	9	1	-	16	4	TCCGC
+chr11	3707	3709	0.16	25	4	B	G	+	9	1	-	16	3	GCCGC
+chr11	3717	3719	0.286	21	6	B	G	+	7	2	-	14	4	TCCGC
+chr11	3960	3962	0.167	18	3	B	G	+	7	2	-	11	1	CCCGA
+chr11	3991	3993	0.167	24	4	B	G	+	9	4	-	15	0	GACGT
+chr11	4056	4058	0.292	24	7	B	G	+	13	4	-	11	3	GACGC
+chr11	4163	4165	0.0333	30	1	B	G	+	17	0	-	13	1	CACGT
+chr11	4195	4197	0.231	26	6	B	G	+	16	4	-	10	2	GGCGG
+chr11	4202	4204	0.0667	30	2	B	G	+	19	1	-	11	1	GTCGC
+chr11	4207	4209	0.194	31	6	B	G	+	19	4	-	12	2	TCCGG
+chr11	4243	4245	0.111	27	3	B	G	+	17	2	-	10	1	GTCGG
+chr11	4248	4250	0.0323	31	1	B	G	+	19	1	-	12	0	GGCGC
+chr11	4284	4286	0.226	31	7	B	G	+	17	5	-	14	2	TCCGC
+chr11	4289	4291	0.0741	27	2	B	G	+	14	1	-	13	1	TCCGG
+chr11	4302	4304	0	30	0	B	G	+	16	0	-	14	0	ATCGT
+chr11	4307	4309	0.0294	34	1	B	G	+	20	0	-	14	1	CACGA
+chr11	4314	4316	0.0606	33	2	B	G	+	19	2	-	14	0	AACGG
+chr11	4320	4322	0.0882	34	3	B	G	+	20	1	-	14	2	ACCGT
+chr11	4330	4332	0.0278	36	1	B	G	+	21	0	-	15	1	AACGA
+chr11	4340	4342	0.0323	31	1	B	G	+	18	1	-	13	0	TCCGT
+chr11	4365	4367	0.0303	33	1	B	G	+	19	1	-	14	0	GACGA
+chr11	4393	4395	0.143	28	4	B	G	+	15	2	-	13	2	TTCGA
+chr11	4403	4405	0	25	0	B	G	+	11	0	-	14	0	GCCGC
+chr11	4405	4407	0	24	0	B	G	+	10	0	-	14	0	CGCGT
+chr11	4420	4422	0.136	22	3	B	G	+	8	1	-	14	2	GCCGA
+chr11	4456	4458	0	15	0	B	G	+	2	0	-	13	0	GGCGC
+chr11	4458	4460	0.0625	16	1	B	G	+	3	0	-	13	1	CGCGG
+chr11	4469	4471	0.154	13	2	B	G	+	2	0	-	11	2	GCCGG
+chr11	4482	4484	0	19	0	B	G	+	5	0	-	14	0	CTCGA
+chr11	4499	4501	0.0417	24	1	B	G	+	9	0	-	15	1	GGCGT
+chr11	4519	4521	0.304	23	7	B	G	+	9	2	-	14	5	TCCGG
+chr11	4523	4525	0.28	25	7	B	G	+	11	3	-	14	4	GCCGG
+chr11	4532	4534	0.174	23	4	B	G	+	10	0	-	13	4	AGCGC
+chr11	4549	4551	0.231	26	6	B	G	+	9	3	-	17	3	TGCGG
+chr11	4558	4560	0.2	25	5	B	G	+	8	1	-	17	4	ACCGC
+chr11	4588	4590	0	15	0	B	G	+	5	0	-	10	0	CACGT
+chr11	4594	4596	0.0556	18	1	B	G	+	5	0	-	13	1	TCCGG
+chr11	4633	4635	0.263	19	5	B	G	+	4	0	-	15	5	TTCGC
+chr11	4692	4694	0.188	16	3	B	G	+	4	1	-	12	2	CCCGG
+chr11	4723	4725	0.235	17	4	B	G	+	7	0	-	10	4	ACCGA
+chr11	4744	4746	0.0952	21	2	B	G	+	8	1	-	13	1	TTCGG
+chr11	4756	4758	0.208	24	5	B	G	+	7	3	-	17	2	CACGG
+chr11	4785	4787	0.182	22	4	B	G	+	5	0	-	17	4	GGCGG
+chr11	4789	4791	0.0833	24	2	B	G	+	7	1	-	17	1	GGCGC
+chr11	4791	4793	0.261	23	6	B	G	+	7	2	-	16	4	CGCGT
+chr11	4799	4801	0	25	0	B	G	+	10	0	-	15	0	GCCGC
+chr11	4845	4847	0	23	0	B	G	+	10	0	-	13	0	AACGT
+chr11	4860	4862	0	17	0	B	G	+	9	0	-	8	0	TCCGC
+chr11	4869	4871	0.238	21	5	B	G	+	11	4	-	10	1	ACCGC
+chr11	4872	4874	0.1	20	2	B	G	+	10	0	-	10	2	GCCGG
+chr11	4894	4896	0.0833	24	2	B	G	+	13	2	-	11	0	GACGC
+chr11	4908	4910	0.0741	27	2	B	G	+	15	2	-	12	0	AGCGC
+chr11	4923	4925	0	28	0	B	G	+	17	0	-	11	0	ATCGA
+chr11	4942	4944	0.136	22	3	B	G	+	12	1	-	10	2	AGCGC
+chr11	5009	5011	0.08	25	2	B	G	+	11	1	-	14	1	GGCGC
+chr11	5015	5017	0.107	28	3	B	G	+	13	1	-	15	2	CTCGC
+chr11	5025	5027	0.192	26	5	B	G	+	12	3	-	14	2	GGCGA
+chr11	5036	5038	0.04	25	1	B	G	+	12	0	-	13	1	TGCGC
+chr11	5047	5049	0	22	0	B	G	+	11	0	-	11	0	CTCGA
+chr11	5051	5053	0.0476	21	1	B	G	+	11	1	-	10	0	ACCGC
+chr11	5053	5055	0	21	0	B	G	+	11	0	-	10	0	CGCGA
+chr11	5097	5099	0	24	0	B	G	+	13	0	-	11	0	TCCGC
+chr11	5099	5101	0.04	25	1	B	G	+	14	0	-	11	1	CGCGC
+chr11	5104	5106	0.0435	23	1	B	G	+	13	1	-	10	0	CTCGC
+chr11	5106	5108	0	22	0	B	G	+	13	0	-	9	0	CGCGC
+chr11	5114	5116	0	19	0	B	G	+	11	0	-	8	0	CCCGG
+chr11	5120	5122	0.1	20	2	B	G	+	12	2	-	8	0	AACGA
+chr11	5124	5126	0	18	0	B	G	+	10	0	-	8	0	AGCGG
+chr11	5140	5142	0.0526	19	1	B	G	+	11	1	-	8	0	GCCGG
+chr11	5144	5146	0	18	0	B	G	+	11	0	-	7	0	GCCGG
+chr11	5149	5151	0	16	0	B	G	+	10	0	-	6	0	CCCGC
+chr11	5151	5153	0	15	0	B	G	+	9	0	-	6	0	CGCGA
+chr11	5157	5159	0.133	15	2	B	G	+	10	0	-	5	2	CCCGC
+chr11	5159	5161	0.0667	15	1	B	G	+	8	0	-	7	1	CGCGC
+chr11	5166	5168	0	16	0	B	G	+	8	0	-	8	0	CTCGG
+chr11	5180	5182	0.0769	13	1	B	G	+	8	0	-	5	1	GCCGG
+chr11	5185	5187	0	12	0	B	G	+	7	0	-	5	0	GCCGC
+chr11	5190	5192	0	10	0	B	G	+	6	0	-	4	0	AGCGA
+chr11	5197	5199	0	10	0	B	G	+	5	0	-	5	0	AACGC
+chr11	5209	5211	0.111	9	1	B	G	+	5	0	-	4	1	TGCGC
+chr11	5211	5213	0	9	0	B	G	+	6	0	-	3	0	CGCGC
+chr11	5227	5229	0.111	9	1	B	G	+	6	1	-	3	0	CCCGC
+chr11	5231	5233	0.0909	11	1	B	G	+	7	0	-	4	1	CCCGA
+chr11	5245	5247	0	13	0	B	G	+	7	0	-	6	0	TGCGG
+chr11	5282	5284	0	12	0	B	G	+	4	0	-	8	0	GCCGC
+chr11	5285	5287	0.167	12	2	B	G	+	4	1	-	8	1	GCCGA
+chr11	5291	5293	0.0769	13	1	B	G	+	4	0	-	9	1	GGCGA
+chr11	5300	5302	0	13	0	B	G	+	4	0	-	9	0	ACCGG
+chr11	5318	5320	0.364	11	4	B	G	+	2	0	-	9	4	TCCGG
+chr11	5342	5344	0.167	12	2	B	G	+	3	0	-	9	2	CCCGC
+chr11	5351	5353	0.429	14	6	B	G	+	4	0	-	10	6	CTCGC
+chr11	5379	5381	0	16	0	B	G	+	6	0	-	10	0	TACGC
+chr11	5381	5383	0	16	0	B	G	+	6	0	-	10	0	CGCGC
+chr11	5394	5396	0.0667	15	1	B	G	+	6	0	-	9	1	CCCGG
+chr11	5410	5412	0	21	0	B	G	+	10	0	-	11	0	GGCGA
+chr11	5414	5416	0.143	21	3	B	G	+	10	2	-	11	1	ACCGG
+chr11	5417	5419	0	20	0	B	G	+	9	0	-	11	0	GGCGT
+chr11	5472	5474	0.174	23	4	B	G	+	13	4	-	10	0	TTCGG
+chr11	5498	5500	0.143	21	3	B	G	+	12	3	-	9	0	CCCGG
+chr11	5532	5534	0.115	26	3	B	G	+	16	2	-	10	1	CCCGA
+chr11	5535	5537	0.231	26	6	B	G	+	17	6	-	9	0	GACGC
+chr11	5538	5540	0.308	26	8	B	G	+	17	7	-	9	1	GCCGT
+chr11	5551	5553	0.348	23	8	B	G	+	14	5	-	9	3	ACCGT
+chr11	5564	5566	0	19	0	B	G	+	12	0	-	7	0	ACCGG
+chr11	5576	5578	0.0476	21	1	B	G	+	15	0	-	6	1	TTCGC
+chr11	5600	5602	0.0833	24	2	B	G	+	18	1	-	6	1	GTCGG
+chr11	5622	5624	0.211	19	4	B	G	+	15	1	-	4	3	GCCGC
+chr11	5624	5626	0.235	17	4	B	G	+	13	4	-	4	0	CGCGT
+chr11	5641	5643	0.0588	17	1	B	G	+	16	1	-	1	0	CACGA
+chr11	5667	5669	0.0417	24	1	B	G	+	22	1	-	2	0	CACGG
+chr11	5720	5722	0.125	24	3	B	G	+	20	2	-	4	1	AGCGC
+chr11	5816	5818	0.238	21	5	B	G	+	11	5	-	10	0	ACCGT
+chr11	5828	5830	0.238	21	5	B	G	+	11	4	-	10	1	CCCGC
+chr11	5830	5832	0.0455	22	1	B	G	+	11	0	-	11	1	CGCGC
+chr11	5843	5845	0.0556	18	1	B	G	+	8	1	-	10	0	AGCGA
+chr11	5905	5907	0	2	0	B	G	+	1	0	-	1	0	CCCGT
+chr11	5969	5971	0.267	15	4	B	G	+	7	2	-	8	2	CCCGC
+chr11	5972	5974	0	15	0	B	G	+	7	0	-	8	0	GCCGG
+chr11	6023	6025	0.333	18	6	B	G	+	8	1	-	10	5	CTCGC
+chr11	6044	6046	0	12	0	B	G	+	6	0	-	6	0	TTCGC
+chr11	6049	6051	0.0833	12	1	B	G	+	6	1	-	6	0	CCCGC
+chr11	6076	6078	0	7	0	B	G	+	1	0	-	6	0	TCCGG
+chr11	6089	6091	0	10	0	B	G	+	3	0	-	7	0	CCCGG
+chr11	6099	6101	0.111	9	1	B	G	+	3	0	-	6	1	GGCGA
+chr11	6148	6150	0	13	0	B	G	+	5	0	-	8	0	ACCGC
+chr11	6156	6158	0	14	0	B	G	+	6	0	-	8	0	TTCGA
+chr11	6164	6166	0	12	0	B	G	+	7	0	-	5	0	CTCGG
+chr11	6168	6170	0	14	0	B	G	+	7	0	-	7	0	GGCGG
+chr11	6208	6210	0.2	10	2	B	G	+	6	1	-	4	1	ATCGA
+chr11	6226	6228	0	10	0	B	G	+	6	0	-	4	0	AACGC
+chr11	6239	6241	0	12	0	B	G	+	7	0	-	5	0	CCCGG
+chr11	6247	6249	0.0909	11	1	B	G	+	6	1	-	5	0	GCCGG
+chr11	6251	6253	0.222	9	2	B	G	+	4	1	-	5	1	GGCGC
+chr11	6253	6255	0.111	9	1	B	G	+	4	0	-	5	1	CGCGG
+chr11	6258	6260	0.0909	11	1	B	G	+	5	0	-	6	1	AGCGG
+chr11	6269	6271	0.0833	12	1	B	G	+	5	0	-	7	1	GCCGC
+chr11	6272	6274	0.0909	11	1	B	G	+	4	0	-	7	1	GCCGC
+chr11	6282	6284	0.1	10	1	B	G	+	5	0	-	5	1	CCCGG
+chr11	6293	6295	0.143	14	2	B	G	+	6	2	-	8	0	GTCGG
+chr11	6297	6299	0	13	0	B	G	+	6	0	-	7	0	GTCGG
+chr11	6300	6302	0.154	13	2	B	G	+	6	1	-	7	1	GGCGG
+chr11	6304	6306	0.0833	12	1	B	G	+	5	1	-	7	0	GCCGG
+chr11	6312	6314	0	12	0	B	G	+	5	0	-	7	0	AGCGC
+chr11	6314	6316	0.167	12	2	B	G	+	5	1	-	7	1	CGCGG
+chr11	6318	6320	0	12	0	B	G	+	5	0	-	7	0	GGCGT
+chr11	6321	6323	0	13	0	B	G	+	5	0	-	8	0	GTCGC
+chr11	6323	6325	0	12	0	B	G	+	5	0	-	7	0	CGCGA
+chr11	6350	6352	0	12	0	B	G	+	5	0	-	7	0	CTCGG
+chr11	6356	6358	0.1	10	1	B	G	+	5	1	-	5	0	TCCGC
+chr11	6358	6360	0	11	0	B	G	+	5	0	-	6	0	CGCGC
+chr11	6365	6367	0	11	0	B	G	+	6	0	-	5	0	CCCGA
+chr11	6371	6373	0.0769	13	1	B	G	+	6	1	-	7	0	CCCGA
+chr11	6430	6432	0.105	19	2	B	G	+	11	0	-	8	2	GCCGC
+chr11	6432	6434	0.105	19	2	B	G	+	10	0	-	9	2	CGCGG
+chr11	6442	6444	0.0909	22	2	B	G	+	12	0	-	10	2	CACGG
+chr11	6449	6451	0.0417	24	1	B	G	+	14	1	-	10	0	CTCGA
+chr11	6491	6493	0.278	18	5	B	G	+	10	2	-	8	3	CTCGG
+chr11	6526	6528	0.115	26	3	B	G	+	12	2	-	14	1	ACCGA
+chr11	6560	6562	0.233	30	7	B	G	+	15	1	-	15	6	TCCGT
+chr11	6596	6598	0.6	30	18	B	G	+	14	7	-	16	11	GGCGA
+chr11	6730	6732	0.516	31	16	B	G	+	9	5	-	22	11	AGCGA
+chr11	6857	6859	0.239	46	11	B	G	+	20	6	-	26	5	CTCGG
+chr11	6900	6902	0.804	46	37	B	G	+	22	19	-	24	18	TGCGC
+chr11	6902	6904	0.804	46	37	B	G	+	22	18	-	24	19	CGCGC
+chr11	6979	6981	0.487	39	19	B	G	+	20	11	-	19	8	TACGT
+chr11	7085	7087	0.704	27	19	B	G	+	14	10	-	13	9	AGCGG
+chr11	7166	7168	0.259	27	7	B	G	+	19	7	-	8	0	TTCGG
+chr11	7336	7338	0.769	13	10	B	G	+	11	9	-	2	1	GCCGG
+chr11	7583	7585	0.176	34	6	B	G	+	7	1	-	27	5	GGCGT
+chr11	7633	7635	0.154	26	4	B	G	+	7	3	-	19	1	CTCGC
+chr11	7655	7657	0.143	21	3	B	G	+	6	0	-	15	3	TACGA
+chr11	7673	7675	0	17	0	B	G	+	4	0	-	13	0	AGCGT
+chr11	7679	7681	0.0476	21	1	B	G	+	7	0	-	14	1	GTCGA
+chr11	7688	7690	0.167	18	3	B	G	+	8	0	-	10	3	GCCGA
+chr11	7695	7697	0	21	0	B	G	+	7	0	-	14	0	GACGG
+chr11	7729	7731	0.19	21	4	B	G	+	7	1	-	14	3	CTCGA
+chr11	7748	7750	0.0435	23	1	B	G	+	8	0	-	15	1	ACCGA
+chr11	7825	7827	0.05	20	1	B	G	+	11	0	-	9	1	GACGC
+chr11	7866	7868	0.0455	22	1	B	G	+	12	0	-	10	1	CCCGA
+chr11	7878	7880	0.0952	21	2	B	G	+	13	1	-	8	1	CCCGT
+chr11	7882	7884	0.3	20	6	B	G	+	12	5	-	8	1	TACGA
+chr11	8027	8029	0.167	18	3	B	G	+	13	2	-	5	1	AACGG
+chr11	8033	8035	0.0625	16	1	B	G	+	11	1	-	5	0	ACCGT
+chr11	8303	8305	0.0435	23	1	B	G	+	10	0	-	13	1	ACCGA
+chr11	8414	8416	0.444	27	12	B	G	+	11	4	-	16	8	CCCGG
+chr11	8583	8585	0.148	27	4	B	G	+	11	1	-	16	3	CACGC
+chr11	8678	8680	0	17	0	B	G	+	8	0	-	9	0	CTCGT
+chr11	8692	8694	0.0526	19	1	B	G	+	8	0	-	11	1	TTCGA
+chr11	8755	8757	0.0417	24	1	B	G	+	3	0	-	21	1	CTCGC
+chr11	8786	8788	0	19	0	B	G	+	1	0	-	18	0	AACGT
+chr11	8812	8814	0	21	0	B	G	+	4	0	-	17	0	CCCGG
+chr11	8829	8831	0.136	22	3	B	G	+	5	0	-	17	3	ATCGA
+chr11	8836	8838	0.0476	21	1	B	G	+	4	0	-	17	1	CCCGG
+chr11	8855	8857	0.0909	22	2	B	G	+	6	1	-	16	1	ACCGG
+chr11	8902	8904	0.333	15	5	B	G	+	4	0	-	11	5	CCCGA
+chr11	8914	8916	0.308	13	4	B	G	+	3	2	-	10	2	TACGG
+chr11	8920	8922	0	12	0	B	G	+	4	0	-	8	0	CTCGC
+chr11	8948	8950	0.0588	17	1	B	G	+	4	0	-	13	1	AACGA
+chr11	8975	8977	0	19	0	B	G	+	4	0	-	15	0	GTCGC
+chr11	9127	9129	0.0968	31	3	B	G	+	13	1	-	18	2	CTCGG
+chr11	9136	9138	0.0606	33	2	B	G	+	15	1	-	18	1	CCCGC
+chr11	9227	9229	0.0526	38	2	B	G	+	17	2	-	21	0	CACGC
+chr11	9270	9272	0.0909	33	3	B	G	+	16	2	-	17	1	ACCGC
+chr11	9301	9303	0.0303	33	1	B	G	+	15	1	-	18	0	CTCGC
+chr11	9371	9373	0	32	0	B	G	+	21	0	-	11	0	TTCGA
+chr11	9380	9382	0	31	0	B	G	+	22	0	-	9	0	AACGT
+chr11	9447	9449	0	30	0	B	G	+	25	0	-	5	0	CACGT
+chr11	9478	9480	0.0476	21	1	B	G	+	15	0	-	6	1	ATCGC
+chr11	9510	9512	0.0435	23	1	B	G	+	14	1	-	9	0	CTCGA
+chr11	9588	9590	0.111	18	2	B	G	+	11	1	-	7	1	GGCGC
+chr11	9677	9679	0.231	13	3	+	G	+	13	3	-	0	0	ATCGA
+chr11	9697	9699	0.2	15	3	B	G	+	14	3	-	1	0	TCCGA
+chr11	9709	9711	0.111	18	2	B	G	+	16	2	-	2	0	GCCGG
+chr11	9729	9731	0	16	0	B	G	+	14	0	-	2	0	GACGA
+chr11	9799	9801	0	7	0	B	G	+	4	0	-	3	0	TTCGG
Binary file test-data/SE_g2_r1.bam has changed
Binary file test-data/SE_g2_r2.bam has changed
Binary file test-data/TKO2PE1.fastq.gz has changed
Binary file test-data/TKO2PE2.fastq.gz has changed
Binary file test-data/WTPE1.fastq.gz has changed
Binary file test-data/WTPE2.fastq.gz has changed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/all_fasta.loc	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,1 @@
+chr11	chr11	chr11	${__HERE__}/chr11.fa
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/chr11.fa	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,601 @@
+>chr11
+CCTGGGGTGTTTCAGGTCCCATCTCTCCTCCCAGACCAGAGCGGCTCCCG
+GCATCCCGCCGTCCTGCTCACAGTATCTTCACCTCACCACTTCCACATAG
+TCGCCACCCACAGAGAGTGAGAAAGTGGGACCCAGAGGGGCGAACGGGGG
+AGGGGCAAGGATAGTGGGGTGGTTTTAGTTGGGAGACCACCACTTCCAAT
+GCCCCATCCCCGACTCGCCCTGCATCTCCCTGCTCACCTTTTCATGTGCA
+CCCTACTCATAGTCACATAATTTGTCAATTACAGCAAAAAGCAACAGGAT
+ACAATTATGTTAACCTCGCAGTCTGTCAGTGCGTTATCTTGTCACACTTC
+TACAAGGGGGCTCTCCCTGCATTCTAATGTCAAATTTGAATAGTGCACCA
+CGCAGCTCGTTCTTATGCACACAGTGCTGCCCACAAACACACACACGCAC
+CCACCAGCCACCTTCATACCGCTCCTGACATACACAACCCTCACATTCCC
+AGGCCAACAAAATGGTTCAACAAACACCACCTACATACAATGTGGTCGAT
+GTGTCCCAGGCCTTCAAATGCAGTCTCACCCATGACATTCCCCAAGCATG
+GAAGCCCTGAGAGGAAATGGTTGGGAGAGTAGGGGACAGGCAAAGGGATG
+CACATATGGAGAGCTGCGTGGATGGCTGCTTGGGTTTTACCTTGCGTAGG
+TTGCCCTGGCACCGAAGTCCCCGGATACCAACCAGGGCGGGTGCCCCAGC
+TTCCGGTCTGCCCGTTCAACATCCTTAGTTTATCATACATGCCGTCTGCG
+CCCATCTGTTGCTTTTCGCTAGCCAGGTTGCGAAGAACTCTGTTTATTGA
+TGACACCTGCAAATCAAACCAAAATCAAACCAAATGGTAGTGTCTCCTGA
+AGACACAGTCACCCATCTCAGCTCACCCACAACCTTAAAAAGCAACTCTC
+AACCCGTTAAAAAGCTCCCAGCCACCCCGGGACAGTGGGTGGATTTGCAG
+ATACACCGTGGAAAAAATGATAGCTATCACTTTGGGCATGGAAATTTAGT
+TGGATAGAGAATTGGGCGGAGGGGGTCCTTCACCATTGCTTGGTGTCTCT
+GGTGACCTAACTCACACGGTTGCTTTGCTTTGTTAAAAGACTGCCCGAAA
+ATAAAAACTATTTAGACTTGTTGTAAAGCTTTCTTTCATTATATTGTTTA
+ATATTAGGATGGAAATGGATGTGTGACACCTGAAGGAAAGTTCTGGCAGG
+TGGATTAGGACGAAGTTACAATCTGGTTCTCATTTTCCTTCTTCATTCAA
+GTGCCTTTCCCTGGTGCCTCCGCCCTCTGTTTTGCAGCATGCAAATGAAG
+TGAATGACTCCCCACACTTGACTCCCATTTTCCAGACAGTCAAAGAAAAG
+TCAGGGCATTCCTCTCTGTTCCCCCAGGTACAAAGGAGACAAATGTGGAG
+CAGGGAGAGGACACAGACTAAGAGACAGTGGAGAGAGAGGGTGGGAGGAG
+GTAAAGAGGAGAGAGCATTGGGCTTAGGGCAGGGAGGGCAGATGTTCTCA
+ATGAACTTACGCTTGGTATGTTATCGTTGGTACAGACCCCCTCGGACAGT
+AATCTGTCTCGGATTTCCCAAGCAAAGATGGACGGGCACTCCCGCTTATA
+CTGGGCTATTTTGCTTACAACTTCTGGAGTCGCTACTCTCGGTTTACTAC
+CACCGATTGCCCTGGGTCTGATGGAGCCAGTCTCGTAATACCTGCCCAGA
+ATTTTACTCACACATCCGTTGGACACCTGCATAGGGGAAGTGGACAGAAA
+ACCACATTATTAATAATTTCAAGACAAAAATAAAATTGTTTAAGTATGCA
+TTAAACAATGACAAGCTTACGTTTTGATTGTCCAGCACTTGGACTTTTGC
+ATCTGCATGGGTCTATAACACAAAAATATACCTTCAATGGTATGAGAACT
+TACTGTAGAGAGCTTTTTTTCTTAAAATTACATTTGTAGCCCTAAAAACT
+ACAAATATGATGATACTTTCAAACAATTTGAACTAAAAAACGAATTTAAA
+AGCTTTTTTTAAAAAAAAAACTTTTCTTAAACACTGCCTGAAGATGCATC
+AAAACGAAGTAGATAAATTTATTTTTGTGCTGACCTTGCTTAAAGTGGCG
+TTATGTTTTAAAGAACAAAGGTAAAAAATAAACTGATTTTTTTACTTCTT
+CTTCTTTAAAAAATACTCAATTGAGTAAATGTGAGTTAGGGAGGAGTAAA
+TAGCATTTGCCTTTTAAAATCAATATATATTATTCATAAACATGAAGAAT
+AAACAAATACCTCAAAAGCTGAGGTAAGCATGGGCTGGGGAGAGCAGAAT
+GAAATTATGCCGGTTTATATAACTCAGCTTATTACGGTTCATAAACTGTT
+CCCACAATGAATGTCTTTACCAAATAAAAAGAGAAAGATTTTTTTAAAAA
+AATCCGCACACAATATTTTTCTGAAGCAGATGAGTTATCTTATGTAACTG
+ACCCAGGTTGAAAGAGATAGGGAAGGATGGTGGAAGGAGAGGGGAAAGTG
+GGAAGGCAGGGGAGTGGGTGGGGGGACTGGGGACTGGGGTGGGTGAGGGT
+GGGGGGGTCCATAATTAGCATCGTTTACAGTAAGAAATGAAGAGAGGGCG
+TTGAGAGTGGAGGGCCGCGGGGGCGGCGAGTGGGGCGGCGCCGGGAGGAT
+CACCTGCAGAATTCGGGAAATGTCGCACGGCCGGGCCCCGCTGTGAGCTA
+GCTCTACAATCTTCTGCCGGGTGGAGTCCGGCAGTGGCCGCCCGTTGACA
+AAGACACCACCGAGCTGATTCACTCCGCTGTGACCTGAGGAAAGGGAGAG
+GAGAGGAAAGGGGAAAAGAAGAGAAGAGAGCAGAGTGAAGAGGAAGAAGA
+GGGAGAAGAAGGAAGAGGAAGATGAGGGAGAACAAGAACAGAAAGGAGAG
+GAGGAAGGGGAAGAGGAAGAAGGAGAGGAGGAGGAGACAACCACAATGCA
+TCCTCATCCCCCTGCCAACCTGTGCCAACCTCGGCGGTCAGTTCCCTGCT
+ATCGATCAAATAAAGGACAACAGGACCACCAGCCATGCCAGACATCGGGC
+AGCCCAGTCCAAGCAACCGTGCCCCTGCCCACTGTACTGCAGCAGCAACA
+CAGAGTCAGACACATGCCAAGGAACAAGTACACATGGACCTCACTACAGC
+CATAACCCCAGCAGTGTCCCCAGCCCACTGCCCCTCTTCATAAACACTGC
+AGCCAGCCTCTGAACCCACAAAGCGCTAAGAACCGGCTCTATAAATCTCC
+CAAAGATCCCCTCTGAGTGTCACCGCAAACCATGGGTTCTACAGATCCCA
+CACTGAACCTTCAAGGTCTCTGCAGTGTTGGGACAGTGATAGCAAACCTT
+ACAGCCCTTGAGTGCTCCCTTAGACGCTAAAATCCAGAAATGTTTCCATC
+ACCGTCCAAACTGGCCAAGGCCTGGCACAAGCATCGAGTTGGGTTTTGTT
+TGGTTAGAATAGAAGAAGGAAGAGGAGAGAACAACCAGAGCAGCCAGTAT
+GTGGCTGGTGAACCCCAGTCCAGGTTTCCGCACCTCCAGTGCAGACGTGA
+GAGTCAGAGCCCGGGCTCGGCTTCAGCCACAGACAAGTCCGGCCACCGCA
+AGTAAGGCTTAGAGCCGGGAGTCGGTTCGCACCAGCCACTGCTCGCTTTA
+TTTCCGCCGCTGAACTCCGCCCTCCCCTTAGCCAAGGATAGATATTATCA
+AAAAGAAAAAACAGAAAAGAAAAATACCACATCATAGTGTTCCCCATCCC
+TCCTTTATCTTTGTATCTCACAACCTCATCCCTGTGGTCCAACCTTAGAG
+AATAGGGAGAGGGAGTTAGCCCCTCTTTGTAGTCCCCCTCCCCAAAAACC
+ACAGACACACAAAGCCAATTGCCAATTTCCACTTCCTTAAAAATCTCCAA
+CCTGCAGCCCCGACTTAGATAGCAGCCCTCCCAGGAAAAGACGTGGCTCC
+CTGTTCAGCTCTCCTCTCCTCTCCTCTCTTCTCCTTGGAGCCTACTGATA
+AATTGACGCCAGGAGCCTGAGCTTCTTAGCAATAAATAAGCTCCAAATAT
+AGAAGAGGGGGGAAAATATGATGAGCCTCTCTGACATTTGTCTTTAAAAT
+AAAACTAGCTGCACGTCAAGTTTGAGCTTAATTTCTGGAAATAGGCGGAA
+GTCGCTCCGGATCATGCATGGAAGGCTAATTGAAAAGATCAGTCGGGGCG
+CTTGTGGCAGCCTTGTCACTTTGTGACAGTGCTCCGCTCCGGGAAATTGC
+ATCGTCACGACAAACGGGACCGTGATAAAACGACCCTTTCCGTCCTTATT
+TGTAGATCACTCAGACGAGATTGAACTGCACTTGTTTCCCCTTCGAGGGG
+AGCCGCGTTTTCAGGGTAGCCGAAGGCTTGGGGCTGAGGGGGGGCCCTCA
+CCAAGGCGCGGGTGGGGGCCGGAGCCTCAACTCGATGAGAAGTGACAGGC
+GTTTGGGGGATCTGGGCTCCGGCCGGGACCAGCGCAAGCAGGGACTTTGC
+GGGGACACCGCTTCTCCAACAGAGCAAGGCCTGGCCCACGTTTCCGGTTT
+CTCCTAACTTCCTTTTATTGCCTTCCTTTGCTTCGCAAGTTCCATCTACC
+CCTCCAGCTACAGAGCCCCACCTCTAGGCACAGGAAGCTTCCCGGAAAAA
+GAAAGGCTGTCCCAGAAAGAGACCGAGAGAGACTTTCCAAACTTCGGGCA
+TAGCCACGGCAATTCCCAGTCTGCTAATGCCAAGGCGGGCGCGTAAGGCC
+GCCTAAATCTAGACCTCCCTCCTCACTCATTTCAAAAAATAACAACGTGC
+CAGCCACCTCCGCAGATACCGCCGGCTGGTGCTTGCCCAGGAGACGCCAG
+GGCCAGAGCGCCACTCCCAGCATCGAAATGGCAGAGAGAAAGCGCAGCTC
+CAAATTCCCCTTCAGAGGTTAAGCCTCAATCATTGTGTCCCTTCCCTAGG
+GACTGCTGGCGCTCTCGCCCACTGGCGATGATTATGCGCCTAGAACTCGA
+CCGCGAAGCAACTAATAGGAAAACATATGGTGTCAATTTGGATGCTCCGC
+GCCTCGCGCACACCCGGGAACGAGCGGCACAAAGCCCTGCCGGCCGGCCC
+GCGACCCCGCGCCCCTCGGGGCCTGCCAGCCGGGCCGCAGCGACAAACGC
+TCAGGGCTGCGCGCCCTGGCTGGGGCCCGCCCGAGAGACAGCCTGCGGCT
+GGGGAGTCTGAGCTCCAAGGGGAGAGCCCAGCCGCCGAAGGCGAGCCTAC
+CGGCCAAGCCCTGGGGTCCGGCAGGTTCTGCACAACTACTCCCGCAAAGC
+TCGCCACCTTTGTGCCCTTTCCTCAGCTACGCGCTTACCAGCCCCGGAAG
+CACCAGGGGGCGACCGGCGTTTCCAAAACAGATGGATAAGGTGCTGGGGA
+CCCTCCTTAAGCCTGTGTGGTTCGGCCAGGGCAGGCTGCTGGAAGACCCG
+GGCCTTTCCACACTTTGAGTTCAGTTGAGCCCCGACGCCGTGTGCCTCAA
+CCGTAACCTCAAACCGGGCCTGGTTTCGCTGGGCCTCCAGGCTCCCAGGT
+CGGAGCTCTAGACTAGCCCTGCCGCGTACCAGCACCCACCACGAACAGGA
+AGGTGGGTGGGACCCCACGGGGTCTGGGCAGAGGAAAGCAGGGAGGTACA
+AGGGCCATAAGTATCAAAAGCGCTCTGGGCCTGGGACTATCAGGACCCAA
+GCAATTGCAGTTGCCAGGGCAGGGATTTTCCAATTCCTGGGCACCTCCTT
+ACCAGACCCAGTAAACCGTGCTCTCTCCCGCGCCCTGTAACAGCGAGTTT
+ACCACCCCTACCCTGCTGCAGGCCTGGGGGCTCCTCTCCCCATCCCTTCT
+AGGCCCGTAGTGGGGTGAGGAGGGGAGGAGCTGCCCTGAGAGTGACTTGG
+GGAGAGAAGCCCCCAACCCCGCCGGGTTGGTGGAGACTTTCAGTCTTGGC
+TAAGCCCAAAGGGGCCTCCTTCTCGCTCAACTAGGAGAAGACTTCGCCCC
+GCCCTGATCAGCCCAGCCCCTCTCTCCGGGGTCAGAGCCCGGGCAGGGGC
+GAGAGGGGGTGTGAGTTACAGGATTTGGGGGGGATGGGGTTTCTCTACCG
+CAGCTTCGAAAACTCGGGCGGGCTGTTCTTAAGGGCCCAGCCCCTGCTTC
+TCCAGTATCGAGAAGAGCCAAGCAAACGCCCTCCCCTCCCGGGCTGCCGG
+GCGCGGAGCGGGGAGCAGCCGCCGCAGGTCCCCGGGCCTCAGTCGGTCGG
+CGGCCGGGCCAGCGCGGGCGTCGCGAGTCCCTGTGTCCTCCCCTGACCCT
+CGGGTCCGCGCACCCCGAGCCCGAAGTCCCAGAAAGACCAGAGGCACTTA
+CTGTTCTGCATGCTGGCTCTGGCTGGGGGCCGCGGGATTCCACGGGGCTC
+GAATATGGGGCTCTGTTAGCAAGAAAATAGGAGTTAATCCTCGGGCAGCT
+GCATCAGGTAGGCCTGAACTCCCAACCGAGGTGGACCTGGGGCTAGGACA
+GGAGAATCTCCGTTCTTTAGGAATTGATCCACTCAACCTCTTGGGGCGAT
+CTGAGGGTCACTCAGCTCTGCCCTGCAGTGCATACAAAATCTACCCTTTG
+GGGCTTGAAACAAAGTACTGGAATATCCCCAATCTGTTTCCCCTACATTC
+TTTAGCAGTTAACCCTTTAATCAGTAGAAGCGATGCCCCAATTTTGATTT
+ATTTTCTTCTATCTGAACTATAGTGAAGAAGCATACTGAAATAGTTTGTT
+TTGACTGTTTGTGGGTTTTGAGCCCAAAGCAGCCACCACAGAACTTGCCT
+GAAATCTCGGATGTCTGTCCACTCTCACAATAAAAGGCCTCACACATCTG
+CGCGCCCCTAGTTAAAGTCTTCCCCCTAAGACAAAACAGAAGGAGAAAAA
+GGTATATCAAGCTGTGGTTCAACAGCATACGTCCATGGGCTCTCTGAGCC
+TCACTCTCAGGAAAAACACATGCATCTAGGGCTCAGTCAAACTGTCATCA
+AACTCTAAGTCCCAGCATTACACCCACATTTTCAGCGGCATCTGCCACAG
+AAGTCTCCATCCTTTAAAAGAGGCTGCTCTGGTCTGGCTTGGGAGAATTG
+GAAAATCTAGCAAGTTCGGATCCCACTTTGCAAAAATATAAAGAACCAGG
+AAAAGGGGAAGGGGGAAACAGGAGGAACAGTGGAGTAAGGAGTAGAGGCC
+AAGAGCTGAGGTTGAAGGGGTGTGAGGGGAGAGGAATACCTGAAGTTTGG
+AGACCAGGGGAAGGGAACCTTGGAACTTCCTGGTGCCGGAGAGCATAGGC
+TGGAGCAGTTTGGGAGGGCTGCAGGGCCACACTCCCTCCTTCCTTAAGGA
+ACTTCACAAGACCCCCAGAACCCAACACTCTGCTTACCCTTTATCTTCTA
+CTCCTCTCCCTACACCAAAAAAACTTCCCAATGGCCCCCTCTCCATCTAC
+CTCCCAGCTTTCTGGGACTTTGAAAGGGGCAAGCCCAACTTTTTATCCCA
+AACCTCACTTATTGTGATTAACTTCTCTTTTGGCGTTTAGATCTTCACCC
+CTTTAACTGTACTGTGACACACCCTGTCTACCTCGCCTCACTGCTTTAAA
+AAGTACGAACAGGAGGGGGTGAGCGTGGTCGACAAAGCCGACAGACGGCA
+AACAATGGCTTCTATTCTCCAGCTTTCCTCGAACCCAAGGAGCAAAACCG
+AGTCTCCCTAACCAGCAACACAGCTTCCAACCAAATGTACCAACAACCCA
+CCAACTCACAAAGGGAACAGAGGGACGCTTCTAGCTTCATAAGCCCTGAC
+CCAGTGCAGGGGTCCCCGAGGCTGCTCCCGTACGAACACCTCCCTACCTT
+TCCCAAAGGAAAAAAAAAAAAAGTCTTTTGGGTGAGTGGAGAAAATAGAA
+AAGTAATTTCTTGTTGTAAATATCTGGAGGTAATTGTGTTACTGGTGTGA
+GTTGTTAGGGGCTGTAAAAAGATACAACGGCACCGTGGGGAGCTCAATTG
+GGAAGAGAGAGTTACAATAAAGAGCCCATTCAAGAGGTCTGCCCCCCAAA
+AAGAGGTAAAATAACAAAACAGTCACAGCATCATTTCTCTCTGTTTTTTT
+TTAAATGCATTATTTTTTCTGTTTTACCTATGGTGATTATTGATCTGCCT
+GCCTGCAAGCAAAGAATTTCAAGAAGCACATTTTACTCTACAAGGCCATT
+CCTGGTCTACTTTTTAACTACTTGGTAGACAGCTTAAATCCATAAAAATA
+AACCGAAAATGAGGACTTTTTCTCAGCATTGGGTAAATACTGGCCACAGC
+AGAGGGTTACTAGTTACCCTTATCTAGGTTGGCTCTCCAATGCTTCCTTC
+CTGTCTCTCAGTCCCGGCACCCATAGCCATGCCCTGGATAGGCAAAATCT
+CCAACTCTAATGCAAAATACCCCACCCCAACAAATGCAGGACCCTGAAGC
+TGCAAAGGAGTAGTAGGTGGGGAGATGGGGAATGGAATTCCCCTAGGAGG
+GGTAAAACTTCCAGAACTCCACTTTGTAACCCACGCTGTGCACAATCACA
+ACAGAATGCTCCCATCCTGGCTTGTAACTTCTTGATACCATTTTTTTTTA
+AATTCAGTAGCTGAGCACAAAGAGGCCTCGTTGATTAGTTTTCGAAGCAG
+CTGGGAAGGTGGTAGTATTTCAAATTCCAGTTCTGCTCCATTCCCACTTC
+CCTCTCGCTAAGAAGCTATCTAATTACTTAAGTAAACGTGAACACTAAGT
+GCAGAAGCAACCCGGGAAGTGGGGCACATCGACACCCGGAAGCAGTCCTG
+TGCACCGGAGAGGTCAGACCCCCAACTGCAGCCTCAGCCCCCCAGCTCAG
+CCCGAGGCTCTGTACGGCCTCGCTGGGGCCATAGAACCACATCCTGAACG
+AAGGAGTCTCCAGACATTTTGGAGTCGCCCCAGCCTTGGACCACAAGAGA
+GTGGTCACTGCCCCCACCTCTCTGTCCCCAACACAAGCCTCAGTTTTCAA
+ATGCTGTTCTTTTTTTAATTACTTAATTCTCAAATTCCATCCCCAAACTC
+TTTAACTGAGGTAAAGGGGAGAGCTCTCGGTCCTCCCGCCCCCAGCCAAA
+AGAGTAATTAGCAATTATGATGACAAAGTTTTACTCTCTAGGGCCAAATA
+AGATAACAACTAATCAAACTGCCTTCACGCTTTTGCAAAGGTCACTGGCT
+ATAAAAAATGTAATGTCAACCGCAACAACTATTATTTCTCTTCAGTAAAC
+TCGCCAGGCTCAAGCTCTTTGCCTTTCATTTCATTTTCTGCTCCCAGCCC
+CTTCTAAAACCTTCCCAAGTTCGAGTTAAACGTCCAGAATGGCTTAGGGA
+GCTTGGGAGAGTATAAGAAAAATAAAAAGGAAGAAAGAAAGAAACCACGT
+TCACTTTGTTTTTAAACTTTCTAGAAATCGCTCATGAATGAAGAGCCAGG
+TTTTAGCACTCGAAAAAGAAAGCATTTCCACTTGAAATTCAACTCCATTT
+CCCCCAATTACACTTGCTAAAATTCAGGGCTATTATGGCGCAGGGAGACT
+GAGAAAATCTGCACAAGGAAAGAGGCAACAACTTTATGAACTTTGAAAAG
+ACACTTTTCCTCCTGGGATCTTAGAATCGAAAGCTCTCAAATTACTCCGA
+CTAATTAGCCGGGAGAAGAAGTGATCAGACGAGGGGGAAGGCCAATGCTT
+TCATACAGATTTATTTTCAAGGCCTTACCAAGAAAAACATAAGGGGATTC
+GGAGTAAAACAAACAAAAAAGCAAACCAACTCTGCACCTCAAAACTTGCT
+GAAATCCCAACACCTGGTGAGACTGGATCAGAGGAAAACAAGGAAAAATA
+CCCACCAAGCGCAGGGGGCCGGAGAACAGGGAAGAGCCCCAGCCTCCGCA
+AGGCAAGGGAGATGAGAGTCGGGGGAAGCCTGGGGAGCAGAGAGTTTTCA
+AGTTCCGAGCTCAGCGCTCCAGTGGAGTTGGGGGCAAGAGCCCGGGGGCT
+GGGTGGAAGCTGACAGTCCTCCCGCTCGGCTCTGAGCTCCCTCGGAGCGC
+GCGTTCCACCGCCTTGGCCTGGTGGCTTCCTTCCCGGCTCGCGCTAGCTC
+GCTCTTCCCTCTACCTCCAGCGTGCGGCTGCCTTGCCGACTGACTCCACT
+CGGGCTGGTGGCGGGCGCGGTGCCCCCCTCGGCGCTCGAGTGCCGGCTGA
+CTCTTCTGCGCCGGAGCCCAGTGCGGCACAAAGCCCGAGGGCCGGCGGGC
+GGGCGGAGACGCGGCGAGCAGGTTCCCTCGGCGCTCCCTCGGCGAGTCCT
+CGGAGCGCCCGCTCTCTGCGCGTTCGCTGGGCCTCGGGCTGGACAAAGCA
+GCCGGCAGCGGTTGCTCTCTGCGGCCGTTCCAGGGCCTTGGGCACCACAG
+AGCCTCTTCGTCTTTTCCCGCGGGTCGGGAGCATTGGGCGACGCTGGGGC
+TAGGAGCCGCGGCTGCTGTCGCGCCTGGAGTCTCCCTTCTTCCCGCGAGC
+TCAGCCCGCGCTGGCGCTGGCGCTGGCGCTGGGGCTGAGCTGGCCGCCCG
+CCCCGAGCCCTGCGCCGCTCCCGCTCGCCGGCAGTCGCCCTCCGACGCGG
+GCTTTCGCTGGAAGTAGAAAGTTTGGGCTCCTGCGTGGAAACTTCTCCGG
+CCCAACTCTCCCGGCGTAGCAGTGGGGGAGGGGACCGGCGAGAGGGGAGG
+AAGGAAGGGGGGAGGAAGGGGGAGACCTGTCTGAATATTGCAATAAAAAT
+AAAGCGAGAAGAAAGAAGCGGACTCACCTTTATGAGGCATCCTTTCTGGT
+TGTCACAGCTTCTGTCAAGAGTTTTGTTTGGTTGGGGTTTTTTGTGCTGC
+TGTTGTTGCTTGAAGACCACAATGGTTTGAAATGACGGTGTTTTAAAGGA
+GTTGCTGGTGAGAGTTTTCTCCACGGATGTTGCTGGGTTGGTGTGTGAGA
+GCAATTCTCAGATTCCTGGGAAGGAGACAGAGATTGACAATAAAATGGGC
+TGTCAGCGGCTGGAGAGTGAGAGATAAAGAGTGTGGGTGAGGGAAGTGGC
+TGCAGCCAGCACACCTATGCTGATTGGTGATGGCTCAAGTGTGTTAATGT
+GTGTGTGCCGGCGCCCGGCCTCGCCTCCACCGCTCCTCACTGGCCCATTA
+GCGAAGCCTGACCTCTGTCATCATCCTCCAGCAAAACACTTCCTCCTGCG
+CCTGAACCAGAGCGGGAAATGAGGCCGAGCCACGGTTCCCTTTTCAAACC
+CACTAATCACTCCGCAACATGCAAATGCACCGCTCGCTCTCACACAAAAT
+ATTGCTTTATGCAAAGCAGCGCCGGGGCCTCGCGCCAGCCGATTGGATGC
+TCCCTCTCCCCGCCTGGTGCAATTGGTCCGTTTGTTTGAATATGAATACA
+CTTGTTTTCGGCTCCGCCGGTGGGTCCCGGCTCCTCTCCCCGGCGCAGGG
+CCCCCGAGCCCGGGCTCTGGCGCTCGCCTGCATCTCCCCGCTCCGCGCAC
+GCCGAGATTCGGCGCGGCCCGCCCTGCTGCCGCGAACTTGAGCGGTCAAG
+TGAAGGTTTCTGGGTTGGGGGTCGCGCCCCCACTCCCCCTAAGCTGGACT
+CGGGACTCCCAGTCCTCGGAGGGTGCAAGTTTCTGTGCTGTCCTTTCTTG
+CCAACTCCAGGACAGCGTGGCCTGCCGCCCCGCCCCCGCCCGCCCCCACC
+CCATCCCAGCCGGGGCCCCAAGACACAGAGCACAGCGGCGGCCGCACCCC
+AGCCCACCTGGCCCCTTGGAAGAAGAGGAAAGGGGCGGGAGCGTCACTGG
+GTCCTGGCAGCTCCATTCCAGCTGCGGAGGCCGCGGCGAAGCCTAGGCCG
+AGAGGAGGTTGCCGGGCCGAGAACTAAACGAGGCCAGAGGCCTCCCCAGC
+CCAAGCCTCCAGGGCCTGTCCAAATCCTCCACCCCGCTCTCCCCGCCCTT
+TCTTCTCTCCTTTCTCCTTTCGGAAGCCGCGGTGCGCAGCGGAGCAGAGG
+CACAGCTCTGGCTGGAGAGGCCCGAGTAAACACGCCACTTACTTTGTGTA
+CAGAGGGTTCTTGTGAAAAGCCCTGAAGAGTCCTCACCAAACACTCACCA
+ACTTCTCCCACGTGCCATTTGTTGACTAAGTGCGGCTGCTTTTCGCAGCG
+CCCGAGAGGAGGGGAAGCCAGGGGAGATAAGAGGGGAGGGGGAGTGGAAG
+CCTGGGTGGGGGGGCGGATGACCAATGCTGGGAGGGATTCTCTCTCTCTC
+TCTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTGTCCCACTGGC
+GATGTCGGTATGTCAGCGTGTGTCGTGTGTGCCACCGGTCATCTGTGTGT
+GGTGCGTGTGTCTCCGCACCACGACTCCCTCTCCTGGGTCTGCTCAGTCC
+ACGGAGGCAGCTCCCCTTCAGCTGTTGCCAGAAAACCGGGGAGGAAAGTA
+AGTAAATGGCTTTTCTCTTCTGACCACACACCAGAAGTCCATTTGTTGAA
+TGCCGCACGATTAGGACACACCTCAGTACACTCCAAAGGCGCCCCTCGGA
+GCCTGCAGTTAAGGTGACAGGAGAGCGAGCGGCTGCCTGGTCCTGCGCTG
+TCCCGAGCTCTCTCGGCGGCTTCGGGCCACGGAGCCGGGAAGGAGAAGGA
+AGGCCTGGGCTGCTGTACACACCCCAGGACAGGGCGAGGGGGAATAGGGC
+TGAAGAGCTGAGGCAGGGAGGCTGAAATGAAGCAGCCGAAGGGTACGCTA
+GACTCTTCTTGACTTAGGGAAAGTGCTGAAAGTCAGAGCTCATCAATAAA
+ATCTGCGCGACCTGAACTCGCTGCTGACTGGCCGCGTCGCCGCAGCCCCT
+GGCAGCCGGGACAGTGAGGACGCTGCAAGCCGCCCTGGGCCGCCGGGCCA
+CATGTGGGCCTCACCCTCTCGTTCCTGCCTGGGACCCAGCAGCCCTGCGA
+CCAAGGAGCCCGGGCAGGGCTCCCGGGAAGCAAAGGCCTCCCCGGCCCTA
+GCCAGAGGTGGGGTCGCATAACCCTCTCCCCTCCCCGGGCTGGGGAGGGC
+TGCCGGCTCTAACCCTTCCCCGGGAGCTCGCAGCGCCGGGAGGAAACCTG
+GGTCAGAGGTGAAGGAGGTGGCGCCGGGTGTCCAGGTCCGCGGCCTCCTG
+GGCGAGGTTTCATTTCGCAGGAGCGCGACTTTCGGAGGATGCCGGCAAGT
+GGGGCCGCGGCCCTGAGTCCCAGACCAGCGAAAGTCGGCTTTCTCTTTCT
+CTCTCCTTTTACTTTCCAAGCGAGCATTAAAGAAGCTTGGACTCAAGACT
+CAAGTCCGCTGTGCATGAGCGAGGGACTTGCGGGGGGGCGGGGGGGGGGG
+AAGTGATTCTGTCCGAACTGGAGGCGAGTATCATTTTTACAACCTCTCAA
+TTAGCGACCGAGTGGATTAGCACCTTCTTAAGCTTTGCCTCGGAAAGAAG
+CAGCCGTAATCTTGTTGCTATTTCCCCCCAGGGCACTGGCTCTTTGAATC
+AGCTCCTTCTCCCACCCCTCCAACCCCCATCTCGCTGCGGAGAAAAGTTT
+TAACAAAAAATCAGAAAAGGCAAGGAGCGAGGAGTGAATGCCACTGACGT
+CATTGGGTGGGGAAGGGGGCTCCGGGAAAGACTCCTGGAAACTCCTCCGC
+ACAAACAAAAAACAAAAAGGACTTGGCTGCTTCCCATGTGTCAATATTTG
+GGGAGGGGAGTTCAACCCCAGAATGAGGCCGAGTTTTCAAAGTTTAAGGA
+GAGAGGGGTGAGAGTGAGATTTCCTCGTTCTCTCTGGACTGGTCTCAAGG
+AAGGTGGGGGTGCTGTTCTACCTCGTGGGGAGCCTCACTCCCTCTGCCAA
+GTACCGGCTCTGGGGGCGAAGCGGCCTAGGGGTGTGGGGGCCGAGACTGA
+GGTGCGGGCGCGGTTGCCAGAGGTGTGTACAAAACCGAGACAACTCGCAG
+GTCGCCATTTTAATTGCGGCTACCAACACGCCCGAGGCCCACTTGCTGCC
+CTCCAAGACGCGGGGGGGCGGCGGGGATTGCGCTTTAGGTTCTGTCTCCT
+TTCCCCGCCTCCCCACAGGCCCCTGTTCCCCCCTACCCTCCGCCCCCTGT
+GGCCAGAGGGAGGACCTCTCAGGCCGGCGGGTGCGTTCCCGCGGGGGCCG
+CCGGGGTCGTCCCACCCGGTAGCGGGGCTGTCTGCGTGGAGGGAGTCCCG
+CTGGGGGACCGTGGCTCGCCAGAGCCTGCCCAGAGCTTCTAGCCAGTATT
+GATGCGTCAACATGAGGTGAGGCAAGCATCTCACCCGCACCGATCCCTGG
+CATTCCCCTGTGCCCACCCCATCAGGCCTTCGGGGAGGCAAAGAAGTGCC
+CACTTCTGGGGGGAAGTGCCTCCAGCTCCCAGCTTCAGAACCAAGAGCCT
+CCCATCCGGAGTTCTGACAGTTCTGACACCTCCAAACTCTCTCTTTCCCC
+GGTCAAGCAGAGGCGCAGCGCCGGCGCGCTGGAAAGGGTGAAGACTGGCA
+TTTTTTGGAGAAGGAAAGTCTTGGCAGACACCTGCGTCACTTTGTAGTGC
+CCCAAACGTCTCTCTTTTTGCAGACCCTGGGAAGGTGGGATCCCTTGTCC
+GCGGAGCTGGCAAACGGGTCCAGTGAGCCCTAAAGGCAAAGGCTCCCATC
+CCGGCCGCCCAGTCCCTGGTGCCTCCCCCACTTCCAGGCCTTTGTCCCTC
+TCCTTTCTCCAGCAACCTGCCGGTTCCCGCTAACCCGGCAGCCTCGCCGA
+GCCCAGAGCCTAGAAAAGCCAAGAAGATCCCAGAGCCGCCCCGTGGCACC
+TTCAGGCTGTGGCCACTGCTCGGCTTGGGGGGAAGCCGCGGTGCCCTCGC
+CCTTCCTCGCCTCCTGCCAAATGGCCAGTGCTAGGGGCTCCCTGGGTTCA
+GTCCTCAGTCCTCTTCCGTTCGCTCTGGGAGGCCAGCCTGTCTCTCTCTC
+TCTCTCTCTCTCTCTCTCTCTCTCTCTCTTTCTCTCTCTGTCTTCTCCCT
+GCCCACCCCCAAAATACCTTTAAAATAACATTGTCCAGTGGAGTCTCCCA
+GCAAAGGGATTCCTAGCAGCCCAGCGGGTCAGACTCCCAAGGCTTTGCCT
+CCCAGTCCCTCCATGCTGTTGTCCAGATAGGAGTTAGGTTTAAAACTCCT
+GACATGTAGGACAACTGAAATCAGCGGCACTCTCCGATTTTTGCAAATAG
+GGTTCCTAGAAAAGTTTGCTTTCCATTTCATCACCCTGTGACGATCCCCC
+CATCTCTCTTTTCCCTGTAACCTCTGCCCCCAAGAGGCAAGAGTTCAGCT
+TCAAAGGAAAAGCAAGCAAGTCAAGGCCACCTCAGTCCCTTCTCCCTGGA
+TATCCCCTCTCTAAACCGCTGCCTATTCTAAAGGCCAACGACCCGGTCCT
+GCGATTTGTCCCGTTGTAGACCTGGGAACAGGCAGGCGGGAACTGGGGGC
+TTTACTGGGGGATTTGAGGCTGGGGAGGGGGAGGGAGCAAATGTCATGGC
+TGGCTCGCTCAAGCATCCAGGGAACCGAAGCTAAGCGCATCCTGACGGGC
+TTTTAAAATGACATTGATTAGGACAAGCTGTTCCCAACCCCAGTAAGAGT
+TAATCTGCCTGTTAATCAAGGCACTAAGGGGCTCAATGCGAGTTTTATTA
+GCGACAGGAGAACAGAGGCGGCAGCCAGGGATCAAAAGCCGGGATCCCAT
+ATATTTGTCAGCGCTATCTCCAAAAAAAAAAAAAAAAGTATTTTCAATTT
+GAAAACAATTCGGCGCTTTTCGTCACTTCCTAACCCAGTCTCACAGAGGG
+TGACTTCCAAACCTGGCTAGCGGGGAAAACCGCTGCCCGGGGGACAGAGG
+GGCTGACAGGAACTGCGGGTTGGCTCAGCCGAATGCGGCCGGGGAGAATT
+TAAGAATTCTCAGCCCGCGCGGCCCGATGCCTCTGATTCCTCACGAGAGG
+AAAGGGAATGAAAAATGAAGCAACAAATGACACCACCCAGGCTGGCAGCC
+CTCGTTCCCGGCCAGACCCCGCTCCTCAGGCCCGGCTCTGGCGCCGGGTG
+GCGTCCAGCCCCTGCACGCGCGGCGCGGCCCGCGGGAAAGTTTGTGCAGC
+GAGAGTGACTGTCCTTCCGCCTCGCGCGCGCTGCCCCCTTCTGCCCCGGA
+GGGGCGTTGGGTTCCCTTCGGTTTTCCTTTCCAATTCTAAAATAAATAAA
+TAAACTCCGGGCCCTGGGGAAATCTGAATAGTAATTCCAGTCATCCAAAT
+TCTATGAGAGGGTCTGAGGAAAAAAAGACAAGCCTATCACGGGCGCCCTC
+CGATCACGAAGGGCACGGGCATGAGAAGGCGACAGAGGAGGTCCTCGGCT
+CGCCCTGGGCGCGGAGCGAGGACTCTGGTCAGAGGTAATTATGTCACCGC
+GTTTTCTCGCCGTGACAGCCGAATAAACACGATCTCCAATAAACATCTCT
+AATGAGGGAGGAGGCCCGAGGATGGCTGGGTTTGATTTATGACTGGAGGA
+GAAGGTCCACTTCCCACTGCGAAGCAGGCGACCTGCTCGCCGCCCAGCTC
+CAGGGAGAGGAACCCGCGGAGGAGAGGATCCCGGCCCAGGTAAGGCGTGC
+GGCCGGACTGCCACTGCGCTCGTCCCCCGTTTCCAGGTGAGGTGAGCAGG
+GGGCGCCACCGCTCGGAGTCGGGCGGAGGTCCCAACACCTGTCACAGGTG
+ACTGAGCTCCGAAGGTGCTCCCAGAAGACCTAGAAAGGGGGCGCAGGAGC
+GCGAGTGGGAGTGGGGGTTGGGGGCTAGCCGGGGCCTGGCCCAGAGGCCG
+CCCCAAACTCCCGGGCGGGGTTCTTGGGTGCCCCGCCTAAGCCGGCCTGG
+GCGTCTCCAGCTGCTTTTAGGGCCGGGCTCCCAGGACTTTCGCCCGGCTG
+AACCACAAGTTAGTCTAAGTTTGGGTCGTCCCTTCCCTGACTCACTAAGG
+GGACACTTCCGGGCCCCGGCTCTGGCCCTCCAAGGGCCGCAGGGTGGTGG
+GCGCTTCTCTTTCCCGGGCTTGGCGGCCCGGCTCCCATGGCGCCCGCGCG
+CTTCCCGCCGAGGCAGCCAGCGCAGCGCCCGCAGTAACCGGAGCCCGCGT
+TCTAGCTCAGTGAGAATCGCTAATTATTAGGTCTTCCGAAGTGAAGCTCA
+AATCACACGCTGGGAGCTTTACGCGCTGCCCAGCTCTGTCGGCCGCATTC
+TTTCGTACTCCGCGTCTTTCTGAGAAACGAACCCCATCTCTAGACTTTGT
+TGCGGTCCAGCCATCTGCACTTTGCAACCGCGGGAAACTTCTGCTGCACC
+TCGCCCGACCCGCGGCAGGGTGTCTTTAAGCGCTCCAGGCGCGGGTACGG
+ACTTTCCCAAGCGTGGATGTCGTGAACTGCGATGCAGAGGAGGTCTGGAT
+GCAAAGTGCGTCCCCACCTCCTGCATTTTTCGCGGTCGCGAGTGGCTTTA
+GCAGACAGCAGCAGCAGACTCCCGCCCCCTACTCCCAGCTCAGGGCCGGA
+CTGGAAGTCCACTCTTAGGAGTGGAACATTTTGCTACTTTTGCAGAAGAG
+CACACTCGGCAAAGACAAGTCCCTCCTCCACCCCCCTTCCTATCCCTCCC
+TCAGTAACTCGCTTCCATCTTTGTATGCCTCCTTTTAAAAGCAAAAATCA
+AGCCCCGGGGCCACTTGGGGCGGGCTTTGATTTGGGGACAGCAGTCGCTG
+GCAGCGGAGAGGCCCCGCGGGCCCCGAGCCGGGCCGGGGCCGGAGCCCGG
+GCGGTTACCTGGGAGCGGAGGGCGGGGAGGCAGGTGAGCCTGGCAGCGCC
+GGTTCCCACACTTCTCACCGCCGCTCGGCAGGGGAAGTGGCAGATCTGAC
+AGCCGCGTTCTACGCGAGGACCTGCCCCAGAGTTTAAATGTCAATGATAA
+GAAAAGAGGGTGCTCAGGCAGGCGCTAACTTTCCTTAATATCCACGCCAG
+CGCCGTCCTCATTGGCTGCCCGGCCCGCGTGACGTCATGGCGGCTAGAGT
+TGGGCACAGCTCTGCGCCGACTAGTTTTCCGGCCGGGCGGGAGCCTGCTT
+CTCCCCACCCAGGGTCCGGCCGGCTCCAACCCCTGCTTTGGCCTTCCTTG
+GCCCGGCCTGTGACTGCTCTAGCCCGGCCCGGTCCAGCCTGGGCCGCCGT
+CCTGTATGGCGTAGAGTCGATCCCCAGGAAACTAGTCTTGCCGAGTGCGC
+TCCTCCGCACCTTCGGGCACTAGGGCCCCAAGATCTCTGCCCTTAGGCTG
+GCCAGCTACCTCTCTGGGTAATCACCCTGTCCCCTCCTCCGATCCCCAAG
+GGAGTCCCTTCATCGGAGAGACCCGAAGGCCGCGCCGGGATGGAGAGGCG
+CTTTTGGCAGGAATCATTTTTAGGAGGAGTTCGAGGCTCGATGAGTACCC
+TCAACCTGGCCCTCGAGCCCGCAAGTCACTAGCTGGGCACAGACGGGGGC
+TCCGGGAAGAGGACATGCGACCCCTTTTTCAGGGCCTAGGCCAAAGGTCT
+CTCCAGGCCTGCTGGGTTCTCAAGAGGGAAGCCCTAAAGGGGGCCTACCC
+GGGCTTCCCTCGCCCTCTGGCTGAAAGGCTCCAACCACATTAGAAAATGT
+AGGTGATTGTATTTATCTTTGACTTTGAATTAACTGTTGACTTTGTGACT
+TAGAGGAGGAAAACAGGAAAATCCCAGATTCTGAAGTCAATACTCATAGA
+AATTGGGGCGGGGTCTGCAGGACACTGACTTGCTGGTCTGTGCCTGTTTT
+TCCCTCACACCTTCTGTGGACTGCGAGATGGCTTAGCGGAGCTAGGCCTG
+GCTATCTTGGAATGGAGATTTAGACCCTCGTTGTCATTCCCGTCCCTCAT
+CCCTCAAAACAGGAAGCTTCTCAACTCCTTCCCCTATGCTTCCATCTGGC
+TCTCCAATTTGAAAGTGACAGTGGTGGGGCCCGGAGGGGCAGGGGTCACT
+TGAGGAGTTTGCGGGAGAGGAAGGGGCGGGGCGTGGAATTTCCTTTAGGC
+AGCTCGTGGAGACCCAGCCTGGGTGGACCCGGGAGTGTGTGTGTAGAAGC
+AGGGAGGGGGACGCCTGGGCCGGAGGAGGGTCCTAGGGCCAGTCCCTGCG
+CGCCTTGATCCGGGCGGCGGCTAAGGGAGAGCCGAGCACAGATGGCGCCG
+AGGCGGAGGTCTAGTGATTCCCGCTCCCGCCCCTGGCACACCTGGACAGG
+TGCCGGGCCTCGCAACCGGCCAGGAGGCAGCTCCAGGCAGCCCGTCTGCC
+CACTCCTCTTAGCCCGGACAGCTCCTAGCTGGACTTGATTTTTTGTCTTG
+CTCAGCCCGTCGAGGACTGGAGGGCAGAGTAGGGGTGGGAGTGGGGGTAG
+GGGTAGAGGTGGGGGTGGCGGTGAAGGTGGGGAGGAGAATGGGAGCTTTG
+GCTGAGAAGTTTCCGCAGGATGTTCTTACTTATTTACAAACATTCTCCAA
+ACTTTCTTTCTTCTTCCTAACGCCCTGGGCATCCTGTGTTAAGGCCTAGG
+AGCGCCCACCTCCCACTGGCCACTCTAGTCTGGGCCTAGACCTTACCATG
+GTCCTGTCGCCCTGCTCGGCTGGGAAGGGATGCTAAGAGAGGCATTTGTT
+CTGGTTTCCCTGTCGCCGTCGGCCATCTCCAGTGACCTAGTTTCTTGGTA
+AGCCCCTAGCTCCCCTCCAGGACCTGGCCCAGCCTCCAGCCTCAGATGTC
+CTGAGGTCTCCCAGACTCCGCCGTCCGCCTTGGCGTGTGGGCTGCGGATA
+GGCCTCGGACCTGTGTAGGCCTTCCTAGCTGGCCAGTTTCTGTCCGGGTC
+TCGCTCCCCACACTGTACCAACCCTAGGCCGCCGAGAGGGCTGTTATGGG
+CTTGTTTGGAGTCCAGGCCATTAAGTGTGCGCAAATCCCTGCAACCCAGC
+AGGGTGCTCGAGAGCGACTCTGGAAACGCCACGGCCCCCGCCTGCTTATC
+TGCCTGTCAATAAAGGGCGCGTTCGGGCTGCGCGGGAGCGCCGGTGGTGC
+AGGCGAGGCCCGAACTAGGCCCTTCCGAGCTCCGGGCCAACGCACGAGCC
+GCTCTCCATTTGTCCCGGGCAAGGAGGCTGGGAGGGAGGGAGCCGGCTCT
+GGGGTCGCCGAGGGAGAAGGTGCGGGTTTAGCCGCTGCCTCAGCTGGATT
+CTTGGAGAAGCGCTCGACCTCGGAGCACTCTGGAGACCCTACCAGAAATC
+TAGTGCTGGGGGCTCCTCTCCGCGCCTGCCACCAGCCTTGGAGGAAGCTG
+GGCAGGGCCCCGTGTGGTGCCCAGCCTACGAGGTCCCACCTCGCCTCTCC
+CTGGCAGAGGCATCCTCTCTAGGCCACACTGGCTCGCATTGCGGGAAGAA
+GGGCACCGCGGGCGGCCCGGGGGCATCATCGACGACTCGTAGAGCCCCAA
+ATGCCGAGTGCTACCTGAGGCAGGCTGTAGCTGGGGTGAAAAGGACCGGA
+AGCGGGAAAGAAACCTCTTTTCTTTCAGGCCTCCCTTCAGCCCGGGGTAG
+CACAGACTGCAAATTTCCGACACCTAGGGCAGCTCCAGGCTCTGATTCTC
+CACCGCACAGTGAGCGAGCTCGGCAGCCGCCCAGTCTGCGTAGGCTCCCA
+GGCCAGCCCCGGCCCCAGAGCTTCCCCATCAGCCAGCAGAGGCTGGGTGG
+AGAGCCTGGATGGGCGTGAGAGGCCTTGAAGTAGGGGAAGTAGTGGAGAG
+AGAGAAAGATAGAGAAATCGGCTCTTGGCTTGAACTGTTGTTCTTGTCGA
+AGGTATGTTTTCCCCTCCCGGCCTCTGTCGCTTTTCCTCTCTCCCCTTAA
+CTACTTTATGCAAATAAAGCCGCTTCCAGACTGGAGGAGCCCCTACCGCT
+GGATCTGTGCTACCAGGTGGGTGGGAATTGATTTAGTATTTGCCTTTTTA
+TTGAGAATAGGCACTGAGTCTCCCACCAATTATAAAATAAAAATGGATTT
+TCCTCTTAAAATTTAATTTTACTCAATCTGGAAGATAATGCATCCTGTTA
+ATTTCTACAATATTTTTAAAAAAGGATTCCAGGGTGAGCGCTAACTTATA
+AAATGGTTTAATAAAGCTTTTCCATATACAAATTAATAGGAGAGGAAAAA
+GCCCCTAAATAATTGATTTAAAATACTGTGTATTATTTGAGTGTCATTTC
+AAATAGCTGTACAAGCCAGTAATCAATATATCATCTAAGGAATCTAATTT
+ACTGGGCAGTACGTGTGGCTCTTGCTCTTCCTTCCCTCTCCACTTCCTCC
+CCTTTAAACAGAAAACTTGTTTTAACAAGATTTTTATTTTAAATCAAAAT
+TGAAGGAAAGTTTTGTGAATTCGAGCGTCAAAAAGGGGACTATTTAGCAA
+AGATAATATTTGCCATCTCTTCACCCCGGTCTCAAGGATTGCATAATGGC
+ATGATTATGTGCAATTCTGTGTTTTCTCTGCTATTAATTTATTTTCCTGC
+TCTCATTTAAATAACTACTTTACTTTTTTCTCTCTTCTATTTTGCCACAC
+CATAAAAAGTTTTCCCGGACCTCTGTGGAAAATCTCCAATCCTTCCAGTT
+AGCTTTGTCACAAACAAACAAACAAACAAAAAACAAACAAACAAACAAAA
+CAAAACAAAAAACTCTACAACAATCAATTTGAGAAAACTCTGCTTTTTCT
+TTCCTTTTAAACATTTTCATATTAAGTCAGTTTGTGGCACCATGTTTGAT
+TTCTTTCTCTCCCTTGTCCGGGTAGAGATGGAGACATTAGCTGAATTCCC
+TGGGGGTGATTGGATGCAATTCACTGAAAGTAGACAGGCACGTCCTGGAT
+GTTAGAAATTCACTTTTTCTCAACGTGACAAGGCCGGCTTCGATCTACAA
+TTGTATTTGAATGGGCAATGAGCGGAAAGATCTCTCATCATTAAGCTTGC
+TTCTACTTTAGCAAAATAAGGAGAGCCTTCATTTGCACACAAATGCTCAT
+AAAACTCGGATTCCCTGCAACCCTTTAATAAAAAAGGATTGCTGCATAAT
+TTAAATTGGAGGCCAGGGGAAGGCTGAGAGAGGGAGATGAGGGGGTTTCA
+ACTGAGTCCCCTCCCTTCCCCCAAAAGGTCCATGTGATCGAGTTGACTTT
+TTTACTTTAAATCTTAAATAAGATGTGAGTTGTAACAAGAGGATTGGCAC
+TTACTCCGTGACCCTTATAATTGGATTAAATATAGACTGAGACGTACAAA
+AACACACCAAGATTGACACCTGCCAACTGGACTTCCTTATGATTGATTGT
+GATGCTTCGTGCTGGCAACAATAATGAAATAAATTGTACGGATAAAACAG
+CACATTTGTCATTTCCGTGCCATGAGTGAGGAAGTCAACGAGTTAGGCAA
+GCTGTTAGCTGACAGGCGCGAGCGGGAGGCAAGCTAAAGAGATAAATAAG
+GATTCCAGGCATTGACAGGACAAATTATAACTTGTCACAACCTTTTTAAA
+TTCAGAAGTTTTCAATTAACATATATATTGCTGTTACATATAAGATAAGG
+GGGCACTGACGTAAACACGCGTGGGGTGCAGGGACCCCCCCCCCCAGTGC
+CAAGGGTTTGGGCATCCGGGCTTGGGCATCCAGGCCTCTCTGGCTCTGTT
+CTCCCATTCTGGCAAGCCGGAGGGCACCCAGGTTTGCATGGCCAGCCTGT
+AGGCCCGGGGCTCCGTAAGGTGGGTGAGAGTGTGATAGGCTAAGTGCAGC
+TGCTTCTTTTGCTAGGAGCCCCAGTTCACAGGTCCAGCACCCTTTCTTCA
+TGGAAGCAAAGCTCCCAGATCAAGGGAGAGGCTGGGCGAGGCCTGGGCCC
+AGGTCTGGAAGGACAAGTCCTGCTCTGTCTCCTGGAGTCCTATGTAAGGA
+GTAAGAGATAGAGGGCAGAGGCTCCCCCTTTTGCGCTCCAATGTCTTGCC
+CTCATAGGAAAATCTGAGCTGCTGGACACAACATATCAACAATGGCCACA
+TCCTGTGGGCAGGCCTAGGCCACTGAAACCTGGTGTGGCCTGAGGCTGCT
+CCTGAGCTCTCAACTCCTTCTCTGACTCACCGGCCCAGTGCCGCCAAGGA
+ACCGGTCCTGTGTCCTACCTAGGAAGAGCTGCTTTTCCCTCCTCAGACTT
+TGTTTCTGGTACAGTTCCCAGGGCTCTTCCTTGACCCATGGGTAGGAACG
+AGTTTGCCAGTGGGTCCTTGATAACCTGAGGCCATTCACCTAGAGCAAGT
+GAGATTTTCAGAGGAGCAAAGCGCTGGAGTCACCCTGTGAGGCTGCAAGT
+TAGGGGAAGCCATGAGATTTCTGGGGAAGTCTAGGTGTTGTGACTGTGAC
+TGGGGACTGCTGTGGGTGAGGTTTGAAGAGTGTGGAGGGACTGTAAACTC
+TGCAACCATCTTCACCCCAAAGGAGGCTGGCTTTTTGGTTTTAAAGCTAA
+CCCCTATCCTAAAGCAGAATCTAGACTTGCAGGGCCTCTCTGAACTCAGG
+TCCCAGCCTCCCCGCAAAGTCCAGCTGGAAAGGAAAGTTTGGGCGGGAAG
+CGAGCGTCCCCGGGTCTGAGGCTGCTGGGTGGAAAGCTGGCCTAGGAGGG
+CGAGTAGGGAGCAAAATAGGGCCTGGTCCTTGCGGGGCCGGGGAGCTAGA
+GGCCGCAGGAAGCTGGGAGGGCCCTTTCCAGGCAGGCAGGGCAGCCTCCT
+GGCGCGGCCCAAGGCGCGAGGCAGCCGCTGGCTAGGGAGGTCCAGTCTCC
+TGCAGGCCAGGGAGGGCTTCCTCTGGGTGCAGCCGTCTCCCTCCGCTCCT
+CAGCTTGGCTACTTGGTGCTCGGACCCCGCATCCTAGTCCTCGTCCTTGT
+GGGACAGAGGCCAACCTTGCTGTCGAGGCCCCAATATGTCACCACAAACA
+CTCCCAGAACAAGAGCATCAACCCCAAGCTCCAGAACTTCTGGTCTTAGC
+CTACTGGAGCCTGCCACCAAATTAAACACCGCGAAGGCTTTGAATCCCTG
+CAAAGACGGAGCCTGGAAGAAAGGGAGGGAGGAAAAGGGAGGAAAAAAAA
+CAAGTTACATGGAATTAGTTGGAAGCGGGCGCATTCCCCAGCGACCGCAG
+CCATTAACATTATAATCTCAAACGCCATTTTTTATGTAACACTATTGTTA
+GTAAATCAACTCCTCAGGTCCTCAGAGACGCTGACCCCTGGCAATTCAGT
+GAAAGTGTAATTATCTCCCCGAATCTTGGACAGATTTAAGAGATTACATT
+TCAACTAAATCTATATCAATGCTAATTCCCTCTCTCCCCCCAGCCTGCAG
+GTGCTAGGCACGGGAAATGGCGTCAGGCCCCAGCAGAACGCTTCACAAAG
+AGGAGAGGGCTAGCTGTGTACAGGGCCAAAAGAAATGAAGCCCCACATAA
+AATTACACCTTTTTTTTTTTGCTACAAAACCTGGTTCTTATATAGCAAAG
+TTAAACAAAGCAGCTATGCTAAATTAGCGTTACATAAAGGCGAGTTTAAT
+GCTCGGTTACAGTCAACTGGTGTATCTCTGATGGCTAATTGTCTATCAAA
+ATTAAAGAGAAAAGAACAAATGATTCGTTATAAGAAGGTGTTAAAAGCCT
+CAACAAAAGGACATTTCACAAATATAATATCGGCAGCTCCCGGATTGCCT
+GCGGGGTTGTTCAGGCCCCGGGCCTGCACAGCTATGAGAAGTGGGGTGCA
+GGGTGCTATTTGCCCAGCAAAAGCTCATTTCCTCCCAAATGAGTATAATG
+AAAACCCGAGAGGCGATATTAGTGTAAAGATCGGGATTAGTCCAAATAGG
+AGCCGCGGACTTGGGCTTCAAAGGCGACCCCTACTCTTTTTGCCCAGTCT
+GGGGGCGGAGGGCGCTCGAAGGGATACTTCAAGGACACCACATCCAACTT
+GTCGCAGCCCCACGTTCAGGCACAGGAGGGAGCGAACATTGGGCAGAGGG
+TAATTTCAACAAAGGACGGCAGGCTCCGGGGATCCCGGCTTTTCCCTGAG
+TTCTCTGGGTTGGAAGGAGAGGGTCCGGGCTCTGTCGGAGGAGTGGGTAC
+AGCCCAGCACATGCTTTTCTGACGCGCAGCGCTGGCCTGGGGAGTGGGTG
+CAGCTGCTGCAGCCAGCCCCCACCCCGCCCCAGGCTGGGAAGGACGCAGA
+GTATAATACGACGGCCCCCTTCCCCAAGGGGACACTAGATCCAGAGCTCC
+AGGCATATTTGTTGTTCAAAGTTAATTATACCTGAAAGCAACCAAAATGC
+GCACGCCAGGCCCACGGAGGCCGATTGCCGACCAGCGTAGCCGGCCCGCA
+GCGCCAGAGAAGCCGCTCACTCCGGCGGAGTCTCTGCAAGAGAATCTGAG
+GCTGGGGAATAAAAGCCGAAACTAAATCAGGCAGCAAGAGTAACAAACCA
+TCTGGAAAATGGACCCACGGAGCCGCAAGAATGTACTTTTTTTTTTTATC
+CTTGACAACTTGTGTCCTAAAGAAAGTGTATTGAAGATAATTGAAATGAT
+TACAATTCATTGCTGCCTGACTTTGAATATGAACAGTCACATAATGACAT
+CCCCGCCACACAACAACAGTTTAATTTCCATTATCCGACTTCGGGATCCT
+AGCGCCAGCCGCCGGGGCAACGTGGCGCTAATGCAGAAATCTTGGCTCCT
+GCCTGCATTATTCATTTCGGGGACGTGGGGATGGGATAGGAGTGGAGAGT
+GGGTTTCTTCCCCTCCCCCGAACTTATTCTGAAAATCCCGGGGACCAATG
+GACACCGAGGAACCGAAAGCCGAGAGAATAGCTGTCAAGGGGCCACTTTG
+ACCGACGGGGCTGCGGTCGACGGCGCTCCCCGGCTGTCTTTTTGTCCCAG
+GGGTTCCCAGGCCCGCAGGGCTAGGATCCCCCAGCTCCCACCCTAGTCAG
+TCCCAGTCAGGGAGGACCTAGGCCAGGGCCTCCTCCCAAGGCTGTCATCT
+GGCGGGGCAAACGCTCTTCTCTCATTCCTTCGTCTCCCTTCCAGCCCTCG
+CATCCTTTCCTTCTTTCTGCTTCAGAGCTCTCAGGCTCGCGTGCAGGGCG
+GGATGGGCCCGGATCGGGGGCCCCTTCGCAGATTTGAAAGTTGCCTCGGC
+CTTGATCTTCTCCCAGCCTCCCTGCTGGGGACTCACCCCCACATCCCGCT
+ACCCCCCGGGGAAAAGCGAAGAAAGCCAGAGTCTCTTTGCACGGCGGGGA
+CGCTGGCGCTCACCTCTGGGCGCGGCCTCGGTGGAGCCTCAGGGCCAGGG
+CGGGGTGTGGCGGGGTGTGGCCCGGTGGGGTGAGGGAGGGACCCGGGGTG
+CAAGCCAGGCCAAGGCCTCCTCAGTGCGCAGGCGCCGTTCGGGCCCGCCG
+CCCCCTCCCAAGGCGTGGGGTCCCGGGGACAGACCCCGAGCCCCTCTCTG
+CGCAACTACCGCCTCTAAAAAAAAGACCATCTGGTAATTGTCTTTCTTTA
+GATCTTCCCTTAATTATGCTTTTAATTAAAGATGGTTCAAGGCACCGGAC
+CACATGAAGGGCAGTAATTACAAATTACATTACAAACCCGACAGACGTAG
+CAGGATCTGCTGGGAAAAAGGTACCGGAATAAAAATTTGACTAAAGTTTG
+GGATAAATTAATCGGTGGCTTGAAGAACAGTAGCACATTTCTGATTATGT
+AGGCTCCAGCCAGGATTTAAGTTTTGTGTCATGATTGTTTTCATATTTCT
+TTGACATCTATTGATTTGGAGGCAGCTGCAAGCCTTTGTCACAGGAGAGA
+AGAGCAGTCACAATATATTAGCTTGATGGCATAACTCGAACCGGAGAAGC
+GAGGGGAGCCTGAGTGACTGGGGCGCACATGACCCCCCTCCTGGGGCCTT
+TGGGCTCAACCCTTACTTTTAACTTGTTCTAACAGGAATAATGGAGTAGT
+ACTACTAATGCATGGTCAGAAGGGGCTTTTCTGTCTAGAACGCAGGCGGC
+GGGAGTTGCACCGCCACAGCCCGGGAGACCAAAGAGTCTCCACTCCCCGC
+TCCCGGCGCCTCCTCGCCGACCGGAGCCCCGAGCCCAGTGAGCCCAGCCG
+GTCCCACGTCGAGAGGCGCTGGGTTGGGAGAAGTTCTGCTGGTCGTGGGC
+TCGGCCCCCAGGCGCCAGGCCGAGTGCCCACCTCGGCTTCTTTAGGGGGC
+CTACAGCGGGCAGCCGAGCGGCTGGCGGACCCGGAGCTTGGGAGGCGACC
+GCCAGGCTGGTGCCCGGCCTGACCCGGCGTTCGCGGCCGCCCGCTCGCCC
+GCCGCGGGCCGGGAGCGTACAGGAGTGTGACGCAGATTGTGAAAACAGAA
+GGGAGGGAGTTGGGTCATTTCCTTCGCTAATCATCGCCCCCTCCTGCTCC
+TCCTCCCCACCCCCATCCCATCCTCGCCGCCACCCAGCCTCGAACTCTCT
+TCCTTGCCAGAGTCTGCCTCACTTTGAGCCAGTGGCTGCACTTTATGGAG
+AAAATATATGTGCCCATCACTCACCCATCTATTTGCAGCGTCCATGGAGC
+GCCTCGGGGAACACGCATGGACAATTAAACACACCCACTTACGGCGGGAT
+GATGAGGCTTTCTGTTAACCACATGTGTTCAGAAATGGATACAGCCCTAA
+TTTTGACACGGAATGCAGCTAATGCACCTTGCTATAGATGTACGTCCGTT
+GGAAAAGTGTGTGATGGGCACCTGATCTACACCAAGAGAACAGGGCTGTC
+TGCAGAAGACCCAGCGGCCTGACACCTTTTCTTGTTGCCGGGGAAGGAAA
+GTTTTTCTTTTGATTTGTCCTCTCACGATTTTTCCCTTTTCCTTTGCATT
+ACAGTTTCAGAGAGGAATGTAAATATTTTGACTAGACCTCTGCAGGGGCC
+TTTTAAACACCAGGCACAATCATCTGTGGAATACTGTAGTTCATACTTTC
+AAAACTTGGGTGTTAAAGAGTCTGTCTTGCTGAGGAAACCCATAGAACTA
+GTGGCGACCTTTCAACTCCCTATTTGTGTGATGAAATGAAGTATGCCCTA
+GTGATTGTCGAATTTTAAGGACTTGTTCTGTCACTGCTACTTTCCGATCG
+TTTGCAGGGAAGAAAAATCAATCAAACAGACCTTTGATGGGGTTGTGGCA
+AGGTCACATTTCACTGTGCTCCAGGCTTACCTCTTCTTAACGCCCCTCCC
+TGCTGTGGAGCGTGTGCAGGTAGTATGCAAAACTCTAGAAATAATATTTG
+AAAAAATTAACACCTTGTTTATCTCTTCTGTGTATTTGATATAGGGGTTT
+CTTTCTTTTCCTGTTTTGGTGTTTAATACAAGGCTGATATATTGGAACAG
+AATGCAGACTTCAAAGAAACTAAAATTAGGAGGCAAATAGCGTCAAGGTT
+TCCTCGGTTTCTGTATGTAGTACAAATAGTATAGCGTCTTTTCATCCAAA
+AAGAAGCAAAAATTGCCAGGAGGCATAAAGCATCTCTGAATATCCGTTTT
+TTTTAAAGAAACAAACATTTCCCAAATATTGATGTACGCATGGTTTTATA
+GTCTGAAAATATTTCTAGGCACAAATATGGCATACAAACATTCGTTTCAA
+CTCTGTCGTTGTTGTTAAAATGACAGTAAATTATAGTTTATTTTAAGATA
+CCTCCCTGTTAATGTCAGTTGAGGTTATTTTAGGGCCTGGATGCTCTAGG
+AAGACCACCTAAAATGTGTTTATTTCTCTGCCCAATTCACCTATAGTGTA
+AGGAGCATTTTATTTTATTTTATTTTATTTTATTTTATTTTATTATACTT
+TAAGTTTTAGGGTACATGTGCACAATGTGCAGGTTTGTTACATGTGTATA
+CTTGTGCCATGTTGGTGTGATGCACCCATTAACTCGTCATTTAGCATTAG
+GTATATCTCCTAATGCTATCCCTCCCCTCCCCCCATCCCACAACAGTCCC
+CAGAGTGTCATGTAAGGAGCATTTTAAATAAAATATTTTGGGCATCCTTA
+AGTCCTTTTTTAAAAATCCAACATTCTACAATTGTGACAAATAGGCTGAA
+AAATTTTGTTTTGTTTATACCAAATTCATTTTAATAAAGCATTTGCTGAG
+ATTTTAATACAACAAATGCCAACAAAATGCATCAGAGTAGAAAAAAGAAA
+ATATTTTCAGAATTGCCTATAATATTGCTTGTGATCCTAAATTTTTTGAT
+GCAAGGAACTTGAATTTTTCAGCTAATATTATTCTGTTTTCTTATTCTGT
+GTTTCTATTAGAACACTCACCTATGCTCTGCTTCCCTTGCTGTCTTCCCA
+TAGCTTATGTGTGCACTATGAATTTCAATAAATATAAATTACTGCTTAAA
+AGAGTAAAGAATGGCTACTTTATGACTCACTGTTCTCTCCTTACAATCTG
+CTCAAAAATCTTTTTTACAAATCCCATTGTAGCTACTTCCCCTCAATTCC
+ACGATGAAAGTGAATGGCAAATGCCATTCAATTTAGAGGAACTCAAGCAC
+AGAGAATATAAAGATGTTGTTCAAAGCCTCCCTCAAATCTATCTGTAGGC
+TTACAAACAAAATGAGACTGGAGGCAACATCAGTACACATCTCCGACTCC
+ATATGAGGTTACTTGAAGCTTTGGGGTTCACTCTGGAAACAGAAGACATG
+CATCAAGGTATGTGTTCCTTAGAGACTGACAAGCAACAGTGGACCTCAGG
+AGAAATAAAACAACTCATGCTCTCAACTGCCTCTGATTTGCAGACAAGGA
+AGAGTTTTCTCTAATTTTTTTTTTTTTTTTTTTTTTTTGCCTGCATTACT
+TTTTCTGTGCTTGGTGGCTAAACCAGAGAGCAGTCTGCAACATTTTGCTT
+TTGGGGATTTCCTGTGGCTTTGCCACATTCCCTGGAGAAACTCATAAAGT
+CCTGGAGTCTATGGAGGGATTACAGAGAGCCACAAACATCATTCTTCCTC
+TAATCCATCTGGGGCAGATACAGTATTGATAACAACCTTTCTGTTACTTG
+ACAGTTTCAGACCATTCAGATGAGAACAAATAATAGGAATTTGGTAGTGG
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/seg.fa	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,140 @@
+>seg
+CGAGGCTCTCTCAGGAAAGGAAACGCCACGCAGTGGTGAGGACAGGAAGTGCGCGGTTCT
+GCTGCGGGCCTCTCTGTGGTCAGGACCCGGAAGCTCTCTGGGCCTCAGTTTCTTCAGCTA
+CTCAGTGCCAATAATGCCAGCTCCCACTGCTCCGCGACCAGGGCCTCTGAGCCTCTGTCA
+CGACCTGCACATTAGCTGAAAATTGGAAATTCCACAGCGTGGCCGCGCTGCCGCCTCCCC
+TCGCCCACCTCGACGCGACGTCATCCATTGTCCCTCATTATTCATTATTGTTAGAGGGCT
+CTTAAGTCATCCTTTCTCCTTCTGAACTAGCCGTTGCTCATCTGCTCATGTCTGTGGCTG
+CAGCAATTAATGCCAGTTTCAGTCAGTACAAACCAATCATTCCCTCGGTTTGCTTCTTTG
+GCTGCTATAAACAATTTTTTTTTTTTTTTTTTTTTTGAGACAGGGTCTCGCTTTGTCACC
+CAGGCTGGAGCGCAGTGGTTTGATCTCAGCTCACTGCAGCCTCGCGGTATGGGAGATGTC
+CATTTAGTCACGCCCCCACCGGCGCCTGCCTGCCCCTGTCCCGCCTGGCCCCCACCGGCC
+CCATTCGATGTTCACCGCGCGAAGAGCAGCCGCACCAGCGCCGGCCGTCCCTGGGCCACC
+AGGGCCCGGGCCAGCCCCAACACCTCGTCCGTGTGCCCGCGCCACACGCGCTGCAGCTCC
+TGGCGCAGGCGCGCGCCTCGGACACAAAGCGCCCGGCCTGGGGAATCTTTTCTACAGAAA
+TAAAAGTACTACTTCAAAAGGATGTTTATTTTTATTGTTGCGCTATTTGTCGGGGCTGAG
+AACGGAAAACAATCTAAACAATCTAAACTATTCAGGTGTCAAAAATACAAGTTATTGGTT
+TTTTGTTTTTTTGTTTTTTTTGGTTTTTTTGAGACGCGTTGGTAGGGCTCCATGGTGCTC
+GCTGCGCCAGGATGCAAGCCAGCTGCTGTTTCAGATCCTCGCGATGGACATTTATATCAC
+GTCTGCTTTTTTGGCCGCTGAAACAAGCAGCCGCGGGGTACAGGTGGCTTCGTCCACCGA
+GGTCCCCTCACCCACGCTGAGGCGTCGGAAGCTGCGGACACTGCTCGCGGCCCCTCCTGC
+CGGACCCCTCGCGGCTCACCAGGCTTGCACTCCGGAATTACGGCCTGGTAATTGGTTCCA
+ACGCGGATCATGCTGTCTGTGGAGCCAGGGGGAGGCAGTCAGGACTCCAGGGAGGGCGCG
+CTCCGTTTCGCGCCGTTGTTCCTCCCAGATTTTATAACCCAGGCGTGTAACAGGCATATT
+TACACTTGTGCCACGCACAGGATCGTGCACTTGAAGCAGCAAGGCTGTGGCAGAAATCTG
+CACAGCTGTAGCCATCAGAAACAGGCTCGCGTCGGCTCTGGCCGCCCTCGAGGAGCCCTT
+CAACCCGCCGCTGCGCTGTGGGGGCCCCTCTCTGGGGCTGGCTGAGGCCGAAGCCGGCTC
+CCTCTGCTGGCGGGGAGGTGTGGAGTGGGAGGCGCGAGCGGGAGCCGGAGCTGCGCGCGG
+AGATGTGAGTGCTAGGATGTCAGGGCTTCCCCAGCCCAGCAGTTGCGAGTCTGGAATCGC
+GTGTTCCGCAAATTTACACAATCAAAGCCACCTCACTGGTTTCCGTTCACTTGCCGCGAT
+GGGTACCACCAAGGGTTTGGTGAGGGAATTGGGGGATGGGGGCAGTGCCAGGCCAACAGG
+TCCCCCATCCACTCTTCAGTCTGAGCAATGCCCATAAACTCGACTGACCTTTTTGGACAA
+AAGCTTGCGCGGCGTGCCCGCGCGGCAGCGTGACAGTTGACGGCAGGGGGCGCTCATGCG
+GTCGCGGCGAAAAGTCTCAGGGCACCCTCCGTTTCATTGGCCACAGCTGAGCCACGCGCG
+CGCCGCGGGAACGTGGCCAGCAGGTAATGCGCGTACGCTTGGTGGTCGTGCACGCGCAGG
+ATGGGCCGAGCCCCGCACGCCCCGAGGCAGCACCCCGCACGCAGTAGCCGGGCGTCGAGG
+CGTTCCCAGCCCACACCAGGCGCGGGCTTGGTGCCAGCAAAGGGACCTGCGCGGTGCCCT
+GGGGCCGGGGAGGCCCGCGCAATGTTTTTGAGGCTCATGGTGTGTGGTGGGGGTACCCGC
+GCGCATCCTCCCTGGGCCCAAGATTCCTGGAGCTGGAGAACCCCCTTCCTCCCTGAGGCA
+TGGAGACTGGGCTGGTCCTGAGCCCACTCCAGGCACAGGGGACACACCTGCTTCCTTCTT
+AACTGAGTAGAGCAGAGTTATTTTGAAGATCGGTTTTCAGCAGGACAGTAAAAAGGCGAA
+CTGGGAGATTGTTCATGCAGCAGCTTCTAAAACCAGTTTTTAGATAAATGCACCTTCGTC
+TGTGGGGGGTGGCCGCACCCAGCCTCGCTTCAGGTGCCTGACGGGGTACGTGCTGCGCGT
+TTTCTCATGTGGGAAAGATGCCGGCCCCCGCTGTTTGGGGTGGAAATCCTCTAACCCGCT
+GGGCCTGTGCGCCAGGGCTGAGCAACAGCGCGGAGCCGTGCATCTGTCGGGGGCGGGACG
+CGGCAGCAGGGCTGTGAGGGCTGCGTTTTCTCTGCCAGCCTGATCCGCATCCCCTGCCCC
+TCGGTCGTCTGCCAGTGTCCGTGTATCTTTTGCCCCAGCTGGTGCGGGTTTAGGGATGTC
+AGCCCCAGGAGTCACAGGTGAGCTGGTTCCAATGGGAACCTCGGCATCACCTGGGAACTG
+AAATATGCATTGTCAGCCCCACGGCAGACCTCTGAATCAGAAACTCCGCGCCATGGGCCG
+CTACCGCATCCGCGTGGCCACCGGGGCCTGGCTCTTCTCCGCGCCTGCGCACCTTTGACG
+CGATGCATACACCACATGAAATATCCTATCATCTGTAGGCTCATTCATTTCCCTAACAGC
+AGTAATATTAATAATTTTCATGATTTGAGAAGCCTTCGCTTCAAAGCGAAAAGTCCTAAT
+AATAGAAGAACCCTCCATAAACCTGGAGTGACTATATGGATGCGCCCCACCCTACCACAC
+ATTCGAAGAACCCGCGCCATCGGGTCCTCTGAAGTGCGCGGAAAACCAAGGTCGCGGGAG
+ACGATCCAAAGAGGATGGGCTACGGTCAGACTTCGCGGGTTCAAGTGATTCTCCTGCCTC
+AGCCTCCCCAGTAGCTGGGACTACAGGCGTGCACCACCACATCTGGCTAATTTTTGTACT
+TTTAGTAAAGACATAGAGACGCGCCGGCCCGGGAGCGCCGCGCAGCCCGGGCTCCCGCCC
+GTGCCTCCGCCTCCACACCTCCCCGCGAGCAGAGGGAGCCGCGCCTCCTGTCAGATCAGC
+GGAGGCATTAGATTCTCATGCCAGAGAGAGCCCTATTCTGAACTGCGCATGCGAGGGATC
+TAGGTAGCGCGCTCGTATGAGAATCTAATTCTTGATGATCTGAGGTGGAACAGTTTCATT
+CTGAAACCATCTGCCCTTCCTCCCCCAAACGCCCCCCACCCCCACCCCACCACTGTCCGC
+GCCCGCGCCCGGCAGCGAAGCAGTTAAGCGCGACTCCTCTACTTTGTAGGGGGGCTCCTT
+GCACACTAAAAAGAAAAGGAAGAAACACACCCATCAGCCTCCTGCCTGTTTCAGGCCCAG
+CTGCCAGCGGCAGTACCGCACCCCTCAACCCCACCACCTCCTCACAGGGGCTTGGCTGGC
+CTCTCGGGGGTGGCTTTTCAGGAGGTTAAGCCCCCGACAGGATTGACTTCCACGCGGCCG
+AGTGCACGAGGTCCAGCACCGCGGACGGCCGCGCCGCCTCCAGGCGCTTGAGCAGGTCCC
+GGGCGTGGCCGCGGAAGGCCTCGGGCGCGTAGACGCACTCCTCGCGAAGAGCGCCTCGGT
+CCCAGCGGGGCTCTGGAAACGGTGCACGCGGCGCCCGCCCCTGTCCTCGCGTCACAGTGG
+TCCTCCGCCGGCTACGCGCGGAGACCCGGGAGCGGCCTAGGACCGGAGGCGCGCGGTCAT
+CAACCCGCTCCTGTACACGGTGTCCTCGCAGCAGTTTCGGCGGGTGTTCGTGCAGGTGCT
+GTGCTGCCGCCTGTCGCTGCAGCACGCCAACCACGCGTTTGATCTACAATCGAGAAAGAA
+CTCAAGGGCCAATAACGCGTGTAACACTCACCGCGAGGGTCCACGCGGCACCTGGAAAAT
+CGGGTCACTCCCACCCGAATACTGCGCTTTTCCGACCGGCTTAAAAAACGGCGCGCGAGC
+GACGGTGCCCGGCCGCGCGGCGTCGCCCTTGACAGCTGCCTTGTGGATCTTCTGCAGTTC
+CGAGTACCGGATTCGGTACCCGGAACCCGCGACCCCCAGCTCGGCTCTCCCCACCCTCCG
+GCTCGGCTCTCTCTGACCCCCAGCTCAGCTCTTCCCCACCTGCCAGGACGATTTCCCCAA
+AGCCAGCGCGCAGAGCGCGTGCCGCGGCACAAATACCGCGCCGCCGCAGTGGCCCGCGCC
+ACCAACGGCAGCGAGCCGAGCCATGCTGCTGCAGGCTTTCAAGCCTCCAGCGTGGGACAG
+GTGACTCCTCCTCTAACCTGGGCATGGAGCAGATGGGCTGGCAGACGCAGAAGTGCCTGG
+AATGGCAAGAACCCCTCCCCCCACAGACCGAGGTTATGCAAGTAAGGTTTCTGGACTACA
+TGTTCAGATTGGATGAGAGAAAAACCTCTAGGTCTACTCTGATTGTGCTTTATTTTCACG
+TTCTGACTGGATGACAGCAAGTCTTAGGATAACCAATCAGAACAAGATAACAAAGCCCGA
+TCAGAGTAGGCCTAGAGGTTTCTCTCATCCGATCAGACTATGTAGTCGAGAAACTTCATT
+TGCATTAAACTCAGTATATAAAGGATGCCAAGAGGGGTGGCTTGCCCTTCCAGGCTCATC
+TGCGTCGGCCTGCGTAGCTGCTCTGTGCCCGGCTTAGAGGAGAAGGAGGAGGAGCCACTT
+GCCGCGCATGCCAGAGGCTAGAGTCATCTACCGCGCATGCTGGAGGCTGGAGCCTACAGT
+GGCTCACCTCGCTGCTGCGGTTGGTGGCCGCGAGGGTGACTACAGCGCGGAGTCTGATCT
+GTACTATAAAGAGGTACAAAGTTCTGAGCGCTGTTTAGGAGCAGCTTTCCATTTGTGTAC
+CGCGCCACATATTACATTTTTACGGATTTCTGCGGAACTCCCGCGCCTCTTCACGGGTAG
+GTCAATTTCACTGGTTGAAAGTAAGAGACAGCTGAACCCTCGCGTAGGCATAGGAGATCA
+TCCAAAAGGTGTGTGGATTTCAAACCAACCCAAAGAGAAGGTAAGATTTCTTCCGCTAAG
+TGATGGCCACCAGATGGCTGTGCATCTCAAGGAATGGAGGGAGGACGTCTCTTCTGGACC
+AGAACTGCCAAAGTCCACCATGAGCACTACCTCTCTGCTCTAAATCAAAGAGACAGGGAT
+GAAGATCTGAAAGCGCGCCGGGCTTCATTACGGGCTCGTTGGCCACACCAGCGGCCCCGC
+CTCTCGCCACAGCCTCGCGGTCCCCACAAGTATGAAATGGAGCCCCACAAAAATCTACTG
+AAGCCTGCGTCCATTCTTGTTGCCTCTGACTTTGGTAGTGGAAAGGGTGTCCTGCAAAAG
+CTGGGACCATTTATCACAGAGCTAAACGCGTCAGGTCAGCAGAGCTGAGCCACCCAAGGT
+AGGAGAGAACTGCTCCCCACAGCTTTGACCACATGGGATTCTGAGGGTACTTTGTGGGAG
+AAGCTGGTGGCAGGTTTCAAAGTGGTGGAAGTGGAGAATGCCTCAATCCAGGAAGCCCAG
+GCAAAGCACAGGCTTCCTGGCTGAGTCTGGGGGTTGCCGGGCATTACTCCGCCAGACCCA
+GTGGTACGGGTCTGGGCGCGCCTTGCAGAGCCCCACGGCCTCCTTCGGTGAGGACAAGAA
+GTGGTAGTTGGTCGCGTTTTAGGGAAAATGACAATTAAGCTAAACTGCTAAAGCTGGAGA
+TGATGGATAAGACAAGGATGAGAAAAAGAAGAAGACAGAAGAATTTCTAGGAAATTCAGA
+TTCCTTCTGAGAAAAGTGAGGAGGCAGGGGACAGGACAATTTAAAGATGCCAAAGAATCC
+TTTATTCTTCAAAATTTCAAAATCTTTCAATGCTACTTACTAAAGTAGAACCACTTTAAT
+GGCGTTGGTGTTTACCTCATTTCCTCTTCACTCATTAACAGCAGTATTAGTGGGTCTTTT
+TAGTCTATTTCATCTTCTGTGTCTATCACACCGCGACCTATGCAGCCTCTGCCCGAGGGA
+GCCCCGCAGCCCAGAACACCTAACAAAAGCAACGCAGGCAGCCACCAGTGATTACAGGAG
+GCTCCCCCAAGACCTAGGAGCTGGTAAAGGGCTATTGCTCTCCCGTCTCACCGCAAAGCA
+CGCGCCCGGCCAGGGTGAACGCCGCGCGTGTCTGGAGTTTGTTCCTTCTGATGTTCGGAT
+ATGTTCGCGCCCACCATCACGCCCAGCTAATTATTTTGTATTTTTAGTAGAGACGGGGTT
+TCATCGTGTTAAGCAGGATGGTCTCCATCTCCTGACCTCGTGATCCGCCCGCGAGGGACG
+TGGCTCAATTCATCCACGAATTTTCTTTCGCGGCGGCAGGCGGAGAGGAGTGCTCGGCTG
+GCCTCTGCCTCCGGGAGCCAGGCCCCCTGCACAGCTGCCCTGGAGCTCCCCAGGGCCCAG
+CTCATGGACAGCCGGGAGAGCCCCAGTGGGAGGAGGGAGTAGCAGAGAATAGGCCGTCAA
+GAGAAAGTGCGCAGCGTGGAATTCAAAGATAGTGTGAAGCTGGCATGGGGCAGACAGGAC
+GGAGGCGTTCTGGAACGCAGGCGTGTGCTCCCAACTGCCTAAATGGCTGTGTGGTGAGGT
+GCGGAGGTGCGGCCACACCCCTAATGGCGGTGTGGTAAGGCGCGGAGGTGCGGCCGCAAC
+ATCCTCATGTGTCCTTTGGTTTGCGTGTGCCATGCAGCCCTGCAGACGCAGACAGCACCC
+CAGTGAGCTCCTCCCCAGGGGGAGCGCGTTAACATCGGTCCTGGGGAAACCCACTGTGCC
+TGCTTCCAAATGGAAATCAGGGCGGCGCGCTCTTTTCCCCCAGTCTGGAGTGCAATGGTG
+CAATCTCGGCTCGGTGCAATCTCGGCTCACTGCAAGCTCCGTCTCCCGCGCTGAGGAAAC
+AAATCGTTGGCTCAATAGTTTTCCTTTTGTGTAATATACTGTCCCAGTGTGTTTAGGGAG
+TTACGCGCCGCCGCCACCACCATCGCATAATTCACCGGCATTGATGAAGGCTGCCGCGGA
+GGGGAAGAGGGCAGACCTCAGCGCTGCGGGGGCCGGCGGCGCGCTCGGGAAGCAGCGCGG
+CCCCAGGAGGAAGGCGGCGCGGCCGAGGCCAGAGCCGGCGGCTACTGCGACCTTCCGGCT
+GGCGCGCAGTGGCAGGCACCTGTAATCTCAGCTACTCGGGAGGCTGAGGGAGGAGATTCA
+CCTGAACCTGGGAGGAAGAATTTGCAGTGATCAGAGATGTCGCCATCGTACTCCAGCCTG
+GGCAACAAGAGTGAAACTCCATTTCAAAAAAAGAAAAAGGATGACCAGGAGCAAATTAAT
+CACATGAATTCCATCTTGCTCAAAAGATATTACAATCAGGAAAAACATGAAAAAAAATCT
+GAGCAAAAAGACAAAAAGGATAATTACATGTCTTTGCTATTTAGTTCTACTCCTCTTTTC
+CTTGTATCTTTTTGTATATCCACCTCTGCAGTCTACTGGTACCTCTATGTTCAGATCCAA
+TATCGCGTGAGGTAATAAGTGGAAAAGAGCTGATGAATTTTCAGGTGATTTTCATTTACA
+CTGAGGAAATGTAACCAAGTTGGGAGGAAAAAAATCAATGATCTCTATGTCCTAATGATA
+CTGAAAATTTTAAATTTTCAATATAGTATTATACACACCTAGAACTCAGGGAAAAAAGAA
+GAGGCTTCGTCCAATTATGATATTCTTCTCTCCCTTTCTGTTTCAGTTCTGGCCCTAGAG
+CAAAAATGTCTGCGAGGGCGTTCTGCCCTTTGTGATTTGTAAAAGGCACCGAGGCGCAGC
+CCGGAGGAGACAGCGCGCCTGGCACTGGCGTCACTACAAAAAATCAGTTCCAAGATGACT
+CACGCGCACAAAGAAGGCTAACGCCGAGGGTGCGCCCGGCTGGCGGAGGGCGCGGGTCTC
+CAGAGACGCGGGAGGGCGGGCCGTGCGGCGGGTGAGGCCGCGGCCAAGGCTATAAATTCG
+CAGGCCGCGGCCGGGCCCCACAGGAGCAGCCGCCCGGGGCACCGGAGCTGCG
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_data_table_conf.xml.sample	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,8 @@
+<!-- Use the file tool_data_table_conf.xml.oldlocstyle if you don't want to update your loc files as changed in revision 4550:535d276c92bc-->
+<tables>
+    <!-- Locations of all fasta files under genome directory -->
+    <table name="all_fasta" comment_char="#" allow_duplicate_entries="False">
+        <columns>value, dbkey, name, path</columns>
+        <file path="tool-data/all_fasta.loc" />
+    </table>
+</tables>
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tool_data_table_conf.xml.test	Fri Sep 06 09:54:27 2019 -0400
@@ -0,0 +1,8 @@
+<!-- Use the file tool_data_table_conf.xml.oldlocstyle if you don't want to update your loc files as changed in revision 4550:535d276c92bc-->
+<tables>
+    <!-- Locations of all fasta files under genome directory -->
+    <table name="all_fasta" comment_char="#" allow_duplicate_entries="False">
+        <columns>value, dbkey, name, path</columns>
+        <file path="${__HERE__}/test-data/all_fasta.loc" />
+    </table>
+</tables>