[med-svn] [pbtestdata] 01/02: Imported Upstream version 0.0.13
Afif Elghraoui
afif at moszumanska.debian.org
Sun Jul 24 03:33:15 UTC 2016
This is an automated email from the git hooks/post-receive script.
afif pushed a commit to branch master
in repository pbtestdata.
commit bc5276c324d909898dc774994f5d85342106a58a
Author: Afif Elghraoui <afif at ghraoui.name>
Date: Sat Jul 23 16:21:28 2016 -0700
Imported Upstream version 0.0.13
---
.gitattributes | 5 +
Makefile | 7 +
README.md | 44 ++
...00001823085912221377_s1_X0.aligned_subreads.bam | 3 +
...1823085912221377_s1_X0.aligned_subreads.bam.bai | 3 +
...1823085912221377_s1_X0.aligned_subreads.bam.pbi | 3 +
...550000001823085912221377_s1_X0.alignmentset.xml | 15 +
.../m54006_160504_020705.alignmentset.xml | 51 ++
...4006_160504_020705.tiny_mapped.alignmentset.xml | 15 +
.../m54006_160504_020705.tiny_mapped.subreads.bam | 3 +
...4006_160504_020705.tiny_mapped.subreads.bam.bai | 3 +
...4006_160504_020705.tiny_mapped.subreads.bam.pbi | 3 +
...4006_160504_020705.tiny_mapped.alignmentset.xml | 15 +
.../m54006_160504_020705.tiny_mapped.subreads.bam | 3 +
...4006_160504_020705.tiny_mapped.subreads.bam.bai | 3 +
...4006_160504_020705.tiny_mapped.subreads.bam.pbi | 3 +
data/BarcodeSet/example_barcodes.barcodeset.xml | 15 +
data/BarcodeSet/example_barcodes.fasta | 6 +
data/BarcodeSet/example_barcodes.fasta.fai | 3 +
.../mapped.consensusalignmentset.bam | 3 +
.../mapped.consensusalignmentset.bam.bai | 3 +
.../mapped.consensusalignmentset.bam.pbi | 3 +
.../mapped.consensusalignmentset.xml | 15 +
data/ConsensusReadSet/barcoded.ccs.bam | 3 +
data/ConsensusReadSet/barcoded.ccs.bam.pbi | 3 +
.../ConsensusReadSet/barcoded.consensusreadset.xml | 17 +
data/ConsensusReadSet/ccs_mixed.bam | 3 +
data/ConsensusReadSet/ccs_mixed.bam.bai | 3 +
data/ConsensusReadSet/ccs_mixed.bam.pbi | 3 +
.../ccs_mixed.consensusreadset.xml | 15 +
...100569412550000001823090301191423_s1_p0.ccs.bam | 3 +
...69412550000001823090301191423_s1_p0.ccs.bam.bai | 3 +
...69412550000001823090301191423_s1_p0.ccs.bam.pbi | 3 +
...00001823090301191423_s1_p0.consensusreadset.xml | 15 +
data/ContigSet/simple.contigset.xml | 20 +
data/ContigSet/simple.fasta | 4 +
data/ContigSet/simple.fasta.fai | 2 +
data/GFF/alignment_summary.gff | 491 +++++++++++++
data/GFF/alignment_summary_consensus.gff | 106 +++
data/GFF/variants.gff | 14 +
...00702482550000001823141103261590_s1_p0.1.bax.h5 | 3 +
...00702482550000001823141103261590_s1_p0.2.bax.h5 | 3 +
...00702482550000001823141103261590_s1_p0.3.bax.h5 | 3 +
...0001823141103261590_s1_p0.all.hdfsubreadset.xml | 12 +
...c100702482550000001823141103261590_s1_p0.bas.h5 | 3 +
data/ReferenceSet/lambdaNEB/reference.info.xml | 24 +
data/ReferenceSet/lambdaNEB/referenceset.xml | 22 +
.../lambdaNEB/sequence/lambdaNEB.fasta | 810 +++++++++++++++++++++
.../sequence/lambdaNEB.fasta.contig.index | 2 +
.../lambdaNEB/sequence/lambdaNEB.fasta.fai | 1 +
.../lambdaNEB/sequence/lambdaNEB.fasta.index | 27 +
.../lambdaNEB/sequence/lambdaNEB.fasta.sa | Bin 0 -> 718324 bytes
data/SubreadSet/barcoded.scraps.bam | 3 +
data/SubreadSet/barcoded.scraps.bam.pbi | 3 +
data/SubreadSet/barcoded.subreads.bam | 3 +
data/SubreadSet/barcoded.subreads.bam.bai | 3 +
data/SubreadSet/barcoded.subreads.bam.pbi | 3 +
data/SubreadSet/barcoded.subreadset.xml | 22 +
...4852550000001823085912221377_s1_X0.subreads.bam | 3 +
...550000001823085912221377_s1_X0.subreads.bam.bai | 3 +
...550000001823085912221377_s1_X0.subreads.bam.pbi | 3 +
...52550000001823085912221377_s1_X0.subreadset.xml | 15 +
...01823085912221377_s1_X0_filtered.subreadset.xml | 22 +
data/SubreadSet/m54006_160504_020705.sts.xml | 1 +
.../m54006_160504_020705.tiny.scraps.bam | 3 +
.../m54006_160504_020705.tiny.scraps.bam.pbi | 3 +
.../m54006_160504_020705.tiny.subreads.bam | 3 +
.../m54006_160504_020705.tiny.subreads.bam.pbi | 3 +
.../m54006_160504_020705.tiny.subreadset.xml | 97 +++
data/files.json | 158 ++++
...01823085912221377_s1_X0.aligned_subreads.cmp.h5 | 3 +
src/python/Makefile | 15 +
src/python/pbtestdata/__init__.py | 1 +
src/python/pbtestdata/core.py | 118 +++
src/python/setup.py | 42 ++
version.json | 3 +
76 files changed, 2378 insertions(+)
diff --git a/.gitattributes b/.gitattributes
new file mode 100644
index 0000000..1708c81
--- /dev/null
+++ b/.gitattributes
@@ -0,0 +1,5 @@
+*.bam filter=lfs diff=lfs merge=lfs -text
+*.pbi filter=lfs diff=lfs merge=lfs -text
+*.bai filter=lfs diff=lfs merge=lfs -text
+*.gz filter=lfs diff=lfs merge=lfs -text
+*.h5 filter=lfs diff=lfs merge=lfs -text
diff --git a/Makefile b/Makefile
new file mode 100644
index 0000000..93b2464
--- /dev/null
+++ b/Makefile
@@ -0,0 +1,7 @@
+.PHONY: python
+
+python:
+ cd src/python && make install
+
+clean:
+ cd src/python && make clean
diff --git a/README.md b/README.md
new file mode 100644
index 0000000..ae71a18
--- /dev/null
+++ b/README.md
@@ -0,0 +1,44 @@
+# PacBioTestData
+Small datasets for testing
+
+Usage
+-----
+
+This repository requires the Git Large File Support extension
+(https://git-lfs.github.com/) to fetch binary data files. To download and
+install as a Python module/cmdline tool:
+
+```
+ $ git clone https://github.com/PacificBiosciences/PacBioTestData.git
+ $ cd PacBioTestData
+ $ git lfs pull
+ $ make python
+```
+
+You can then access the `pbtestdata` Python module programatically:
+
+```
+ $ python
+ >>> import pbtestdata
+ >>> pbtestdata.get_file("subreads_bam")
+ '/path/to/movie.subreads.bam'
+```
+
+Or use the command-line tool:
+
+```
+ $ pbdata --help
+ usage: pbdata [-h] [-v] {show,get,validate} ...
+
+ Utility and API for accessing representative PacBio data files for testing.
+ Run 'pbdata show [--verbose]' to display a list of files sorted by type.
+
+ positional arguments:
+ {show,get,validate}
+
+ optional arguments:
+ -h, --help show this help message and exit
+ -v, --version show program's version number and exit
+ $ pbtestdata get subreads-bam
+ /path/to/movie.subreads.bam
+```
diff --git a/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam
new file mode 100644
index 0000000..a177a36
--- /dev/null
+++ b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:235c12691bff133e429804f7407ef462064894a4bdf7dc0324dc932b2e7ed4a7
+size 201960
diff --git a/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.bai b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.bai
new file mode 100644
index 0000000..e4bee5e
--- /dev/null
+++ b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c55e8603ac09b5b97b8f5b0cc31edea78f73bff15d349274f42935603d7c6436
+size 112
diff --git a/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.pbi b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.pbi
new file mode 100644
index 0000000..374e87c
--- /dev/null
+++ b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d3de1e1222033ebe8fd6d62fa50571a11652b770495382a42536c976a75d75f4
+size 2609
diff --git a/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.alignmentset.xml b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.alignmentset.xml
new file mode 100644
index 0000000..58bafd4
--- /dev/null
+++ b/data/AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.alignmentset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<AlignmentSet CreatedAt="2016-05-31T13:28:07" MetaType="PacBio.DataSet.AlignmentSet" Name="rsii-aligned" Tags="" TimeStampedName="pacbio_dataset_alignmentset-160531_202807168" UniqueId="65a1916e-fd46-37aa-c224-5b6aae10f85d" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.AlignmentFile.AlignmentBamFile" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam" TimeStampedName="pacbio_alignmentfile_alignmentbamfile-160531_202807167" UniqueId="b5084853-4668-40a8-bbfc-f76b8dc615e5">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160531_202807167" UniqueId="e6801e5c-0284-47f2-9103-be1f7d186ece"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160531_202807167" UniqueId="24997237-2841-47d5-a4e8-49b8e9ad9ddc"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>60467</TotalLength>
+ <NumRecords>112</NumRecords>
+ </DataSetMetadata>
+</AlignmentSet>
diff --git a/data/AlignmentSet/m54006_160504_020705.alignmentset.xml b/data/AlignmentSet/m54006_160504_020705.alignmentset.xml
new file mode 100644
index 0000000..cb04cfd
--- /dev/null
+++ b/data/AlignmentSet/m54006_160504_020705.alignmentset.xml
@@ -0,0 +1,51 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<AlignmentSet CreatedAt="2016-07-15T08:22:36" MetaType="PacBio.DataSet.AlignmentSet" Name="aligned-ds-2" Tags="" TimeStampedName="pacbio_dataset_alignmentset-160715_152236743" UniqueId="d18bbef7-d6db-a25f-09a8-b1e302a1f85b" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.AlignmentFile.AlignmentBamFile" ResourceId="mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam" TimeStampedName="pacbio_alignmentfile_alignmentbamfile-160715_152155892" UniqueId="7383dfb3-67b9-40c8-8f1a-db29db84ebb9">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160715_152155892" UniqueId="934014e7-2dbd-421d-9c16-0dba123b9f78"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160715_152155892" UniqueId="160059c9-dc51-495f-b77e-237ecd96ef00"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ <pbbase:ExternalResource MetaType="PacBio.AlignmentFile.AlignmentBamFile" ResourceId="mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam" TimeStampedName="pacbio_alignmentfile_alignmentbamfile-160715_152236742" UniqueId="f04ec688-2679-48aa-8872-31d490b81939">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160715_152236743" UniqueId="0fbdb364-0a48-40ec-b9ae-76deea54fb02"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160715_152236743" UniqueId="b9411aad-e869-43f0-9f47-c5c4fc49ef3b"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSets>
+ <AlignmentSet CreatedAt="2016-07-15T08:22:36" MetaType="PacBio.DataSet.AlignmentSet" Name="aligned-split-2" Tags="" TimeStampedName="pacbio_dataset_alignmentset-160715_152236743" UniqueId="c03e7d00-ed7d-3220-75c7-80b6cbd0fe1d" Version="3.0.1" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.AlignmentFile.AlignmentBamFile" ResourceId="mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam" TimeStampedName="pacbio_alignmentfile_alignmentbamfile-160715_152155892" UniqueId="7383dfb3-67b9-40c8-8f1a-db29db84ebb9">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160715_152155892" UniqueId="934014e7-2dbd-421d-9c16-0dba123b9f78"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160715_152155892" UniqueId="160059c9-dc51-495f-b77e-237ecd96ef00"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>113886</TotalLength>
+ <NumRecords>11</NumRecords>
+ </DataSetMetadata>
+ </AlignmentSet>
+ <AlignmentSet CreatedAt="2016-07-15T08:22:36" MetaType="PacBio.DataSet.AlignmentSet" Name="aligned-split-2" Tags="" TimeStampedName="pacbio_dataset_alignmentset-160715_152236743" UniqueId="ef50dbf9-3ef3-bf59-afd8-72278108c95b" Version="3.0.1" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.AlignmentFile.AlignmentBamFile" ResourceId="mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam" TimeStampedName="pacbio_alignmentfile_alignmentbamfile-160715_152236742" UniqueId="f04ec688-2679-48aa-8872-31d490b81939">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160715_152236743" UniqueId="0fbdb364-0a48-40ec-b9ae-76deea54fb02"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160715_152236743" UniqueId="b9411aad-e869-43f0-9f47-c5c4fc49ef3b"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>99864</TotalLength>
+ <NumRecords>10</NumRecords>
+ </DataSetMetadata>
+ </AlignmentSet>
+ </DataSets>
+ <DataSetMetadata>
+ <TotalLength>213750</TotalLength>
+ <NumRecords>21</NumRecords>
+ </DataSetMetadata>
+</AlignmentSet>
diff --git a/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.alignmentset.xml b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.alignmentset.xml
new file mode 100644
index 0000000..1f7a8ca
--- /dev/null
+++ b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.alignmentset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<AlignmentSet CreatedAt="2016-07-15T08:21:55" MetaType="PacBio.DataSet.AlignmentSet" Name="aligned-split-1" Tags="" TimeStampedName="pacbio_dataset_alignmentset-160715_152155893" UniqueId="fce14729-d663-01b5-9c06-c451e2790a11" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.AlignmentFile.AlignmentBamFile" ResourceId="m54006_160504_020705.tiny_mapped.subreads.bam" TimeStampedName="pacbio_alignmentfile_alignmentbamfile-160715_152155892" UniqueId="7383dfb3-67b9-40c8-8f1a-db29db84ebb9">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m54006_160504_020705.tiny_mapped.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160715_152155892" UniqueId="934014e7-2dbd-421d-9c16-0dba123b9f78"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="m54006_160504_020705.tiny_mapped.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160715_152155892" UniqueId="160059c9-dc51-495f-b77e-237ecd96ef00"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>113886</TotalLength>
+ <NumRecords>11</NumRecords>
+ </DataSetMetadata>
+</AlignmentSet>
diff --git a/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam
new file mode 100644
index 0000000..8e15e4c
--- /dev/null
+++ b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:78d193926de13ea596e8890e4662bc08fcf1efe33547eb71a5eaa96f66d29232
+size 242878
diff --git a/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.bai b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.bai
new file mode 100644
index 0000000..f578c08
--- /dev/null
+++ b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:a55ec8b34946d943e17fc84a4686b4d4d9a0e73cf7d97ab3e21beb1118d19774
+size 168
diff --git a/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.pbi b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.pbi
new file mode 100644
index 0000000..ed3d7b2
--- /dev/null
+++ b/data/AlignmentSet/mapped-1/m54006_160504_020705.tiny_mapped.subreads.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8072f4c3fb898f7dcc5262d2de320a734f40193f3d175f7f1ed75f53d56744f6
+size 448
diff --git a/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.alignmentset.xml b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.alignmentset.xml
new file mode 100644
index 0000000..1dbeda5
--- /dev/null
+++ b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.alignmentset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<AlignmentSet CreatedAt="2016-07-15T08:22:36" MetaType="PacBio.DataSet.AlignmentSet" Name="aligned-split-2" Tags="" TimeStampedName="pacbio_dataset_alignmentset-160715_152236743" UniqueId="d18bbef7-d6db-a25f-09a8-b1e302a1f85b" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.AlignmentFile.AlignmentBamFile" ResourceId="m54006_160504_020705.tiny_mapped.subreads.bam" TimeStampedName="pacbio_alignmentfile_alignmentbamfile-160715_152236742" UniqueId="f04ec688-2679-48aa-8872-31d490b81939">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m54006_160504_020705.tiny_mapped.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160715_152236743" UniqueId="0fbdb364-0a48-40ec-b9ae-76deea54fb02"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="m54006_160504_020705.tiny_mapped.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160715_152236743" UniqueId="b9411aad-e869-43f0-9f47-c5c4fc49ef3b"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>99864</TotalLength>
+ <NumRecords>10</NumRecords>
+ </DataSetMetadata>
+</AlignmentSet>
diff --git a/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam
new file mode 100644
index 0000000..590615b
--- /dev/null
+++ b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:31f062c514bbf98c3a920d1d1e265b14b318d44031f73c9daa772e67187f290e
+size 215752
diff --git a/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.bai b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.bai
new file mode 100644
index 0000000..f7c5352
--- /dev/null
+++ b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:56553f958d1614e2bdd9ded9fb864f1642961380575545eb9a03ad4ad231f82a
+size 168
diff --git a/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.pbi b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.pbi
new file mode 100644
index 0000000..525edef
--- /dev/null
+++ b/data/AlignmentSet/mapped-2/m54006_160504_020705.tiny_mapped.subreads.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:6c134d1262c29ddf65c4b7212da542374eb821c2654a44c0bc2a5f6c6cafa38d
+size 425
diff --git a/data/BarcodeSet/example_barcodes.barcodeset.xml b/data/BarcodeSet/example_barcodes.barcodeset.xml
new file mode 100644
index 0000000..f88d4e8
--- /dev/null
+++ b/data/BarcodeSet/example_barcodes.barcodeset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<BarcodeSet CreatedAt="2016-06-17T10:07:57" MetaType="PacBio.DataSet.BarcodeSet" Name="pacbio_example_barcodes" Tags="" TimeStampedName="pacbio_dataset_barcodeset-160617_170757950" UniqueId="afa1b6fd-6a44-50d8-838f-cb9d485ce9da" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.BarcodeFile.BarcodeFastaFile" ResourceId="example_barcodes.fasta" TimeStampedName="pacbio_barcodefile_barcodefastafile-160617_170757950" UniqueId="2826be05-8d49-41b0-a057-66af9980c0ee">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.SamIndex" ResourceId="example_barcodes.fasta.fai" TimeStampedName="pacbio_index_samindex-160617_170757950" UniqueId="e6247b21-9204-4f32-9992-5fa5931bce35"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>48</TotalLength>
+ <NumRecords>3</NumRecords>
+ <BarcodeConstruction/>
+ </DataSetMetadata>
+</BarcodeSet>
diff --git a/data/BarcodeSet/example_barcodes.fasta b/data/BarcodeSet/example_barcodes.fasta
new file mode 100644
index 0000000..2b02060
--- /dev/null
+++ b/data/BarcodeSet/example_barcodes.fasta
@@ -0,0 +1,6 @@
+>lbc1
+TCAGACGATGCGTCAT
+>lbc2
+CTATACATGACTCTGC
+>lbc3
+TACTAGAGTAGCACTC
diff --git a/data/BarcodeSet/example_barcodes.fasta.fai b/data/BarcodeSet/example_barcodes.fasta.fai
new file mode 100644
index 0000000..c47060b
--- /dev/null
+++ b/data/BarcodeSet/example_barcodes.fasta.fai
@@ -0,0 +1,3 @@
+lbc1 16 6 16 17
+lbc2 16 29 16 17
+lbc3 16 52 16 17
diff --git a/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam
new file mode 100644
index 0000000..3cafe04
--- /dev/null
+++ b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8536b592ed37d7fb97d822c6bf6c5654662fd7715607fef94331936594eeff98
+size 9281
diff --git a/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam.bai b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam.bai
new file mode 100644
index 0000000..0a15c2d
--- /dev/null
+++ b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:bc670ff5f0d5bfae0750a9253b916d1ba6b044f8deb7264209e0f5c0541e5595
+size 96
diff --git a/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam.pbi b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam.pbi
new file mode 100644
index 0000000..b532504
--- /dev/null
+++ b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b161b4ea6efbfc769ac577a6cfb468288fa363bbeb88acc8083f36f940849fac
+size 306
diff --git a/data/ConsensusAlignmentSet/mapped.consensusalignmentset.xml b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.xml
new file mode 100644
index 0000000..108baf2
--- /dev/null
+++ b/data/ConsensusAlignmentSet/mapped.consensusalignmentset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ConsensusAlignmentSet CreatedAt="2016-06-17T10:28:24" MetaType="PacBio.DataSet.ConsensusAlignmentSet" Name="rsii-ccs-aligned" Tags="" TimeStampedName="pacbio_dataset_consensusalignmentset-160617_172824317" UniqueId="aac61a94-6046-6719-efc6-773b0f56d792" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbioscie [...]
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.ConsensusReadFile.ConsensusReadBamFile" ResourceId="mapped.consensusalignmentset.bam" TimeStampedName="pacbio_consensusreadfile_consensusreadbamfile-160617_172824316" UniqueId="70c9875b-a799-44e3-890c-8f52e9887572">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="mapped.consensusalignmentset.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160617_172824316" UniqueId="ae544058-5c00-4399-8857-893e3988002e"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="mapped.consensusalignmentset.bam.bai" TimeStampedName="pacbio_index_bamindex-160617_172824316" UniqueId="196feb28-3c69-45f0-a65f-01a558ef8ab9"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>11229</TotalLength>
+ <NumRecords>10</NumRecords>
+ </DataSetMetadata>
+</ConsensusAlignmentSet>
diff --git a/data/ConsensusReadSet/barcoded.ccs.bam b/data/ConsensusReadSet/barcoded.ccs.bam
new file mode 100644
index 0000000..c9e29fc
--- /dev/null
+++ b/data/ConsensusReadSet/barcoded.ccs.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:cea7d6d485df2705478b56fe2a0d8d5be924c5d9816a817d0345223670787e7f
+size 5800
diff --git a/data/ConsensusReadSet/barcoded.ccs.bam.pbi b/data/ConsensusReadSet/barcoded.ccs.bam.pbi
new file mode 100644
index 0000000..f1e2d15
--- /dev/null
+++ b/data/ConsensusReadSet/barcoded.ccs.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4962d99a2bbd5cd966de852dde0dfd543457eafc8403a67541414c309c86678d
+size 110
diff --git a/data/ConsensusReadSet/barcoded.consensusreadset.xml b/data/ConsensusReadSet/barcoded.consensusreadset.xml
new file mode 100644
index 0000000..139aadf
--- /dev/null
+++ b/data/ConsensusReadSet/barcoded.consensusreadset.xml
@@ -0,0 +1,17 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ConsensusReadSet CreatedAt="2016-06-29T05:36:27" MetaType="PacBio.DataSet.ConsensusReadSet" Name="ccs-barcoded" Tags="" TimeStampedName="pacbio_dataset_consensusreadset-160629_123627937" UniqueId="6bdaaeb8-2d8f-aa27-a4b0-82db7d227d9d" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioData [...]
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.ConsensusReadFile.ConsensusReadBamFile" ResourceId="barcoded.ccs.bam" TimeStampedName="pacbio_consensusreadfile_consensusreadbamfile-160629_123636057" UniqueId="10a378d7-634a-4212-bdd8-9154483025dd">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="barcoded.ccs.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160629_123636057" UniqueId="486f486b-7b8b-4dd4-9923-cc5d62e803c1"/>
+ </pbbase:FileIndices>
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.DataSet.BarcodeSet" ResourceId="../BarcodeSet/example_barcodes.barcodeset.xml" TimeStampedName="PacBio.DataSet.BarcodeSet-160629_123636" UniqueId="1b186027-27c0-4dc6-9652-f35b92d85561"/>
+ </pbbase:ExternalResources>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>0</TotalLength>
+ <NumRecords>2</NumRecords>
+ </DataSetMetadata>
+</ConsensusReadSet>
diff --git a/data/ConsensusReadSet/ccs_mixed.bam b/data/ConsensusReadSet/ccs_mixed.bam
new file mode 100644
index 0000000..3c5c148
--- /dev/null
+++ b/data/ConsensusReadSet/ccs_mixed.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:09239ef1e6997213709171fcb425275bcb3bd219d36db3a74f425ba9e857227b
+size 1659
diff --git a/data/ConsensusReadSet/ccs_mixed.bam.bai b/data/ConsensusReadSet/ccs_mixed.bam.bai
new file mode 100644
index 0000000..dfe6a78
--- /dev/null
+++ b/data/ConsensusReadSet/ccs_mixed.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8b27e414ac5ca899c577eb45c26167178cc5cf844944b6212a297b751a387554
+size 16
diff --git a/data/ConsensusReadSet/ccs_mixed.bam.pbi b/data/ConsensusReadSet/ccs_mixed.bam.pbi
new file mode 100644
index 0000000..1d5325b
--- /dev/null
+++ b/data/ConsensusReadSet/ccs_mixed.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:51b0d1cdc06a25cdf13602920b3bd41ded6411aec6977db19f243988be158741
+size 102
diff --git a/data/ConsensusReadSet/ccs_mixed.consensusreadset.xml b/data/ConsensusReadSet/ccs_mixed.consensusreadset.xml
new file mode 100644
index 0000000..4591cc8
--- /dev/null
+++ b/data/ConsensusReadSet/ccs_mixed.consensusreadset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ConsensusReadSet CreatedAt="2016-06-17T10:31:43" MetaType="PacBio.DataSet.ConsensusReadSet" Name="rsii-ccs-multi-cell" Tags="" TimeStampedName="pacbio_dataset_consensusreadset-160617_173143383" UniqueId="11a80fb3-b369-8924-5408-29dff906b7bf" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/Pac [...]
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.ConsensusReadFile.ConsensusReadBamFile" ResourceId="ccs_mixed.bam" TimeStampedName="pacbio_consensusreadfile_consensusreadbamfile-160617_173143382" UniqueId="215e028c-03aa-4b23-a7a1-1f0a1a4044a2">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="ccs_mixed.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160617_173143383" UniqueId="1330fc73-8756-4230-a0e7-fa079dc5fe36"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="ccs_mixed.bam.bai" TimeStampedName="pacbio_index_bamindex-160617_173143383" UniqueId="6d6886fb-2544-4c47-9c48-07bf5f305ccc"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>0</TotalLength>
+ <NumRecords>2</NumRecords>
+ </DataSetMetadata>
+</ConsensusReadSet>
diff --git a/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam
new file mode 100644
index 0000000..2dd39fb
--- /dev/null
+++ b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:52fe41fc9707641735d36fb181efefd1e4ee5e18721f49f03e7ef2b0c97ac55e
+size 278396
diff --git a/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.bai b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.bai
new file mode 100644
index 0000000..35a32a5
--- /dev/null
+++ b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:130eb68ccb54362037a69b02eb25063f743f6ad6d6392058cbf8613bfeecb7fe
+size 16
diff --git a/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.pbi b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.pbi
new file mode 100644
index 0000000..ba991fd
--- /dev/null
+++ b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:19244d22b361f3de649630dc0492e9624d2ee542998ecf96bfff01f9014cbec5
+size 301
diff --git a/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.consensusreadset.xml b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.consensusreadset.xml
new file mode 100644
index 0000000..b37e21c
--- /dev/null
+++ b/data/ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.consensusreadset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ConsensusReadSet CreatedAt="2016-06-17T10:31:23" MetaType="PacBio.DataSet.ConsensusReadSet" Name="rsii-ccs" Tags="" TimeStampedName="pacbio_dataset_consensusreadset-160617_173123347" UniqueId="7a70efd2-18ae-9861-3368-05f35562e3be" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.ConsensusReadFile.ConsensusReadBamFile" ResourceId="m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam" TimeStampedName="pacbio_consensusreadfile_consensusreadbamfile-160617_173123346" UniqueId="3a901d8f-20cc-4075-a3ea-0c74e232d6e4">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160617_173123346" UniqueId="56038397-ee3c-4920-8307-ca30f332d2fd"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam.bai" TimeStampedName="pacbio_index_bamindex-160617_173123346" UniqueId="9fefd001-6e68-4a3e-8ac0-07e50e7b218b"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>0</TotalLength>
+ <NumRecords>18</NumRecords>
+ </DataSetMetadata>
+</ConsensusReadSet>
diff --git a/data/ContigSet/simple.contigset.xml b/data/ContigSet/simple.contigset.xml
new file mode 100644
index 0000000..1e3cbce
--- /dev/null
+++ b/data/ContigSet/simple.contigset.xml
@@ -0,0 +1,20 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ContigSet CreatedAt="2016-06-24T17:52:57" MetaType="PacBio.DataSet.ContigSet" Name="simple_contigs" Tags="" TimeStampedName="pacbio_dataset_contigset-160625_005257667" UniqueId="0ef5ddf7-6a61-45b5-81e8-939e5ae06996" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.ContigFile.ContigFastaFile" ResourceId="simple.fasta" TimeStampedName="pacbio_contigfile_contigfastafile-160625_005257666" UniqueId="49d53821-9aac-4b54-b384-7bfb686dd821">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.SamIndex" ResourceId="simple.fasta.fai" TimeStampedName="pacbio_index_samindex-160625_005257666" UniqueId="92718de9-80ce-48b5-9319-0bb76e36f927"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>56</TotalLength>
+ <NumRecords>2</NumRecords>
+ <Organism/>
+ <Ploidy/>
+ <Contigs>
+ <Contig Description="very minimal example of a multi-contig FASTA file" Digest="DEPRECATED" Length="28" Name="chr1"/>
+ <Contig Description="note that GenomicConsensus may output mixed case" Digest="DEPRECATED" Length="28" Name="chr2"/>
+ </Contigs>
+ </DataSetMetadata>
+</ContigSet>
diff --git a/data/ContigSet/simple.fasta b/data/ContigSet/simple.fasta
new file mode 100644
index 0000000..a79352c
--- /dev/null
+++ b/data/ContigSet/simple.fasta
@@ -0,0 +1,4 @@
+>chr1 very minimal example of a multi-contig FASTA file
+ACGTACGTACGTACGTACGTACGTACGT
+>chr2 note that GenomicConsensus may output mixed case
+AAACCCGGGTTTaaacccGGGTTTAAAT
diff --git a/data/ContigSet/simple.fasta.fai b/data/ContigSet/simple.fasta.fai
new file mode 100644
index 0000000..7b976fc
--- /dev/null
+++ b/data/ContigSet/simple.fasta.fai
@@ -0,0 +1,2 @@
+chr1 28 56 28 29
+chr2 28 140 28 29
diff --git a/data/GFF/alignment_summary.gff b/data/GFF/alignment_summary.gff
new file mode 100644
index 0000000..69e6d9e
--- /dev/null
+++ b/data/GFF/alignment_summary.gff
@@ -0,0 +1,491 @@
+##gff-version 3
+##date 2016-01-11T09:14:15
+##source PACBIO_AlignmentSummary 1.0
+##source-commandline summarize_coverage.py m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam lambdaNEB.fa alignment_summary.gff --num_regions=500
+##sequence-header lambda_NEB3011 lambda_NEB3011
+##sequence-region lambda_NEB3011 1 48502
+lambda_NEB3011 . region 1 100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 101 200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 201 300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 301 400 0.00 + . cov=0,1,1;cov2=0.980,0.140;gaps=1,2
+lambda_NEB3011 . region 401 500 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 501 600 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 601 700 0.00 + . cov=1,1,2;cov2=1.250,0.433;gaps=0,0
+lambda_NEB3011 . region 701 800 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 801 900 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 901 1000 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 1001 1100 0.00 + . cov=0,0,2;cov2=0.450,0.792;gaps=1,74
+lambda_NEB3011 . region 1101 1200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1201 1300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1301 1400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1401 1500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1501 1600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1601 1700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1701 1800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1801 1900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 1901 2000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 2001 2100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 2101 2200 0.00 + . cov=0,0,2;cov2=0.550,0.865;gaps=1,70
+lambda_NEB3011 . region 2201 2300 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 2301 2400 0.00 + . cov=1,1,2;cov2=1.260,0.439;gaps=0,0
+lambda_NEB3011 . region 2401 2500 0.00 + . cov=0,0,1;cov2=0.130,0.336;gaps=1,87
+lambda_NEB3011 . region 2501 2600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 2601 2700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 2701 2800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 2801 2900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 2901 3000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 3001 3100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 3101 3200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 3201 3300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 3301 3400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 3401 3500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 3501 3600 0.00 + . cov=0,0,1;cov2=0.280,0.449;gaps=1,72
+lambda_NEB3011 . region 3601 3700 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 3701 3800 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 3801 3900 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 3901 4000 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 4001 4100 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 4101 4200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 4201 4300 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 4301 4400 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 4401 4500 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 4501 4600 0.00 + . cov=1,3,4;cov2=2.950,0.572;gaps=0,0
+lambda_NEB3011 . region 4601 4700 0.00 + . cov=4,4,5;cov2=4.420,0.494;gaps=0,0
+lambda_NEB3011 . region 4701 4800 0.00 + . cov=5,5,5;cov2=5.000,0.000;gaps=0,0
+lambda_NEB3011 . region 4801 4900 0.00 + . cov=5,6,6;cov2=5.580,0.494;gaps=0,0
+lambda_NEB3011 . region 4901 5000 0.00 + . cov=6,7,7;cov2=6.970,0.171;gaps=0,0
+lambda_NEB3011 . region 5001 5100 0.00 + . cov=5,5,7;cov2=5.260,0.642;gaps=0,0
+lambda_NEB3011 . region 5101 5200 0.00 + . cov=5,6,6;cov2=5.650,0.477;gaps=0,0
+lambda_NEB3011 . region 5201 5300 0.00 + . cov=4,4,6;cov2=4.300,0.714;gaps=0,0
+lambda_NEB3011 . region 5301 5400 0.00 + . cov=2,4,4;cov2=3.770,0.630;gaps=0,0
+lambda_NEB3011 . region 5401 5500 0.00 + . cov=2,3,3;cov2=2.740,0.439;gaps=0,0
+lambda_NEB3011 . region 5501 5600 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 5601 5700 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 5701 5800 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 5801 5900 0.00 + . cov=2,3,3;cov2=2.600,0.490;gaps=0,0
+lambda_NEB3011 . region 5901 6000 0.00 + . cov=2,3,3;cov2=2.580,0.494;gaps=0,0
+lambda_NEB3011 . region 6001 6100 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 6101 6200 0.00 + . cov=1,1,3;cov2=1.540,0.805;gaps=0,0
+lambda_NEB3011 . region 6201 6300 0.00 + . cov=1,1,2;cov2=1.430,0.495;gaps=0,0
+lambda_NEB3011 . region 6301 6400 0.00 + . cov=2,3,3;cov2=2.890,0.313;gaps=0,0
+lambda_NEB3011 . region 6401 6500 0.00 + . cov=3,3,5;cov2=3.280,0.694;gaps=0,0
+lambda_NEB3011 . region 6501 6600 0.00 + . cov=3,5,5;cov2=4.140,0.990;gaps=0,0
+lambda_NEB3011 . region 6601 6700 0.00 + . cov=3,4,4;cov2=3.690,0.462;gaps=0,0
+lambda_NEB3011 . region 6701 6800 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 6801 6900 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 6901 7000 0.00 + . cov=3,4,4;cov2=3.650,0.477;gaps=0,0
+lambda_NEB3011 . region 7001 7100 0.00 + . cov=3,3,5;cov2=3.690,0.902;gaps=0,0
+lambda_NEB3011 . region 7101 7200 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 7201 7300 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 7301 7400 0.00 + . cov=1,2,3;cov2=2.010,0.985;gaps=0,0
+lambda_NEB3011 . region 7401 7500 0.00 + . cov=0,1,1;cov2=0.960,0.196;gaps=1,4
+lambda_NEB3011 . region 7501 7600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 7601 7700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 7701 7800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 7801 7900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 7901 8000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 8001 8100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 8101 8200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 8201 8300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 8301 8400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 8401 8500 0.00 + . cov=0,0,2;cov2=0.920,0.987;gaps=1,53
+lambda_NEB3011 . region 8501 8600 0.00 + . cov=2,4,4;cov2=3.560,0.791;gaps=0,0
+lambda_NEB3011 . region 8601 8700 0.00 + . cov=4,5,5;cov2=4.920,0.271;gaps=0,0
+lambda_NEB3011 . region 8701 8800 0.00 + . cov=5,6,6;cov2=5.690,0.462;gaps=0,0
+lambda_NEB3011 . region 8801 8900 0.00 + . cov=6,6,6;cov2=6.000,0.000;gaps=0,0
+lambda_NEB3011 . region 8901 9000 0.00 + . cov=5,6,7;cov2=5.910,0.680;gaps=0,0
+lambda_NEB3011 . region 9001 9100 0.00 + . cov=7,7,7;cov2=7.000,0.000;gaps=0,0
+lambda_NEB3011 . region 9101 9200 0.00 + . cov=5,7,7;cov2=6.220,0.965;gaps=0,0
+lambda_NEB3011 . region 9201 9300 0.00 + . cov=5,6,7;cov2=5.930,0.495;gaps=0,0
+lambda_NEB3011 . region 9301 9400 0.00 + . cov=6,6,7;cov2=6.160,0.367;gaps=0,0
+lambda_NEB3011 . region 9401 9500 0.00 + . cov=6,6,6;cov2=6.000,0.000;gaps=0,0
+lambda_NEB3011 . region 9501 9600 0.00 + . cov=6,7,7;cov2=6.930,0.255;gaps=0,0
+lambda_NEB3011 . region 9601 9700 0.00 + . cov=6,6,7;cov2=6.310,0.462;gaps=0,0
+lambda_NEB3011 . region 9701 9800 0.00 + . cov=6,6,6;cov2=6.000,0.000;gaps=0,0
+lambda_NEB3011 . region 9801 9900 0.00 + . cov=5,6,6;cov2=5.930,0.255;gaps=0,0
+lambda_NEB3011 . region 9901 10000 0.00 + . cov=0,0,4;cov2=0.230,0.773;gaps=1,91
+lambda_NEB3011 . region 10001 10100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 10101 10200 0.00 + . cov=0,0,1;cov2=0.120,0.325;gaps=1,88
+lambda_NEB3011 . region 10201 10300 0.00 + . cov=1,1,3;cov2=1.350,0.753;gaps=0,0
+lambda_NEB3011 . region 10301 10400 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 10401 10500 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 10501 10600 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 10601 10700 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 10701 10800 0.00 + . cov=3,3,3;cov2=3.000,0.000;gaps=0,0
+lambda_NEB3011 . region 10801 10900 0.00 + . cov=1,1,3;cov2=1.770,0.882;gaps=0,0
+lambda_NEB3011 . region 10901 11000 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 11001 11100 0.00 + . cov=1,2,3;cov2=2.290,0.621;gaps=0,0
+lambda_NEB3011 . region 11101 11200 0.00 + . cov=3,4,4;cov2=3.850,0.357;gaps=0,0
+lambda_NEB3011 . region 11201 11300 0.00 + . cov=2,4,4;cov2=3.720,0.694;gaps=0,0
+lambda_NEB3011 . region 11301 11400 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 11401 11500 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 11501 11600 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 11601 11700 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 11701 11800 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 11801 11900 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 11901 12000 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 12001 12100 0.00 + . cov=0,2,2;cov2=1.560,0.766;gaps=1,17
+lambda_NEB3011 . region 12101 12200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12201 12300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12301 12400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12401 12500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12501 12600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12601 12700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12701 12800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12801 12900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 12901 13000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 13001 13100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 13101 13200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 13201 13300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 13301 13400 0.00 + . cov=0,0,6;cov2=1.460,2.484;gaps=1,74
+lambda_NEB3011 . region 13401 13500 0.00 + . cov=6,8,8;cov2=7.600,0.648;gaps=0,0
+lambda_NEB3011 . region 13501 13600 0.00 + . cov=7,7,7;cov2=7.000,0.000;gaps=0,0
+lambda_NEB3011 . region 13601 13700 0.00 + . cov=7,7,9;cov2=7.630,0.924;gaps=0,0
+lambda_NEB3011 . region 13701 13800 0.00 + . cov=9,10,10;cov2=9.860,0.347;gaps=0,0
+lambda_NEB3011 . region 13801 13900 0.00 + . cov=10,10,10;cov2=10.000,0.000;gaps=0,0
+lambda_NEB3011 . region 13901 14000 0.00 + . cov=10,10,10;cov2=10.000,0.000;gaps=0,0
+lambda_NEB3011 . region 14001 14100 0.00 + . cov=9,9,10;cov2=9.390,0.488;gaps=0,0
+lambda_NEB3011 . region 14101 14200 0.00 + . cov=1,1,9;cov2=2.960,3.209;gaps=0,0
+lambda_NEB3011 . region 14201 14300 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 14301 14400 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 14401 14500 0.00 + . cov=0,1,1;cov2=0.940,0.237;gaps=1,6
+lambda_NEB3011 . region 14501 14600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 14601 14700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 14701 14800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 14801 14900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 14901 15000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15001 15100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15101 15200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15201 15300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15301 15400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15401 15500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15501 15600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15601 15700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15701 15800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15801 15900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 15901 16000 0.00 + . cov=0,0,1;cov2=0.400,0.490;gaps=1,60
+lambda_NEB3011 . region 16001 16100 0.00 + . cov=1,1,2;cov2=1.350,0.477;gaps=0,0
+lambda_NEB3011 . region 16101 16200 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 16201 16300 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 16301 16400 0.00 + . cov=2,2,3;cov2=2.310,0.462;gaps=0,0
+lambda_NEB3011 . region 16401 16500 0.00 + . cov=0,2,2;cov2=1.670,0.649;gaps=1,10
+lambda_NEB3011 . region 16501 16600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 16601 16700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 16701 16800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 16801 16900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 16901 17000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17001 17100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17101 17200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17201 17300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17301 17400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17401 17500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17501 17600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17601 17700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17701 17800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17801 17900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 17901 18000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 18001 18100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 18101 18200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 18201 18300 0.00 + . cov=0,1,1;cov2=0.700,0.458;gaps=1,30
+lambda_NEB3011 . region 18301 18400 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 18401 18500 0.00 + . cov=0,1,1;cov2=0.930,0.255;gaps=1,7
+lambda_NEB3011 . region 18501 18600 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 18601 18700 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 18701 18800 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 18801 18900 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 18901 19000 0.00 + . cov=0,1,1;cov2=0.890,0.313;gaps=1,11
+lambda_NEB3011 . region 19001 19100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19101 19200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19201 19300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19301 19400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19401 19500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19501 19600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19601 19700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19701 19800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 19801 19900 0.00 + . cov=0,0,1;cov2=0.190,0.392;gaps=1,81
+lambda_NEB3011 . region 19901 20000 0.00 + . cov=1,1,2;cov2=1.020,0.140;gaps=0,0
+lambda_NEB3011 . region 20001 20100 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 20101 20200 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 20201 20300 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 20301 20400 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 20401 20500 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 20501 20600 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 20601 20700 0.00 + . cov=0,2,2;cov2=1.180,0.984;gaps=1,41
+lambda_NEB3011 . region 20701 20800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 20801 20900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 20901 21000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21001 21100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21101 21200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21201 21300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21301 21400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21401 21500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21501 21600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21601 21700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 21701 21800 0.00 + . cov=0,2,2;cov2=1.460,0.877;gaps=1,26
+lambda_NEB3011 . region 21801 21900 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 21901 22000 0.00 + . cov=1,1,2;cov2=1.440,0.496;gaps=0,0
+lambda_NEB3011 . region 22001 22100 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 22101 22200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 22201 22300 0.00 + . cov=1,2,2;cov2=1.650,0.477;gaps=0,0
+lambda_NEB3011 . region 22301 22400 0.00 + . cov=0,2,2;cov2=1.720,0.694;gaps=1,14
+lambda_NEB3011 . region 22401 22500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 22501 22600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 22601 22700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 22701 22800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 22801 22900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 22901 23000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23001 23100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23101 23200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23201 23300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23301 23400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23401 23500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23501 23600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23601 23700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23701 23800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23801 23900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 23901 24000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24001 24100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24101 24200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24201 24300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24301 24400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24401 24500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24501 24600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24601 24700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24701 24800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24801 24900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 24901 25000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 25001 25100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 25101 25200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 25201 25300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 25301 25400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 25401 25500 0.00 + . cov=0,1,1;cov2=0.990,0.099;gaps=1,1
+lambda_NEB3011 . region 25501 25600 0.00 + . cov=1,2,2;cov2=1.650,0.477;gaps=0,0
+lambda_NEB3011 . region 25601 25700 0.00 + . cov=1,1,2;cov2=1.270,0.444;gaps=0,0
+lambda_NEB3011 . region 25701 25800 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 25801 25900 0.00 + . cov=0,0,1;cov2=0.200,0.400;gaps=1,80
+lambda_NEB3011 . region 25901 26000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26001 26100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26101 26200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26201 26300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26301 26400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26401 26500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26501 26600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26601 26700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26701 26800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26801 26900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 26901 27000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 27001 27100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 27101 27200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 27201 27300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 27301 27400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 27401 27500 0.00 + . cov=0,0,1;cov2=0.470,0.499;gaps=1,53
+lambda_NEB3011 . region 27501 27600 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 27601 27700 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 27701 27800 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 27801 27900 0.00 + . cov=1,2,2;cov2=1.970,0.171;gaps=0,0
+lambda_NEB3011 . region 27901 28000 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 28001 28100 0.00 + . cov=0,2,2;cov2=1.540,0.842;gaps=1,23
+lambda_NEB3011 . region 28101 28200 0.00 + . cov=0,0,1;cov2=0.200,0.400;gaps=1,80
+lambda_NEB3011 . region 28201 28300 0.00 + . cov=1,2,2;cov2=1.660,0.474;gaps=0,0
+lambda_NEB3011 . region 28301 28400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 28401 28500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 28501 28600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 28601 28700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 28701 28800 0.00 + . cov=0,0,1;cov2=0.020,0.140;gaps=1,98
+lambda_NEB3011 . region 28801 28900 0.00 + . cov=1,2,2;cov2=1.650,0.477;gaps=0,0
+lambda_NEB3011 . region 28901 29000 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 29001 29100 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 29101 29200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 29201 29300 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 29301 29400 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 29401 29500 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 29501 29600 0.00 + . cov=1,1,2;cov2=1.300,0.458;gaps=0,0
+lambda_NEB3011 . region 29601 29700 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 29701 29800 0.00 + . cov=2,2,3;cov2=2.440,0.496;gaps=0,0
+lambda_NEB3011 . region 29801 29900 0.00 + . cov=2,2,3;cov2=2.230,0.421;gaps=0,0
+lambda_NEB3011 . region 29901 30000 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 30001 30100 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 30101 30200 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 30201 30300 0.00 + . cov=0,0,2;cov2=0.770,0.968;gaps=1,61
+lambda_NEB3011 . region 30301 30400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 30401 30500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 30501 30600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 30601 30700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 30701 30800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 30801 30900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 30901 31000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31001 31100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31101 31200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31201 31300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31301 31400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31401 31500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31501 31600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31601 31700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 31701 31800 0.00 + . cov=0,0,1;cov2=0.430,0.495;gaps=1,57
+lambda_NEB3011 . region 31801 31900 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 31901 32000 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 32001 32100 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 32101 32200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 32201 32300 0.00 + . cov=0,1,1;cov2=0.900,0.300;gaps=1,10
+lambda_NEB3011 . region 32301 32400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 32401 32500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 32501 32600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 32601 32700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 32701 32800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 32801 32900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 32901 33000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 33001 33100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 33101 33200 0.00 + . cov=0,0,1;cov2=0.390,0.488;gaps=1,61
+lambda_NEB3011 . region 33201 33300 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 33301 33400 0.00 + . cov=1,9,9;cov2=5.390,3.957;gaps=0,0
+lambda_NEB3011 . region 33401 33500 0.00 + . cov=9,9,10;cov2=9.110,0.313;gaps=0,0
+lambda_NEB3011 . region 33501 33600 0.00 + . cov=9,9,10;cov2=9.390,0.488;gaps=0,0
+lambda_NEB3011 . region 33601 33700 0.00 + . cov=9,9,9;cov2=9.000,0.000;gaps=0,0
+lambda_NEB3011 . region 33701 33800 0.00 + . cov=9,9,9;cov2=9.000,0.000;gaps=0,0
+lambda_NEB3011 . region 33801 33900 0.00 + . cov=9,10,10;cov2=9.930,0.255;gaps=0,0
+lambda_NEB3011 . region 33901 34000 0.00 + . cov=1,3,11;cov2=5.880,4.068;gaps=0,0
+lambda_NEB3011 . region 34001 34100 0.00 + . cov=3,5,5;cov2=4.200,0.949;gaps=0,0
+lambda_NEB3011 . region 34101 34200 0.00 + . cov=4,5,5;cov2=4.700,0.458;gaps=0,0
+lambda_NEB3011 . region 34201 34300 0.00 + . cov=4,4,4;cov2=4.000,0.000;gaps=0,0
+lambda_NEB3011 . region 34301 34400 0.00 + . cov=4,4,5;cov2=4.450,0.497;gaps=0,0
+lambda_NEB3011 . region 34401 34500 0.00 + . cov=4,5,6;cov2=4.890,0.706;gaps=0,0
+lambda_NEB3011 . region 34501 34600 0.00 + . cov=3,3,5;cov2=3.370,0.702;gaps=0,0
+lambda_NEB3011 . region 34601 34700 0.00 + . cov=0,3,3;cov2=1.570,1.485;gaps=1,47
+lambda_NEB3011 . region 34701 34800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 34801 34900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 34901 35000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35001 35100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35101 35200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35201 35300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35301 35400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35401 35500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35501 35600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35601 35700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35701 35800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35801 35900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 35901 36000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36001 36100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36101 36200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36201 36300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36301 36400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36401 36500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36501 36600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36601 36700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36701 36800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36801 36900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 36901 37000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37001 37100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37101 37200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37201 37300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37301 37400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37401 37500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37501 37600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37601 37700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37701 37800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37801 37900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 37901 38000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38001 38100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38101 38200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38201 38300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38301 38400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38401 38500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38501 38600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38601 38700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38701 38800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38801 38900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 38901 39000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39001 39100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39101 39200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39201 39300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39301 39400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39401 39500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39501 39600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39601 39700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39701 39800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39801 39900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 39901 40000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 40001 40100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 40101 40200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 40201 40300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 40301 40400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 40401 40500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 40501 40600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 40601 40700 0.00 + . cov=0,0,2;cov2=0.960,0.999;gaps=1,52
+lambda_NEB3011 . region 40701 40800 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 40801 40900 0.00 + . cov=1,2,2;cov2=1.670,0.470;gaps=0,0
+lambda_NEB3011 . region 40901 41000 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 41001 41100 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 41101 41200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 41201 41300 0.00 + . cov=1,14,14;cov2=10.340,5.615;gaps=0,0
+lambda_NEB3011 . region 41301 41400 0.00 + . cov=14,14,15;cov2=14.300,0.458;gaps=0,0
+lambda_NEB3011 . region 41401 41500 0.00 + . cov=15,15,15;cov2=15.000,0.000;gaps=0,0
+lambda_NEB3011 . region 41501 41600 0.00 + . cov=15,15,15;cov2=15.000,0.000;gaps=0,0
+lambda_NEB3011 . region 41601 41700 0.00 + . cov=1,15,15;cov2=10.170,6.447;gaps=0,0
+lambda_NEB3011 . region 41701 41800 0.00 + . cov=0,0,1;cov2=0.290,0.454;gaps=1,71
+lambda_NEB3011 . region 41801 41900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 41901 42000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 42001 42100 0.00 + . cov=0,0,1;cov2=0.460,0.498;gaps=1,54
+lambda_NEB3011 . region 42101 42200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 42201 42300 0.00 + . cov=0,0,1;cov2=0.450,0.497;gaps=1,55
+lambda_NEB3011 . region 42301 42400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 42401 42500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 42501 42600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 42601 42700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 42701 42800 0.00 + . cov=0,1,1;cov2=0.950,0.218;gaps=1,5
+lambda_NEB3011 . region 42801 42900 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 42901 43000 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 43001 43100 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 43101 43200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 43201 43300 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 43301 43400 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 43401 43500 0.00 + . cov=1,2,3;cov2=1.950,0.497;gaps=0,0
+lambda_NEB3011 . region 43501 43600 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 43601 43700 0.00 + . cov=2,2,2;cov2=2.000,0.000;gaps=0,0
+lambda_NEB3011 . region 43701 43800 0.00 + . cov=1,1,2;cov2=1.070,0.255;gaps=0,0
+lambda_NEB3011 . region 43801 43900 0.00 + . cov=0,1,1;cov2=0.520,0.500;gaps=1,48
+lambda_NEB3011 . region 43901 44000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44001 44100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44101 44200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44201 44300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44301 44400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44401 44500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44501 44600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44601 44700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44701 44800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44801 44900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 44901 45000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 45001 45100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 45101 45200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 45201 45300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 45301 45400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 45401 45500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 45501 45600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 45601 45700 0.00 + . cov=0,0,1;cov2=0.370,0.483;gaps=1,63
+lambda_NEB3011 . region 45701 45800 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 45801 45900 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 45901 46000 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 46001 46100 0.00 + . cov=0,0,1;cov2=0.210,0.407;gaps=1,79
+lambda_NEB3011 . region 46101 46200 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46201 46300 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46301 46400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46401 46500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46501 46600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46601 46700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46701 46800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46801 46900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 46901 47000 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47001 47100 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47101 47200 0.00 + . cov=0,0,1;cov2=0.020,0.140;gaps=1,98
+lambda_NEB3011 . region 47201 47300 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 47301 47400 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47401 47500 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47501 47600 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47601 47700 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47701 47800 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47801 47900 0.00 + . cov=0,0,0;cov2=0.000,0.000;gaps=1,100
+lambda_NEB3011 . region 47901 48000 0.00 + . cov=0,1,1;cov2=0.880,0.325;gaps=1,12
+lambda_NEB3011 . region 48001 48100 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 48101 48200 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 48201 48300 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 48301 48400 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
+lambda_NEB3011 . region 48401 48502 0.00 + . cov=1,1,1;cov2=1.000,0.000;gaps=0,0
diff --git a/data/GFF/alignment_summary_consensus.gff b/data/GFF/alignment_summary_consensus.gff
new file mode 100644
index 0000000..e91b1f7
--- /dev/null
+++ b/data/GFF/alignment_summary_consensus.gff
@@ -0,0 +1,106 @@
+##gff-version 3
+##date 2016-02-10T14:51:19
+##source PACBIO_AlignmentSummary 1.0
+##source-commandline summarize_coverage.py --resolved-tool-contract resolved-tool-contract.json
+##sequence-header lambda_NEB3011 lambda_NEB3011
+##sequence-region lambda_NEB3011 1 48502
+##source GenomicConsensus 2.0.0
+##pacbio-alignment-summary-version 0.6
+##source-commandline /pbi/dept/secondary/builds/3.0.4/smrtanalysis/prebuilt/linux/x86-64/libc-2.5/gcc-4.8.4/ubuntu-1404/173065/smrtanalysis_3.0.4.173065.installdir/private/pacbio/pythonpkgs/GenomicConsensus/binwrap/../../../../../private/pacbio/pythonpkgs/GenomicConsensus/bin/summarizeConsensus --resolved-tool-contract resolved-tool-contract.json
+lambda_NEB3011 . region 1 500 0.00 + . cov=49,91,101;cov2=90.502,6.470;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 501 1000 0.00 + . cov=100,127,151;cov2=122.438,17.046;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 1001 1500 0.00 + . cov=127,160,209;cov2=171.850,23.623;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 1501 2000 0.00 + . cov=185,265,350;cov2=274.568,38.682;gaps=0,0;cQv=20,20,20;del=0;ins=1;sub=0
+lambda_NEB3011 . region 2001 2500 0.00 + . cov=308,364,455;cov2=372.406,37.221;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 2501 3000 0.00 + . cov=320,353,446;cov2=360.066,35.763;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 3001 3500 0.00 + . cov=180,246,330;cov2=253.032,38.331;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 3501 4000 0.00 + . cov=195,225,316;cov2=233.552,33.991;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 4001 4500 0.00 + . cov=157,216,319;cov2=217.900,36.644;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 4501 5000 0.00 + . cov=169,258,317;cov2=256.412,50.398;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 5001 5500 0.00 + . cov=161,213,280;cov2=224.480,33.049;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 5501 6000 0.00 + . cov=268,328,385;cov2=328.710,33.056;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 6001 6500 0.00 + . cov=286,402,435;cov2=386.484,43.509;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 6501 7000 0.00 + . cov=217,269,321;cov2=270.580,29.099;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 7001 7500 0.00 + . cov=245,312,380;cov2=315.726,32.046;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 7501 8000 0.00 + . cov=176,218,271;cov2=229.116,31.685;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 8001 8500 0.00 + . cov=151,202,256;cov2=202.208,28.555;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 8501 9000 0.00 + . cov=204,236,268;cov2=231.834,16.947;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 9001 9500 0.00 + . cov=269,338,365;cov2=327.368,25.019;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=1
+lambda_NEB3011 . region 9501 10000 0.00 + . cov=243,322,380;cov2=316.668,28.057;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 10001 10500 0.00 + . cov=179,246,317;cov2=240.464,46.679;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 10501 11000 0.00 + . cov=146,198,243;cov2=191.940,21.456;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 11001 11500 0.00 + . cov=218,252,279;cov2=249.698,16.401;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 11501 12000 0.00 + . cov=242,272,295;cov2=274.188,14.782;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 12001 12500 0.00 + . cov=227,259,311;cov2=257.820,19.488;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 12501 13000 0.00 + . cov=225,260,297;cov2=260.748,20.330;gaps=0,0;cQv=20,20,20;del=0;ins=1;sub=0
+lambda_NEB3011 . region 13001 13500 0.00 + . cov=171,240,297;cov2=241.808,36.791;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 13501 14000 0.00 + . cov=180,233,273;cov2=233.436,24.094;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 14001 14500 0.00 + . cov=208,252,333;cov2=255.078,23.409;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 14501 15000 0.00 + . cov=196,289,341;cov2=279.358,42.049;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 15001 15500 0.00 + . cov=197,297,326;cov2=276.706,43.330;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 15501 16000 0.00 + . cov=182,240,290;cov2=242.270,35.844;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 16001 16500 0.00 + . cov=190,224,271;cov2=227.482,19.057;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 16501 17000 0.00 + . cov=267,314,372;cov2=311.484,26.776;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 17001 17500 0.00 + . cov=291,332,390;cov2=337.172,27.122;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 17501 18000 0.00 + . cov=268,346,382;cov2=336.360,30.927;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 18001 18500 0.00 + . cov=291,322,361;cov2=326.836,19.506;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 18501 19000 0.00 + . cov=261,289,323;cov2=284.620,17.287;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 19001 19500 0.00 + . cov=181,227,332;cov2=247.282,47.914;gaps=0,0;cQv=20,20,20;del=1;ins=0;sub=0
+lambda_NEB3011 . region 19501 20000 0.00 + . cov=219,281,304;cov2=271.088,22.547;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 20001 20500 0.00 + . cov=280,339,394;cov2=330.854,34.915;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 20501 21000 0.00 + . cov=342,372,432;cov2=383.126,30.093;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 21001 21500 0.00 + . cov=168,284,424;cov2=276.126,73.598;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 21501 22000 0.00 + . cov=128,150,202;cov2=154.580,19.151;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 22001 22500 0.00 + . cov=134,170,177;cov2=163.336,13.015;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 22501 23000 0.00 + . cov=169,229,304;cov2=233.472,43.275;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 23001 23500 0.00 + . cov=283,364,388;cov2=344.762,39.062;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 23501 24000 0.00 + . cov=205,246,330;cov2=253.078,42.092;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 24001 24500 0.00 + . cov=249,322,368;cov2=323.258,26.928;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 24501 25000 0.00 + . cov=279,326,370;cov2=331.720,21.054;gaps=0,0;cQv=20,20,20;del=0;ins=1;sub=0
+lambda_NEB3011 . region 25001 25500 0.00 + . cov=220,279,462;cov2=323.392,79.100;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 25501 26000 0.00 + . cov=369,488,552;cov2=475.902,52.870;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 26001 26500 0.00 + . cov=221,268,368;cov2=268.544,28.232;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 26501 27000 0.00 + . cov=247,304,506;cov2=343.504,80.233;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 27001 27500 0.00 + . cov=485,577,622;cov2=566.380,32.134;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 27501 28000 0.00 + . cov=427,480,618;cov2=508.008,60.657;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 28001 28500 0.00 + . cov=317,359,444;cov2=363.154,23.458;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 28501 29000 0.00 + . cov=230,357,389;cov2=330.050,53.604;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 29001 29500 0.00 + . cov=241,273,351;cov2=273.834,23.984;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 29501 30000 0.00 + . cov=213,248,352;cov2=260.496,37.641;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 30001 30500 0.00 + . cov=256,312,353;cov2=305.746,30.544;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 30501 31000 0.00 + . cov=203,280,344;cov2=278.232,28.766;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 31001 31500 0.00 + . cov=186,236,243;cov2=226.550,17.669;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 31501 32000 0.00 + . cov=151,262,287;cov2=239.110,42.423;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 32001 32500 0.00 + . cov=127,186,198;cov2=176.916,20.022;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 32501 33000 0.00 + . cov=191,322,347;cov2=288.392,55.792;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 33001 33500 0.00 + . cov=141,281,328;cov2=247.040,64.494;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 33501 34000 0.00 + . cov=176,256,363;cov2=255.372,61.858;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 34001 34500 0.00 + . cov=254,335,401;cov2=324.832,48.078;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 34501 35000 0.00 + . cov=174,241,296;cov2=242.956,25.408;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 35001 35500 0.00 + . cov=295,402,464;cov2=390.042,53.516;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 35501 36000 0.00 + . cov=404,482,555;cov2=481.176,51.538;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 36001 36500 0.00 + . cov=295,366,418;cov2=363.634,28.498;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 36501 37000 0.00 + . cov=307,380,425;cov2=374.140,32.156;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 37001 37500 0.00 + . cov=187,261,344;cov2=249.096,48.607;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 37501 38000 0.00 + . cov=169,209,337;cov2=226.880,43.702;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 38001 38500 0.00 + . cov=238,321,339;cov2=307.618,27.202;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 38501 39000 0.00 + . cov=230,276,306;cov2=275.444,19.108;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 39001 39500 0.00 + . cov=193,267,278;cov2=255.792,22.974;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 39501 40000 0.00 + . cov=169,204,246;cov2=202.452,22.566;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 40001 40500 0.00 + . cov=193,262,289;cov2=247.608,35.968;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 40501 41000 0.00 + . cov=273,298,368;cov2=312.362,29.436;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 41001 41500 0.00 + . cov=275,320,354;cov2=314.668,28.350;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 41501 42000 0.00 + . cov=289,326,399;cov2=336.002,29.788;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 42001 42500 0.00 + . cov=368,410,516;cov2=431.598,48.413;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 42501 43000 0.00 + . cov=333,423,557;cov2=439.252,61.237;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 43001 43500 0.00 + . cov=253,293,372;cov2=297.768,24.497;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 43501 44000 0.00 + . cov=234,256,349;cov2=267.674,34.289;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 44001 44500 0.00 + . cov=130,207,242;cov2=199.206,31.665;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 44501 45000 0.00 + . cov=163,203,264;cov2=206.454,23.094;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 45001 45500 0.00 + . cov=177,199,235;cov2=205.556,16.965;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 45501 46000 0.00 + . cov=168,219,283;cov2=230.132,25.942;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 46001 46500 0.00 + . cov=251,335,411;cov2=331.218,53.726;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 46501 47000 0.00 + . cov=267,376,422;cov2=363.474,38.282;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 47001 47500 0.00 + . cov=240,268,354;cov2=274.806,30.188;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 47501 48000 0.00 + . cov=88,145,243;cov2=173.026,54.269;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
+lambda_NEB3011 . region 48001 48502 0.00 + . cov=20,54,88;cov2=58.958,10.589;gaps=0,0;cQv=20,20,20;del=0;ins=0;sub=0
diff --git a/data/GFF/variants.gff b/data/GFF/variants.gff
new file mode 100644
index 0000000..e84e372
--- /dev/null
+++ b/data/GFF/variants.gff
@@ -0,0 +1,14 @@
+##gff-version 3
+##pacbio-variant-version 2.1
+##date Wed Feb 10 07:19:31 2016
+##feature-ontology http://song.cvs.sourceforge.net/*checkout*/song/ontology/sofa.obo?revision=1.12
+##source GenomicConsensus 2.0.0
+##source-commandline variantCaller --resolved-tool-contract resolved-tool-contract.json
+##source-alignment-file aligned.subreads.alignmentset.xml
+##source-reference-file references/lambdaNEB/lambdaNEB.referenceset.xml
+##sequence-region lambda_NEB3011 1 48502
+lambda_NEB3011 . insertion 1536 1536 . . . reference=.;variantSeq=C;coverage=18;confidence=47
+lambda_NEB3011 . substitution 9251 9251 . . . reference=G;variantSeq=A;coverage=21;confidence=43
+lambda_NEB3011 . insertion 12534 12534 . . . reference=.;variantSeq=T;coverage=20;confidence=44
+lambda_NEB3011 . deletion 19183 19183 . . . reference=T;variantSeq=.;coverage=20;confidence=48
+lambda_NEB3011 . insertion 24872 24872 . . . reference=.;variantSeq=T;coverage=17;confidence=41
diff --git a/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.1.bax.h5 b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.1.bax.h5
new file mode 100644
index 0000000..76ac690
--- /dev/null
+++ b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.1.bax.h5
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:70396e9e03dc16506135faa801e282299aff5f0e058746c44d146594ba5a9214
+size 856704
diff --git a/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.2.bax.h5 b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.2.bax.h5
new file mode 100644
index 0000000..c178d96
--- /dev/null
+++ b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.2.bax.h5
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3e52e0545846b93b155746e2bae6eed7edce1055cf1c8fc0dde47fe788c1ce2f
+size 404288
diff --git a/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.3.bax.h5 b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.3.bax.h5
new file mode 100644
index 0000000..5926d41
--- /dev/null
+++ b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.3.bax.h5
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:48539d67fe9a3410654d647c688ec615cfb1e4f67d889348356524bda0afc294
+size 610688
diff --git a/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.all.hdfsubreadset.xml b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.all.hdfsubreadset.xml
new file mode 100644
index 0000000..ac31616
--- /dev/null
+++ b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.all.hdfsubreadset.xml
@@ -0,0 +1,12 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<HdfSubreadSet CreatedAt="2016-06-17T11:36:54" MetaType="PacBio.DataSet.HdfSubreadSet" Name="hdfsubreads" Tags="" TimeStampedName="pacbio_dataset_hdfsubreadset-160617_183654177" UniqueId="43c51ade-cb34-f699-9b36-f39744bc421d" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.BaxFile" ResourceId="m140912_020930_00114_c100702482550000001823141103261590_s1_p0.1.bax.h5" TimeStampedName="pacbio_subreadfile_baxfile-160617_183653985" UniqueId="1226d3bb-d249-45d2-ba90-8e4dded376ed"/>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.BaxFile" ResourceId="m140912_020930_00114_c100702482550000001823141103261590_s1_p0.2.bax.h5" TimeStampedName="pacbio_subreadfile_baxfile-160617_183653985" UniqueId="0d9f8a76-9d1a-410f-a700-174df17972ec"/>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.BaxFile" ResourceId="m140912_020930_00114_c100702482550000001823141103261590_s1_p0.3.bax.h5" TimeStampedName="pacbio_subreadfile_baxfile-160617_183653985" UniqueId="c4521ac9-ba73-478a-939b-9a808fb10d5d"/>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>58468</TotalLength>
+ <NumRecords>9</NumRecords>
+ </DataSetMetadata>
+</HdfSubreadSet>
diff --git a/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.bas.h5 b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.bas.h5
new file mode 100644
index 0000000..f947909
--- /dev/null
+++ b/data/HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.bas.h5
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:2395ba67ec75fc808ebba3104db6a37132ff608db2a5c4a43a96b692507a9521
+size 1318480
diff --git a/data/ReferenceSet/lambdaNEB/reference.info.xml b/data/ReferenceSet/lambdaNEB/reference.info.xml
new file mode 100644
index 0000000..06d4a2c
--- /dev/null
+++ b/data/ReferenceSet/lambdaNEB/reference.info.xml
@@ -0,0 +1,24 @@
+<?xml version="1.0" encoding="UTF-8" standalone="yes"?>
+<reference_info xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" id="lambdaNEB" last_modified="2016-05-31T13:28:56-0700" xsi:schemaLocation="http://www.w3.org/2001/XMLSchema-instance">
+ <organism>
+ <name>lambda</name>
+ <ploidy>haploid</ploidy>
+ </organism>
+ <reference>
+ <description>Fasta converter v0.4.0 to PacBio Reference.info.xml and ReferenceSet XML</description>
+ <file format="text/fasta">./sequence/lambdaNEB.fasta</file>
+ <index_file type="indexer">./sequence/lambdaNEB.fasta.index</index_file>
+ <index_file type="sawriter">./sequence/lambdaNEB.fasta.sa</index_file>
+ <index_file type="sam_idx">./sequence/lambdaNEB.fasta.fai</index_file>
+ <index_file type="fasta_contig_index">./sequence/lambdaNEB.fasta.contig.index</index_file>
+ <max_contig_length>48502</max_contig_length>
+ <num_contigs>1</num_contigs>
+ <type>sample</type>
+ </reference>
+ <contigs>
+ <contig displayName="lambda_NEB3011" id="ref000001" length="48502">
+ <digest type="md5">a1319ff90e994c8190a4fe6569d0822a</digest>
+ <header>lambda_NEB3011</header>
+ </contig>
+ </contigs>
+</reference_info>
diff --git a/data/ReferenceSet/lambdaNEB/referenceset.xml b/data/ReferenceSet/lambdaNEB/referenceset.xml
new file mode 100644
index 0000000..d548f98
--- /dev/null
+++ b/data/ReferenceSet/lambdaNEB/referenceset.xml
@@ -0,0 +1,22 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ReferenceSet xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:ns0="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:ns4="http://pacificbiosciences.com/PacBioCollectionMetadata.xsd" xmlns:ns1="http://pacificbiosciences.com/PacBioSampleInfo.xsd" xmlns:ns6="http://pacificbiosciences.com/PacBioReagentKit.xsd" Name="lambdaNEB" Description="Converted Reference lambdaNEB" Tags="converted, reference" Version="3.0.1" CreatedAt="2016-05-31T13:28:56.786-07:00" ModifiedAt= [...]
+ <ns0:ExternalResources>
+ <ns0:ExternalResource Name="Fasta lambdaNEB" Description="Converted with v0.5.0 Reference Name lambdaNEB" Tags="converted, reference" ResourceId="./sequence/lambdaNEB.fasta" CreatedAt="2016-05-31T13:28:56.786-07:00" ModifiedAt="2016-05-31T13:28:56.786-07:00" UniqueId="7cf17215-9c3a-4230-aba3-eb8548d8e960" MetaType="PacBio.ReferenceFile.ReferenceFastaFile" TimeStampedName="fasta_160531_132856">
+ <ns0:FileIndices>
+ <ns0:FileIndex ResourceId="./sequence/lambdaNEB.fasta.index" UniqueId="8936c5e7-c2af-45e8-a0ff-bc21913d3b8f" MetaType="PacBio.Index.Indexer" TimeStampedName="index_160531_132856"/>
+ <ns0:FileIndex ResourceId="./sequence/lambdaNEB.fasta.sa" UniqueId="9b480a1c-e461-4ef9-b3f6-8ce0ab1025fa" MetaType="PacBio.Index.SaWriterIndex" TimeStampedName="index_160531_132856"/>
+ <ns0:FileIndex ResourceId="./sequence/lambdaNEB.fasta.fai" UniqueId="2d0b8ed7-717c-4daa-a7e6-42464cf2c02f" MetaType="PacBio.Index.SamIndex" TimeStampedName="index_160531_132856"/>
+ <ns0:FileIndex ResourceId="./sequence/lambdaNEB.fasta.contig.index" UniqueId="c986b2ff-77dc-4164-aa6d-bae6163bfe27" MetaType="PacBio.Index.FastaContigIndex" TimeStampedName="index_160531_132856"/>
+ </ns0:FileIndices>
+ </ns0:ExternalResource>
+ </ns0:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>48502</TotalLength>
+ <NumRecords>1</NumRecords>
+ <Organism>lambda</Organism>
+ <Ploidy>haploid</Ploidy>
+ <Contigs>
+ <Contig Name="ref000001" Description="lambda_NEB3011" Length="48502" Digest="a1319ff90e994c8190a4fe6569d0822a"/>
+ </Contigs>
+ </DataSetMetadata>
+</ReferenceSet>
diff --git a/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta
new file mode 100644
index 0000000..8cf0bce
--- /dev/null
+++ b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta
@@ -0,0 +1,810 @@
+>lambda_NEB3011
+GGGCGGCGACCTCGCGGGTTTTCGCTATTTATGAAAATTTTCCGGTTTAAGGCGTTTCCG
+TTCTTCTTCGTCATAACTTAATGTTTTTATTTAAAATACCCTCTGAAAAGAAAGGAAACG
+ACAGGTGCTGAAAGCGAGCTTTTTGGCCTCTGTCGTTTCCTTTCTCTGTTTTTGTCCGTG
+GAATGAACAATGGAAGTCAACAAAAAGCAGCTGGCTGACATTTTCGGTGCGAGTATCCGT
+ACCATTCAGAACTGGCAGGAACAGGGAATGCCCGTTCTGCGAGGCGGTGGCAAGGGTAAT
+GAGGTGCTTTATGACTCTGCCGCCGTCATAAAATGGTATGCCGAAAGGGATGCTGAAATT
+GAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCGGCAGGCCAGCGAGGCAGATCTC
+CAGCCAGGAACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGACGCACAG
+GAACTGAAGAATGCCAGAGACTCCGCTGAAGTGGTGGAAACCGCATTCTGTACTTTCGTG
+CTGTCGCGGATCGCAGGTGAAATTGCCAGTATTCTCGACGGGCTCCCCCTGTCGGTGCAG
+CGGCGTTTTCCGGAACTGGAAAACCGACATGTTGATTTCCTGAAACGGGATATCATCAAA
+GCCATGAACAAAGCAGCCGCGCTGGATGAACTGATACCGGGGTTGCTGAGTGAATATATC
+GAACAGTCAGGTTAACAGGCTGCGGCATTTTGTCCGCGCCGGGCTTCGCTCACTGTTCAG
+GCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATC
+CGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGAT
+GGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAA
+AATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATCTG
+GTTGCCGACGGATGGTGATGCCGAGAACTTTATGAAAACCCACGTTGAGCCGACTATTCG
+TGATATTCCGTCGCTGCTGGCGCTGGCCCCGTGGTATGGCAAAAAGCACCGGGATAACAC
+GCTCACCATGAAGCGTTTCACTAATGGGCGTGGCTTCTGGTGCCTGGGCGGTAAAGCGGC
+AAAAAACTACCGTGAAAAGTCGGTGGATGTGGCGGGTTATGATGAACTTGCTGCTTTTGA
+TGATGATATTGAACAGGAAGGCTCTCCGACGTTCCTGGGTGACAAGCGTATTGAAGGCTC
+GGTCTGGCCAAAGTCCATCCGTGGCTCCACGCCAAAAGTGAGAGGCACCTGTCAGATTGA
+GCGTGCAGCCAGTGAATCCCCGCATTTTATGCGTTTTCATGTTGCCTGCCCGCATTGCGG
+GGAGGAGCAGTATCTTAAATTTGGCGACAAAGAGACGCCGTTTGGCCTCAAATGGACGCC
+GGATGACCCCTCCAGCGTGTTTTATCTCTGCGAGCATAATGCCTGCGTCATCCGCCAGCA
+GGAGCTGGACTTTACTGATGCCCGTTATATCTGCGAAAAGACCGGGATCTGGACCCGTGA
+TGGCATTCTCTGGTTTTCGTCATCCGGTGAAGAGATTGAGCCACCTGACAGTGTGACCTT
+TCACATCTGGACAGCGTACAGCCCGTTCACCACCTGGGTGCAGATTGTCAAAGACTGGAT
+GAAAACGAAAGGGGATACGGGAAAACGTAAAACCTTCGTAAACACCACGCTCGGTGAGAC
+GTGGGAGGCGAAAATTGGCGAACGTCCGGATGCTGAAGTGATGGCAGAGCGGAAAGAGCA
+TTATTCAGCGCCCGTTCCTGACCGTGTGGCTTACCTGACCGCCGGTATCGACTCCCAGCT
+GGACCGCTACGAAATGCGCGTATGGGGATGGGGGCCGGGTGAGGAAAGCTGGCTGATTGA
+CCGGCAGATTATTATGGGCCGCCACGACGATGAACAGACGCTGCTGCGTGTGGATGAGGC
+CATCAATAAAACCTATACCCGCCGGAATGGTGCAGAAATGTCGATATCCCGTATCTGCTG
+GGATACTGGCGGGATTGACCCGACCATTGTGTATGAACGCTCGAAAAAACATGGGCTGTT
+CCGGGTGATCCCCATTAAAGGGGCATCCGTCTACGGAAAGCCGGTGGCCAGCATGCCACG
+TAAGCGAAACAAAAACGGGGTTTACCTTACCGAAATCGGTACGGATACCGCGAAAGAGCA
+GATTTATAACCGCTTCACACTGACGCCGGAAGGGGATGAACCGCTTCCCGGTGCCGTTCA
+CTTCCCGAATAACCCGGATATTTTTGATCTGACCGAAGCGCAGCAGCTGACTGCTGAAGA
+GCAGGTCGAAAAATGGGTGGATGGCAGGAAAAAAATACTGTGGGACAGCAAAAAGCGACG
+CAATGAGGCACTCGACTGCTTCGTTTATGCGCTGGCGGCGCTGCGCATCAGTATTTCCCG
+CTGGCAGCTGGATCTCAGTGCGCTGCTGGCGAGCCTGCAGGAAGAGGATGGTGCAGCAAC
+CAACAAGAAAACACTGGCAGATTACGCCCGTGCCTTATCCGGAGAGGATGAATGACGCGA
+CAGGAAGAACTTGCCGCTGCCCGTGCGGCACTGCATGACCTGATGACAGGTAAACGGGTG
+GCAACAGTACAGAAAGACGGACGAAGGGTGGAGTTTACGGCCACTTCCGTGTCTGACCTG
+AAAAAATATATTGCAGAGCTGGAAGTGCAGACCGGCATGACACAGCGACGCAGGGGACCT
+GCAGGATTTTATGTATGAAAACGCCCACCATTCCCACCCTTCTGGGGCCGGACGGCATGA
+CATCGCTGCGCGAATATGCCGGTTATCACGGCGGTGGCAGCGGATTTGGAGGGCAGTTGC
+GGTCGTGGAACCCACCGAGTGAAAGTGTGGATGCAGCCCTGTTGCCCAACTTTACCCGTG
+GCAATGCCCGCGCAGACGATCTGGTACGCAATAACGGCTATGCCGCCAACGCCATCCAGC
+TGCATCAGGATCATATCGTCGGGTCTTTTTTCCGGCTCAGTCATCGCCCAAGCTGGCGCT
+ATCTGGGCATCGGGGAGGAAGAAGCCCGTGCCTTTTCCCGCGAGGTTGAAGCGGCATGGA
+AAGAGTTTGCCGAGGATGACTGCTGCTGCATTGACGTTGAGCGAAAACGCACGTTTACCA
+TGATGATTCGGGAAGGTGTGGCCATGCACGCCTTTAACGGTGAACTGTTCGTTCAGGCCA
+CCTGGGATACCAGTTCGTCGCGGCTTTTCCGGACACAGTTCCGGATGGTCAGCCCGAAGC
+GCATCAGCAACCCGAACAATACCGGCGACAGCCGGAACTGCCGTGCCGGTGTGCAGATTA
+ATGACAGCGGTGCGGCGCTGGGATATTACGTCAGCGAGGACGGGTATCCTGGCTGGATGC
+CGCAGAAATGGACATGGATACCCCGTGAGTTACCCGGCGGGCGCGCCTCGTTCATTCACG
+TTTTTGAACCCGTGGAGGACGGGCAGACTCGCGGTGCAAATGTGTTTTACAGCGTGATGG
+AGCAGATGAAGATGCTCGACACGCTGCAGAACACGCAGCTGCAGAGCGCCATTGTGAAGG
+CGATGTATGCCGCCACCATTGAGAGTGAGCTGGATACGCAGTCAGCGATGGATTTTATTC
+TGGGCGCGAACAGTCAGGAGCAGCGGGAAAGGCTGACCGGCTGGATTGGTGAAATTGCCG
+CGTATTACGCCGCAGCGCCGGTCCGGCTGGGAGGCGCAAAAGTACCGCACCTGATGCCGG
+GTGACTCACTGAACCTGCAGACGGCTCAGGATACGGATAACGGCTACTCCGTGTTTGAGC
+AGTCACTGCTGCGGTATATCGCTGCCGGGCTGGGTGTCTCGTATGAGCAGCTTTCCCGGA
+ATTACGCCCAGATGAGCTACTCCACGGCACGGGCCAGTGCGAACGAGTCGTGGGCGTACT
+TTATGGGGCGGCGAAAATTCGTCGCATCCCGTCAGGCGAGCCAGATGTTTCTGTGCTGGC
+TGGAAGAGGCCATCGTTCGCCGCGTGGTGACGTTACCTTCAAAAGCGCGCTTCAGTTTTC
+AGGAAGCCCGCAGTGCCTGGGGGAACTGCGACTGGATAGGCTCCGGTCGTATGGCCATCG
+ATGGTCTGAAAGAAGTTCAGGAAGCGGTGATGCTGATAGAAGCCGGACTGAGTACCTACG
+AGAAAGAGTGCGCAAAACGCGGTGACGACTATCAGGAAATTTTTGCCCAGCAGGTCCGTG
+AAACGATGGAGCGCCGTGCAGCCGGTCTTAAACCGCCCGCCTGGGCGGCTGCAGCATTTG
+AATCCGGGCTGCGACAATCAACAGAGGAGGAGAAGAGTGACAGCAGAGCTGCGTAATCTC
+CCGCATATTGCCAGCATGGCCTTTAATGAGCCGCTGATGCTTGAACCCGCCTATGCGCGG
+GTTTTCTTTTGTGCGCTTGCAGGCCAGCTTGGGATCAGCAGCCTGACGGATGCGGTGTCC
+GGCGACAGCCTGACTGCCCAGGAGGCACTCGCGACGCTGGCATTATCCGGTGATGATGAC
+GGACCACGACAGGCCCGCAGTTATCAGGTCATGAACGGCATCGCCGTGCTGCCGGTGTCC
+GGCACGCTGGTCAGCCGGACGCGGGCGCTGCAGCCGTACTCGGGGATGACCGGTTACAAC
+GGCATTATCGCCCGTCTGCAACAGGCTGCCAGCGATCCGATGGTGGACGGCATTCTGCTC
+GATATGGACACGCCCGGCGGGATGGTGGCGGGGGCATTTGACTGCGCTGACATCATCGCC
+CGTGTGCGTGACATAAAACCGGTATGGGCGCTTGCCAACGACATGAACTGCAGTGCAGGT
+CAGTTGCTTGCCAGTGCCGCCTCCCGGCGTCTGGTCACGCAGACCGCCCGGACAGGCTCC
+ATCGGCGTCATGATGGCTCACAGTAATTACGGTGCTGCGCTGGAGAAACAGGGTGTGGAA
+ATCACGCTGATTTACAGCGGCAGCCATAAGGTGGATGGCAACCCCTACAGCCATCTTCCG
+GATGACGTCCGGGAGACACTGCAGTCCCGGATGGACGCAACCCGCCAGATGTTTGCGCAG
+AAGGTGTCGGCATATACCGGCCTGTCCGTGCAGGTTGTGCTGGATACCGAGGCTGCAGTG
+TACAGCGGTCAGGAGGCCATTGATGCCGGACTGGCTGATGAACTTGTTAACAGCACCGAT
+GCGATCACCGTCATGCGTGATGCACTGGATGCACGTAAATCCCGTCTCTCAGGAGGGCGA
+ATGACCAAAGAGACTCAATCAACAACTGTTTCAGCCACTGCTTCGCAGGCTGACGTTACT
+GACGTGGTGCCAGCGACGGAGGGCGAGAACGCCAGCGCGGCGCAGCCGGACGTGAACGCG
+CAGATCACCGCAGCGGTTGCGGCAGAAAACAGCCGCATTATGGGGATCCTCAACTGTGAG
+GAGGCTCACGGACGCGAAGAACAGGCACGCGTGCTGGCAGAAACCCCCGGTATGACCGTG
+AAAACGGCCCGCCGCATTCTGGCCGCAGCACCACAGAGTGCACAGGCGCGCAGTGACACT
+GCGCTGGATCGTCTGATGCAGGGGGCACCGGCACCGCTGGCTGCAGGTAACCCGGCATCT
+GATGCCGTTAACGATTTGCTGAACACACCAGTGTAAGGGATGTTTATGACGAGCAAAGAA
+ACCTTTACCCATTACCAGCCGCAGGGCAACAGTGACCCGGCTCATACCGCAACCGCGCCC
+GGCGGATTGAGTGCGAAAGCGCCTGCAATGACCCCGCTGATGCTGGACACCTCCAGCCGT
+AAGCTGGTTGCGTGGGATGGCACCACCGACGGTGCTGCCGTTGGCATTCTTGCGGTTGCT
+GCTGACCAGACCAGCACCACGCTGACGTTCTACAAGTCCGGCACGTTCCGTTATGAGGAT
+GTGCTCTGGCCGGAGGCTGCCAGCGACGAGACGAAAAAACGGACCGCGTTTGCCGGAACG
+GCAATCAGCATCGTTTAACTTTACCCTTCATCACTAAAGGCCGCCTGTGCGGCTTTTTTT
+ACGGGATTTTTTTATGTCGATGTACACAACCGCCCAACTGCTGGCGGCAAATGAGCAGAA
+ATTTAAGTTTGATCCGCTGTTTCTGCGTCTCTTTTTCCGTGAGAGCTATCCCTTCACCAC
+GGAGAAAGTCTATCTCTCACAAATTCCGGGACTGGTAAACATGGCGCTGTACGTTTCGCC
+GATTGTTTCCGGTGAGGTTATCCGTTCCCGTGGCGGCTCCACCTCTGAATTTACGCCGGG
+ATATGTCAAGCCGAAGCATGAAGTGAATCCGCAGATGACCCTGCGTCGCCTGCCGGATGA
+AGATCCGCAGAATCTGGCGGACCCGGCTTACCGCCGCCGTCGCATCATCATGCAGAACAT
+GCGTGACGAAGAGCTGGCCATTGCTCAGGTCGAAGAGATGCAGGCAGTTTCTGCCGTGCT
+TAAGGGCAAATACACCATGACCGGTGAAGCCTTCGATCCGGTTGAGGTGGATATGGGCCG
+CAGTGAGGAGAATAACATCACGCAGTCCGGCGGCACGGAGTGGAGCAAGCGTGACAAGTC
+CACGTATGACCCGACCGACGATATCGAAGCCTACGCGCTGAACGCCAGCGGTGTGGTGAA
+TATCATCGTGTTCGATCCGAAAGGCTGGGCGCTGTTCCGTTCCTTCAAAGCCGTCAAGGA
+GAAGCTGGATACCCGTCGTGGCTCTAATTCCGAGCTGGAGACAGCGGTGAAAGACCTGGG
+CAAAGCGGTGTCCTATAAGGGGATGTATGGCGATGTGGCCATCGTCGTGTATTCCGGACA
+GTACGTGGAAAACGGCGTCAAAAAGAACTTCCTGCCGGACAACACGATGGTGCTGGGGAA
+CACTCAGGCACGCGGTCTGCGCACCTATGGCTGCATTCAGGATGCGGACGCACAGCGCGA
+AGGCATTAACGCCTCTGCCCGTTACCCGAAAAACTGGGTGACCACCGGCGATCCGGCGCG
+TGAGTTCACCATGATTCAGTCAGCACCGCTGATGCTGCTGGCTGACCCTGATGAGTTCGT
+GTCCGTACAACTGGCGTAATCATGGCCCTTCGGGGCCATTGTTTCTCTGTGGAGGAGTCC
+ATGACGAAAGATGAACTGATTGCCCGTCTCCGCTCGCTGGGTGAACAACTGAACCGTGAT
+GTCAGCCTGACGGGGACGAAAGAAGAACTGGCGCTCCGTGTGGCAGAGCTGAAAGAGGAG
+CTTGATGACACGGATGAAACTGCCGGTCAGGACACCCCTCTCAGCCGGGAAAATGTGCTG
+ACCGGACATGAAAATGAGGTGGGATCAGCGCAGCCGGATACCGTGATTCTGGATACGTCT
+GAACTGGTCACGGTCGTGGCACTGGTGAAGCTGCATACTGATGCACTTCACGCCACGCGG
+GATGAACCTGTGGCATTTGTGCTGCCGGGAACGGCGTTTCGTGTCTCTGCCGGTGTGGCA
+GCCGAAATGACAGAGCGCGGCCTGGCCAGAATGCAATAACGGGAGGCGCTGTGGCTGATT
+TCGATAACCTGTTCGATGCTGCCATTGCCCGCGCCGATGAAACGATACGCGGGTACATGG
+GAACGTCAGCCACCATTACATCCGGTGAGCAGTCAGGTGCGGTGATACGTGGTGTTTTTG
+ATGACCCTGAAAATATCAGCTATGCCGGACAGGGCGTGCGCGTTGAAGGCTCCAGCCCGT
+CCCTGTTTGTCCGGACTGATGAGGTGCGGCAGCTGCGGCGTGGAGACACGCTGACCATCG
+GTGAGGAAAATTTCTGGGTAGATCGGGTTTCGCCGGATGATGGCGGAAGTTGTCATCTCT
+GGCTTGGACGGGGCGTACCGCCTGCCGTTAACCGTCGCCGCTGAAAGGGGGATGTATGGC
+CATAAAAGGTCTTGAGCAGGCCGTTGAAAACCTCAGCCGTATCAGCAAAACGGCGGTGCC
+TGGTGCCGCCGCAATGGCCATTAACCGCGTTGCTTCATCCGCGATATCGCAGTCGGCGTC
+ACAGGTTGCCCGTGAGACAAAGGTACGCCGGAAACTGGTAAAGGAAAGGGCCAGGCTGAA
+AAGGGCCACGGTCAAAAATCCGCAGGCCAGAATCAAAGTTAACCGGGGGGATTTGCCCGT
+AATCAAGCTGGGTAATGCGCGGGTTGTCCTTTCGCGCCGCAGGCGTCGTAAAAAGGGGCA
+GCGTTCATCCCTGAAAGGTGGCGGCAGCGTGCTTGTGGTGGGTAACCGTCGTATTCCCGG
+CGCGTTTATTCAGCAACTGAAAAATGGCCGGTGGCATGTCATGCAGCGTGTGGCTGGGAA
+AAACCGTTACCCCATTGATGTGGTGAAAATCCCGATGGCGGTGCCGCTGACCACGGCGTT
+TAAACAAAATATTGAGCGGATACGGCGTGAACGTCTTCCGAAAGAGCTGGGCTATGCGCT
+GCAGCATCAACTGAGGATGGTAATAAAGCGATGAAACATACTGAACTCCGTGCAGCCGTA
+CTGGATGCACTGGAGAAGCATGACACCGGGGCGACGTTTTTTGATGGTCGCCCCGCTGTT
+TTTGATGAGGCGGATTTTCCGGCAGTTGCCGTTTATCTCACCGGCGCTGAATACACGGGC
+GAAGAGCTGGACAGCGATACCTGGCAGGCGGAGCTGCATATCGAAGTTTTCCTGCCTGCT
+CAGGTGCCGGATTCAGAGCTGGATGCGTGGATGGAGTCCCGGATTTATCCGGTGATGAGC
+GATATCCCGGCACTGTCAGATTTGATCACCAGTATGGTGGCCAGCGGCTATGACTACCGG
+CGCGACGATGATGCGGGCTTGTGGAGTTCAGCCGATCTGACTTATGTCATTACCTATGAA
+ATGTGAGGACGCTATGCCTGTACCAAATCCTACAATGCCGGTGAAAGGTGCCGGGACCAC
+CCTGTGGGTTTATAAGGGGAGCGGTGACCCTTACGCGAATCCGCTTTCAGACGTTGACTG
+GTCGCGTCTGGCAAAAGTTAAAGACCTGACGCCCGGCGAACTGACCGCTGAGTCCTATGA
+CGACAGCTATCTCGATGATGAAGATGCAGACTGGACTGCGACCGGGCAGGGGCAGAAATC
+TGCCGGAGATACCAGCTTCACGCTGGCGTGGATGCCCGGAGAGCAGGGGCAGCAGGCGCT
+GCTGGCGTGGTTTAATGAAGGCGATACCCGTGCCTATAAAATCCGCTTCCCGAACGGCAC
+GGTCGATGTGTTCCGTGGCTGGGTCAGCAGTATCGGTAAGGCGGTGACGGCGAAGGAAGT
+GATCACCCGCACGGTGAAAGTCACCAATGTGGGACGTCCGTCGATGGCAGAAGATCGCAG
+CACGGTAACAGCGGCAACCGGCATGACCGTGACGCCTGCCAGCACCTCGGTGGTGAAAGG
+GCAGAGCACCACGCTGACCGTGGCCTTCCAGCCGGAGGGCGTAACCGACAAGAGCTTTCG
+TGCGGTGTCTGCGGATAAAACAAAAGCCACCGTGTCGGTCAGTGGTATGACCATCACCGT
+GAACGGCGTTGCTGCAGGCAAGGTCAACATTCCGGTTGTATCCGGTAATGGTGAGTTTGC
+TGCGGTTGCAGAAATTACCGTCACCGCCAGTTAATCCGGAGAGTCAGCGATGTTCCTGAA
+AACCGAATCATTTGAACATAACGGTGTGACCGTCACGCTTTCTGAACTGTCAGCCCTGCA
+GCGCATTGAGCATCTCGCCCTGATGAAACGGCAGGCAGAACAGGCGGAGTCAGACAGCAA
+CCGGAAGTTTACTGTGGAAGACGCCATCAGAACCGGCGCGTTTCTGGTGGCGATGTCCCT
+GTGGCATAACCATCCGCAGAAGACGCAGATGCCGTCCATGAATGAAGCCGTTAAACAGAT
+TGAGCAGGAAGTGCTTACCACCTGGCCCACGGAGGCAATTTCTCATGCTGAAAACGTGGT
+GTACCGGCTGTCTGGTATGTATGAGTTTGTGGTGAATAATGCCCCTGAACAGACAGAGGA
+CGCCGGGCCCGCAGAGCCTGTTTCTGCGGGAAAGTGTTCGACGGTGAGCTGAGTTTTGCC
+CTGAAACTGGCGCGTGAGATGGGGCGACCCGACTGGCGTGCCATGCTTGCCGGGATGTCA
+TCCACGGAGTATGCCGACTGGCACCGCTTTTACAGTACCCATTATTTTCATGATGTTCTG
+CTGGATATGCACTTTTCCGGGCTGACGTACACCGTGCTCAGCCTGTTTTTCAGCGATCCG
+GATATGCATCCGCTGGATTTCAGTCTGCTGAACCGGCGCGAGGCTGACGAAGAGCCTGAA
+GATGATGTGCTGATGCAGAAAGCGGCAGGGCTTGCCGGAGGTGTCCGCTTTGGCCCGGAC
+GGGAATGAAGTTATCCCCGCTTCCCCGGATGTGGCGGACATGACGGAGGATGACGTAATG
+CTGATGACAGTATCAGAAGGGATCGCAGGAGGAGTCCGGTATGGCTGAACCGGTAGGCGA
+TCTGGTCGTTGATTTGAGTCTGGATGCGGCCAGATTTGACGAGCAGATGGCCAGAGTCAG
+GCGTCATTTTTCTGGTACGGAAAGTGATGCGAAAAAAACAGCGGCAGTCGTTGAACAGTC
+GCTGAGCCGACAGGCGCTGGCTGCACAGAAAGCGGGGATTTCCGTCGGGCAGTATAAAGC
+CGCCATGCGTATGCTGCCTGCACAGTTCACCGACGTGGCCACGCAGCTTGCAGGCGGGCA
+AAGTCCGTGGCTGATCCTGCTGCAACAGGGGGGGCAGGTGAAGGACTCCTTCGGCGGGAT
+GATCCCCATGTTCAGGGGGCTTGCCGGTGCGATCACCCTGCCGATGGTGGGGGCCACCTC
+GCTGGCGGTGGCGACCGGTGCGCTGGCGTATGCCTGGTATCAGGGCAACTCAACCCTGTC
+CGATTTCAACAAAACGCTGGTCCTTTCCGGCAATCAGGCGGGACTGACGGCAGATCGTAT
+GCTGGTCCTGTCCAGAGCCGGGCAGGCGGCAGGGCTGACGTTTAACCAGACCAGCGAGTC
+ACTCAGCGCACTGGTTAAGGCGGGGGTAAGCGGTGAGGCTCAGATTGCGTCCATCAGCCA
+GAGTGTGGCGCGTTTCTCCTCTGCATCCGGCGTGGAGGTGGACAAGGTCGCTGAAGCCTT
+CGGGAAGCTGACCACAGACCCGACGTCGGGGCTGACGGCGATGGCTCGCCAGTTCCATAA
+CGTGTCGGCGGAGCAGATTGCGTATGTTGCTCAGTTGCAGCGTTCCGGCGATGAAGCCGG
+GGCATTGCAGGCGGCGAACGAGGCCGCAACGAAAGGGTTTGATGACCAGACCCGCCGCCT
+GAAAGAGAACATGGGCACGCTGGAGACCTGGGCAGACAGGACTGCGCGGGCATTCAAATC
+CATGTGGGATGCGGTGCTGGATATTGGTCGTCCTGATACCGCGCAGGAGATGCTGATTAA
+GGCAGAGGCTGCGTATAAGAAAGCAGACGACATCTGGAATCTGCGCAAGGATGATTATTT
+TGTTAACGATGAAGCGCGGGCGCGTTACTGGGATGATCGTGAAAAGGCCCGTCTTGCGCT
+TGAAGCCGCCCGAAAGAAGGCTGAGCAGCAGACTCAACAGGACAAAAATGCGCAGCAGCA
+GAGCGATACCGAAGCGTCACGGCTGAAATATACCGAAGAGGCGCAGAAGGCTTACGAACG
+GCTGCAGACGCCGCTGGAGAAATATACCGCCCGTCAGGAAGAACTGAACAAGGCACTGAA
+AGACGGGAAAATCCTGCAGGCGGATTACAACACGCTGATGGCGGCGGCGAAAAAGGATTA
+TGAAGCGACGCTGAAAAAGCCGAAACAGTCCAGCGTGAAGGTGTCTGCGGGCGATCGTCA
+GGAAGACAGTGCTCATGCTGCCCTGCTGACGCTTCAGGCAGAACTCCGGACGCTGGAGAA
+GCATGCCGGAGCAAATGAGAAAATCAGCCAGCAGCGCCGGGATTTGTGGAAGGCGGAGAG
+TCAGTTCGCGGTACTGGAGGAGGCGGCGCAACGTCGCCAGCTGTCTGCACAGGAGAAATC
+CCTGCTGGCGCATAAAGATGAGACGCTGGAGTACAAACGCCAGCTGGCTGCACTTGGCGA
+CAAGGTTACGTATCAGGAGCGCCTGAACGCGCTGGCGCAGCAGGCGGATAAATTCGCACA
+GCAGCAACGGGCAAAACGGGCCGCCATTGATGCGAAAAGCCGGGGGCTGACTGACCGGCA
+GGCAGAACGGGAAGCCACGGAACAGCGCCTGAAGGAACAGTATGGCGATAATCCGCTGGC
+GCTGAATAACGTCATGTCAGAGCAGAAAAAGACCTGGGCGGCTGAAGACCAGCTTCGCGG
+GAACTGGATGGCAGGCCTGAAGTCCGGCTGGAGTGAGTGGGAAGAGAGCGCCACGGACAG
+TATGTCGCAGGTAAAAAGTGCAGCCACGCAGACCTTTGATGGTATTGCACAGAATATGGC
+GGCGATGCTGACCGGCAGTGAGCAGAACTGGCGCAGCTTCACCCGTTCCGTGCTGTCCAT
+GATGACAGAAATTCTGCTTAAGCAGGCAATGGTGGGGATTGTCGGGAGTATCGGCAGCGC
+CATTGGCGGGGCTGTTGGTGGCGGCGCATCCGCGTCAGGCGGTACAGCCATTCAGGCCGC
+TGCGGCGAAATTCCATTTTGCAACCGGAGGATTTACGGGAACCGGCGGCAAATATGAGCC
+AGCGGGGATTGTTCACCGTGGTGAGTTTGTCTTCACGAAGGAGGCAACCAGCCGGATTGG
+CGTGGGGAATCTTTACCGGCTGATGCGCGGCTATGCCACCGGCGGTTATGTCGGTACACC
+GGGCAGCATGGCAGACAGCCGGTCGCAGGCGTCCGGGACGTTTGAGCAGAATAACCATGT
+GGTGATTAACAACGACGGCACGAACGGGCAGATAGGTCCGGCTGCTCTGAAGGCGGTGTA
+TGACATGGCCCGCAAGGGTGCCCGTGATGAAATTCAGACACAGATGCGTGATGGTGGCCT
+GTTCTCCGGAGGTGGACGATGAAGACCTTCCGCTGGAAAGTGAAACCCGGTATGGATGTG
+GCTTCGGTCCCTTCTGTAAGAAAGGTGCGCTTTGGTGATGGCTATTCTCAGCGAGCGCCT
+GCCGGGCTGAATGCCAACCTGAAAACGTACAGCGTGACGCTTTCTGTCCCCCGTGAGGAG
+GCCACGGTACTGGAGTCGTTTCTGGAAGAGCACGGGGGCTGGAAATCCTTTCTGTGGACG
+CCGCCTTATGAGTGGCGGCAGATAAAGGTGACCTGCGCAAAATGGTCGTCGCGGGTCAGT
+ATGCTGCGTGTTGAGTTCAGCGCAGAGTTTGAACAGGTGGTGAACTGATGCAGGATATCC
+GGCAGGAAACACTGAATGAATGCACCCGTGCGGAGCAGTCGGCCAGCGTGGTGCTCTGGG
+AAATCGACCTGACAGAGGTCGGTGGAGAACGTTATTTTTTCTGTAATGAGCAGAACGAAA
+AAGGTGAGCCGGTCACCTGGCAGGGGCGACAGTATCAGCCGTATCCCATTCAGGGGAGCG
+GTTTTGAACTGAATGGCAAAGGCACCAGTACGCGCCCCACGCTGACGGTTTCTAACCTGT
+ACGGTATGGTCACCGGGATGGCGGAAGATATGCAGAGTCTGGTCGGCGGAACGGTGGTCC
+GGCGTAAGGTTTACGCCCGTTTTCTGGATGCGGTGAACTTCGTCAACGGAAACAGTTACG
+CCGATCCGGAGCAGGAGGTGATCAGCCGCTGGCGCATTGAGCAGTGCAGCGAACTGAGCG
+CGGTGAGTGCCTCCTTTGTACTGTCCACGCCGACGGAAACGGATGGCGCTGTTTTTCCGG
+GACGTATCATGCTGGCCAACACCTGCACCTGGACCTATCGCGGTGACGAGTGCGGTTATA
+GCGGTCCGGCTGTCGCGGATGAATATGACCAGCCAACGTCCGATATCACGAAGGATAAAT
+GCAGCAAATGCCTGAGCGGTTGTAAGTTCCGCAATAACGTCGGCAACTTTGGCGGCTTCC
+TTTCCATTAACAAACTTTCGCAGTAAATCCCATGACACAGACAGAATCAGCGATTCTGGC
+GCACGCCCGGCGATGTGCGCCAGCGGAGTCGTGCGGCTTCGTGGTAAGCACGCCGGAGGG
+GGAAAGATATTTCCCCTGCGTGAATATCTCCGGTGAGCCGGAGGCGTATTTCCGTATGTC
+GCCGGAAGACTGGCTGCAGGCAGAAATGCAGGGTGAGATTGTGGCGCTGGTCCACAGCCA
+CCCCGGTGGTCTGCCCTGGCTGAGTGAGGCCGACCGGCGGCTGCAGGTGCAGAGTGATTT
+GCCGTGGTGGCTGGTCTGCCGGGGGACGATTCATAAGTTCCGCTGTGTGCCGCATCTCAC
+CGGGCGGCGCTTTGAGCACGGTGTGACGGACTGTTACACACTGTTCCGGGATGCTTATCA
+TCTGGCGGGGATTGAGATGCCGGACTTTCATCGTGAGGATGACTGGTGGCGTAACGGCCA
+GAATCTCTATCTGGATAATCTGGAGGCGACGGGGCTGTATCAGGTGCCGTTGTCAGCGGC
+ACAGCCGGGCGATGTGCTGCTGTGCTGTTTTGGTTCATCAGTGCCGAATCACGCCGCAAT
+TTACTGCGGCGACGGCGAGCTGCTGCACCATATTCCTGAACAACTGAGCAAACGAGAGAG
+GTACACCGACAAATGGCAGCGACGCACACACTCCCTCTGGCGTCACCGGGCATGGCGCGC
+ATCTGCCTTTACGGGGATTTACAACGATTTGGTCGCCGCATCGACCTTCGTGTGAAAACG
+GGGGCTGAAGCCATCCGGGCACTGGCCACACAGCTCCCGGCGTTTCGTCAGAAACTGAGC
+GACGGCTGGTATCAGGTACGGATTGCCGGGCGGGACGTCAGCACGTCCGGGTTAACGGCG
+CAGTTACATGAGACTCTGCCTGATGGCGCTGTAATTCATATTGTTCCCAGAGTCGCCGGG
+GCCAAGTCAGGTGGCGTATTCCAGATTGTCCTGGGGGCTGCCGCCATTGCCGGATCATTC
+TTTACCGCCGGAGCCACCCTTGCAGCATGGGGGGCAGCCATTGGGGCCGGTGGTATGACC
+GGCATCCTGTTTTCTCTCGGTGCCAGTATGGTGCTCGGTGGTGTGGCGCAGATGCTGGCA
+CCGAAAGCCAGAACTCCCCGTATACAGACAACGGATAACGGTAAGCAGAACACCTATTTC
+TCCTCACTGGATAACATGGTTGCCCAGGGCAATGTTCTGCCTGTTCTGTACGGGGAAATG
+CGCGTGGGGTCACGCGTGGTTTCTCAGGAGATCAGCACGGCAGACGAAGGGGACGGTGGT
+CAGGTTGTGGTGATTGGTCGCTGATGCAAAATGTTTTATGTGAAACCGCCTGCGGGCGGT
+TTTGTCATTTATGGAGCGTGAGGAATGGGTAAAGGAAGCAGTAAGGGGCATACCCCGCGC
+GAAGCGAAGGACAACCTGAAGTCCACGCAGTTGCTGAGTGTGATCGATGCCATCAGCGAA
+GGGCCGATTGAAGGTCCGGTGGATGGCTTAAAAAGCGTGCTGCTGAACAGTACGCCGGTG
+CTGGACACTGAGGGGAATACCAACATATCCGGTGTCACGGTGGTGTTCCGGGCTGGTGAG
+CAGGAGCAGACTCCGCCGGAGGGATTTGAATCCTCCGGCTCCGAGACGGTGCTGGGTACG
+GAAGTGAAATATGACACGCCGATCACCCGCACCATTACGTCTGCAAACATCGACCGTCTG
+CGCTTTACCTTCGGTGTACAGGCACTGGTGGAAACCACCTCAAAGGGTGACAGGAATCCG
+TCGGAAGTCCGCCTGCTGGTTCAGATACAACGTAACGGTGGCTGGGTGACGGAAAAAGAC
+ATCACCATTAAGGGCAAAACCACCTCGCAGTATCTGGCCTCGGTGGTGATGGGTAACCTG
+CCGCCGCGCCCGTTTAATATCCGGATGCGCAGGATGACGCCGGACAGCACCACAGACCAG
+CTGCAGAACAAAACGCTCTGGTCGTCATACACTGAAATCATCGATGTGAAACAGTGCTAC
+CCGAACACGGCACTGGTCGGCGTGCAGGTGGACTCGGAGCAGTTCGGCAGCCAGCAGGTG
+AGCCGTAATTATCATCTGCGCGGGCGTATTCTGCAGGTGCCGTCGAACTATAACCCGCAG
+ACGCGGCAATACAGCGGTATCTGGGACGGAACGTTTAAACCGGCATACAGCAACAACATG
+GCCTGGTGTCTGTGGGATATGCTGACCCATCCGCGCTACGGCATGGGGAAACGTCTTGGT
+GCGGCGGATGTGGATAAATGGGCGCTGTATGTCATCGGCCAGTACTGCGACCAGTCAGTG
+CCGGACGGCTTTGGCGGCACGGAGCCGCGCATCACCTGTAATGCGTACCTGACCACACAG
+CGTAAGGCGTGGGATGTGCTCAGCGATTTCTGCTCGGCGATGCGCTGTATGCCGGTATGG
+AACGGGCAGACGCTGACGTTCGTGCAGGACCGACCGTCGGATAAGACGTGGACCTATAAC
+CGCAGTAATGTGGTGATGCCGGATGATGGCGCGCCGTTCCGCTACAGCTTCAGCGCCCTG
+AAGGACCGCCATAATGCCGTTGAGGTGAACTGGATTGACCCGAACAACGGCTGGGAGACG
+GCGACAGAGCTTGTTGAAGATACGCAGGCCATTGCCCGTTACGGTCGTAATGTTACGAAG
+ATGGATGCCTTTGGCTGTACCAGCCGGGGGCAGGCACACCGCGCCGGGCTGTGGCTGATT
+AAAACAGAACTGCTGGAAACGCAGACCGTGGATTTCAGCGTCGGCGCAGAAGGGCTTCGC
+CATGTACCGGGCGATGTTATTGAAATCTGCGATGATGACTATGCCGGTATCAGCACCGGT
+GGTCGTGTGCTGGCGGTGAACAGCCAGACCCGGACGCTGACGCTCGACCGTGAAATCACG
+CTGCCATCCTCCGGTACCGCGCTGATAAGCCTGGTTGACGGAAGTGGCAATCCGGTCAGC
+GTGGAGGTTCAGTCCGTCACCGACGGCGTGAAGGTAAAAGTGAGCCGTGTTCCTGACGGT
+GTTGCTGAATACAGCGTATGGGAGCTGAAGCTGCCGACGCTGCGCCAGCGACTGTTCCGC
+TGCGTGAGTATCCGTGAGAACGACGACGGCACGTATGCCATCACCGCCGTGCAGCATGTG
+CCGGAAAAAGAGGCCATCGTGGATAACGGGGCGCACTTTGACGGCGAACAGAGTGGCACG
+GTGAATGGTGTCACGCCGCCAGCGGTGCAGCACCTGACCGCAGAAGTCACTGCAGACAGC
+GGGGAATATCAGGTGCTGGCGCGATGGGACACACCGAAGGTGGTGAAGGGCGTGAGTTTC
+CTGCTCCGTCTGACCGTAACAGCGGACGACGGCAGTGAGCGGCTGGTCAGCACGGCCCGG
+ACGACGGAAACCACATACCGCTTCACGCAACTGGCGCTGGGGAACTACAGGCTGACAGTC
+CGGGCGGTAAATGCGTGGGGGCAGCAGGGCGATCCGGCGTCGGTATCGTTCCGGATTGCC
+GCACCGGCAGCACCGTCGAGGATTGAGCTGACGCCGGGCTATTTTCAGATAACCGCCACG
+CCGCATCTTGCCGTTTATGACCCGACGGTACAGTTTGAGTTCTGGTTCTCGGAAAAGCAG
+ATTGCGGATATCAGACAGGTTGAAACCAGCACGCGTTATCTTGGTACGGCGCTGTACTGG
+ATAGCCGCCAGTATCAATATCAAACCGGGCCATGATTATTACTTTTATATCCGCAGTGTG
+AACACCGTTGGCAAATCGGCATTCGTGGAGGCCGTCGGTCGGGCGAGCGATGATGCGGAA
+GGTTACCTGGATTTTTTCAAAGGCAAGATAACCGAATCCCATCTCGGCAAGGAGCTGCTG
+GAAAAAGTCGAGCTGACGGAGGATAACGCCAGCAGACTGGAGGAGTTTTCGAAAGAGTGG
+AAGGATGCCAGTGATAAGTGGAATGCCATGTGGGCTGTCAAAATTGAGCAGACCAAAGAC
+GGCAAACATTATGTCGCGGGTATTGGCCTCAGCATGGAGGACACGGAGGAAGGCAAACTG
+AGCCAGTTTCTGGTTGCCGCCAATCGTATCGCATTTATTGACCCGGCAAACGGGAATGAA
+ACGCCGATGTTTGTGGCGCAGGGCAACCAGATATTCATGAACGACGTGTTCCTGAAGCGC
+CTGACGGCCCCCACCATTACCAGCGGCGGCAATCCTCCGGCCTTTTCCCTGACACCGGAC
+GGAAAGCTGACCGCTAAAAATGCGGATATCAGTGGCAGTGTGAATGCGAACTCCGGGACG
+CTCAGTAATGTGACGATAGCTGAAAACTGTACGATAAACGGTACGCTGAGGGCGGAAAAA
+ATCGTCGGGGACATTGTAAAGGCGGCGAGCGCGGCTTTTCCGCGCCAGCGTGAAAGCAGT
+GTGGACTGGCCGTCAGGTACCCGTACTGTCACCGTGACCGATGACCATCCTTTTGATCGC
+CAGATAGTGGTGCTTCCGCTGACGTTTCGCGGAAGTAAGCGTACTGTCAGCGGCAGGACA
+ACGTATTCGATGTGTTATCTGAAAGTACTGATGAACGGTGCGGTGATTTATGATGGCGCG
+GCGAACGAGGCGGTACAGGTGTTCTCCCGTATTGTTGACATGCCAGCGGGTCGGGGAAAC
+GTGATCCTGACGTTCACGCTTACGTCCACACGGCATTCGGCAGATATTCCGCCGTATACG
+TTTGCCAGCGATGTGCAGGTTATGGTGATTAAGAAACAGGCGCTGGGCATCAGCGTGGTC
+TGAGTGTGTTACAGAGGTTCGTCCGGGAACGGGCGTTTTATTATAAAACAGTGAGAGGTG
+AACGATGCGTAATGTGTGTATTGCCGTTGCTGTCTTTGCCGCACTTGCGGTGACAGTCAC
+TCCGGCCCGTGCGGAAGGTGGACATGGTACGTTTACGGTGGGCTATTTTCAAGTGAAACC
+GGGTACATTGCCGTCGTTGTCGGGCGGGGATACCGGTGTGAGTCATCTGAAAGGGATTAA
+CGTGAAGTACCGTTATGAGCTGACGGACAGTGTGGGGGTGATGGCTTCCCTGGGGTTCGC
+CGCGTCGAAAAAGAGCAGCACAGTGATGACCGGGGAGGATACGTTTCACTATGAGAGCCT
+GCGTGGACGTTATGTGAGCGTGATGGCCGGACCGGTTTTACAAATCAGTAAGCAGGTCAG
+TGCGTACGCCATGGCCGGAGTGGCTCACAGTCGGTGGTCCGGCAGTACAATGGATTACCG
+TAAGACGGAAATCACTCCCGGGTATATGAAAGAGACGACCACTGCCAGGGACGAAAGTGC
+AATGCGGCATACCTCAGTGGCGTGGAGTGCAGGTATACAGATTAATCCGGCAGCGTCCGT
+CGTTGTTGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCAT
+CGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCG
+GAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAG
+ACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCA
+CCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCA
+TGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGC
+ACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCT
+GTGCCATGACGGAGGATGATGCCCGGCCGGAGGTGCTGCGTCGTCTTGAACTGATGGTGG
+AAGAGGTGGCGCGTAACGCGTCCGTGGTGGCACAGAGTACGGCAGACGCGAAGAAATCAG
+CCGGCGATGCCAGTGCATCAGCTGCTCAGGTCGCGGCCCTTGTGACTGATGCAACTGACT
+CAGCACGCGCCGCCAGCACGTCCGCCGGACAGGCTGCATCGTCAGCTCAGGAAGCGTCCT
+CCGGCGCAGAAGCGGCATCAGCAAAGGCCACTGAAGCGGAAAAAAGTGCCGCAGCCGCAG
+AGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATG
+CTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAG
+AGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAATCATCAGAAA
+CGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTG
+CCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGA
+GCGCCTCTGCCGCGGCAGACGCAAAAACAGCGGCGGCGGGGAGTGCGTCAACGGCATCCA
+CGAAGGCGACAGAGGCTGCGGGAAGTGCGGTATCAGCATCGCAGAGCAAAAGTGCGGCAG
+AAGCGGCGGCAATACGTGCAAAAAATTCGGCAAAACGTGCAGAAGATATAGCTTCAGCTG
+TCGCGCTTGAGGATGCGGACACAACGAGAAAGGGGATAGTGCAGCTCAGCAGTGCAACCA
+ACAGCACGTCTGAAACGCTTGCTGCAACGCCAAAGGCGGTTAAGGTGGTAATGGATGAAA
+CGAACAGAAAAGCCCACTGGACAGTCCGGCACTGACCGGAACGCCAACAGCACCAACCGC
+GCTCAGGGGAACAAACAATACCCAGATTGCGAACACCGCTTTTGTACTGGCCGCGATTGC
+AGATGTTATCGACGCGTCACCTGACGCACTGAATACGCTGAATGAACTGGCCGCAGCGCT
+CGGGAATGATCCAGATTTTGCTACCACCATGACTAACGCGCTTGCGGGTAAACAACCGAA
+GAATGCGACACTGACGGCGCTGGCAGGGCTTTCCACGGCGAAAAATAAATTACCGTATTT
+TGCGGAAAATGATGCCGCCAGCCTGACTGAACTGACTCAGGTTGGCAGGGATATTCTGGC
+AAAAAATTCCGTTGCAGATGTTCTTGAATACCTTGGGGCCGGTGAGAATTCGGCCTTTCC
+GGCAGGTGCGCCGATCCCGTGGCCATCAGATATCGTTCCGTCTGGCTACGTCCTGATGCA
+GGGGCAGGCGTTTGACAAATCAGCCTACCCAAAACTTGCTGTCGCGTATCCATCGGGTGT
+GCTTCCTGATATGCGAGGCTGGACAATCAAGGGGAAACCCGCCAGCGGTCGTGCTGTATT
+GTCTCAGGAACAGGATGGAATTAAGTCGCACACCCACAGTGCCAGTGCATCCGGTACGGA
+TTTGGGGACGAAAACCACATCGTCGTTTGATTACGGGACGAAAACAACAGGCAGTTTCGA
+TTACGGCACCAAATCGACGAATAACACGGGGGCTCATGCTCACAGTCTGAGCGGTTCAAC
+AGGGGCCGCGGGTGCTCATGCCCACACAAGTGGTTTAAGGATGAACAGTTCTGGCTGGAG
+TCAGTATGGAACAGCAACCATTACAGGAAGTTTATCCACAGTTAAAGGAACCAGCACACA
+GGGTATTGCTTATTTATCGAAAACGGACAGTCAGGGCAGCCACAGTCACTCATTGTCCGG
+TACAGCCGTGAGTGCCGGTGCACATGCGCATACAGTTGGTATTGGTGCGCACCAGCATCC
+GGTTGTTATCGGTGCTCATGCCCATTCTTTCAGTATTGGTTCACACGGACACACCATCAC
+CGTTAACGCTGCGGGTAACGCGGAAAACACCGTCAAAAACATTGCATTTAACTATATTGT
+GAGGCTTGCATAATGGCATTCAGAATGAGTGAACAACCACGGACCATAAAAATTTATAAT
+CTGCTGGCCGGAACTAATGAATTTATTGGTGAAGGTGACGCATATATTCCGCCTCATACC
+GGTCTGCCTGCAAACAGTACCGATATTGCACCGCCAGATATTCCGGCTGGCTTTGTGGCT
+GTTTTCAACAGTGATGAGGCATCGTGGCATCTCGTTGAAGACCATCGGGGTAAAACCGTC
+TATGACGTGGCTTCCGGCGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAAT
+TTTACCTGGTTATCGCCGGGAGGGGAATATCAGAAGTGGAACGGCACAGCCTGGGTGAAG
+GATACGGAAGCAGAAAAACTGTTCCGGATCCGGGAGGCGGAAGAAACAAAAAAAAGCCTG
+ATGCAGGTAGCCAGTGAGCATATTGCGCCGCTTCAGGATGCTGCAGATCTGGAAATTGCA
+ACGAAGGAAGAAACCTCGTTGCTGGAAGCCTGGAAGAAGTATCGGGTGTTGCTGAACCGT
+GTTGATACATCAACTGCACCTGATATTGAGTGGCCTGCTGTCCCTGTTATGGAGTAATCG
+TTTTGTGATATGCCGCAGAAACGTTGTATGAAATAACGTTCTGCGGTTAGTTAGTATATT
+GTAAAGCTGAGTATTGGTTTATTTGGCGATTATTATCTTCAGGAGAATAATGGAAGTTCT
+ATGACTCAATTGTTCATAGTGTTTACATCACCGCCAATTGCTTTTAAGACTGAACGCATG
+AAATATGGTTTTTCGTCATGTTTTGAGTCTGCTGTTGATATTTCTAAAGTCGGTTTTTTT
+TCTTCGTTTTCTCTAACTATTTTCCATGAAATACATTTTTGATTATTATTTGAATCAATT
+CCAATTACCTGAAGTCTTTCATCTATAATTGGCATTGTATGTATTGGTTTATTGGAGTAG
+ATGCTTGCTTTTCTGAGCCATAGCTCTGATATCCAAATGAAGCCATAGGCATTTGTTATT
+TTGGCTCTGTCAGCTGCATAACGCCAAAAAATATATTTATCTGCTTGATCTTCAAATGTT
+GTATTGATTAAATCAATTGGATGGAATTGTTTATCATAAAAAATTAATGTTTGAATGTGA
+TAACCGTCCTTTAAAAAAGTCGTTTCTGCAAGCTTGGCTGTATAGTCAACTAACTCTTCT
+GTCGAAGTGATATTTTTAGGCTTATCTACCAGTTTTAGACGCTCTTTAATATCTTCAGGA
+ATTATTTTATTGTCATATTGTATCATGCTAAATGACAATTTGCTTATGGAGTAATCTTTT
+AATTTTAAATAAGTTATTCTCCTGGCTTCATCAAATAAAGAGTCGAATGATGTTGGCGAA
+ATCACATCGTCACCCATTGGATTGTTTATTTGTATGCCAAGAGAGTTACAGCAGTTATAC
+ATTCTGCCATAGATTATAGCTAAGGCATGTAATAATTCGTAATCTTTTAGCGTATTAGCG
+ACCCATCGTCTTTCTGATTTAATAATAGATGATTCAGTTAAATATGAAGGTAATTTCTTT
+TGTGCAAGTCTGACTAACTTTTTTATACCAATGTTTAACATACTTTCATTTGTAATAAAC
+TCAATGTCATTTTCTTCAATGTAAGATGAAATAAGAGTAGCCTTTGCCTCGCTATACATT
+TCTAAATCGCCTTGTTTTTCTATCGTATTGCGAGAATTTTTAGCCCAAGCCATTAATGGA
+TCATTTTTCCATTTTTCAATAACATTATTGTTATACCAAATGTCATATCCTATAATCTGG
+TTTTTGTTTTTTTGAATAATAAATGTTACTGTTCTTGCGGTTTGGAGGAATTGATTCAAA
+TTCAAGCGAAATAATTCAGGGTCAAAATATGTATCAATGCAGCATTTGAGCAAGTGCGAT
+AAATCTTTAAGTCTTCTTTCCCATGGTTTTTTAGTCATAAAACTCTCCATTTTGATAGGT
+TGCATGCTAGATGCTGATATATTTTAGAGGTGATAAAATTAACTGCTTAACTGTCAATGT
+AATACAAGTTGTTTGATCTTTGCAATGATTCTTATCAGAAACCATATAGTAAATTAGTTA
+CACAGGAAATTTTTAATATTATTATTATCATTCATTATGTATTAAAATTAGAGTTGTGGC
+TTGGCTCTGCTAACACGTTGCTCATAGGAGATATGGTAGAGCCGCAGACACGTCGTATGC
+AGGAACGTGCTGCGGCTGGCTGGTGAACTTCCGATAGTGCGGGTGTTGAATGATTTCCAG
+TTGCTACCGATTTTACATATTTTTTGCATGAGAGAATTTGTACCACCTCCCACCGACCAT
+CTATGACTGTACGCCACTGTCCCTAGGACTGCTATGTGCCGGAGCGGACATTACAAACGT
+CCTTCTCGGTGCATGCCACTGTTGCCAATGACCTGCCTAGGAATTGGTTAGCAAGTTACT
+ACCGGATTTTGTAAAAACAGCCCTCCTCATATAAAAAGTATTCGTTCACTTCCGATAAGC
+GTCGTAATTTTCTATCTTTCATCATATTCTAGATCCCTCTGAAAAAATCTTCCGAGTTTG
+CTAGGCACTGATACATAACTCTTTTCCAATAATTGGGGAAGTCATTCAAATCTATAATAG
+GTTTCAGATTTGCTTCAATAAATTCTGACTGTAGCTGCTGAAACGTTGCGGTTGAACTAT
+ATTTCCTTATAACTTTTACGAAAGAGTTTCTTTGAGTAATCACTTCACTCAAGTGCTTCC
+CTGCCTCCAAACGATACCTGTTAGCAATATTTAATAGCTTGAAATGATGAAGAGCTCTGT
+GTTTGTCTTCCTGCCTCCAGTTCGCCGGGCATTCAACATAAAAACTGATAGCACCCGGAG
+TTCCGGAAACGAAATTTGCATATACCCATTGCTCACGAAAAAAAATGTCCTTGTCGATAT
+AGGGATGAATCGCTTGGTGTACCTCATCTACTGCGAAAACTTGACCTTTCTCTCCCATAT
+TGCAGTCGCGGCACGATGGAACTAAATTAATAGGCATCACCGAAAATTCAGGATAATGTG
+CAATAGGAAGAAAATGATCTATATTTTTTGTCTGTCCTATATCACCACAAAATGGACATT
+TTTCACCTGATGAAACAAGCATGTCATCGTAATATGTTCTAGCGGGTTTGTTTTTATCTC
+GGAGATTATTTTCATAAAGCTTTTCTAATTTAACCTTTGTCAGGTTACCAACTACTAAGG
+TTGTAGGCTCAAGAGGGTGTGTCCTGTCGTAGGTAAATAACTGACCTGTCGAGCTTAATA
+TTCTATATTGTTGTTCTTTCTGCAAAAAAGTGGGGAAGTGAGTAATGAAATTATTTCTAA
+CATTTATCTGCATCATACCTTCCGAGCATTTATTAAGCATTTCGCTATAAGTTCTCGCTG
+GAAGAGGTAGTTTTTTCATTGTACTTTACCTTCATCTCTGTTCATTATCATCGCTTTTAA
+AACGGTTCGACCTTCTAATCCTATCTGACCATTATAATTTTTTAGAATGGTTTCATAAGA
+AAGCTCTGAATCAACGGACTGCGATAATAAGTGGTGGTATCCAGAATTTGTCACTTCAAG
+TAAAAACACCTCACGAGTTAAAACACCTAAGTTCTCACCGAATGTCTCAATATCCGGACG
+GATAATATTTATTGCTTCTCTTGACCGTAGGACTTTCCACATGCAGGATTTTGGAACCTC
+TTGCAGTACTACTGGGGAATGAGTTGCAATTATTGCTACACCATTGCGTGCATCGAGTAA
+GTCGCTTAATGTTCGTAAAAAAGCAGAGAGCAAAGGTGGATGCAGATGAACCTCTGGTTC
+ATCGAATAAAACTAATGACTTTTCGCCAACGACATCTACTAATCTTGTGATAGTAAATAA
+AACAATTGCATGTCCAGAGCTCATTCGAAGCAGATATTTCTGGATATTGTCATAAAACAA
+TTTAGTGAATTTATCATCGTCCACTTGAATCTGTGGTTCATTACGTCTTAACTCTTCATA
+TTTAGAAATGAGGCTGATGAGTTCCATATTTGAAAAGTTTTCATCACTACTTAGTTTTTT
+GATAGCTTCAAGCCAGAGTTGTCTTTTTCTATCTACTCTCATACAACCAATAAATGCTGA
+AATGAATTCTAAGCGGAGATCGCCTAGTGATTTTAAACTATTGCTGGCAGCATTCTTGAG
+TCCAATATAAAAGTATTGTGTACCTTTTGCTGGGTCAGGTTGTTCTTTAGGAGGAGTAAA
+AGGATCAAATGCACTAAACGAAACTGAAACAAGCGATCGAAAATATCCCTTTGGGATTCT
+TGACTCGATAAGTCTATTATTTTCAGAGAAAAAATATTCATTGTTTTCTGGGTTGGTGAT
+TGCACCAATCATTCCATTCAAAATTGTTGTTTTACCACACCCATTCCGCCCGATAAAAGC
+ATGAATGTTCGTGCTGGGCATAGAATTAACCGTCACCTCAAAAGGTATAGTTAAATCACT
+GAATCCGGGAGCACTTTTTCTATTAAATGAAAAGTGGAAATCTGACAATTCTGGCAAACC
+ATTTAACACACGTGCGAACTGTCCATGAATTTCTGAAAGAGTTACCCCTCTAAGTAATGA
+GGTGTTAAGGACGCTTTCATTTTCAATGTCGGCTAATCGATTTGGCCATACTACTAAATC
+CTGAATAGCTTTAAGAAGGTTATGTTTAAAACCATCGCTTAATTTGCTGAGATTAACATA
+GTAGTCAATGCTTTCACCTAAGGAAAAAAACATTTCAGGGAGTTGACTGAATTTTTTATC
+TATTAATGAATAAGTGCTTACTTCTTCTTTTTGACCTACAAAACCAATTTTAACATTTCC
+GATATCGCATTTTTCACCATGCTCATCAAAGACAGTAAGATAAAACATTGTAACAAAGGA
+ATAGTCATTCCAACCATCTGCTCGTAGGAATGCCTTATTTTTTTCTACTGCAGGAATATA
+CCCGCCTCTTTCAATAACACTAAACTCCAACATATAGTAACCCTTAATTTTATTAAAATA
+ACCGCAATTTATTTGGCGGCAACACAGGATCTCTCTTTTAAGTTACTCTCTATTACATAC
+GTTTTCCATCTAAAAATTAGTAGTATTGAACTTAACGGGGCATCGTATTGTAGTTTTCCA
+TATTTAGCTTTCTGCTTCCTTTTGGATAACCCACTGTTATTCATGTTGCATGGTGCACTG
+TTTATACCAACGATATAGTCTATTAATGCATATATAGTATCGCCGAACGATTAGCTCTTC
+AGGCTTCTGAAGAAGCGTTTCAAGTACTAATAAGCCGATAGATAGCCACGGACTTCGTAG
+CCATTTTTCATAAGTGTTAACTTCCGCTCCTCGCTCATAACAGACATTCACTACAGTTAT
+GGCGGAAAGGTATGCATGCTGGGTGTGGGGAAGTCGTGAAAGAAAAGAAGTCAGCTGCGT
+CGTTTGACATCACTGCTATCTTCTTACTGGTTATGCAGGTCGTAGTGGGTGGCACACAAA
+GCTTTGCACTGGATTGCGAGGCTTTGTGCTTCTCTGGAGTGCGACAGGTTTGATGACAAA
+AAATTAGCGCAAGAAGACAAAAATCACCTTGCGCTAATGCTCTGTTACAGGTCACTAATA
+CCATCTAAGTAGTTGATTCATAGTGACTGCATATGTTGTGTTTTACAGTATTATGTAGTC
+TGTTTTTTATGCAAAATCTAATTTAATATATTGATATTTATATCATTTTACGTTTCTCGT
+TCAGCTTTTTTATACTAAGTTGGCATTATAAAAAAGCATTGCTTATCAATTTGTTGCAAC
+GAACAGGTCACTATCAGTCAAAATAAAATCATTATTTGATTTCAATTTTGTCCCACTCCC
+TGCCTCTGTCATCACGATACTGTGATGCCATGGTGTCCGACTTATGCCCGAGAAGATGTT
+GAGCAAACTTATCGCTTATCTGCTTCTCATAGAGTCTTGCAGACAAACTGCGCAACTCGT
+GAAAGGTAGGCGGATCCCCTTCGAAGGAAAGACCTGATGCTTTTCGTGCGCGCATAAAAT
+ACCTTGATACTGTGCCGGATGAAAGCGGTTCGCGACGAGTAGATGCAATTATGGTTTCTC
+CGCCAAGAATCTCTTTGCATTTATCAAGTGTTTCCTTCATTGATATTCCGAGAGCATCAA
+TATGCAATGCTGTTGGGATGGCAATTTTTACGCCTGTTTTGCTTTGCTCGACATAAAGAT
+ATCCATCTACGATATCAGACCACTTCATTTCGCATAAATCACCAACTCGTTGCCCGGTAA
+CAACAGCCAGTTCCATTGCAAGTCTGAGCCAACATGGTGATGATTCTGCTGCTTGATAAA
+TTTTCAGGTATTCGTCAGCCGTAAGTCTTGATCTCCTTACCTCTGATTTTGCTGCGCGAG
+TGGCAGCGACATGGTTTGTTGTTATATGGCCTTCAGCTATTGCCTCTCGGAATGCATCGC
+TCAGTGTTGATCTGATTAACTTGGCTGACGCCGCCTTGCCCTCGTCTATGTATCCATTGA
+GCATTGCCGCAATTTCTTTTGTGGTGATGTCTTCAAGTGGAGCATCAGGCAGACCCCTCC
+TTATTGCTTTAATTTTGCTCATGTAATTTATGAGTGTCTTCTGCTTGATTCCTCTGCTGG
+CCAGGATTTTTTCGTAGCGATCAAGCCATGAATGTAACGTAACGGAATTATCACTGTTGA
+TTCTCGCTGTCAGAGGCTTGTGTTTGTGTCCTGAAAATAACTCAATGTTGGCCTGTATAG
+CTTCAGTGATTGCGATTCGCCTGTCTCTGCCTAATCCAAACTCTTTACCCGTCCTTGGGT
+CCCTGTAGCAGTAATATCCATTGTTTCTTATATAAAGGTTAGGGGGTAAATCCCGGCGCT
+CATGACTTCGCCTTCTTCCCATTTCTGATCCTCTTCAAAAGGCCACCTGTTACTGGTCGA
+TTTAAGTCAACCTTTACCGCTGATTCGTGGAACAGATACTCTCTTCCATCCTTAACCGGA
+GGTGGGAATATCCTGCATTCCCGAACCCATCGACGAACTGTTTCAAGGCTTCTTGGACGT
+CGCTGGCGTGCGTTCCACTCCTGAAGTGTCAAGTACATCGCAAAGTCTCCGCAATTACAC
+GCAAGAAAAAACCGCCATCAGGCGGCTTGGTGTTCTTTCAGTTCTTCAATTCGAATATTG
+GTTACGTCTGCATGTGCTATCTGCGCCCATATCATCCAGTGGTCGTAGCAGTCGTTGATG
+TTCTCCGCTTCGATAACTCTGTTGAATGGCTCTCCATTCCATTCTCCTGTGACTCGGAAG
+TGCATTTATCATCTCCATAAAACAAAACCCGCCGTAGCGAGTTCAGATAAAATAAATCCC
+CGCGAGTGCGAGGATTGTTATGTAATATTGGGTTTAATCATCTATATGTTTTGTACAGAG
+AGGGCAAGTATCGTTTCCACCGTACTCGTGATAATAATTTTGCACGGTATCAGTCATTTC
+TCGCACATTGCAGAATGGGGATTTGTCTTCATTAGACTTATAAACCTTCATGGAATATTT
+GTATGCCGACTCTATATCTATACCTTCATCTACATAAACACCTTCGTGATGTCTGCATGG
+AGACAAGACACCGGATCTGCACAACATTGATAACGCCCAATCTTTTTGCTCAGACTCTAA
+CTCATTGATACTCATTTATAAACTCCTTGCAATGTATGTCGTTTCAGCTAAACGGTATCA
+GCAATGTTTATGTAAAGAAACAGTAAGATAATACTCAACCCGATGTTTGAGTACGGTCAT
+CATCTGACACTACAGACTCTGGCATCGCTGTGAAGACGACGCGAAATTCAGCATTTTCAC
+AAGCGTTATCTTTTACAAAACCGATCTCACTCTCCTTTGATGCGAATGCCAGCGTCAGAC
+ATCATATGCAGATACTCACCTGCATCCTGAACCCATTGACCTCCAACCCCGTAATAGCGA
+TGCGTAATGATGTCGATAGTTACTAACGGGTCTTGTTCGATTAACTGCCGCAGAAACTCT
+TCCAGGTCACCAGTGCAGTGCTTGATAACAGGAGTCTTCCCAGGATGGCGAACAACAAGA
+AACTGGTTTCCGTCTTCACGGACTTCGTTGCTTTCCAGTTTAGCAATACGCTTACTCCCA
+TCCGAGATAACACCTTCGTAATACTCACGCTGCTCGTTGAGTTTTGATTTTGCTGTTTCA
+AGCTCAACACGCAGTTTCCCTACTGTTAGCGCAATATCCTCGTTCTCCTGGTCGCGGCGT
+TTGATGTATTGCTGGTTTCTTTCCCGTTCATCCAGCAGTTCCAGCACAATCGATGGTGTT
+ACCAATTCATGGAAAAGGTCTGCGTCAAATCCCCAGTCGTCATGCATTGCCTGCTCTGCC
+GCTTCACGCAGTGCCTGAGAGTTAATTTCGCTCACTTCGAACCTCTCTGTTTACTGATAA
+GTTCCAGATCCTCCTGGCAACTTGCACAAGTCCGACAACCCTGAACGACCAGGCGTCTTC
+GTTCATCTATCGGATCGCCACACTCACAACAATGAGTGGCAGATATAGCCTGGTGGTTCA
+GGCGGCGCATTTTTATTGCTGTGTTGCGCTGTAATTCTTCTATTTCTGATGCTGAATCAA
+TGATGTCTGCCATCTTTCATTAATCCCTGAACTGTTGGTTAATACGCTTGAGGGTGAATG
+CGAATAATAAAAAAGGAGCCTGTAGCTCCCTGATGATTTTGCTTTTCATGTTCATCGTTC
+CTTAAAGACGCCGTTTAACATGCCGATTGCCAGGCTTAAATGAGTCGGTGTGAATCCCAT
+CAGCGTTACCGTTTCGCGGTGCTTCTTCAGTACGCTACGGCAAATGTCATCGACGTTTTT
+ATCCGGAAACTGCTGTCTGGCTTTTTTTGATTTCAGAATTAGCCTGACGGGCAATGCTGC
+GAAGGGCGTTTTCCTGCTGAGGTGTCATTGAACAAGTCCCATGTCGGCAAGCATAAGCAC
+ACAGAATATGAAGCCCGCTGCCAGAAAAATGCATTCCGTGGTTGTCATACCTGGTTTCTC
+TCATCTGCTTCTGCTTTCGCCACCATCATTTCCAGCTTTTGTGAAAGGGATGCGGCTAAC
+GTATGAAATTCTTCGTCTGTTTCTACTGGTATTGGCACAAACCTGATTCCAATTTGAGCA
+AGGCTATGTGCCATCTCGATACTCGTTCTTAACTCAACAGAAGATGCTTTGTGCATACAG
+CCCCTCGTTTATTATTTATCTCCTCAGCCAGCCGCTGTGCTTTCAGTGGATTTCGGATAA
+CAGAAAGGCCGGGAAATACCCAGCCTCGCTTTGTAACGGAGTAGACGAAAGTGATTGCGC
+CTACCCGGATATTATCGTGAGGATGCGTCATCGCCATTGCTCCCCAAATACAAAACCAAT
+TTCAGCCAGTGCCTCGTCCATTTTTTCGATGAACTCCGGCACGATCTCGTCAAAACTCGC
+CATGTACTTTTCATCCCGCTCAATCACGACATAATGCAGGCCTTCACGCTTCATACGCGG
+GTCATAGTTGGCAAAGTACCAGGCATTTTTTCGCGTCACCCACATGCTGTACTGCACCTG
+GGCCATGTAAGCTGACTTTATGGCCTCGAAACCACCGAGCCGGAACTTCATGAAATCCCG
+GGAGGTAAACGGGCATTTCAGTTCAAGGCCGTTGCCGTCACTGCATAAACCATCGGGAGA
+GCAGGCGGTACGCATACTTTCGTCGCGATAGATGATCGGGGATTCAGTAACATTCACGCC
+GGAAGTGAATTCAAACAGGGTTCTGGCGTCGTTCTCGTACTGTTTTCCCCAGGCCAGTGC
+TTTAGCGTTAACTTCCGGAGCCACACCGGTGCAAACCTCAGCAAGCAGGGTGTGGAAGTA
+GGACATTTTCATGTCAGGCCACTTCTTTCCGGAGCGGGGTTTTGCTATCACGTTGTGAAC
+TTCTGAAGCGGTGATGACGCCGAGCCGTAATTTGTGCCACGCATCATCCCCCTGTTCGAC
+AGCTCTCACATCGATCCCGGTACGCTGCAGGATAATGTCCGGTGTCATGCTGCCACCTTC
+TGCTCTGCGGCTTTCTGTTTCAGGAATCCAAGAGCTTTTACTGCTTCGGCCTGTGTCAGT
+TCTGACGATGCACGAATGTCGCGGCGAAATATCTGGGAACAGAGCGGCAATAAGTCGTCA
+TCCCATGTTTTATCCAGGGCGATCAGCAGAGTGTTAATCTCCTGCATGGTTTCATCGTTA
+ACCGGAGTGATGTCGCGTTCCGGCTGACGTTCTGCAGTGTATGCAGTATTTTCGACAATG
+CGCTCGGCTTCATCCTTGTCATAGATACCAGCAAATCCGAAGGCCAGACGGGCACACTGA
+ATCATGGCTTTATGACGTAACATCCGTTTGGGATGCGACTGCCACGGCCCCGTGATTTCT
+CTGCCTTCGCGAGTTTTGAATGGTTCGCGGCGGCATTCATCCATCCATTCGGTAACGCAG
+ATCGGATGATTACGGTCCTTGCGGTAAATCCGGCATGTACAGGATTCATTGTCCTGCTCA
+AAGTCCATGCCATCAAACTGCTGGTTTTCATTGATGATGCGGGACCAGCCATCAACGCCC
+ACCACCGGAACGATGCCATTCTGCTTATCAGGAAAGGCGTAAATTTCTTTCGTCCACGGA
+TTAAGGCCGTACTGGTTGGCAACGATCAGTAATGCGATGAACTGCGCATCGCTGGCATCA
+CCTTTAAATGCCGTCTGGCGAAGAGTGGTGATCAGTTCCTGTGGGTCGACAGAATCCATG
+CCGACACGTTCAGCCAGCTTCCCAGCCAGCGTTGCGAGTGCAGTACTCATTCGTTTTATA
+CCTCTGAATCAATATCAACCTGGTGGTGAGCAATGGTTTCAACCATGTACCGGATGTGTT
+CTGCCATGCGCTCCTGAAACTCAACATCGTCATCAAACGCACGGGTAATGGATTTTTTGC
+TGGCCCCGTGGCGTTGCAAATGATCGATGCATAGCGATTCAAACAGGTGCTGGGGCAGGC
+CTTTTTCCATGTCGTCTGCCAGTTCTGCCTCTTTCTCTTCACGGGCGAGCTGCTGGTAGT
+GACGCGCCCAGCTCTGAGCCTCAAGACGATCCTGAATGTAATAAGCGTTCATGGCTGAAC
+TCCTGAAATAGCTGTGAAAATATCGCCCGCGAAATGCCGGGCTGATTAGGAAAACAGGAA
+AGGGGGTTAGTGAATGCTTTTGCTTGATCTCAGTTTCAGTATTAATATCCATTTTTTATA
+AGCGTCGACGGCTTCACGAAACATCTTTTCATCGCCAATAAAAGTGGCGATAGTGAATTT
+AGTCTGGATAGCCATAAGTGTTTGATCCATTCTTTGGGACTCCTGGCTGATTAAGTATGT
+CGATAAGGCGTTTCCATCCGTCACGTAATTTACGGGTGATTCGTTCAAGTAAAGATTCGG
+AAGGGCAGCCAGCAACAGGCCACCCTGCAATGGCATATTGCATGGTGTGCTCCTTATTTA
+TACATAACGAAAAACGCCTCGAGTGAAGCGTTATTGGTATGCGGTAAAACCGCACTCAGG
+CGGCCTTGATAGTCATATCATCTGAATCAAATATTCCTGATGTATCGATATCGGTAATTC
+TTATTCCTTCGCTACCATCCATTGGAGGCCATCCTTCCTGACCATTTCCATCATTCCAGT
+CGAACTCACACACAACACCATATGCATTTAAGTCGCTTGAAATTGCTATAAGCAGAGCAT
+GTTGCGCCAGCATGATTAATACAGCATTTAATACAGAGCCGTGTTTATTGAGTCGGTATT
+CAGAGTCTGACCAGAAATTATTAATCTGGTGAAGTTTTTCCTCTGTCATTACGTCATGGT
+CGATTTCAATTTCTATTGATGCTTTCCAGTCGTAATCAATGATGTATTTTTTGATGTTTG
+ACATCTGTTCATATCCTCACAGATAAAAAATCGCCCTCACACTGGAGGGCAAAGAAGATT
+TCCAATAATCAGAACAAGTCGGCTCCTGTTTAGTTACGAGCGACATTGCTCCGTGTATTC
+ACTCGTTGGAATGAATACACAGTGCAGTGTTTATTCTGTTATTTATGCCAAAAATAAAGG
+CCACTATCAGGCAGCTTTGTTGTTCTGTTTACCAAGTTCTCTGGCAATCATTGCCGTCGT
+TCGTATTGCCCATTTATCGACATATTTCCCATCTTCCATTACAGGAAACATTTCTTCAGG
+CTTAACCATGCATTCCGATTGCAGCTTGCATCCATTGCATCGCTTGAATTGTCCACACCA
+TTGATTTTTATCAATAGTCGTAGTCATACGGATAGTCCTGGTATTGTTCCATCACATCCT
+GAGGATGCTCTTCGAACTCTTCAAATTCTTCTTCCATATATCACCTTAAATAGTGGATTG
+CGGTAGTAAAGATTGTGCCTGTCTTTTAACCACATCAGGCTCGGTGGTTCTCGTGTACCC
+CTACAGCGAGAAATCGGATAAACTATTACAACCCCTACAGTTTGATGAGTATAGAAATGG
+ATCCACTCGTTATTCTCGGACGAGTGTTCAGTAATGAACCTCTGGAGAGAACCATGTATA
+TGATCGTTATCTGGGTTGGACTTCTGCTTTTAAGCCCAGATAACTGGCCTGAATATGTTA
+ATGAGAGAATCGGTATTCCTCATGTGTGGCATGTTTTCGTCTTTGCTCTTGCATTTTCGC
+TAGCAATTAATGTGCATCGATTATCAGCTATTGCCAGCGCCAGATATAAGCGATTTAAGC
+TAAGAAAACGCATTAAGATGCAAAACGATAAAGTGCGATCAGTAATTCAAAACCTTACAG
+AAGAGCAATCTATGGTTTTGTGCGCAGCCCTTAATGAAGGCAGGAAGTATGTGGTTACAT
+CAAAACAATTCCCATACATTAGTGAGTTGATTGAGCTTGGTGTGTTGAACAAAACTTTTT
+CCCGATGGAATGGAAAGCATATATTATTCCCTATTGAGGATATTTACTGGACTGAATTAG
+TTGCCAGCTATGATCCATATAATATTGAGATAAAGCCAAGGCCAATATCTAAGTAACTAG
+ATAAGAGGAATCGATTTTCCCTTAATTTTCTGGCGTCCACTGCATGTTATGCCGCGTTCG
+CCAGGCTTGCTGTACCATGTGCGCTGATTCTTGCGCTCAATACGTTGCAGGTTGCTTTCA
+ATCTGTTTGTGGTATTCAGCCAGCACTGTAAGGTCTATCGGATTTAGTGCGCTTTCTACT
+CGTGATTTCGGTTTGCGATTCAGCGAGAGAATAGGGCGGTTAACTGGTTTTGCGCTTACC
+CCAACCAACAGGGGATTTGCTGCTTTCCATTGAGCCTGTTTCTCTGCGCGACGTTCGCGG
+CGGCGTGTTTGTGCATCCATCTGGATTCTCCTGTCAGTTAGCTTTGGTGGTGTGTGGCAG
+TTGTAGTCCTGAACGAAAACCCCCCGCGATTGGCACATTGGCAGCTAATCCGGAATCGCA
+CTTACGGCCAATGCTTCGTTTCGTATCACACACCCCAAAGCCTTCTGCTTTGAATGCTGC
+CCTTCTTCAGGGCTTAATTTTTAAGAGCGTCACCTTCATGGTGGTCAGTGCGTCCTGCTG
+ATGTGCTCAGTATCACCGCCAGTGGTATTTATGTCAACACCGCCAGAGATAATTTATCAC
+CGCAGATGGTTATCTGTATGTTTTTTATATGAATTTATTTTTTGCAGGGGGGCATTGTTT
+GGTAGGTGAGAGATCTGAATTGCTATGTTTAGTGAGTTGTATCTATTTATTTTTCAATAA
+ATACAATTGGTTATGTGTTTTGGGGGCGATCGTGAGGCAAAGAAAACCCGGCGCTGAGGC
+CGGGTTATTCTTGTTCTCTGGTCAAATTATATAGTTGGAAAACAAGGATGCATATATGAA
+TGAACGATGCAGAGGCAATGCCGATGGCGATAGTGGGTATCATGTAGCCGCTTATGCTGG
+AAAGAAGCAATAACCCGCAGAAAAACAAAGCTCCAAGCTCAACAAAACTAAGGGCATAGA
+CAATAACTACCGATGTCATATACCCATACTCTCTAATCTTGGCCAGTCGGCGCGTTCTGC
+TTCCGATTAGAAACGTCAAGGCAGCAATCAGGATTGCAATCATGGTTCCTGCATATGATG
+ACAATGTCGCCCCAAGACCATCTCTATGAGCTGAAAAAGAAACACCAGGAATGTAGTGGC
+GGAAAAGGAGATAGCAAATGCTTACGATAACGTAAGGAATTATTACTATGTAAACACCAG
+GCATGATTCTGTTCCGCATAATTACTCCTGATAATTAATCCTTAACTTTGCCCACCTGCC
+TTTTAAAACATTCCAGTATATCACTTTTCATTCTTGCGTAGCAATATGCCATCTCTTCAG
+CTATCTCAGCATTGGTGACCTTGTTCAGAGGCGCTGAGAGATGGCCTTTTTCTGATAGAT
+AATGTTCTGTTAAAATATCTCCGGCCTCATCTTTTGCCCGCAGGCTAATGTCTGAAAATT
+GAGGTGACGGGTTAAAAATAATATCCTTGGCAACCTTTTTTATATCCCTTTTAAATTTTG
+GCTTAATGACTATATCCAATGAGTCAAAAAGCTCCCCTTCAATATCTGTTGCCCCTAAGA
+CCTTTAATATATCGCCAAATACAGGTAGCTTGGCTTCTACCTTCACCGTTGTTCGGCCGA
+TGAAATGCATATGCATAACATCGTCTTTGGTGGTTCCCCTCATCAGTGGCTCTATCTGAA
+CGCGCTCTCCACTGCTTAATGACATTCCTTTCCCGATTAAAAAATCTGTCAGATCGGATG
+TGGTCGGCCCGAAAACAGTTCTGGCAAAACCAATGGTGTCGCCTTCAACAAACAAAAAAG
+ATGGGAATCCCAATGATTCGTCATCTGCGAGGCTGTTCTTAATATCTTCAACTGAAGCTT
+TAGAGCGATTTATCTTCTGAACCAGACTCTTGTCATTTGTTTTGGTAAAGAGAAAAGTTT
+TTCCATCGATTTTATGAATATACAAATAATTGGAGCCAACCTGCAGGTGATGATTATCAG
+CCAGCAGAGAATTAAGGAAAACAGACAGGTTTATTGAGCGCTTATCTTTCCCTTTATTTT
+TGCTGCGGTAAGTCGCATAAAAACCATTCTTCATAATTCAATCCATTTACTATGTTATGT
+TCTGAGGGGAGTGAAAATTCCCCTAATTCGATGAAGATTCTTGCTCAATTGTTATCAGCT
+ATGCGCCGACCAGAACACCTTGCCGATCAGCCAAACGTCTCTTCAGGCCACTGACTAGCG
+ATAACTTTCCCCACAACGGAACAACTCTCATTGCATGGGATCATTGGGTACTGTGGGTTT
+AGTGGTTGTAAAAACACCTGACCGCTATCCCTGATCAGTTTCTTGAAGGTAAACTCATCA
+CCCCCAAGTCTGGCTATGCAGAAATCACCTGGCTCAACAGCCTGCTCAGGGTCAACGAGA
+ATTAACATTCCGTCAGGAAAGCTTGGCTTGGAGCCTGTTGGTGCGGTCATGGAATTACCT
+TCAACCTCAAGCCAGAATGCAGAATCACTGGCTTTTTTGGTTGTGCTTACCCATCTCTCC
+GCATCACCTTTGGTAAAGGTTCTAAGCTTAGGTGAGAACATCCCTGCCTGAACATGAGAA
+AAAACAGGGTACTCATACTCACTTCTAAGTGACGGCTGCATACTAACCGCTTCATACATC
+TCGTAGATTTCTCTGGCGATTGAAGGGCTAAATTCTTCAACGCTAACTTTGAGAATTTTT
+GTAAGCAATGCGGCGTTATAAGCATTTAATGCATTGATGCCATTAAATAAAGCACCAACG
+CCTGACTGCCCCATCCCCATCTTGTCTGCGACAGATTCCTGGGATAAGCCAAGTTCATTT
+TTCTTTTTTTCATAAATTGCTTTAAGGCGACGTGCGTCCTCAAGCTGCTCTTGTGTTAAT
+GGTTTCTTTTTTGTGCTCATACGTTAAATCTATCACCGCAAGGGATAAATATCTAACACC
+GTGCGTGTTGACTATTTTACCTCTGGCGGTGATAATGGTTGCATGTACTAAGGAGGTTGT
+ATGGAACAACGCATAACCCTGAAAGATTATGCAATGCGCTTTGGGCAAACCAAGACAGCT
+AAAGATCTCGGCGTATATCAAAGCGCGATCAACAAGGCCATTCATGCAGGCCGAAAGATT
+TTTTTAACTATAAACGCTGATGGAAGCGTTTATGCGGAAGAGGTAAAGCCCTTCCCGAGT
+AACAAAAAAACAACAGCATAAATAACCCCGCTCTTACACATTCCAGCCCTGAAAAAGGGC
+ATCAAATTAAACCACACCTATGGTGTATGCATTTATTTGCATACATTCAATCAATTGTTA
+TCTAAGGAAATACTTACATATGGTTCGTGCAAACAAACGCAACGAGGCTCTACGAATCGA
+GAGTGCGTTGCTTAACAAAATCGCAATGCTTGGAACTGAGAAGACAGCGGAAGCTGTGGG
+CGTTGATAAGTCGCAGATCAGCAGGTGGAAGAGGGACTGGATTCCAAAGTTCTCAATGCT
+GCTTGCTGTTCTTGAATGGGGGGTCGTTGACGACGACATGGCTCGATTGGCGCGACAAGT
+TGCTGCGATTCTCACCAATAAAAAACGCCCGGCGGCAACCGAGCGTTCTGAACAAATCCA
+GATGGAGTTCTGAGGTCATTACTGGATCTATCAACAGGAGTCATTATGACAAATACAGCA
+AAAATACTCAACTTCGGCAGAGGTAACTTTGCCGGACAGGAGCGTAATGTGGCAGATCTC
+GATGATGGTTACGCCAGACTATCAAATATGCTGCTTGAGGCTTATTCGGGCGCAGATCTG
+ACCAAGCGACAGTTTAAAGTGCTGCTTGCCATTCTGCGTAAAACCTATGGGTGGAATAAA
+CCAATGGACAGAATCACCGATTCTCAACTTAGCGAGATTACAAAGTTACCTGTCAAACGG
+TGCAATGAAGCCAAGTTAGAACTCGTCAGAATGAATATTATCAAGCAGCAAGGCGGCATG
+TTTGGACCAAATAAAAACATCTCAGAATGGTGCATCCCTCAAAACGAGGGAAAATCCCCT
+AAAACGAGGGATAAAACATCCCTCAAATTGGGGGATTGCTATCCCTCAAAACAGGGGGAC
+ACAAAAGACACTATTACAAAAGAAAAAAGAAAAGATTATTCGTCAGAGAATTCTGGCGAA
+TCCTCTGACCAGCCAGAAAACGACCTTTCTGTGGTGAAACCGGATGCTGCAATTCAGAGC
+GGCAGCAAGTGGGGGACAGCAGAAGACCTGACCGCCGCAGAGTGGATGTTTGACATGGTG
+AAGACTATCGCACCATCAGCCAGAAAACCGAATTTTGCTGGGTGGGCTAACGATATCCGC
+CTGATGCGTGAACGTGACGGACGTAACCACCGCGACATGTGTGTGCTGTTCCGCTGGGCA
+TGCCAGGACAACTTCTGGTCCGGTAACGTGCTGAGCCCGGCCAAACTCCGCGATAAGTGG
+ACCCAACTCGAAATCAACCGTAACAAGCAACAGGCAGGCGTGACAGCCAGCAAACCAAAA
+CTCGACCTGACAAACACAGACTGGATTTACGGGGTGGATCTATGAAAAACATCGCCGCAC
+AGATGGTTAACTTTGACCGTGAGCAGATGCGTCGGATCGCCAACAACATGCCGGAACAGT
+ACGACGAAAAGCCGCAGGTACAGCAGGTAGCGCAGATCATCAACGGTGTGTTCAGCCAGT
+TACTGGCAACTTTCCCGGCGAGCCTGGCTAACCGTGACCAGAACGAAGTGAACGAAATCC
+GTCGCCAGTGGGTTCTGGCTTTTCGGGAAAACGGGATCACCACGATGGAACAGGTTAACG
+CAGGAATGCGCGTAGCCCGTCGGCAGAATCGACCATTTCTGCCATCACCCGGGCAGTTTG
+TTGCATGGTGCCGGGAAGAAGCATCCGTTACCGCCGGACTGCCAAACGTCAGCGAGCTGG
+TTGATATGGTTTACGAGTATTGCCGGAAGCGAGGCCTGTATCCGGATGCGGAGTCTTATC
+CGTGGAAATCAAACGCGCACTACTGGCTGGTTACCAACCTGTATCAGAACATGCGGGCCA
+ATGCGCTTACTGATGCGGAATTACGCCGTAAGGCCGCAGATGAGCTTGTCCATATGACTG
+CGAGAATTAACCGTGGTGAGGCGATCCCTGAACCAGTAAAACAACTTCCTGTCATGGGCG
+GTAGACCTCTAAATCGTGCACAGGCTCTGGCGAAGATCGCAGAAATCAAAGCTAAGTTCG
+GACTGAAAGGAGCAAGTGTATGACGGGCAAAGAGGCAATTATTCATTACCTGGGGACGCA
+TAATAGCTTCTGTGCGCCGGACGTTGCCGCGCTAACAGGCGCAACAGTAACCAGCATAAA
+TCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCG
+AACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCT
+GGTTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGGCGGTATATGGAGT
+TAAAAGATGACCATCTACATTACTGAGCTAATAACAGGCCTGCTGGTAATCGCAGGCCTT
+TTTATTTGGGGGAGAGGGAAGTCATGAAAAAACTAACCTTTGAAATTCGATCTCCAGCAC
+ATCAGCAAAACGCTATTCACGCAGTACAGCAAATCCTTCCAGACCCAACCAAACCAATCG
+TAGTAACCATTCAGGAACGCAACCGCAGCTTAGACCAAAACAGGAAGCTATGGGCCTGCT
+TAGGTGACGTCTCTCGTCAGGTTGAATGGCATGGTCGCTGGCTGGATGCAGAAAGCTGGA
+AGTGTGTGTTTACCGCAGCATTAAAGCAGCAGGATGTTGTTCCTAACCTTGCCGGGAATG
+GCTTTGTGGTAATAGGCCAGTCAACCAGCAGGATGCGTGTAGGCGAATTTGCGGAGCTAT
+TAGAGCTTATACAGGCATTCGGTACAGAGCGTGGCGTTAAGTGGTCAGACGAAGCGAGAC
+TGGCTCTGGAGTGGAAAGCGAGATGGGGAGACAGGGCTGCATGATAAATGTCGTTAGTTT
+CTCCGGTGGCAGGACGTCAGCATATTTGCTCTGGCTAATGGAGCAAAAGCGACGGGCAGG
+TAAAGACGTGCATTACGTTTTCATGGATACAGGTTGTGAACATCCAATGACATATCGGTT
+TGTCAGGGAAGTTGTGAAGTTCTGGGATATACCGCTCACCGTATTGCAGGTTGATATCAA
+CCCGGAGCTTGGACAGCCAAATGGTTATACGGTATGGGAACCAAAGGATATTCAGACGCG
+AATGCCTGTTCTGAAGCCATTTATCGATATGGTAAAGAAATATGGCACTCCATACGTCGG
+CGGCGCGTTCTGCACTGACAGATTAAAACTCGTTCCCTTCACCAAATACTGTGATGACCA
+TTTCGGGCGAGGGAATTACACCACGTGGATTGGCATCAGAGCTGATGAACCGAAGCGGCT
+AAAGCCAAAGCCTGGAATCAGATATCTTGCTGAACTGTCAGACTTTGAGAAGGAAGATAT
+CCTCGCATGGTGGAAGCAACAACCATTCGATTTGCAAATACCGGAACATCTCGGTAACTG
+CATATTCTGCATTAAAAAATCAACGCAAAAAATCGGACTTGCCTGCAAAGATGAGGAGGG
+ATTGCAGCGTGTTTTTAATGAGGTCATCACGGGATCCCATGTGCGTGACGGACATCGGGA
+AACGCCAAAGGAGATTATGTACCGAGGAAGAATGTCGCTGGACGGTATCGCGAAAATGTA
+TTCAGAAAATGATTATCAAGCCCTGTATCAGGACATGGTACGAGCTAAAAGATTCGATAC
+CGGCTCTTGTTCTGAGTCATGCGAAATATTTGGAGGGCAGCTTGATTTCGACTTCGGGAG
+GGAAGCTGCATGATGCGATGTTATCGGTGCGGTGAATGCAAAGAAGATAACCGCTTCCGA
+CCAAATCAACCTTACTGGAATCGATGGTGTCTCCGGTGTGAAAGAACACCAACAGGGGTG
+TTACCACTACCGCAGGAAAAGGAGGACGTGTGGCGAGACAGCGACGAAGTATCACCGACA
+TAATCTGCGAAAACTGCAAATACCTTCCAACGAAACGCACCAGAAATAAACCCAAGCCAA
+TCCCAAAAGAATCTGACGTAAAAACCTTCAACTACACGGCTCACCTGTGGGATATCCGGT
+GGCTAAGACGTCGTGCGAGGAAAACAAGGTGATTGACCAAAATCGAAGTTACGAACAAGA
+AAGCGTCGAGCGAGCTTTAACGTGCGCTAACTGCGGTCAGAAGCTGCATGTGCTGGAAGT
+TCACGTGTGTGAGCACTGCTGCGCAGAACTGATGAGCGATCCGAATAGCTCGATGCACGA
+GGAAGAAGATGATGGCTAAACCAGCGCGAAGACGATGTAAAAACGATGAATGCCGGGAAT
+GGTTTCACCCTGCATTCGCTAATCAGTGGTGGTGCTCTCCAGAGTGTGGAACCAAGATAG
+CACTCGAACGACGAAGTAAAGAACGCGAAAAAGCGGAAAAAGCAGCAGAGAAGAAACGAC
+GACGAGAGGAGCAGAAACAGAAAGATAAACTTAAGATTCGAAAACTCGCCTTAAAGCCCC
+GCAGTTACTGGATTAAACAAGCCCAACAAGCCGTAAACGCCTTCATCAGAGAAAGAGACC
+GCGACTTACCATGTATCTCGTGCGGAACGCTCACGTCTGCTCAGTGGGATGCCGGACATT
+ACCGGACAACTGCTGCGGCACCTCAACTCCGATTTAATGAACGCAATATTCACAAGCAAT
+GCGTGGTGTGCAACCAGCACAAAAGCGGAAATCTCGTTCCGTATCGCGTCGAACTGATTA
+GCCGCATCGGGCAGGAAGCAGTAGACGAAATCGAATCAAACCATAACCGCCATCGCTGGA
+CTATCGAAGAGTGCAAGGCGATCAAGGCAGAGTACCAACAGAAACTCAAAGACCTGCGAA
+ATAGCAGAAGTGAGGCCGCATGACGTTCTCAGTAAAAACCATTCCAGACATGCTCGTTGA
+AACATACGGAAATCAGACAGAAGTAGCACGCAGACTGAAATGTAGTCGCGGTACGGTCAG
+AAAATACGTTGATGATAAAGACGGGAAAATGCACGCCATCGTCAACGACGTTCTCATGGT
+TCATCGCGGATGGAGTGAAAGAGATGCGCTATTACGAAAAAATTGATGGCAGCAAATACC
+GAAATATTTGGGTAGTTGGCGATCTGCACGGATGCTACACGAACCTGATGAACAAACTGG
+ATACGATTGGATTCGACAACAAAAAAGACCTGCTTATCTCGGTGGGCGATTTGGTTGATC
+GTGGTGCAGAGAACGTTGAATGCCTGGAATTAATCACATTCCCCTGGTTCAGAGCTGTAC
+GTGGAAACCATGAGCAAATGATGATTGATGGCTTATCAGAGCGTGGAAACGTTAATCACT
+GGCTGCTTAATGGCGGTGGCTGGTTCTTTAATCTCGATTACGACAAAGAAATTCTGGCTA
+AAGCTCTTGCCCATAAAGCAGATGAACTTCCGTTAATCATCGAACTGGTGAGCAAAGATA
+AAAAATATGTTATCTGCCACGCCGATTATCCCTTTGACGAATACGAGTTTGGAAAGCCAG
+TTGATCATCAGCAGGTAATCTGGAACCGCGAACGAATCAGCAACTCACAAAACGGGATCG
+TGAAAGAAATCAAAGGCGCGGACACGTTCATCTTTGGTCATACGCCAGCAGTGAAACCAC
+TCAAGTTTGCCAACCAAATGTATATCGATACCGGCGCAGTGTTCTGCGGAAACCTAACAT
+TGATTCAGGTACAGGGAGAAGGCGCATGAGACTCGAAAGCGTAGCTAAATTTCATTCGCC
+AAAAAGCCCGATGATGAGCGACTCACCACGGGCCACGGCTTCTGACTCTCTTTCCGGTAC
+TGATGTGATGGCTGCTATGGGGATGGCGCAATCACAAGCCGGATTCGGTATGGCTGCATT
+CTGCGGTAAGCACGAACTCAGCCAGAACGACAAACAAAAGGCTATCAACTATCTGATGCA
+ATTTGCACACAAGGTATCGGGGAAATACCGTGGTGTGGCAAAGCTTGAAGGAAATACTAA
+GGCAAAGGTACTGCAAGTGCTCGCAACATTCGCTTATGCGGATTATTGCCGTAGTGCCGC
+GACGCCGGGGGCAAGATGCAGAGATTGCCATGGTACAGGCCGTGCGGTTGATATTGCCAA
+AACAGAGCTGTGGGGGAGAGTTGTCGAGAAAGAGTGCGGAAGATGCAAAGGCGTCGGCTA
+TTCAAGGATGCCAGCAAGCGCAGCATATCGCGCTGTGACGATGCTAATCCCAAACCTTAC
+CCAACCCACCTGGTCACGCACTGTTAAGCCGCTGTATGACGCTCTGGTGGTGCAATGCCA
+CAAAGAAGAGTCAATCGCAGACAACATTTTGAATGCGGTCACACGTTAGCAGCATGATTG
+CCACGGATGGCAACATATTAACGGCATGATATTGACTTATTGAATAAAATTGGGTAAATT
+TGACTCAACGATGGGTTAATTCGCTCGTTGTGGTAGTGAGATGAAAAGAGGCGGCGCTTA
+CTACCGATTCCGCCTAGTTGGTCACTTCGACGTATCGTCTGGAACTCCAACCATCGCAGG
+CAGAGAGGTCTGCAAAATGCAATCCCGAAACAGTTCGCAGGTAATAGTTAGAGCCTGCAT
+AACGGTTTCGGGATTTTTTATATCTGCACAACAGGTAAGAGCATTGAGTCGATAATCGTG
+AAGAGTCGGCGAGCCTGGTTAGCCAGTGCTCTTTCCGTTGTGCTGAATTAAGCGAATACC
+GGAAGCAGAACCGGATCACCAAATGCGTACAGGCGTCATCGCCGCCCAGCAACAGCACAA
+CCCAAACTGAGCCGTAGCCACTGTCTGTCCTGAATTCATTAGTAATAGTTACGCTGCGGC
+CTTTTACACATGACCTTCGTGAAAGCGGGTGGCAGGAGGTCGCGCTAACAACCTCCTGCC
+GTTTTGCCCGTGCATATCGGTCACGAACAAATCTGATTACTAAACACAGTAGCCTGGATT
+TGTTCTATCAGTAATCGACCTTATTCCTAATTAAATAGAGCAAATCCCCTTATTGGGGGT
+AAGACATGAAGATGCCAGAAAAACATGACCTGTTGGCCGCCATTCTCGCGGCAAAGGAAC
+AAGGCATCGGGGCAATCCTTGCGTTTGCAATGGCGTACCTTCGCGGCAGATATAATGGCG
+GTGCGTTTACAAAAACAGTAATCGACGCAACGATGTGCGCCATTATCGCCTAGTTCATTC
+GTGACCTTCTCGACTTCGCCGGACTAAGTAGCAATCTCGCTTATATAACGAGCGTGTTTA
+TCGGCTACATCGGTACTGACTCGATTGGTTCGCTTATCAAACGCTTCGCTGCTAAAAAAG
+CCGGAGTAGAAGATGGTAGAAATCAATAATCAACGTAAGGCGTTCCTCGATATGCTGGCG
+TGGTCGGAGGGAACTGATAACGGACGTCAGAAAACCAGAAATCATGGTTATGACGTCATT
+GTAGGCGGAGAGCTATTTACTGATTACTCCGATCACCCTCGCAAACTTGTCACGCTAAAC
+CCAAAACTCAAATCAACAGGCGCCGGACGCTACCAGCTTCTTTCCCGTTGGTGGGATGCC
+TACCGCAAGCAGCTTGGCCTGAAAGACTTCTCTCCGAAAAGTCAGGACGCTGTGGCATTG
+CAGCAGATTAAGGAGCGTGGCGCTTTACCTATGATTGATCGTGGTGATATCCGTCAGGCA
+ATCGACCGTTGCAGCAATATCTGGGCTTCACTGCCGGGCGCTGGTTATGGTCAGTTCGAG
+CATAAGGCTGACAGCCTGATTGCAAAATTCAAAGAAGCGGGCGGAACGGTCAGAGAGATT
+GATGTATGAGCAGAGTCACCGCGATTATCTCCGCTCTGGTTATCTGCATCATCGTCTGCC
+TGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACA
+AAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGC
+GTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAA
+ATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCT
+GTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGAC
+TGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAA
+AACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGCCCATATC
+GATGGGCAACTCATGCAATTATTGTGAGCAATACACACGCGCTTCCAGCGGAGTATAAAT
+GCCTAAAGTAATAAAACCGAGCAATCCATTTACGAATGTTTGCTGGGTTTCTGTTTTAAC
+AACATTTTCTGCGCCGCCACAAATTTTGGCTGCATCGACAGTTTTCTTCTGCCCAATTCC
+AGAAACGAAGAAATGATGGGTGATGGTTTCCTTTGGTGCTACTGCTGCCGGTTTGTTTTG
+AACAGTAAACGTCTGTTGAGCACATCCTGTAATAAGCAGGGCCAGCGCAGTAGCGAGTAG
+CATTTTTTTCATGGTGTTATTCCCGATGCTTTTTGAAGTTCGCAGAATCGTATGTGTAGA
+AAATTAAACAAACCCTAAACAATGAGTTGAAATTTCATATTGTTAATATTTATTAATGTA
+TGTCAGGTGCGATGAATCGTCATTGTATTCCCGGATTAACTATGTCCACAGCCCTGACGG
+GGAACTTCTCTGCGGGAGTGTCCGGGAATAATTAAAACGATGCACACAGGGTTTAGCGCG
+TACACGTATTGCATTATGCCAACGCCCCGGTGCTGACACGGAAGAAACCGGACGTTATGA
+TTTAGCGTGGAAAGATTTGTGTAGTGTTCTGAATGCTCTCAGTAAATAGTAATGAATTAT
+CAAAGGTATAGTAATATCTTTTATGTTCATGGATATTTGTAACCCATCGGAAAACTCCTG
+CTTTAGCAAGATTTTCCCTGTATTGCTGAAATGTGATTTCTCTTGATTTCAACCTATCAT
+AGGACGTTTCTATAAGATGCGTGTTTCTTGAGAATTTAACATTTACAACCTTTTTAAGTC
+CTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGAT
+AGTAAATATAATGTGAGACGTTGTGACGTTTTAGTTCAGAATAAAACAATTCACAGTCTA
+AATCTTTTCGCACTTGATCGAATATTTCTTTAAAAATGGCAACCTGAGCCATTGGTAAAA
+CCTTCCATGTGATACGAGGGCGCGTAGTTTGCATTATCGTTTTTATCGTTTCAATCTGGT
+CTGACCTCCTTGTGTTTTGTTGATGATTTATGTCAAATATTAGGAATGTTTTCACTTAAT
+AGTATTGGTTGCGTAACAAAGTGCGGTCCTGCTGGCATTCTGGAGGGAAATACAACCGAC
+AGATGTATGTAAGGCCAACGTGCTCAAATCTTCATACAGAAAGATTTGAAGTAATATTTT
+AACCGCTAGATGAAGAGCAAGCGCATGGAGCGACAAAATGAATAAAGAACAATCTGCTGA
+TGATCCCTCCGTGGATCTGATTCGTGTAAAAAATATGCTTAATAGCACCATTTCTATGAG
+TTACCCTGATGTTGTAATTGCATGTATAGAACATAAGGTGTCTCTGGAAGCATTCAGAGC
+AATTGAGGCAGCGTTGGTGAAGCACGATAATAATATGAAGGATTATTCCCTGGTGGTTGA
+CTGATCACCATAACTGCTAATCATTCAAACTATTTAGTCTGTGACAGAGCCAACACGCAG
+TCTGTCACTGTCAGGAAAGTGGTAAAACTGCAACTCAATTACTGCAATGCCCTCGTAATT
+AAGTGAATTTACAATATCGTCCTGTTCGGAGGGAAGAACGCGGGATGTTCATTCTTCATC
+ACTTTTAATTGATGTATATGCTCTCTTTTCTGACGTTAGTCTCCGACGGCAGGCTTCAAT
+GACCCAGGCTGAGAAATTCCCGGACCCTTTTTGCTCAAGAGCGATGTTAATTTGTTCAAT
+CATTTGGTTAGGAAAGCGGATGTTGCGGGTTGTTGTTCTGCGGGTTCTGTTCTTCGTTGA
+CATGAGGTTGCCCCGTATTCAGTGTCGCTGATTTGTATTGTCTGAAGTTGTTTTTACGTT
+AAGTTGATGCAGATCAATTAATACGATACCTGCGTCATAATTGATTATTTGACGTGGTTT
+GATGGCCTCCACGCACGTTGTGATATGTAGATGATAATCATTATCACTTTACGGGTCCTT
+TCCGGTGATCCGACAGGTTACG
diff --git a/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.contig.index b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.contig.index
new file mode 100644
index 0000000..a8732a3
--- /dev/null
+++ b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.contig.index
@@ -0,0 +1,2 @@
+>lambda_NEB3011
+ref000001 0
diff --git a/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.fai b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.fai
new file mode 100644
index 0000000..7bc16a4
--- /dev/null
+++ b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.fai
@@ -0,0 +1 @@
+lambda_NEB3011 48502 16 60 61
diff --git a/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.index b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.index
new file mode 100644
index 0000000..c21d4f7
--- /dev/null
+++ b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.index
@@ -0,0 +1,27 @@
+>lambda_NEB3011
+1 16
+2000 2048
+4000 4081
+6000 6114
+8000 8148
+10000 10181
+12000 12214
+14000 14248
+16000 16281
+18000 18314
+20000 20348
+22000 22381
+24000 24414
+26000 26448
+28000 28481
+30000 30514
+32000 32548
+34000 34581
+36000 36614
+38000 38648
+40000 40681
+42000 42714
+44000 44748
+46000 46781
+48000 48814
+48502 49325
\ No newline at end of file
diff --git a/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.sa b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.sa
new file mode 100644
index 0000000..df41c12
Binary files /dev/null and b/data/ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta.sa differ
diff --git a/data/SubreadSet/barcoded.scraps.bam b/data/SubreadSet/barcoded.scraps.bam
new file mode 100644
index 0000000..36c03d2
--- /dev/null
+++ b/data/SubreadSet/barcoded.scraps.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c8197573c7d3f20fcfb37015ce08171d998b48c8d24109c2a84e9b26a56456c7
+size 10332
diff --git a/data/SubreadSet/barcoded.scraps.bam.pbi b/data/SubreadSet/barcoded.scraps.bam.pbi
new file mode 100644
index 0000000..99a3619
--- /dev/null
+++ b/data/SubreadSet/barcoded.scraps.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:1fbc05ab649c467e6e02ffbf6bc6d84ae60bc5a3e77b88a49875703467995ccc
+size 113
diff --git a/data/SubreadSet/barcoded.subreads.bam b/data/SubreadSet/barcoded.subreads.bam
new file mode 100644
index 0000000..202102b
--- /dev/null
+++ b/data/SubreadSet/barcoded.subreads.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:2cc171d81c1da3bafe7ab4c868cb59f5df82ef2c7a6ffff5e6cc173d25c4f243
+size 4025
diff --git a/data/SubreadSet/barcoded.subreads.bam.bai b/data/SubreadSet/barcoded.subreads.bam.bai
new file mode 100644
index 0000000..aedcb14
--- /dev/null
+++ b/data/SubreadSet/barcoded.subreads.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:a6ff93b0d4fab5efadb2301aad792598ab4d1ebb497c337f61ca11d8cd9059ae
+size 16
diff --git a/data/SubreadSet/barcoded.subreads.bam.pbi b/data/SubreadSet/barcoded.subreads.bam.pbi
new file mode 100644
index 0000000..b7064cc
--- /dev/null
+++ b/data/SubreadSet/barcoded.subreads.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:29cd41b01b9244731872d07696348e340ab8b32ef7afee4463179d98e8d83d65
+size 121
diff --git a/data/SubreadSet/barcoded.subreadset.xml b/data/SubreadSet/barcoded.subreadset.xml
new file mode 100644
index 0000000..d1b7c09
--- /dev/null
+++ b/data/SubreadSet/barcoded.subreadset.xml
@@ -0,0 +1,22 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<SubreadSet CreatedAt="2016-03-31T18:18:39" MetaType="PacBio.DataSet.SubreadSet" Name="barcoded-subreadset" Tags="" TimeStampedName="pacbio_dataset_subreadset-160401_011839830" UniqueId="ceec5b21-69cc-0c4b-91db-d00835210816" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.SubreadBamFile" ResourceId="barcoded.subreads.bam" TimeStampedName="pacbio_subreadfile_subreadbamfile-160401_011839830" UniqueId="1c4a4fa3-2c6a-4d0a-867d-9bbb932af555">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="barcoded.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160401_011839830" UniqueId="c17c4401-0f56-40b5-85c6-46f7771105c7"/>
+ </pbbase:FileIndices>
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource Description="Points to the scraps bam file." MetaType="PacBio.SubreadFile.ScrapsBamFile" Name="scraps bam" ResourceId="barcoded.scraps.bam" TimeStampedName="pacbio_subreadfile_scrapsbamfile-160219_060028519" UniqueId="c7fc0c63-234b-4b93-9976-d7d82cb03e71" Version="3.0.1">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="barcoded.scraps.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160219_060028519" UniqueId="77ecadb3-e2b0-48b4-80ca-3de1e87b339d" Version="3.0.1" />
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ <pbbase:ExternalResource Description="BarcodeSet" MetaType="PacBio.DataSet.BarcodeSet" Name="barcodeset" ResourceId="../BarcodeSet/example_barcodes.barcodeset.xml" UniqueId="afa1b6fd-6a44-50d8-838f-cb9d485ce9da"/>
+ </pbbase:ExternalResources>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>1640</TotalLength>
+ <NumRecords>2</NumRecords>
+ </DataSetMetadata>
+</SubreadSet>
diff --git a/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam
new file mode 100644
index 0000000..b6842cd
--- /dev/null
+++ b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:a97b7b1adaf67172132d81b3886ad9d73e10228ffe3d92b946137309840a3d05
+size 204543
diff --git a/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.bai b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.bai
new file mode 100644
index 0000000..9432789
--- /dev/null
+++ b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.bai
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:551ebf0bf906f881c0079d0ebf35dffae5644369db72fe720ddd60fd0fba15c5
+size 16
diff --git a/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.pbi b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.pbi
new file mode 100644
index 0000000..6a1b1b7
--- /dev/null
+++ b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:659dd0c0e8cd03a3de37c931d3f10ba17f1de7f6788f3268567e3d1c4572c0fc
+size 1517
diff --git a/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreadset.xml b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreadset.xml
new file mode 100644
index 0000000..2ba54da
--- /dev/null
+++ b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreadset.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<SubreadSet CreatedAt="2016-05-31T15:23:03" MetaType="PacBio.DataSet.SubreadSet" Name="subreads-xml" Tags="" TimeStampedName="pacbio_dataset_subreadset-160531_222303338" UniqueId="af916621-ab06-46cc-94f2-bc922c48b277" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.SubreadBamFile" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam" TimeStampedName="pacbio_subreadfile_subreadbamfile-160531_222303338" UniqueId="51549ac1-f5d4-4130-9699-e2713079654f">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160531_222303338" UniqueId="843290b6-23bb-422c-88a8-0bb4c55378e9"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160531_222303338" UniqueId="09b3b42a-73ea-4b84-847e-b87a321abae5"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>62340</TotalLength>
+ <NumRecords>117</NumRecords>
+ </DataSetMetadata>
+</SubreadSet>
diff --git a/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0_filtered.subreadset.xml b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0_filtered.subreadset.xml
new file mode 100644
index 0000000..e36ae36
--- /dev/null
+++ b/data/SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0_filtered.subreadset.xml
@@ -0,0 +1,22 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<SubreadSet CreatedAt="2016-05-31T15:23:03" MetaType="PacBio.DataSet.SubreadSet" Name="subreads-filtered" Tags="" TimeStampedName="pacbio_dataset_subreadset-160531_222303338" UniqueId="af916621-ab06-46cc-94f2-bc922c48b277" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://pacificbiosciences.com/PacBioDataModel.xsd">
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.SubreadBamFile" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam" TimeStampedName="pacbio_subreadfile_subreadbamfile-160531_222303338" UniqueId="51549ac1-f5d4-4130-9699-e2713079654f">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160531_222303338" UniqueId="843290b6-23bb-422c-88a8-0bb4c55378e9"/>
+ <pbbase:FileIndex MetaType="PacBio.Index.BamIndex" ResourceId="m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam.bai" TimeStampedName="pacbio_index_bamindex-160531_222303338" UniqueId="09b3b42a-73ea-4b84-847e-b87a321abae5"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <Filters>
+ <Filter>
+ <pbbase:Properties>
+ <pbbase:Property Name="rq" Operator=">" Value="0.901"/>
+ </pbbase:Properties>
+ </Filter>
+ </Filters>
+ <DataSetMetadata>
+ <TotalLength>40910</TotalLength>
+ <NumRecords>78</NumRecords>
+ </DataSetMetadata>
+</SubreadSet>
diff --git a/data/SubreadSet/m54006_160504_020705.sts.xml b/data/SubreadSet/m54006_160504_020705.sts.xml
new file mode 100644
index 0000000..5e93a78
--- /dev/null
+++ b/data/SubreadSet/m54006_160504_020705.sts.xml
@@ -0,0 +1 @@
+<?xml version="1.0" encoding="utf-8"?><PipeStats xmlns="http://pacificbiosciences.com/PacBioPipelineStats.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ns="http://pacificbiosciences.com/PacBioBaseDataModel.xsd"><MovieName>/data/pb/m54006_160504_020705.baz</MovieName><MovieLength>120</MovieLength><NumSequencingZmws>1035426</NumSequencingZmws><TotalBaseFractionPerChannel Channel="A"><TotalBaseFractionValue>0.280779</TotalBaseFractionValue></TotalBaseFractionPerChannel><T [...]
\ No newline at end of file
diff --git a/data/SubreadSet/m54006_160504_020705.tiny.scraps.bam b/data/SubreadSet/m54006_160504_020705.tiny.scraps.bam
new file mode 100644
index 0000000..b93b267
--- /dev/null
+++ b/data/SubreadSet/m54006_160504_020705.tiny.scraps.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3bc87753ff0341227db7b8a708015243c835ef56e2660663de54e50ac8305291
+size 5375
diff --git a/data/SubreadSet/m54006_160504_020705.tiny.scraps.bam.pbi b/data/SubreadSet/m54006_160504_020705.tiny.scraps.bam.pbi
new file mode 100644
index 0000000..32a80e9
--- /dev/null
+++ b/data/SubreadSet/m54006_160504_020705.tiny.scraps.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:841b3f7f0b91c149d72994355cdafdf72acfa8ada260c0b4d4ae9108a90b609a
+size 329
diff --git a/data/SubreadSet/m54006_160504_020705.tiny.subreads.bam b/data/SubreadSet/m54006_160504_020705.tiny.subreads.bam
new file mode 100644
index 0000000..2fc5585
--- /dev/null
+++ b/data/SubreadSet/m54006_160504_020705.tiny.subreads.bam
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0a19d7c18bb5f82442f0407d61f78943f49ae9188b670ac063ec6a2abcb732e1
+size 37731
diff --git a/data/SubreadSet/m54006_160504_020705.tiny.subreads.bam.pbi b/data/SubreadSet/m54006_160504_020705.tiny.subreads.bam.pbi
new file mode 100644
index 0000000..5ac53c8
--- /dev/null
+++ b/data/SubreadSet/m54006_160504_020705.tiny.subreads.bam.pbi
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:29fbb5baf60b45bdda3690ae8e99814a4954ba4f98334202b8f466294413530e
+size 315
diff --git a/data/SubreadSet/m54006_160504_020705.tiny.subreadset.xml b/data/SubreadSet/m54006_160504_020705.tiny.subreadset.xml
new file mode 100644
index 0000000..9868be0
--- /dev/null
+++ b/data/SubreadSet/m54006_160504_020705.tiny.subreadset.xml
@@ -0,0 +1,97 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<SubreadSet CreatedAt="2016-06-13T16:26:20" MetaType="PacBio.DataSet.SubreadSet" Name="subreads-sequel" Tags="" TimeStampedName="pacbio_dataset_subreadset-160613_232620602" UniqueId="c9074f64-37e0-d795-9b8e-7190ac6f87a5" Version="3.0.1" xmlns="http://pacificbiosciences.com/PacBioDatasets.xsd" xmlns:ns4="http://pacificbiosciences.com/PacBioReagentKit.xsd" xmlns:pbbase="http://pacificbiosciences.com/PacBioBaseDataModel.xsd" xmlns:pbmeta="http://pacificbiosciences.com/PacBioCollectionMetada [...]
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.SubreadBamFile" ResourceId="m54006_160504_020705.tiny.subreads.bam" TimeStampedName="pacbio_subreadfile_subreadbamfile-160613_232620601" UniqueId="b5412e74-6eb4-4a63-bd82-731c87fe0122">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m54006_160504_020705.tiny.subreads.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160613_232620601" UniqueId="d356c436-892d-4199-b946-0f4fbdb834da"/>
+ </pbbase:FileIndices>
+ <pbbase:ExternalResources>
+ <pbbase:ExternalResource MetaType="PacBio.SubreadFile.ScrapsBamFile" ResourceId="m54006_160504_020705.tiny.scraps.bam" TimeStampedName="pacbio_subreadfile_scrapsbamfile-160613_232620601" UniqueId="824fa041-e354-4fc2-9669-ae6ae54a80fc">
+ <pbbase:FileIndices>
+ <pbbase:FileIndex MetaType="PacBio.Index.PacBioIndex" ResourceId="m54006_160504_020705.tiny.scraps.bam.pbi" TimeStampedName="pacbio_index_pacbioindex-160504_063727354" UniqueId="cc0d00fe-07da-4a6b-bf42-746d3fde43ea" Version="3.0.1"/>
+ </pbbase:FileIndices>
+ </pbbase:ExternalResource>
+ <pbbase:ExternalResource Description="Points to the summary sts.xml file." MetaType="PacBio.SubreadFile.ChipStatsFile" Name="Chipstats XML" ResourceId="m54006_160504_020705.sts.xml" TimeStampedName="pacbio_subreadfile_chipstatsfile-160504_063727355" UniqueId="60051785-dbc1-4760-a117-b67d50f105bd" Version="3.0.1"/>
+ </pbbase:ExternalResources>
+ </pbbase:ExternalResource>
+ </pbbase:ExternalResources>
+ <DataSetMetadata>
+ <TotalLength>19014</TotalLength>
+ <NumRecords>20</NumRecords>
+ <pbmeta:Collections>
+ <pbmeta:CollectionMetadata Context="m54006_160504_020705" CreatedAt="2016-05-03T23:37:42.775Z" InstrumentId="54006" InstrumentName="Alpha6" MetaType="CollectionMetadata" ModifiedAt="0001-01-01T00:00:00" Status="Ready" TimeStampedName="54006-CollectionMetadata-2016-37-03T23:37:56.212Z" UniqueId="7fac56f8-92bf-44f1-a474-7d3a96a8d347">
+ <pbmeta:InstCtrlVer>3.1.0.179041</pbmeta:InstCtrlVer>
+ <pbmeta:SigProcVer>3.1.0-179029</pbmeta:SigProcVer>
+ <pbmeta:RunDetails>
+ <pbmeta:TimeStampedName>r54006_20160503_233756</pbmeta:TimeStampedName>
+ <pbmeta:Name>lambda</pbmeta:Name>
+ <pbmeta:WhenCreated>2016-05-03T23:37:42.775Z</pbmeta:WhenCreated>
+ <pbmeta:WhenStarted>2016-05-03T23:41:59.065365Z</pbmeta:WhenStarted>
+ </pbmeta:RunDetails>
+ <pbmeta:WellSample CreatedAt="2016-05-03T23:37:42.775Z" ModifiedAt="0001-01-01T00:00:00" Name="2k_lambda_Mag_2hr">
+ <pbmeta:WellName>B01</pbmeta:WellName>
+ <pbmeta:Concentration>0.0</pbmeta:Concentration>
+ <pbmeta:InsertSize>10000</pbmeta:InsertSize>
+ <pbmeta:SampleReuseEnabled>false</pbmeta:SampleReuseEnabled>
+ <pbmeta:StageHotstartEnabled>false</pbmeta:StageHotstartEnabled>
+ <pbmeta:SizeSelectionEnabled>false</pbmeta:SizeSelectionEnabled>
+ <pbmeta:UseCount>0</pbmeta:UseCount>
+ <pbmeta:DNAControlComplex>Sequel™ SMRT®Cell Control Complex 1.0</pbmeta:DNAControlComplex>
+ </pbmeta:WellSample>
+ <pbmeta:Automation Name="Workflow_Magbead.py">
+ <pbbase:AutomationParameters>
+ <pbbase:AutomationParameter CreatedAt="2016-05-03T23:37:42.775Z" ModifiedAt="0001-01-01T00:00:00" Name="MovieLength" SimpleValue="120" ValueDataType="Double"/>
+ <pbbase:AutomationParameter CreatedAt="2016-05-03T23:37:42.775Z" ModifiedAt="0001-01-01T00:00:00" Name="ImmobilizationTime" SimpleValue="120" ValueDataType="Double"/>
+ <pbbase:AutomationParameter Name="CollectionNumber" SimpleValue="1" ValueDataType="Int32"/>
+ <pbbase:AutomationParameter Name="UseStageHotStart" SimpleValue="False" ValueDataType="String"/>
+ <pbbase:AutomationParameter Name="InsertSize" SimpleValue="10000" ValueDataType="Int32"/>
+ </pbbase:AutomationParameters>
+ </pbmeta:Automation>
+ <pbmeta:CollectionNumber>1</pbmeta:CollectionNumber>
+ <pbmeta:CellIndex>1</pbmeta:CellIndex>
+ <pbmeta:CellPac Barcode="BA010713635979244259265850" ExpirationDate="2016-11-04" LotNumber="320296" PartNumber="100-512-700"/>
+ <pbmeta:TemplatePrepKit Barcode="DM1117100259100111716" ExpirationDate="2016-11-17" LotNumber="DM1117" MaxInsertSize="20000" MinInsertSize="500" Name="SMRTbell™ Template Prep Kit" PartNumber="100-259-100">
+ <pbbase:LeftAdaptorSequence>ATCTCTCTCAACAACAACAACGGAGGAGGAGGAAAAGAGAGAGAT</pbbase:LeftAdaptorSequence>
+ <pbbase:RightAdaptorSequence>ATCTCTCTCAACAACAACAACGGAGGAGGAGGAAAAGAGAGAGAT</pbbase:RightAdaptorSequence>
+ </pbmeta:TemplatePrepKit>
+ <pbmeta:BindingKit Barcode="DM1117100619300111716" ExpirationDate="2016-11-17" LotNumber="DM1117" Name="Sequel™ Binding Kit 1.0" PartNumber="100-619-300"/>
+ <pbmeta:SequencingKitPlate Barcode="160318100620000031817" ExpirationDate="2017-03-18" LotNumber="160318" Name="Sequel™ Sequencing Plate 1.0" PartNumber="100-620-000">
+ <ns4:ReagentTubes ExpirationDate="2018-05-12" LotNumber="706520" PartNumber="100-619-700"/>
+ <ns4:ReagentTubes ExpirationDate="2020-02-28" LotNumber="007045" PartNumber="100-619-600"/>
+ </pbmeta:SequencingKitPlate>
+ <pbmeta:Primary>
+ <pbmeta:AutomationName>Default</pbmeta:AutomationName>
+ <pbmeta:ConfigFileName>SqlPoC_SubCrf_2C2A-t2.xml</pbmeta:ConfigFileName>
+ <pbmeta:SequencingCondition>DefaultPrimarySequencingCondition</pbmeta:SequencingCondition>
+ <pbmeta:OutputOptions>
+ <pbmeta:ResultsFolder>315/3150198/r54006_20160503_233756/2_B01/</pbmeta:ResultsFolder>
+ <pbmeta:CollectionPathUri>/pbi/collections/315/3150198/r54006_20160503_233756/2_B01/</pbmeta:CollectionPathUri>
+ <pbmeta:CopyFiles>
+ <pbmeta:CollectionFileCopy>Fasta</pbmeta:CollectionFileCopy>
+ <pbmeta:CollectionFileCopy>Bam</pbmeta:CollectionFileCopy>
+ </pbmeta:CopyFiles>
+ <pbmeta:Readout>Bases_Without_QVs</pbmeta:Readout>
+ <pbmeta:MetricsVerbosity>Minimal</pbmeta:MetricsVerbosity>
+ <pbmeta:TransferResource>
+ <pbmeta:Id>rsync-pbi-collections</pbmeta:Id>
+ <pbmeta:TransferScheme>RSYNC</pbmeta:TransferScheme>
+ <pbmeta:Name>PBI Collections Transfer test using Rsync Server</pbmeta:Name>
+ <pbmeta:Description>Location for writing Transfer services to write to. For testing, Internal tools (PA SIM and ICS) tests must explicitly set the relative path prefix to 'xfer-test'</pbmeta:Description>
+ <pbmeta:DestPath>/pbi/collections</pbmeta:DestPath>
+ </pbmeta:TransferResource>
+ </pbmeta:OutputOptions>
+ </pbmeta:Primary>
+ <pbmeta:Secondary>
+ <pbmeta:AutomationName>DefaultSecondaryAutomationName</pbmeta:AutomationName>
+ <pbmeta:AutomationParameters>
+ <pbmeta:AutomationParameter CreatedAt="0001-01-01T00:00:00" ModifiedAt="0001-01-01T00:00:00" Name="Reference" SimpleValue="DefaultSecondaryAnalysisReferenceName" ValueDataType="String"/>
+ </pbmeta:AutomationParameters>
+ <pbmeta:CellCountInJob>0</pbmeta:CellCountInJob>
+ </pbmeta:Secondary>
+ <pbmeta:UserDefinedFields>
+ <pbbase:DataEntities Name=" LIMS_IMPORT " SimpleValue="DefaultUserDefinedFieldLIMS" ValueDataType="String"/>
+ </pbmeta:UserDefinedFields>
+ </pbmeta:CollectionMetadata>
+ </pbmeta:Collections>
+ </DataSetMetadata>
+</SubreadSet>
diff --git a/data/files.json b/data/files.json
new file mode 100644
index 0000000..c73cc7e
--- /dev/null
+++ b/data/files.json
@@ -0,0 +1,158 @@
+[
+ {
+ "id": "subreads-xml",
+ "description": "SubreadSet with RSII lambda data",
+ "path": "SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreadset.xml",
+ "filetype": "PacBio.DataSet.SubreadSet"
+ },
+ {
+ "id": "subreads-sequel",
+ "description": "SubreadSet with Sequel data and sts.xml",
+ "path": "SubreadSet/m54006_160504_020705.tiny.subreadset.xml",
+ "filetype": "PacBio.DataSet.SubreadSet"
+ },
+ {
+ "id": "subreads-filtered",
+ "description": "SubreadSet with readQual filter",
+ "path": "SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0_filtered.subreadset.xml",
+ "filetype": "PacBio.DataSet.SubreadSet"
+ },
+ {
+ "id": "stats-xml",
+ "description": "XML of primary analysis statistics",
+ "path": "SubreadSet/m54006_160504_020705.sts.xml",
+ "filetype": "PacBio.SubreadFile.ChipStatsFile"
+ },
+ {
+ "id": "aligned-xml",
+ "description": "AlignmentSet with RSII lambda data",
+ "path": "AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.alignmentset.xml",
+ "filetype": "PacBio.DataSet.AlignmentSet"
+ },
+ {
+ "id": "lambdaNEB",
+ "description": "Lambda virus ReferenceSet",
+ "path": "ReferenceSet/lambdaNEB/referenceset.xml",
+ "filetype": "PacBio.DataSet.ReferenceSet"
+ },
+ {
+ "id": "subreads-bam",
+ "description": "BAM file with unmapped RSII subreads",
+ "path": "SubreadSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.subreads.bam",
+ "filetype": "PacBio.SubreadFile.SubreadBamFile"
+ },
+ {
+ "id": "aligned-bam",
+ "description": "RSII BAM file mapped to lambda reference",
+ "path": "AlignmentSet/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.bam",
+ "filetype": "PacBio.AlignmentFile.AlignmentBamFile"
+ },
+ {
+ "id": "lambda-fasta",
+ "description": "Lambda virus genome FASTA",
+ "path": "ReferenceSet/lambdaNEB/sequence/lambdaNEB.fasta",
+ "filetype": "PacBio.FileTypes.Fasta"
+ },
+ {
+ "id": "cmph5",
+ "description": "Aligned RSII sequences in legacy cmp.h5 format",
+ "path": "legacy/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.cmp.h5",
+ "filetype": "PacBio.FileTypes.alignment_cmp_h5"
+ },
+ {
+ "id": "barcoded-subreadset",
+ "description": "Barcoded SubreadSet",
+ "path": "SubreadSet/barcoded.subreadset.xml",
+ "filetype": "PacBio.DataSet.SubreadSet"
+ },
+ {
+ "id": "barcoded-subreads-bam",
+ "description": "Barcoded subreads BAM file",
+ "path": "SubreadSet/barcoded.subreads.bam",
+ "filetype": "PacBio.SubreadFile.SubreadBamFile"
+ },
+ {
+ "id": "barcodeset",
+ "description": "BarcodeSet XML",
+ "path": "BarcodeSet/example_barcodes.barcodeset.xml",
+ "filetype": "PacBio.DataSet.BarcodeSet"
+ },
+ {
+ "id": "alignment-summary-gff",
+ "description": "Alignment summary GFF from resequencing reports",
+ "path": "GFF/alignment_summary.gff",
+ "filetype": "PacBio.FileTypes.gff"
+ },
+ {
+ "id": "variants-gff",
+ "description": "Variants GFF from GenomicConsensus/Quiver",
+ "path": "GFF/variants.gff",
+ "filetype": "PacBio.FileTypes.gff"
+ },
+ {
+ "id": "consensus-summary-gff",
+ "description": "Expanded alignment summary GFF with consensus calling annotations",
+ "path": "GFF/alignment_summary_consensus.gff",
+ "filetype": "PacBio.FileTypes.gff"
+ },
+ {
+ "id": "rsii-ccs-aligned",
+ "description": "ConsensusAlignmentSet of RSII lambda data",
+ "path": "ConsensusAlignmentSet/mapped.consensusalignmentset.xml",
+ "filetype": "PacBio.DataSet.ConsensusAlignmentSet"
+ },
+ {
+ "id": "ccs-bam-aligned",
+ "description": "Mapped CCS BAM file",
+ "path": "ConsensusAlignmentSet/mapped.consensusalignmentset.bam",
+ "filetype": "PacBio.AlignmentFile.ConsensusAlignmentBamFile"
+ },
+ {
+ "id": "rsii-ccs",
+ "description": "ConsensusReadSet of RSII lambda data",
+ "path": "ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.consensusreadset.xml",
+ "filetype": "PacBio.DataSet.ConsensusReadSet"
+ },
+ {
+ "id": "ccs-bam",
+ "description": "CCS BAM file",
+ "path": "ConsensusReadSet/m130727_114215_42211_c100569412550000001823090301191423_s1_p0.ccs.bam",
+ "filetype": "PacBio.ConsensusReadFile.ConsensusReadBamFile"
+ },
+ {
+ "id": "rsii-ccs-multi-cell",
+ "description": "ConsensusReadSet from multiple RSII cells",
+ "path": "ConsensusReadSet/ccs_mixed.consensusreadset.xml",
+ "filetype": "PacBio.DataSet.ConsensusReadSet"
+ },
+ {
+ "id": "hdfsubreads",
+ "description": "HdfSubreadSet of RSII data",
+ "path": "HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.all.hdfsubreadset.xml",
+ "filetype": "PacBio.DataSet.HdfSubreadSet"
+ },
+ {
+ "id": "rsii-bax-h5",
+ "description": "RSII bax.h5 file",
+ "path": "HdfSubreadSet/m140912_020930_00114_c100702482550000001823141103261590_s1_p0.1.bax.h5",
+ "filetype": "PacBio.SubreadFile.BaxFile"
+ },
+ {
+ "id": "contigset",
+ "description": "Minimal mock ContigSet",
+ "path": "ContigSet/simple.contigset.xml",
+ "filetype": "PacBio.DataSet.ContigSet"
+ },
+ {
+ "id": "ccs-barcoded",
+ "description": "Barcoded ConsensusReadSet",
+ "path": "ConsensusReadSet/barcoded.consensusreadset.xml",
+ "filetype": "PacBio.DataSet.ConsensusReadSet"
+ },
+ {
+ "id": "aligned-ds-2",
+ "description": "Chunked AlignmentSet",
+ "path": "AlignmentSet/m54006_160504_020705.alignmentset.xml",
+ "filetype": "PacBio.DataSet.AlignmentSet"
+ }
+]
diff --git a/data/legacy/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.cmp.h5 b/data/legacy/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.cmp.h5
new file mode 100644
index 0000000..0bcb00d
--- /dev/null
+++ b/data/legacy/m140905_042212_sidney_c100564852550000001823085912221377_s1_X0.aligned_subreads.cmp.h5
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:22166497cccce5da73beaa65f55e45376a2f0b4cfbe3efbaf034a12f8de90e8a
+size 297430
diff --git a/src/python/Makefile b/src/python/Makefile
new file mode 100644
index 0000000..16ee088
--- /dev/null
+++ b/src/python/Makefile
@@ -0,0 +1,15 @@
+.PHONY: all clean install
+
+all: install
+
+install:
+ python setup.py install
+
+clean:
+ rm -rf build dist
+ rm -rf pbtestdata/data
+ rm -rf pbtestdata/version.py
+ find . -name "*.egg-info" | xargs rm -rf
+ find . -name "*.pyc" | xargs rm -f
+ find . -name "*.err" | xargs rm -f
+ find . -name "*.log" | xargs rm -f
diff --git a/src/python/pbtestdata/__init__.py b/src/python/pbtestdata/__init__.py
new file mode 100644
index 0000000..108c0c3
--- /dev/null
+++ b/src/python/pbtestdata/__init__.py
@@ -0,0 +1 @@
+from .core import get_file
diff --git a/src/python/pbtestdata/core.py b/src/python/pbtestdata/core.py
new file mode 100644
index 0000000..1041916
--- /dev/null
+++ b/src/python/pbtestdata/core.py
@@ -0,0 +1,118 @@
+
+"""
+Utility and API for accessing representative PacBio data files for testing.
+Run 'pbdata show [--verbose]' to display a list of files sorted by type.
+"""
+
+from collections import defaultdict
+import pkg_resources
+import argparse
+import json
+import re
+import os.path as op
+import sys
+
+import version
+
+RAW_DICTS = None
+DATA_DICT = None
+
+
+def _load_data():
+ global DATA_DICT
+ global RAW_DICTS
+ if DATA_DICT is None:
+ RAW_DICTS = json.loads(pkg_resources.resource_string("pbtestdata",
+ "data/files.json"))
+ DATA_DICT = {f["id"]:f for f in RAW_DICTS}
+ return DATA_DICT
+
+
+def get_file(id_):
+ d = _load_data()
+ try:
+ file_name = pkg_resources.resource_filename("pbtestdata",
+ op.join("data", d[id_]['path']))
+ except KeyError:
+ raise KeyError("No data file with ID '{i}' is available".format(i=id_))
+ else:
+ if not op.isfile(file_name):
+ raise IOError("The data file {f} does not exist.".format(f=file_name))
+ return file_name
+
+
+def validate():
+ global RAW_DICTS
+ d = _load_data()
+ have_keys = set()
+ have_files = set()
+ for f in RAW_DICTS:
+ file_path = pkg_resources.resource_filename("pbtestdata",
+ op.join("data", f['path']))
+ if f['id'] in have_keys:
+ raise KeyError("ID '{i}' occurs more than once".format(i=f['id']))
+ elif re.search("[^0-9A-Za-z_-]{1,}", f['id']):
+ raise ValueError("ID '{i}' contains non-alphanumeric characters".format(i=f['id']))
+ elif f['path'] in have_files:
+ raise ValueError("Path '{p}' occurs more than once".format(
+ p=f['path']))
+ elif not op.exists(file_path):
+ raise OSError("The path {p} does not exist".format(p=f['path']))
+ elif not "description" in f or not "filetype" in f:
+ raise ValueError(
+ "File {i} is missing description or filetype".format(i=f['id']))
+ have_keys.add(f['id'])
+ have_files.add(f['path'])
+ return 0
+
+
+def show_all(verbose=False):
+ global RAW_DICTS
+ d = _load_data()
+ by_type = defaultdict(list)
+ for f in RAW_DICTS:
+ by_type[f['filetype']].append(f)
+ filetypes = sorted(by_type.keys())
+ def _show_files_by_type(file_type):
+ file_ids = sorted([f['id'] for f in by_type[file_type]])
+ if not verbose:
+ print "{t}:".format(t=file_type)
+ print " {i}".format(i=", ".join(file_ids))
+ else:
+ print "{t}:".format(t=file_type)
+ for file_id in file_ids:
+ print " {i} - {d}".format(i=file_id,
+ d=d[file_id]['description'])
+ print " data/{p}".format(p=d[file_id]['path'])
+ # DataSet XML first
+ for file_type in sorted(by_type.keys()):
+ if file_type.startswith("PacBio.DataSet"):
+ _show_files_by_type(file_type)
+ # other files
+ for file_type in sorted(by_type.keys()):
+ if not file_type.startswith("PacBio.DataSet"):
+ _show_files_by_type(file_type)
+ return 0
+
+
+def main(argv=sys.argv):
+ p = argparse.ArgumentParser(description=__doc__, version=version.VERSION)
+ sp = p.add_subparsers(dest="command")
+ p1 = sp.add_parser("show")
+ p1.add_argument("--verbose", action="store_true")
+ p2 = sp.add_parser("get")
+ p2.add_argument("file_id")
+ p3 = sp.add_parser("validate")
+ args = p.parse_args(argv[1:])
+ if args.command == "show":
+ return show_all(args.verbose)
+ elif args.command == "validate":
+ return validate()
+ elif args.command == "get":
+ print get_file(args.file_id)
+ return 0
+ else:
+ raise ValueError("Unrecognized command '{c}'".format(c=args.command))
+
+if __name__ == "__main__":
+ sys.exit(main(sys.argv))
diff --git a/src/python/setup.py b/src/python/setup.py
new file mode 100644
index 0000000..2a9aba0
--- /dev/null
+++ b/src/python/setup.py
@@ -0,0 +1,42 @@
+
+import shutil
+import os.path
+import json
+import sys
+
+try:
+ from setuptools import setup, find_packages
+except ImportError:
+ from distutils.core import setup
+
+version = json.load(open(os.path.join("..", "..", "version.json")))["version"]
+with open(os.path.join("pbtestdata", "version.py"), "w") as version_py:
+ version_py.write("VERSION = \"{v}\"".format(v=version))
+
+DATA = os.path.join("pbtestdata", "data")
+
+if sys.argv[1] == "install":
+ if os.path.exists(DATA):
+ print "Removing old pbtestdata/data..."
+ shutil.rmtree(DATA)
+ print "Copying over pbtestdata/data..."
+ shutil.copytree(os.path.abspath(os.path.join("..", "..", "data")), DATA)
+package_data = []
+for dirname, dirnames, filenames in os.walk(DATA):
+ for file_name in filenames:
+ if not file_name.startswith("."):
+ path_name = os.path.relpath(dirname, "pbtestdata")
+ package_data.append(os.path.join(path_name, file_name))
+
+setup(
+ name="pbtestdata",
+ version=version,
+ license="BSD",
+ author="natechols",
+ url="https://github.com/PacificBiosciences/PacBioTestData",
+ keywords="pacbio".split(),
+ packages=find_packages(),
+ package_data={"pbtestdata": package_data},
+ include_package_data=True,
+ entry_points={"console_scripts": ["pbdata = pbtestdata.core:main"]},
+ zip_safe=False)
diff --git a/version.json b/version.json
new file mode 100644
index 0000000..93db032
--- /dev/null
+++ b/version.json
@@ -0,0 +1,3 @@
+{
+ "version": "0.0.13"
+}
--
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-med/pbtestdata.git
More information about the debian-med-commit
mailing list