Repository 'gstf_preparation'
hg clone https://toolshed.g2.bx.psu.edu/repos/earlhaminst/gstf_preparation

Changeset 6:56bbdbfe3eaa (2018-04-25)
Previous changeset 5:b3ba0c84667c (2018-04-16) Next changeset 7:9ef7661e8e9c (2018-04-25)
Commit message:
planemo upload for repository https://github.com/TGAC/earlham-galaxytools/tree/master/tools/gstf_preparation commit fa875eea77a9471acada2b7b8882a0467994c960
modified:
gstf_preparation.py
gstf_preparation.xml
schema/gstf.mwb
schema/gstf.png
schema/gstf.svg
test-data/test1.sqlite
added:
test-data/test1.ns.fasta
test-data/test4.fasta
test-data/test4.ns.fasta
test-data/test4.sqlite
test-data/test5.ns.fasta
test-data/test5.sqlite
test-data/test5_filtered.fasta
removed:
test-data/test2.sqlite
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa gstf_preparation.py
--- a/gstf_preparation.py Mon Apr 16 14:05:09 2018 -0400
+++ b/gstf_preparation.py Wed Apr 25 11:00:33 2018 -0400
[
@@ -6,7 +6,7 @@
 import sqlite3
 import sys
 
-version = "0.3.0"
+version = "0.4.0"
 gene_count = 0
 
 Sequence = collections.namedtuple('Sequence', ['header', 'sequence'])
@@ -41,6 +41,10 @@
     cur.execute('''CREATE TABLE gene (
         gene_id VARCHAR PRIMARY KEY NOT NULL,
         gene_symbol VARCHAR,
+        seq_region_name VARCHAR NOT NULL,
+        seq_region_start INTEGER NOT NULL,
+        seq_region_end INTEGER NOT NULL,
+        seq_region_strand INTEGER NOT NULL,
         species VARCHAR NOT NULL,
         gene_json VARCHAR NOT NULL)''')
     cur.execute('CREATE INDEX gene_symbol_index ON gene (gene_symbol)')
@@ -52,7 +56,7 @@
         gene_id VARCHAR NOT NULL REFERENCES gene(gene_id))''')
 
     cur.execute('''CREATE VIEW transcript_species AS
-        SELECT transcript_id, species
+        SELECT transcript_id, species, seq_region_name
         FROM transcript JOIN gene
         ON transcript.gene_id = gene.gene_id''')
 
@@ -225,8 +229,8 @@
             # This can happen when loading a JSON file from Ensembl
             continue
         gene_id = gene['id']
-        cur.execute('INSERT INTO gene (gene_id, gene_symbol, species, gene_json) VALUES (?, ?, ?, ?)',
-                    (gene_id, gene.get('display_name', None), gene['species'], json.dumps(gene)))
+        cur.execute('INSERT INTO gene (gene_id, gene_symbol, seq_region_name, seq_region_start, seq_region_end, seq_region_strand, species, gene_json) VALUES (?, ?, ?, ?, ?, ?, ?, ?)',
+                    (gene_id, gene.get('display_name', None), gene['seq_region_name'], gene['start'], gene['end'], gene['strand'], gene['species'], json.dumps(gene)))
 
         if "Transcript" in gene:
             for transcript in gene["Transcript"]:
@@ -241,15 +245,15 @@
     conn.commit()
 
 
-def fetch_species_for_transcript(conn, transcript_id):
+def fetch_species_and_seq_region_for_transcript(conn, transcript_id):
     cur = conn.cursor()
 
-    cur.execute('SELECT species FROM transcript_species WHERE transcript_id=?',
+    cur.execute('SELECT species, seq_region_name FROM transcript_species WHERE transcript_id=?',
                 (transcript_id, ))
     results = cur.fetchone()
     if not results:
         return None
-    return results[0]
+    return results
 
 
 def fetch_gene_id_for_transcript(conn, transcript_id):
@@ -280,8 +284,11 @@
     parser.add_option('--fasta', action='append', default=[], help='Path of the input FASTA files')
     parser.add_option('-l', action='store_true', default=False, dest='longestCDS', help='Keep only the longest CDS per gene')
     parser.add_option('--headers', action='store_true', default=False, help='Change the header line of the FASTA sequences to the >TranscriptId_species format')
+    parser.add_option('--regions', default="", help='Comma-separated list of region IDs for which FASTA sequences should be filtered')
     parser.add_option('-o', '--output', help='Path of the output SQLite file')
     parser.add_option('--of', help='Path of the output FASTA file')
+    parser.add_option('--ff', help='Path of the filtered sequences output FASTA file')
+
     options, args = parser.parse_args()
     if args:
         raise Exception('Use options to provide inputs')
@@ -368,14 +375,15 @@
         # first one to appear in the FASTA file is selected
         selected_transcript_ids = [max(transcript_id_lengths, key=lambda _: _[1])[0] for transcript_id_lengths in gene_transcripts_dict.values()]
 
-    with open(options.of, 'w') as output_fasta_file:
+    regions = [_.strip().lower() for _ in options.regions.split(",")]
+    with open(options.of, 'w') as output_fasta_file, open(options.ff, 'w') as filtered_fasta_file:
         for fasta_arg in options.fasta:
             for entry in FASTAReader_gen(fasta_arg):
                 transcript_id = remove_id_version(entry.header[1:].lstrip().split(' ')[0])
                 if options.longestCDS and transcript_id not in selected_transcript_ids:
                     continue
 
-                species_for_transcript = fetch_species_for_transcript(conn, transcript_id)
+                species_for_transcript, seq_region_for_transcript = fetch_species_and_seq_region_for_transcript(conn, transcript_id)
                 if not species_for_transcript:
                     print("Transcript '%s' in file '%s' not found in the gene feature information" % (transcript_id, fasta_arg), file=sys.stderr)
                     continue
@@ -387,7 +395,10 @@
                 else:
                     header = entry.header
 
-                output_fasta_file.write("%s\n%s\n" % (header, entry.sequence))
+                if seq_region_for_transcript.lower() in regions:
+                    filtered_fasta_file.write("%s\n%s\n" % (header, entry.sequence))
+                else:
+                    output_fasta_file.write("%s\n%s\n" % (header, entry.sequence))
 
     conn.close()
 
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa gstf_preparation.xml
--- a/gstf_preparation.xml Mon Apr 16 14:05:09 2018 -0400
+++ b/gstf_preparation.xml Wed Apr 25 11:00:33 2018 -0400
b
@@ -20,8 +20,12 @@
 #if $longestCDS
     -l
 #end if
+#if $regions
+    --regions '$regions'
+#end if
 -o '$output_db'
 --of '$output_fasta'
+--ff '$filtered_fasta'
 ]]>
     </command>
 
@@ -36,11 +40,13 @@
         <param name="fasta_inputs" type="data" format="fasta" multiple="true" label="Corresponding FASTA datasets" help="Each FASTA header line should start with a transcript id" />
         <param name="longestCDS" type="boolean" checked="false" label="Keep only the longest CDS per gene" />
         <param name="headers" type="boolean" checked="true" label="Change the header line of the FASTA sequences to the &gt;TranscriptId_species format" help="As required by TreeBest, part of the GeneSeqToFamily workflow" />
+        <param name="regions" type="text" optional="true" label="Comma-separated list of region IDs (e.g. chromosomes or scaffolds) for which FASTA sequences should be filtered" help="Region IDs are in the `seqid` column for GFF3 and in the `seq_region_name` field in JSON. This is typically used to filter chromosomes with a non-standard genetic code, like mitochondria, to be analysed separately" />
     </inputs>
 
     <outputs>
          <data name="output_db" format="sqlite" label="${tool.name} on ${on_string}: SQLite" />
          <data name="output_fasta" format="fasta" label="${tool.name} on ${on_string}: FASTA" />
+         <data name="filtered_fasta" format="fasta" label="${tool.name} on ${on_string}: filtered sequences" />
     </outputs>
 
     <tests>
@@ -53,6 +59,7 @@
 
             <output name="output_db" file="test1.sqlite" compare="sim_size" />
             <output name="output_fasta" file="test1.fasta" />
+            <output name="filtered_fasta" file="test1.ns.fasta" />
         </test>
         <test>
             <param name="fasta_inputs" ftype="fasta" value="Caenorhabditis_elegans.WBcel235.cds.all.shortened.fa" />
@@ -63,6 +70,7 @@
 
             <output name="output_db" file="test1.sqlite" compare="sim_size" />
             <output name="output_fasta" file="test1_longest.fasta" />
+            <output name="filtered_fasta" file="test1.ns.fasta" />
         </test>
         <test>
             <param name="fasta_inputs" ftype="fasta" value="Caenorhabditis_elegans.WBcel235.cds.all.shortened.fa" />
@@ -73,6 +81,7 @@
 
             <output name="output_db" file="test1.sqlite" compare="sim_size" />
             <output name="output_fasta" file="Caenorhabditis_elegans.WBcel235.cds.all.shortened.fa" />
+            <output name="filtered_fasta" file="test1.ns.fasta" />
         </test>
         <test>
             <param name="fasta_inputs" ftype="fasta" value="CDS.fasta" />
@@ -80,8 +89,20 @@
             <param name="longestCDS" value="false" />
             <param name="headers" value="true" />
 
-            <output name="output_db" file="test2.sqlite" compare="sim_size" />
-            <output name="output_fasta" file="test2.fasta" />
+            <output name="output_db" file="test4.sqlite" compare="sim_size" />
+            <output name="output_fasta" file="test4.fasta" />
+            <output name="filtered_fasta" file="test4.ns.fasta" />
+        </test>
+        <test>
+            <param name="fasta_inputs" ftype="fasta" value="CDS.fasta" />
+            <param name="json" ftype="json" value="gene.json" />
+            <param name="longestCDS" value="false" />
+            <param name="headers" value="true" />
+            <param name="regions" value="X" />
+
+            <output name="output_db" file="test5.sqlite" compare="sim_size" />
+            <output name="output_fasta" file="test5_filtered.fasta" />
+            <output name="filtered_fasta" file="test5.ns.fasta" />
         </test>
     </tests>
     <help>
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa schema/gstf.mwb
b
Binary file schema/gstf.mwb has changed
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa schema/gstf.png
b
Binary file schema/gstf.png has changed
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa schema/gstf.svg
--- a/schema/gstf.svg Mon Apr 16 14:05:09 2018 -0400
+++ b/schema/gstf.svg Wed Apr 25 11:00:33 2018 -0400
b
b'@@ -9,7 +9,7 @@\n <path style="stroke:none;" d="M 3.8125 -0.03125 L 3.8125 -3.984375 L 2.640625 -3.984375 L 2.640625 -3.265625 L 2.875 -3.375 C 2.578125 -3.8125 2.171875 -4.046875 1.75 -4.046875 C 0.890625 -4.046875 0.09375 -3.078125 0.09375 -1.953125 C 0.09375 -0.8125 0.828125 0.0625 1.734375 0.0625 C 2.15625 0.0625 2.515625 -0.109375 2.640625 -0.234375 L 2.640625 -0.03125 C 2.640625 0.40625 2.4375 0.546875 1.96875 0.546875 C 1.609375 0.546875 1.5 0.5625 1.375 0.125 L 0.15625 0.125 C 0.1875 0.890625 0.921875 1.484375 1.9375 1.484375 C 3.03125 1.484375 3.8125 0.828125 3.8125 -0.03125 Z M 2.65625 -1.96875 C 2.65625 -1.296875 2.4375 -1.046875 1.9375 -1.046875 C 1.5 -1.046875 1.3125 -1.296875 1.3125 -1.96875 C 1.3125 -2.65625 1.5 -2.9375 1.953125 -2.9375 C 2.4375 -2.9375 2.65625 -2.640625 2.65625 -1.96875 Z M 2.65625 -1.96875 "/>\n </symbol>\n <symbol overflow="visible" id="glyph0-2">\n-<path style="stroke:none;" d="M 3.703125 -1.859375 C 3.703125 -3.125 2.921875 -4.046875 1.84375 -4.046875 C 0.796875 -4.046875 0.015625 -3.171875 0.015625 -1.9375 C 0.015625 -0.765625 0.78125 0.0625 1.828125 0.0625 C 2.65625 0.0625 3.453125 -0.46875 3.703125 -1.359375 L 2.515625 -1.359375 C 2.375 -0.90625 2.203125 -0.984375 1.859375 -0.984375 C 1.40625 -0.984375 1.265625 -1.03125 1.234375 -1.578125 L 3.6875 -1.578125 Z M 2.578125 -2.53125 L 1.265625 -2.53125 C 1.3125 -2.875 1.40625 -3 1.828125 -3 C 2.25 -3 2.390625 -2.890625 2.421875 -2.53125 Z M 2.578125 -2.53125 "/>\n+<path style="stroke:none;" d="M 3.703125 -1.859375 C 3.703125 -3.125 2.921875 -4.046875 1.84375 -4.046875 C 0.796875 -4.046875 0.015625 -3.171875 0.015625 -1.9375 C 0.015625 -0.765625 0.78125 0.0625 1.828125 0.0625 C 2.65625 0.0625 3.453125 -0.46875 3.703125 -1.359375 L 2.5 -1.359375 C 2.375 -0.921875 2.203125 -0.984375 1.859375 -0.984375 C 1.40625 -0.984375 1.265625 -1.03125 1.234375 -1.578125 L 3.6875 -1.578125 Z M 2.578125 -2.53125 L 1.265625 -2.53125 C 1.3125 -2.875 1.40625 -3 1.828125 -3 C 2.25 -3 2.390625 -2.890625 2.421875 -2.53125 Z M 2.578125 -2.53125 "/>\n </symbol>\n <symbol overflow="visible" id="glyph0-3">\n <path style="stroke:none;" d="M 3.84375 -0.15625 L 3.84375 -2.609375 C 3.84375 -3.421875 3.265625 -4.046875 2.484375 -4.046875 C 1.984375 -4.046875 1.5625 -3.8125 1.28125 -3.40625 L 1.515625 -3.296875 L 1.515625 -3.984375 L 0.296875 -3.984375 L 0.296875 0 L 1.515625 0 L 1.515625 -2.359375 C 1.515625 -2.796875 1.6875 -2.921875 2.140625 -2.921875 C 2.5625 -2.921875 2.625 -2.859375 2.625 -2.421875 L 2.625 0 L 3.84375 0 Z M 3.84375 -0.15625 "/>\n@@ -18,16 +18,16 @@\n <path style="stroke:none;" d="M 2.1875 -0.15625 L 2.1875 -1 C 1.953125 -0.96875 1.890625 -0.96875 1.828125 -0.96875 C 1.578125 -0.96875 1.65625 -0.875 1.65625 -1.203125 L 1.65625 -2.96875 L 2.1875 -2.96875 L 2.1875 -3.921875 L 1.65625 -3.921875 L 1.65625 -4.90625 L 0.421875 -4.90625 L 0.421875 -3.921875 L -0.046875 -3.921875 L -0.046875 -2.96875 L 0.421875 -2.96875 L 0.421875 -0.9375 C 0.421875 -0.359375 0.875 0.03125 1.5 0.03125 C 1.703125 0.03125 1.875 0 2.1875 0 Z M 2.1875 -0.15625 "/>\n </symbol>\n <symbol overflow="visible" id="glyph0-5">\n-<path style="stroke:none;" d="M 2.65625 -2.921875 L 2.65625 -4.03125 C 2.453125 -4.046875 2.421875 -4.046875 2.40625 -4.046875 C 1.96875 -4.046875 1.453125 -3.578125 1.25 -3.078125 L 1.515625 -3.109375 L 1.515625 -3.984375 L 0.296875 -3.984375 L 0.296875 0 L 1.515625 0 L 1.515625 -2.109375 C 1.515625 -2.671875 1.65625 -2.78125 2.21875 -2.78125 C 2.3125 -2.78125 2.390625 -2.78125 2.65625 -2.734375 Z M 2.65625 -2.921875 "/>\n+<path style="stroke:none;" d="M 2.65625 -2.921875 L 2.65625 -4.03125 C 2.453125 -4.046875 2.421875 -4.046875 2.40625 -4.046875 C 1.96875 -4.046875 1.453125 -3.609375 1.25 -3.109375 L 1.515625 -3.109375 L 1.515625 -3.984375 L 0.296875 -3.984375 L 0.296875 0 L 1.515625 0 L 1.515625 -2.109375 C 1.515625 -2.671875 1.65625 -2.78125 2.21875 -2.78125 C 2.3125 -2.78125 2.390625 -2.78125 2.65625 -2.734375 Z M 2.65625 -2.921875 "/>\n </symbol>\n <symbol overflow="visibl'..b'(59.607843%,74.901961%,85.490196%);fill-opacity:1;stroke-width:1;stroke-linecap:butt;stroke-linejoin:miter;stroke:rgb(70%,70%,70%);stroke-opacity:1;stroke-miterlimit:10;" d="M 0.499837 8.501845 C 0.499837 4.078342 4.082737 0.502333 8.499349 0.502333 L 125.501682 0.502333 C 129.918294 0.502333 133.501194 4.078342 133.501194 8.501845 L 133.501194 26.499023 L 0.499837 26.499023 Z M 0.499837 8.501845 " transform="matrix(0.566929,0,0,0.566929,61.228346,174.047244)"/>\n-<use xlink:href="#image19033" transform="matrix(0.566929,0,0,0.566929,64.629921,176.88189)"/>\n+<use xlink:href="#image15744" transform="matrix(0.566929,0,0,0.566929,64.629921,176.88189)"/>\n <g style="fill:rgb(0%,0%,0%);fill-opacity:1;">\n   <use xlink:href="#glyph0-11" x="76.535156" y="184.199219"/>\n   <use xlink:href="#glyph0-2" x="82.577545" y="184.199219"/>\n@@ -547,11 +667,11 @@\n   <use xlink:href="#glyph0-6" x="88.619934" y="184.199219"/>\n </g>\n <path style=" stroke:none;fill-rule:nonzero;fill:rgb(40%,40%,40%);fill-opacity:1;" d="M 127.84375 179.433594 L 132.945312 179.433594 L 130.394531 183.96875 Z M 127.84375 179.433594 "/>\n-<use xlink:href="#image19034" transform="matrix(0.566929,0,0,0.566929,63.496063,192.472441)"/>\n+<use xlink:href="#image15745" transform="matrix(0.566929,0,0,0.566929,63.496063,192.472441)"/>\n <g style="fill:rgb(0%,0%,0%);fill-opacity:1;">\n-  <use xlink:href="#glyph1-43" x="70.867188" y="198.726562"/>\n+  <use xlink:href="#glyph1-45" x="70.867188" y="198.726562"/>\n   <use xlink:href="#glyph1-2" x="73.98439" y="198.726562"/>\n-  <use xlink:href="#glyph1-33" x="77.450714" y="198.726562"/>\n+  <use xlink:href="#glyph1-24" x="77.450714" y="198.726562"/>\n   <use xlink:href="#glyph1-17" x="79.526779" y="198.726562"/>\n   <use xlink:href="#glyph1-5" x="82.643982" y="198.726562"/>\n   <use xlink:href="#glyph1-21" x="84.028015" y="198.726562"/>\n@@ -565,7 +685,7 @@\n   <use xlink:href="#glyph1-9" x="114.514252" y="198.726562"/>\n   <use xlink:href="#glyph1-10" x="118.672607" y="198.726562"/>\n   <use xlink:href="#glyph1-13" x="123.173843" y="198.726562"/>\n-  <use xlink:href="#glyph1-25" x="125.249908" y="198.726562"/>\n+  <use xlink:href="#glyph1-26" x="125.249908" y="198.726562"/>\n   <use xlink:href="#glyph1-15" x="128.716232" y="198.726562"/>\n   <use xlink:href="#glyph1-16" x="132.182556" y="198.726562"/>\n </g>\n@@ -582,7 +702,7 @@\n <path style=" stroke:none;fill-rule:nonzero;fill:rgb(40%,40%,40%);fill-opacity:1;" d="M 127.84375 206.644531 L 132.378906 209.195312 L 127.84375 211.746094 Z M 127.84375 206.644531 "/>\n <path style="fill-rule:nonzero;fill:rgb(59.607843%,74.901961%,85.490196%);fill-opacity:1;stroke-width:1;stroke-linecap:butt;stroke-linejoin:miter;stroke:rgb(70%,70%,70%);stroke-opacity:1;stroke-miterlimit:10;" d="M 133.501194 48.502984 C 133.501194 52.919596 129.918294 56.502496 125.501682 56.502496 L 8.499349 56.502496 C 4.082737 56.502496 0.499837 52.919596 0.499837 48.502984 Z M 133.501194 48.502984 " transform="matrix(0.566929,0,0,0.566929,61.228346,188.787402)"/>\n <path style="fill-rule:nonzero;fill:rgb(99.607843%,87.058824%,34.509804%);fill-opacity:1;stroke-width:1;stroke-linecap:butt;stroke-linejoin:miter;stroke:rgb(70%,70%,70%);stroke-opacity:1;stroke-miterlimit:10;" d="M 0.499295 8.497179 C 0.499295 4.080566 4.082194 0.497667 8.498806 0.497667 L 132.501573 0.497667 C 136.918186 0.497667 140.501085 4.080566 140.501085 8.497179 L 140.501085 18.501736 C 140.501085 22.918349 136.918186 26.501248 132.501573 26.501248 L 8.498806 26.501248 C 4.082194 26.501248 0.499295 22.918349 0.499295 18.501736 Z M 0.499295 8.497179 " transform="matrix(0.566929,0,0,0.566929,218.267717,185.385827)"/>\n-<use xlink:href="#image19038" transform="matrix(0.566929,0,0,0.566929,221.669291,188.220472)"/>\n+<use xlink:href="#image15749" transform="matrix(0.566929,0,0,0.566929,221.669291,188.220472)"/>\n <g style="fill:rgb(0%,0%,0%);fill-opacity:1;">\n   <use xlink:href="#glyph0-4" x="233.574219" y="195.539062"/>\n   <use xlink:href="#glyph0-5" x="235.83757" y="195.539062"/>\n'
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa test-data/test1.sqlite
b
Binary file test-data/test1.sqlite has changed
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa test-data/test2.sqlite
b
Binary file test-data/test2.sqlite has changed
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa test-data/test4.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test4.fasta Wed Apr 25 11:00:33 2018 -0400
b
b'@@ -0,0 +1,1265 @@\n+>ENST00000338702_homosapiens\n+ATGGAGAATCAAGAGAAGGCGAGTATCGCGGGCCACATGTTCGACGTAGTCGTGATCGGA\n+GGTGGCATTTCAGGACTATCTGCTGCCAAACTCTTGACTGAATATGGCGTTAGTGTTTTG\n+GTTTTAGAAGCTCGGGACAGGGTTGGAGGAAGAACATATACTATAAGGAATGAGCATGTT\n+GATTACGTAGATGTTGGTGGAGCTTATGTGGGACCAACCCAAAACAGAATCTTACGCTTG\n+TCTAAGGAGCTGGGCATAGAGACTTACAAAGTGAATGTCAGTGAGCGTCTCGTTCAATAT\n+GTCAAGGGGAAAACATATCCATTTCGGGGCGCCTTTCCACCAGTATGGAATCCCATTGCA\n+TATTTGGATTACAATAATCTGTGGAGGACAATAGATAACATGGGGAAGGAGATTCCAACT\n+GATGCACCCTGGGAGGCTCAACATGCTGACAAATGGGACAAAATGACCATGAAAGAGCTC\n+ATTGACAAAATCTGCTGGACAAAGACTGCTAGGCGGTTTGCTTATCTTTTTGTGAATATC\n+AATGTGACCTCTGAGCCTCACGAAGTGTCTGCCCTGTGGTTCTTGTGGTATGTGAAGCAG\n+TGCGGGGGCACCACTCGGATATTCTCTGTCACCAATGGTGGCCAGGAACGGAAGTTTGTA\n+GGTGGATCTGGTCAAGTGAGCGAACGGATAATGGACCTCCTCGGAGACCAAGTGAAGCTG\n+AACCATCCTGTCACTCACGTTGACCAGTCAAGTGACAACATCATCATAGAGACGCTGAAC\n+CATGAACATTATGAGTGCAAATACGTAATTAATGCGATCCCTCCGACCTTGACTGCCAAG\n+ATTCACTTCAGACCAGAGCTTCCAGCAGAGAGAAACCAGTTAATTCAGCGGCTTCCAATG\n+GGAGCTGTCATTAAGTGCATGATGTATTACAAGGAGGCCTTCTGGAAGAAGAAGGATTAC\n+TGTGGCTGCATGATCATTGAAGATGAAGATGCTCCAATTTCAATAACCTTGGATGACACC\n+AAGCCAGATGGGTCACTGCCTGCCATCATGGGCTTCATTCTTGCCCGGAAAGCTGATCGA\n+CTTGCTAAGCTACATAAGGAAATAAGGAAGAAGAAAATCTGTGAGCTCTATGCCAAAGTG\n+CTGGGATCCCAAGAAGCTTTACATCCAGTGCATTATGAAGAGAAGAACTGGTGTGAGGAG\n+CAGTACTCTGGGGGCTGCTACACGGCCTACTTCCCTCCTGGGATCATGACTCAATATGGA\n+AGGGTGATTCGTCAACCCGTGGGCAGGATTTTCTTTGCGGGCACAGAGACTGCCACAAAG\n+TGGAGCGGCTACATGGAAGGGGCAGTTGAGGCTGGAGAACGAGCAGCTAGGGAGGTCTTA\n+AATGGTCTCGGGAAGGTGACCGAGAAAGATATCTGGGTACAAGAACCTGAATCAAAGGAC\n+GTTCCAGCGGTAGAAATCACCCACACCTTCTGGGAAAGGAACCTGCCCTCTGTTTCTGGC\n+CTGCTGAAGATCATTGGATTTTCCACATCAGTAACTGCCCTGGGGTTTGTGCTGTACAAA\n+TACAAGCTCCTGCCACGGTCTTGA\n+>ENST00000542639_homosapiens\n+ATGGGGAAGGAGATTCCAACTGATGCACCCTGGGAGGCTCAACATGCTGACAAATGGGAC\n+AAAATGACCATGAAAGAGCTCATTGACAAAATCTGCTGGACAAAGACTGCTAGGCGGTTT\n+GCTTATCTTTTTGTGAATATCAATGTGACCTCTGAGCCTCACGAAGTGTCTGCCCTGTGG\n+TTCTTGTGGTATGTGAAGCAGTGCGGGGGCACCACTCGGATATTCTCTGTCACCAATGGT\n+GGCCAGGAACGGAAGTTTGTAGGTGGATCTGGTCAAGTGAGCGAACGGATAATGGACCTC\n+CTCGGAGACCAAGTGAAGCTGAACCATCCTGTCACTCACGTTGACCAGTCAAGTGACAAC\n+ATCATCATAGAGACGCTGAACCATGAACATTATGAGTGCAAATACGTAATTAATGCGATC\n+CCTCCGACCTTGACTGCCAAGATTCACTTCAGACCAGAGCTTCCAGCAGAGAGAAACCAG\n+TTAATTCAGCGGCTTCCAATGGGAGCTGTCATTAAGTGCATGATGTATTACAAGGAGGCC\n+TTCTGGAAGAAGAAGGATTACTGTGGCTGCATGATCATTGAAGATGAAGATGCTCCAATT\n+TCAATAACCTTGGATGACACCAAGCCAGATGGGTCACTGCCTGCCATCATGGGCTTCATT\n+CTTGCCCGGAAAGCTGATCGACTTGCTAAGCTACATAAGGAAATAAGGAAGAAGAAAATC\n+TGTGAGCTCTATGCCAAAGTGCTGGGATCCCAAGAAGCTTTACATCCAGTGCATTATGAA\n+GAGAAGAACTGGTGTGAGGAGCAGTACTCTGGGGGCTGCTACACGGCCTACTTCCCTCCT\n+GGGATCATGACTCAATATGGAAGGGTGATTCGTCAACCCGTGGGCAGGATTTTCTTTGCG\n+GGCACAGAGACTGCCACAAAGTGGAGCGGCTACATGGAAGGGGCAGTTGAGGCTGGAGAA\n+CGAGCAGCTAGGGAGGTCTTAAATGGTCTCGGGAAGGTGACCGAGAAAGATATCTGGGTA\n+CAAGAACCTGAATCAAAGGACGTTCCAGCGGTAGAAATCACCCACACCTTCTGGGAAAGG\n+AACCTGCCCTCTGTTTCTGGCCTGCTGAAGATCATTGGATTTTCCACATCAGTAACTGCC\n+CTGGGGTTTGTGCTGTACAAATACAAGCTCCTGCCACGGTCTTGA\n+>ENSPTRT00000040520_pantroglodytes\n+ATGGAGAATCAAGAGAAGGCGAGTATCGCGGGCCACATGTTCGACGTAGTCGTGATCGGA\n+GGTGGCATTTCAGGACTATCTGCTGCCAAACTCTTGACTGAATATGGCGTTAGTGTTTTA\n+GTTTTAGAAGCTCGGGACAGGGTTGGAGGAAGAACATATACTATAAGGAATGAGCATGTT\n+GATTACGTAGATGTTGGTGGAGCTTATGTGGGACCAACCCAAAACAGAATCTTACGCTTG\n+TCTAAGGAGCTGGGCATAGAGACTTACAAAGTGAATGTCAGTGAGCGTCTCGTTCAATAT\n+GTCAAGGGGAAAACATATCCATTTCGGGGCGCCTTTCCACCAGTATGGAATCCCATTGCA\n+TATTTGGATTACAATAATCTGTGGCGGACAATAGATAACATGGGGAAGGAGATTCCAAAT\n+GATGCACCCTGGGAGGCTCAACATGCTGACGAATGGGACAAAATGACCATGAAAGAGCTC\n+ATTGACAAAATCTGCTGGACAAAGACTGCTAGGCGGTTTGCTTATCTTTTTGTGAATATC\n+AATGTGACCTCTGAGCCTCACGAAGTGTCTGCCCTGTGGTTCTTGTGGTATGTGAAGCAG\n+TGCGGGGGCACCACTCGGATATTCTCTGTCACCAATGGCGGCCAGGAACGGAAGTTTGTA\n+GGTGGATCTGGTCAAGTGAGCGAACGGATAATGGACCTCCTTGGAGACCAAGTGAAGCTG\n+AACCATCCTGTCACTCATGTTGACCAGTCAAGTGACAACATCATCATAGAGACACTGAAC\n+CATGAACATTATGAGTGCAAATACGTAATTAATGCGATCCCTCCGACCTTGACTGCCAAG\n+ATTCACTTCAGACCAGAGCTTCCAGCAGAGAGAAACCAGTTAATTCAGCGTCTTCCAATG\n+GGAGCTATCATTAAGTGCATGATGTATTACAAGGAGGCCTTCTGGAAGAAGAAGGATTAC\n+TGTGGCTGCATGATCATTGAAGATGAAG'..b'GCACAATCAG\n+AGTGAGTATGACGACTCGGCCAGCGAGTGCTGCTCATGTCCTAAGACTGACTCTCAGATC\n+CTGAAGGAGCTGGAGGAGTCTTCATTCAGGAAGACCTTCGAGGATTACCTGCACAACGTG\n+GTTTTTGTTCCCAGAAAAACCTCTTCAGGCAATGGTGCTGAGGACACTAGGCCATCCCGA\n+AAGCGAAGATCCCTTGAAGAGGTGGGCAATGTGACAGCCACTACACCCACACTTCCAGAT\n+TTTCCCAACATCTCCTCCACCATCGCGCCCACAAGCCACGAAGAGCACAGACCATTTGAG\n+AAAGTAGTAAACAAGGAGTCACTTGTCATCTCTGGCCTGAGACACTTCACTGGGTACCGC\n+ATTGAGCTGCAGGCATGCAATCAGGACTCCCCAGAAGAGAGGTGCAGCGTGGCTGCCTAC\n+GTCAGTGCCCGGACCATGCCTGAAGCTAAGGCAGATGACATCGTTGGCCCTGTGACCCAT\n+GAAATCTTTGAGAACAATGTTGTACACTTAATGTGGCAAGAGCCAAAGGAACCTAATGGT\n+CTGATTGTGCTATATGAAGTGAGCTATCGGCGATATGGTGATGAGGAGCTGCACCTCTGT\n+GTCTCCCGGAAGCATTTTGCCCTGGAGCGGGGCTGCAGGCTTCGAGGGCTCTCTCCAGGA\n+AACTACAGTGTTCGAGTCCGGGCTACCTCTCTGGCAGGAAATGGCTCCTGGACAGAACCC\n+ACCTATTTTTATGTGACTGATTATTTAGATGTCCCATCAAATATTGCCAAAATTATCATC\n+GGGCCCCTCATCTTCGTCTTCCTCTTCAGTGTCGTGATCGGAAGTATTTATCTATTCTTG\n+AGGAAGAGGCAGCCAGATGGGCCAATGGGACCACTGTACGCTTCTTCAAACCCAGAGTAC\n+CTCAGTGCCAGTGATGTCTTTCCATCTTCCGTATACGTTCCGGATGAGTGGGAGGTACCT\n+CGAGAGAAGATCACCCTCCTCCGAGAGCTGGGGCAGGGATCCTTCGGTATGGTGTACGAA\n+GGCAATGCCAAGGATATCATCAAGGGTGAGGTAGAGACCCGTGTTGCGGTGAAGACGGTC\n+AATGAGTCAGCCAGTCTTCGAGAACGGATCGAGTTCCTCAATGAGGCATCAGTCATGAAG\n+GGCTTCACCTGTCATCACGTGGTCCGCCTTCTTGGGGTGGTGTCCAAAGGCCAGCCCACA\n+TTGGTAGTGATGGAACTGATGGCTCATGGAGACCTGAAAAGTCACCTCCGTTCTCTGCGG\n+CCCGATGCTGAGAACAACCCAGGCCGTCCTCCCCCTACCTTGCAAGAAATGATTCAGATG\n+ACAGCAGAAATTGCCGATGGCATGGCATACTTGAACGCCAAGAAGTTTGTGCACCGGGAC\n+CTGGCAGCTCGGAACTGCATGGTTGCCCATGATTTTACTGTCAAAATCGGAGACTTTGGA\n+ATGACGAGAGACATCTACGAGACAGATTACTATCGGAAAGGGGGCAAGGGGTTGCTGCCC\n+GTGAGGTGGATGTCACCCGAGTCCCTGAAGGACGGAGTCTTCACTGCTTCTTCCGACATG\n+TGGTCCTTTGGGGTGGTCCTTTGGGAAATCACCAGCCTGGCTGAGCAACCTTACCAAGGC\n+CTGTCTAATGAACAGGTGTTGAAATTTGTCATGGATGGAGGCTATCTGGATCCCCCTGAT\n+AACTGTCCAGAGAGACTCACTGACCTGATGCGCATGTGCTGGCAGTTCAACCCCAAGATG\n+AGGCCGACCTTCCTGGAAATCGTCAACCTGCTCAAGGACGACCTCCACCCCAGCTTTCCG\n+GAAGTTTCCTTCTTCTACAGCGAGGAGAACAAGGCTCCCGAGAGTGAAGAGCTGGAGATG\n+GAGTTCGAGGACATGGAGAATGTCCCCTTGGATCGTTCCTCTCACTGTCAGAGAGAAGAG\n+GCTGGATGCCGGGAGGGAGGGTCCTCTCTGAGCATCAAACGGACCTATGATGAACACATC\n+CCCTACACCCACATGAACGGGGGCAAGAAGAATGGGCGGGTCCTCACCCTGCCGAGGTCG\n+AACCCTTCCTAA\n+>ENSSSCT00000014817_susscrofa\n+GTGTGCCCAGGGATGGATATCCGGAATAACCTTACACGGCTGCACGAGTTGGCCAACTGC\n+TCGGTCATCGAAGGACATTTGCAGATCCTGTTGATGTTCAAAACGCGGCCCGAGGATTTC\n+CGAGACCTCAGTTTCCCCAAACTCATCATGATCACTGATTACTTGCTGCTCTTCCGGGTC\n+TACGGGCTGGAGAGCCTGAAGGACCTGTTCCCCAACCTCACCGTCATCCGGGGCTCACGC\n+CTCTTCTTTAACTATGCGCTGGTCATCTTTGAGATGGTTCACTTGAAGGAGCTTGGCCTC\n+TACAATTTGATGAACATCACCAGGGGTGCTGTCCGCATCGAGAAGAACAATGAGCTCTGC\n+TACCTGGCGACCATTGACTGGTCGCGCATCCTGGACTCTGTGGAGGATAATTACATTGTG\n+CTGAACAAAGACGACAACGAGGAGTGTGGGGACATTTGCCCAGGCACTGCGAAGGGCAAG\n+ACCAATTGCCCTGCCACCGTCATCAATGGGCAATTTGTCGAGCGGTGTTGGACGCACAGT\n+CACTGCCAGAAAGTGTGCCCGACCATCTGTAAGTCGCACGGCTGCACTGCTGAGGGCCTC\n+TGCTGTCACAGCGAGTGTTTGGGCAACTGCTCTGAGCCAGACGACCCCACCAAGTGCGTG\n+GCCTGCCGCAACTTCTACCTGGACGGCAGATGCGTGGAGACCTGCCCGCCCCCCTACTAC\n+CACTTCCAAGACTGGCGCTGCGTGAACTTCAGCTTCTGCCAGGACCTGCACAACAAATGC\n+AAGAACTCAAGGAGGCAGGGCTGCCACCAGTACGTCATTCACAACAACAAGTGTATCCCT\n+GAGTGCCCCTCAGGGTACACGATGAATTCCAGCAACTTGATGTGCACTCCGTGCCTAGGC\n+CCCTGTCCCAAAGTGTGTCACCTCCTGGAAGGCGAGAAGACCATCGACTCAGTGACATCC\n+GCCCAGGAGCTCCGAGGCTGCACCATTATCAACGGGAGCCTAATCATCAACATTCGAGGA\n+GGCAACAACCTGGCAGCCGAACTAGAGGCCAACCTTGGACTCATTGAGGAGATTTCAGGG\n+TACCTGAAAATCCGCCGATCCTATGCCCTCGTGTCACTTTCCTTCTTCCGGAAGTTGCGT\n+CTGATCCGAGGGGAGACGTTGGAAATTGGGAACTATTCTTTCTATGCCTTGGACAACCAG\n+AACCTAAGGCAACTGTGGGACTGGAGCAAACACAACCTCACCATCACTCAGGGGAAACTC\n+TTCTTCCATTATAATCCCAAACTCTGCTTGTCGGAAATTCACAAGATGGAGGAAGTTTCT\n+GGAACCAAGGGGCGCCAGGAGAGAAATGATATTGCCCTGAAGACCAATGGGGACCAGGCG\n+TCCTGTGAAAATGAGTTACTTAAATTTTCTTACATTCGGACATCTTATGACAAGATCTTG\n+CTGAAGTGGGAGCCGTATTGGCCCCCCGACTTCCGAGACCTCCTGGGGTTCATGCTCTTC\n+TACAAAGAGGCCCCTTATCAGAACGTGACGGAGTTTGACGGGCAGGATGCGTGTGGCTCC\n+AACAGCTGGACGGTGGTGGACATTGACCCGCCTACGAGGTCCAATGACCCCAAGTCCCAG\n+AACCATCCTGGGTGGCTGATGCGTGGTCTCAAGCCCTGGACCCAGTATGCCATCTTTGTC\n+AAGACTTTGGTCACCTTTTCTGATGAACGACGCACCTATGGAGCCAAGAGTGACATCATC\n+TACGTCCAGACAGATGCCACAAGTAAGCATGTC\n+\n'
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa test-data/test4.sqlite
b
Binary file test-data/test4.sqlite has changed
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa test-data/test5.ns.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test5.ns.fasta Wed Apr 25 11:00:33 2018 -0400
b
b'@@ -0,0 +1,423 @@\n+>ENST00000338702_homosapiens\n+ATGGAGAATCAAGAGAAGGCGAGTATCGCGGGCCACATGTTCGACGTAGTCGTGATCGGA\n+GGTGGCATTTCAGGACTATCTGCTGCCAAACTCTTGACTGAATATGGCGTTAGTGTTTTG\n+GTTTTAGAAGCTCGGGACAGGGTTGGAGGAAGAACATATACTATAAGGAATGAGCATGTT\n+GATTACGTAGATGTTGGTGGAGCTTATGTGGGACCAACCCAAAACAGAATCTTACGCTTG\n+TCTAAGGAGCTGGGCATAGAGACTTACAAAGTGAATGTCAGTGAGCGTCTCGTTCAATAT\n+GTCAAGGGGAAAACATATCCATTTCGGGGCGCCTTTCCACCAGTATGGAATCCCATTGCA\n+TATTTGGATTACAATAATCTGTGGAGGACAATAGATAACATGGGGAAGGAGATTCCAACT\n+GATGCACCCTGGGAGGCTCAACATGCTGACAAATGGGACAAAATGACCATGAAAGAGCTC\n+ATTGACAAAATCTGCTGGACAAAGACTGCTAGGCGGTTTGCTTATCTTTTTGTGAATATC\n+AATGTGACCTCTGAGCCTCACGAAGTGTCTGCCCTGTGGTTCTTGTGGTATGTGAAGCAG\n+TGCGGGGGCACCACTCGGATATTCTCTGTCACCAATGGTGGCCAGGAACGGAAGTTTGTA\n+GGTGGATCTGGTCAAGTGAGCGAACGGATAATGGACCTCCTCGGAGACCAAGTGAAGCTG\n+AACCATCCTGTCACTCACGTTGACCAGTCAAGTGACAACATCATCATAGAGACGCTGAAC\n+CATGAACATTATGAGTGCAAATACGTAATTAATGCGATCCCTCCGACCTTGACTGCCAAG\n+ATTCACTTCAGACCAGAGCTTCCAGCAGAGAGAAACCAGTTAATTCAGCGGCTTCCAATG\n+GGAGCTGTCATTAAGTGCATGATGTATTACAAGGAGGCCTTCTGGAAGAAGAAGGATTAC\n+TGTGGCTGCATGATCATTGAAGATGAAGATGCTCCAATTTCAATAACCTTGGATGACACC\n+AAGCCAGATGGGTCACTGCCTGCCATCATGGGCTTCATTCTTGCCCGGAAAGCTGATCGA\n+CTTGCTAAGCTACATAAGGAAATAAGGAAGAAGAAAATCTGTGAGCTCTATGCCAAAGTG\n+CTGGGATCCCAAGAAGCTTTACATCCAGTGCATTATGAAGAGAAGAACTGGTGTGAGGAG\n+CAGTACTCTGGGGGCTGCTACACGGCCTACTTCCCTCCTGGGATCATGACTCAATATGGA\n+AGGGTGATTCGTCAACCCGTGGGCAGGATTTTCTTTGCGGGCACAGAGACTGCCACAAAG\n+TGGAGCGGCTACATGGAAGGGGCAGTTGAGGCTGGAGAACGAGCAGCTAGGGAGGTCTTA\n+AATGGTCTCGGGAAGGTGACCGAGAAAGATATCTGGGTACAAGAACCTGAATCAAAGGAC\n+GTTCCAGCGGTAGAAATCACCCACACCTTCTGGGAAAGGAACCTGCCCTCTGTTTCTGGC\n+CTGCTGAAGATCATTGGATTTTCCACATCAGTAACTGCCCTGGGGTTTGTGCTGTACAAA\n+TACAAGCTCCTGCCACGGTCTTGA\n+>ENST00000542639_homosapiens\n+ATGGGGAAGGAGATTCCAACTGATGCACCCTGGGAGGCTCAACATGCTGACAAATGGGAC\n+AAAATGACCATGAAAGAGCTCATTGACAAAATCTGCTGGACAAAGACTGCTAGGCGGTTT\n+GCTTATCTTTTTGTGAATATCAATGTGACCTCTGAGCCTCACGAAGTGTCTGCCCTGTGG\n+TTCTTGTGGTATGTGAAGCAGTGCGGGGGCACCACTCGGATATTCTCTGTCACCAATGGT\n+GGCCAGGAACGGAAGTTTGTAGGTGGATCTGGTCAAGTGAGCGAACGGATAATGGACCTC\n+CTCGGAGACCAAGTGAAGCTGAACCATCCTGTCACTCACGTTGACCAGTCAAGTGACAAC\n+ATCATCATAGAGACGCTGAACCATGAACATTATGAGTGCAAATACGTAATTAATGCGATC\n+CCTCCGACCTTGACTGCCAAGATTCACTTCAGACCAGAGCTTCCAGCAGAGAGAAACCAG\n+TTAATTCAGCGGCTTCCAATGGGAGCTGTCATTAAGTGCATGATGTATTACAAGGAGGCC\n+TTCTGGAAGAAGAAGGATTACTGTGGCTGCATGATCATTGAAGATGAAGATGCTCCAATT\n+TCAATAACCTTGGATGACACCAAGCCAGATGGGTCACTGCCTGCCATCATGGGCTTCATT\n+CTTGCCCGGAAAGCTGATCGACTTGCTAAGCTACATAAGGAAATAAGGAAGAAGAAAATC\n+TGTGAGCTCTATGCCAAAGTGCTGGGATCCCAAGAAGCTTTACATCCAGTGCATTATGAA\n+GAGAAGAACTGGTGTGAGGAGCAGTACTCTGGGGGCTGCTACACGGCCTACTTCCCTCCT\n+GGGATCATGACTCAATATGGAAGGGTGATTCGTCAACCCGTGGGCAGGATTTTCTTTGCG\n+GGCACAGAGACTGCCACAAAGTGGAGCGGCTACATGGAAGGGGCAGTTGAGGCTGGAGAA\n+CGAGCAGCTAGGGAGGTCTTAAATGGTCTCGGGAAGGTGACCGAGAAAGATATCTGGGTA\n+CAAGAACCTGAATCAAAGGACGTTCCAGCGGTAGAAATCACCCACACCTTCTGGGAAAGG\n+AACCTGCCCTCTGTTTCTGGCCTGCTGAAGATCATTGGATTTTCCACATCAGTAACTGCC\n+CTGGGGTTTGTGCTGTACAAATACAAGCTCCTGCCACGGTCTTGA\n+>ENSPTRT00000040520_pantroglodytes\n+ATGGAGAATCAAGAGAAGGCGAGTATCGCGGGCCACATGTTCGACGTAGTCGTGATCGGA\n+GGTGGCATTTCAGGACTATCTGCTGCCAAACTCTTGACTGAATATGGCGTTAGTGTTTTA\n+GTTTTAGAAGCTCGGGACAGGGTTGGAGGAAGAACATATACTATAAGGAATGAGCATGTT\n+GATTACGTAGATGTTGGTGGAGCTTATGTGGGACCAACCCAAAACAGAATCTTACGCTTG\n+TCTAAGGAGCTGGGCATAGAGACTTACAAAGTGAATGTCAGTGAGCGTCTCGTTCAATAT\n+GTCAAGGGGAAAACATATCCATTTCGGGGCGCCTTTCCACCAGTATGGAATCCCATTGCA\n+TATTTGGATTACAATAATCTGTGGCGGACAATAGATAACATGGGGAAGGAGATTCCAAAT\n+GATGCACCCTGGGAGGCTCAACATGCTGACGAATGGGACAAAATGACCATGAAAGAGCTC\n+ATTGACAAAATCTGCTGGACAAAGACTGCTAGGCGGTTTGCTTATCTTTTTGTGAATATC\n+AATGTGACCTCTGAGCCTCACGAAGTGTCTGCCCTGTGGTTCTTGTGGTATGTGAAGCAG\n+TGCGGGGGCACCACTCGGATATTCTCTGTCACCAATGGCGGCCAGGAACGGAAGTTTGTA\n+GGTGGATCTGGTCAAGTGAGCGAACGGATAATGGACCTCCTTGGAGACCAAGTGAAGCTG\n+AACCATCCTGTCACTCATGTTGACCAGTCAAGTGACAACATCATCATAGAGACACTGAAC\n+CATGAACATTATGAGTGCAAATACGTAATTAATGCGATCCCTCCGACCTTGACTGCCAAG\n+ATTCACTTCAGACCAGAGCTTCCAGCAGAGAGAAACCAGTTAATTCAGCGTCTTCCAATG\n+GGAGCTATCATTAAGTGCATGATGTATTACAAGGAGGCCTTCTGGAAGAAGAAGGATTAC\n+TGTGGCTGCATGATCATTGAAGATGAAGA'..b'TGAGAAACCAGCTCATCACTCGTGTACCTCTGGGCTCTGTCATCAAG\n+TGTATAGTTTATTACAAAGAGCCCTTCTGGAGGCATAAGGATTACTGTGGAAGCATGATT\n+ATTGAAGGAGAGGAAGCTCCAATCGCCTACACGTTGGATGATTCCAAGCCTGATGGCAGC\n+TGTGCCGCCATCATAGGATTTATCCTTGCCCACAAAGCCAGAAAACTGGCCCGTCTTACC\n+AAAGAAGAAAGGCTGAAGAAACTTTGCGACCTCTATGCAAAAGTTCTGGGTTCAAAAGAA\n+GCTTTGAACCCCGTGCACTATGAAGAGAAGAACTGGTGCGAGGAGCAGTACTCGGCGGGC\n+TGCTACACGACCTACTTCCCCCCTGGGATCATGACTCAGTATGGAAGGGTTCTACGCCAG\n+CCAGTCGGCAGGATTTATTTCGCCGGCACGGAGACTGCCACGCACTGGAGTGGCTACATG\n+GAGGGGGCCGTGGAGGCCGGAGAGAGAGCGGCCCGAGAGATCCTGCATGCTATGGGAAAG\n+ATCCCAGAAGATGAAATCTGGCAGTCTGAACCAGAGTCCGTGGATGTGCCTGCGAAGCCC\n+ATTACCACGACCTTCTTGGAGAGACACTTGCCCTCGGTGCCCGGCCTGCTGAGGCTGATT\n+GGATTGACCGCCATCTTTTCAGCCACTGCTCTCGGCTACCTGGCCCACAAAAGGGGGCTA\n+CTCGTGCGGGTCTGA\n+>ENSSSCT00000023183_susscrofa\n+ATGGCAGCGGCCAAACTTCTGCATGACTCTGGCCTGAGTGTGATTGTTCTGGAAGCCCGG\n+GACCGCGTGGGAGGCAGGACTTACACCGTCAGGAACCAACAAGTTAAATATGTGGACCTT\n+GGAGGATCTTATGTTGGGCCAACTCAGAATCGCATCTTAAGATTGTCCAAGGAGCTAGGA\n+TTGGAGACCTACAAAGTGAATGAAGTGGAGCGTCTGATTCACTATATCAAAGGCAAATCC\n+TACCCCTTCAGGGGCCCATTACCACCTGTGAGGAATCCGATTACCTTCCTAGATCTTAAC\n+AACCTTTGGAGGACGGTGGATGACATGGGACGAGAGATTCCCAGTGATGCCCCATGGAAG\n+GCGCCCCTTGCAGAACAGTGGGACCAGATGACAATGAAGGAGCTGTTGGACAAGCTCTGC\n+TGGACTGAATCTTCGAAGCAGCTGGCCACCCTTTTTGTGAACCTGTGTGTCACCGCGGAG\n+ACCCATGAGGTCTCTGCTCTCTGGTTCCTGTGGTATGTGAAGCAGTGTGGAGGCACCACC\n+AGGATCATCTCAACAACTAACGGAGGGCAGGAGAGGAAATTTGTGGGCGGATCTGGTCAA\n+GTGACCGAGCGGATAAAGGACCTCCTTGGAGACCGAGTGAAGCTGGAGAGGCCTGTGGTC\n+CACATTGACCAGACAGGAGAAAATGTCCTCGTGGAGACCCTAAACCACGAGGTGTACGAG\n+GCTAAGTATGTGATTAGCGCCATTCCTCCTGTCCTGGGCATGAAGATTCATTTCAGTCCC\n+CCTCTGCCAATGATGAGAAACCAGCTCATCACTCGTGTACCTCTGGGCTCTGTCATCAAG\n+TGTATAGTTTATTACAAAGAGCCCTTCTGGAGGCATAAGGATTACTGTGGAAGCATGATT\n+ATTGAAGGAGAGGAAGCTCCAATCGCCTACACGTTGGATGATTCCAAGCCTGATGGCAGC\n+TGTGCCGCCATCATAGGATTTATCCTTGCCCACAAAGCCAGAAAACTGGCCCGTCTTACC\n+AAAGAAGAAAGGCTGAAGAAACTTTGCGACCTCTATGCAAAAGTTCTGGGTTCAAAAGAA\n+GCTTTGAACCCCGTGCACTATGAAGAGAAGAACTGGTGCGAGGAGCAGTACTCGGCGGGC\n+TGCTACACGACCTACTTCCCCCCTGGGATCATGACTCAGTATGGAAGGGTTCTACGCCAG\n+CCAGTCGGCAGGATTTATTTCGCCGGCACGGAGACTGCCACGCACTGGAGTGGCTACATG\n+GAGGGGGCCGTGGAGGCCGGAGAGAGAGCGGCCCGAGAGATCCTGCATGCTATGGGAAAG\n+ATCCCAGAAGATGAAATCTGGCAGTCTGAACCAGAGTCCGTGGATGTGCCTGCGAAGCCC\n+ATTACCACGACCTTCTTGGAGAGACACTTGCCCTCGGTGCCCGGCCTGCTGAGGCTGATT\n+GGATTGACCGCCATCTTTTCAGCCACTGCTCTCGGCTACCTGGCCCACAAAAGGGGGCTA\n+CTCGTGCGGGTCTGA\n+>ENSCAFT00000022963_canisfamiliaris\n+ATGAGCGGCAAGTGCGACGTGGTCATGGTGGGGGGCGGCATCTCAGGCATGGCAGCAGCC\n+AAACTTCTGCATGATTTTGGCCTGAACGTGGTTGTTCTGGAGGCCCGGGACCGAGTGGGA\n+GGCAGGACTTACACCATCAGGAACCAAAAGGTTAAATATTTGGACCTTGGAGGATCTTAT\n+GTTGGGCCAACTCAGAATTGTATCTTAAGGTTAGCCAAGGAGCTAGGATTGGAGACCTAC\n+AAAGTGAATGAAGTAGAGCGTCTCATCCATCATGTAAAGGGCAAATCGTACCCCTTCAGG\n+GGCCCCTTCCCACCTGTATGGAACCCAATTGCATACCTAGATCATAACAACCTCTGGAGG\n+ACGATGGACGACATGGGGCGAGAGATTCCCAGTGATGCCCCATGGAAGGCACCTCTTGCG\n+GAGGAGTGGGACCACATGACAATGAAGGAGCTTCTGGACAAGATCTGCTGGACAGAATCT\n+GCCAAGCAACTTGCTACTCTCTTTGTGAATCTGTGTGTCACTGCAGAGACCCATGAGGTC\n+TCTGCTCTCTGGTTCCTGTGGTATGTGAAGCAGTGTGGAGGCACGACCAGGATCATCTCA\n+ACAACCAATGGAGGGCAGGAGAGAAAATTTGTGGGCGGATCTGGTCAAGTGAGTGAACGG\n+ATAATGGACCTCCTTGGAGACCAAGTGAAGCTGGAGAGGCCTGTGACCCACATTGACCAG\n+ACAGGAGAAAATGTCCTTGTGGAGACCCTAAACCATGAGGTGTATGAGGCTAAGTATGTG\n+ATTAGTGCCATTCCTCCTACTCTGGGCATGAAGATTCATTTCAACCCCCCTCTGCCAATG\n+ATGAGAAACCAGCTGATCACTCGTGTGCCTTTGGGCTCAGTCATCAAGTGCATAGTTTAT\n+TATAAAGAGCCCTTCTGGAGGAAAAAGGATTACTGTGGAACCATGATTATCGAAGGAGAG\n+GAAGCTCCAATTGCTTACACGTTGGATGATACCAAACCTGATGGCAACTATGCCGCCATA\n+ATGGGATTTATCCTTGCCCACAAAGCCAGAAAACTGGCACGTCTTACCAAAGATGAAAGG\n+ATGAAGAAACTTTGTGAGCTCTATGCAAAAGTCCTGGGTTCCCAAGAAGCTTTGCAGCCC\n+GTACACTACGAAGAGAAGAACTGGTGCGAGGAGCAGTACTCCGGGGGCTGCTACACCACC\n+TACTTTCCCCCTGGGATCATGACTCAATATGGAAGGGTTCTACGACAGCCTGTGGGCAGG\n+ATTTATTTTGCAGGCACGGAGACTGCCACCCACTGGAGTGGCTACATGGAGGGCGCTGTG\n+GAGGCTGGAGAGAGAGCGGCCCGAGAGATCCTGCATGCCATGGGGAAGATCCCAGAGGAT\n+GAAATCTGGCAGTCAGAACCAGAATCAGTGGATGTCCCTGCGCAGCCCATCACCACGACC\n+TTCTTGGAGAGACATTTGCCCTCTGTGCCAGGCCTGCTGAGGCTGATCGGATTGACCGCC\n+ATCTTTTCTGCAACTGCTCTTGGCGTCCTGGCACACAAAAGGGGTCTGCTTGTGCGGGTC\n+TAA\n'
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa test-data/test5.sqlite
b
Binary file test-data/test5.sqlite has changed
b
diff -r b3ba0c84667c -r 56bbdbfe3eaa test-data/test5_filtered.fasta
--- /dev/null Thu Jan 01 00:00:00 1970 +0000
+++ b/test-data/test5_filtered.fasta Wed Apr 25 11:00:33 2018 -0400
b
b'@@ -0,0 +1,842 @@\n+>ENSPTRT00000034846_pantroglodytes\n+ATGGACCCGGAATGCGCCCAGCTGCTCCCGGCTCTCTGTGCTGTTCTGGTAGATCCCAGG\n+CAGCCGGTGGCAGATGACACCTGTTTGGAGAAGCTCCTGGACTGGTTTAAAACGGTCACT\n+GAAGGAGAGTCCAGTGTCGTGCTGCTGCAGGAGCACCCCTGCCTGGTGGAGCTGCTGTCC\n+CATGTGCTGAAAGTCCAGGACCTGAGTTCTGGGGTCCTCTCCTTCTCACTCCGCCTGGCA\n+GGAACCTTCGCAGCCCAGGAAAACTGCTTCCAGTATCTTCAGCAGGGGGAGTTACTACCA\n+GGGCTCTTTGGGGAGGCAGGACCCCTCGGCCGAGCAGCCTGGGCCGTCCCCACCGTGCGC\n+AGCGGCTGGATCCAGGGCCTGCGCTCCCTGGCACAGCACCCCAGCGCCCTGCGCTTCCTG\n+GCCGACCACGGTGCGGTCGACACCATCTTCTCCCTTCAGGGAGACTCCAGCCTGTTTGTG\n+GCCTCGGCGGCCAGTCAGCTCCTGGTACACGTCCTGGCTTTGTCCATGCAAGGTGGAGCC\n+GAGGGACAGCCCTGCCTGCCGGGGGGTGACTGGCCCGCGTGTGCCCAGAGGATCATGGAT\n+CACGTTGAAGAGTCCTTGTGCTCCGCGGCCACCCCCAAGGTCACTCAGGCCCTGAACGTC\n+CTGACCACAACCTTCGGGCGCTGCCAGAGCCCCTGGACGGAAGCCCTGTGGGTGCGGCTG\n+AGTCCCCGCGTGGCGTGTCTGCTGGAGAGAGACCCCATCCCCGCCGCACACTCGTTCGTG\n+GACCTGCTTCTCTGTGTGGCTCGTTCTCCCGTGTTCAGTTCTTCCGACGGCAGCCTGTGG\n+GAGACAGTGGCGCGGGCTCTGAGCTGCCTGGGTCCCACCCACATGGGACCCCTGGCTTTG\n+GGGATCCTGAAGCTCGAGCACTGTCCACAGGCACTGAGGACCCAGGCCTTCCAGGTCCTT\n+CTCCAGCCCCTGGCCTGTGTCCTGAAGGCCACGGTTCAGGCCCCCGGACCCCCAGGCTTG\n+CTGGACGGGACGGCAGACGATGCCACGACGGTGGACACACTCCTGGCCTCCAAGTCGTCC\n+TGCGCGGGCCTCCTGTGCCGCACCCTGGCTCACCTGGAGGAGCTGCAGCCGCTGCCCCAG\n+CGCCCTTCACCGTGGCCCCAGGCGTCTCTACTGGGGGCTACAGTGACTGTCCTGCGGCTC\n+TGTGACGGCTCAGCTGCCCCTGCCTCCAGTGTGGGGGGCCACCTCTGTGGGACCCTGGCG\n+GGCTGCGTCCGGGTCCAGCGAGCAGCCCTCGACTTCCTGGGGACACTGTCACAGGGGACA\n+GGCCCCCAGGAGCTGGTGACGCAGGCGCTTGCTGTCCTCCTGGAGTGCCTCGAGAGCCCC\n+GGCTCCAGCCCCACGGTTCTGAAGAAGGCCTTCCAGGCCACGCTCAGGTGGCTCCTGAGC\n+TCACCCAAGACCCCCGGCTGCTCTGATCTCGGCCCCCTCATCCCGCAGTTCCTCAGAGAG\n+CTGTTCCCTGTGCTGCAGAAACGCCTGTGCCACCCCTGCTGGGAGGTGAGGGACTCCGCC\n+CTCGAGTTCCTGACCCAGCTGAGCAGGCACTGGGGAGGACAGGCCGACTTCAGATGCGCA\n+CTCTTGGCTTCAGAGGTGCCTGAGCTGGCCCTGCAGCTCCTCCAGGACCCTGAGAGTTAT\n+GTCCGAGCGAGCGCGGTGACCGCCATGGGGCAGCTGTCCAGCCAGGGCCTGCACGCCCCC\n+ACCAGCCCTGAGCATGCAGAGGCCCGGCAGAGCCTGTTCCCGGAGCTCCTGCACATCCTC\n+TCCGTAGACTCGGAGGGCTTCCCACGGCGGGCGGTCATGCAAGTCTTCACTGAGTGGCTG\n+CGGGACGGCCACGCCGACGCGGCCCGGGACACGGAGCAGTTCGTGGCCACTGTGCTGCAG\n+GTGGCGAGCCGGGACCTGGACTGGGAGGTCCGCGCCCAGGGCCTGGAGCTGGCCCTCGTG\n+TTCCTGGGCCAGACTTTGGGGCCGCCGCGTACCCACTGCCCCTATGCCGTGGCCCTACCC\n+GAGGTGGCCCCAGCCCAGCCACTCACCGAGGCACTGAGGGCTCTCTGCCACGTGGGGCTC\n+TTTGACTTCGCCTTTTGTGCCTTGTTTGACTGCGACCGCCCTGTGGCGCAGAAGTCTTGT\n+GACCTCCTTCTCTTCCTGAGGGACAAGATTGCTTCCTACAGCAGCCTGCGGGAGGCCAGG\n+GGCGGCCCCAACACTGCCTCCGCAGAGGCCACCCTGCCGAGGTGGCGGGCGGGTGAGCAG\n+GCCCAGCCCCCAGGGGACCAGGAGCCTGAGGCTGTGCTGGCCATGCTCAGGTCCCTAGAC\n+CTGGAGGGCCTGCGGAGCACACTGGCCGAGAGCAGCGACCACGTGGAGAAGAGTCCCCAG\n+TCCCTCCTGCAGGACATGCTGGCCACGGGAGGCTTCCTGCAGGGGGATGAGGCCGACTGC\n+TACTGA\n+>ENST00000340611_homosapiens\n+ATGGACCCAGAATGCGCCCAGCTGCTCCCGGCTCTCTGTGCTGTTCTGGTAGATCCCAGG\n+CAGCCGGTGGCAGATGACACCTGTTTGGAGAAGCTCCTGGACTGGTTTAAAACGGTCACT\n+GAAGGAGAGTCCAGTGTCGTGCTGCTGCAGGAGCACCCCTGCCTGGTGGAGCTGCTGTCC\n+CATGTGCTGAAAGTCCAGGACCTGAGTTCTGGGGTCCTCTCCTTCTCACTGCGCCTGGCA\n+GGAACCTTCGCAGCCCAGGAAAACTGCTTCCAGTATCTTCAGCAGGGGGAGTTACTACCA\n+GGGCTCTTTGGGGAGCCAGGACCCCTCGGCCGAGCAACCTGGGCCGTCCCCACCGTGCGC\n+AGCGGCTGGATCCAGGGCCTGCGCTCCCTGGCACAGCACCCCAGCGCCCTGCGCTTCCTG\n+GCCGACCATGGTGCGGTCGACACCATCTTCTCCCTGCAGGGAGACTCCAGCCTGTTTGTG\n+GCCTCGGCGGCCAGTCAGCTCCTGGTGCACGTCCTGGCTTTGTCCATGCGAGGTGGAGCC\n+GAGGGGCAGCCCTGCCTGCCGGGGGGTGACTGGCCCGCGTGTGCCCAGAAGATCATGGAT\n+CACGTTGAAGAGTCCTTGTGCTCCGCGGCCACCCCCAAGGTCACTCAGGCCCTGAACGTC\n+CTGACCACGACCTTCGGGCGCTGCCAGAGCCCCTGGACGGAAGCCCTGTGGGTGCGGCTG\n+AGTCCCCGCGTGGCCTGTCTGCTGGAGAGAGACCCCATCCCCGCCGCACACTCGTTCGTG\n+GACCTGCTTCTCTGTGTGGCTCGTTCTCCCGTGTTCAGTTCTTCCGACGGCAGCCTGTGG\n+GAGACAGTGGCGCGGGCTCTGAGCTGCCTGGGTCCCACCCACATGGGACCCCTGGCTTTG\n+GGGATCCTGAAGCTCGAGCACTGTCCACAGGCACTGAGGACCCAGGCCTTCCAGGTCCTT\n+CTCCAGCCCCTGGCCTGTGTCCTGAAGGCCACGGTTCAGGCCCCCGGACCCCCAGGCTTG\n+CTGGACGGGACGGCAGACGATGCCACGACGGTGGACACACTCCTGGCCTCCAAGTCGTCC\n+TGCGCCGGCCTCCTGTGCCGCACCCTGGCTCACCTGGAGGAGCTGCAGCCGCTGCCCCAG\n+CGCCCTTCACCGTGGCCCCAGGCGTCTCTACTGGGGGCTACAGTGACTGTCCTGCGGCTC\n+TGTGACGGCTCGGCTGCCCCTGCCTCCAGTGTGGGGGGCCACCTCTGTGGGACCCTGGCG\n+GGCTGCGTCCGGGTCCAGCGAGCAGCCCTCGACTTCCTGGGGACGCTGTCACAGGGGACA\n+'..b'GCACAATCAG\n+AGTGAGTATGACGACTCGGCCAGCGAGTGCTGCTCATGTCCTAAGACTGACTCTCAGATC\n+CTGAAGGAGCTGGAGGAGTCTTCATTCAGGAAGACCTTCGAGGATTACCTGCACAACGTG\n+GTTTTTGTTCCCAGAAAAACCTCTTCAGGCAATGGTGCTGAGGACACTAGGCCATCCCGA\n+AAGCGAAGATCCCTTGAAGAGGTGGGCAATGTGACAGCCACTACACCCACACTTCCAGAT\n+TTTCCCAACATCTCCTCCACCATCGCGCCCACAAGCCACGAAGAGCACAGACCATTTGAG\n+AAAGTAGTAAACAAGGAGTCACTTGTCATCTCTGGCCTGAGACACTTCACTGGGTACCGC\n+ATTGAGCTGCAGGCATGCAATCAGGACTCCCCAGAAGAGAGGTGCAGCGTGGCTGCCTAC\n+GTCAGTGCCCGGACCATGCCTGAAGCTAAGGCAGATGACATCGTTGGCCCTGTGACCCAT\n+GAAATCTTTGAGAACAATGTTGTACACTTAATGTGGCAAGAGCCAAAGGAACCTAATGGT\n+CTGATTGTGCTATATGAAGTGAGCTATCGGCGATATGGTGATGAGGAGCTGCACCTCTGT\n+GTCTCCCGGAAGCATTTTGCCCTGGAGCGGGGCTGCAGGCTTCGAGGGCTCTCTCCAGGA\n+AACTACAGTGTTCGAGTCCGGGCTACCTCTCTGGCAGGAAATGGCTCCTGGACAGAACCC\n+ACCTATTTTTATGTGACTGATTATTTAGATGTCCCATCAAATATTGCCAAAATTATCATC\n+GGGCCCCTCATCTTCGTCTTCCTCTTCAGTGTCGTGATCGGAAGTATTTATCTATTCTTG\n+AGGAAGAGGCAGCCAGATGGGCCAATGGGACCACTGTACGCTTCTTCAAACCCAGAGTAC\n+CTCAGTGCCAGTGATGTCTTTCCATCTTCCGTATACGTTCCGGATGAGTGGGAGGTACCT\n+CGAGAGAAGATCACCCTCCTCCGAGAGCTGGGGCAGGGATCCTTCGGTATGGTGTACGAA\n+GGCAATGCCAAGGATATCATCAAGGGTGAGGTAGAGACCCGTGTTGCGGTGAAGACGGTC\n+AATGAGTCAGCCAGTCTTCGAGAACGGATCGAGTTCCTCAATGAGGCATCAGTCATGAAG\n+GGCTTCACCTGTCATCACGTGGTCCGCCTTCTTGGGGTGGTGTCCAAAGGCCAGCCCACA\n+TTGGTAGTGATGGAACTGATGGCTCATGGAGACCTGAAAAGTCACCTCCGTTCTCTGCGG\n+CCCGATGCTGAGAACAACCCAGGCCGTCCTCCCCCTACCTTGCAAGAAATGATTCAGATG\n+ACAGCAGAAATTGCCGATGGCATGGCATACTTGAACGCCAAGAAGTTTGTGCACCGGGAC\n+CTGGCAGCTCGGAACTGCATGGTTGCCCATGATTTTACTGTCAAAATCGGAGACTTTGGA\n+ATGACGAGAGACATCTACGAGACAGATTACTATCGGAAAGGGGGCAAGGGGTTGCTGCCC\n+GTGAGGTGGATGTCACCCGAGTCCCTGAAGGACGGAGTCTTCACTGCTTCTTCCGACATG\n+TGGTCCTTTGGGGTGGTCCTTTGGGAAATCACCAGCCTGGCTGAGCAACCTTACCAAGGC\n+CTGTCTAATGAACAGGTGTTGAAATTTGTCATGGATGGAGGCTATCTGGATCCCCCTGAT\n+AACTGTCCAGAGAGACTCACTGACCTGATGCGCATGTGCTGGCAGTTCAACCCCAAGATG\n+AGGCCGACCTTCCTGGAAATCGTCAACCTGCTCAAGGACGACCTCCACCCCAGCTTTCCG\n+GAAGTTTCCTTCTTCTACAGCGAGGAGAACAAGGCTCCCGAGAGTGAAGAGCTGGAGATG\n+GAGTTCGAGGACATGGAGAATGTCCCCTTGGATCGTTCCTCTCACTGTCAGAGAGAAGAG\n+GCTGGATGCCGGGAGGGAGGGTCCTCTCTGAGCATCAAACGGACCTATGATGAACACATC\n+CCCTACACCCACATGAACGGGGGCAAGAAGAATGGGCGGGTCCTCACCCTGCCGAGGTCG\n+AACCCTTCCTAA\n+>ENSSSCT00000014817_susscrofa\n+GTGTGCCCAGGGATGGATATCCGGAATAACCTTACACGGCTGCACGAGTTGGCCAACTGC\n+TCGGTCATCGAAGGACATTTGCAGATCCTGTTGATGTTCAAAACGCGGCCCGAGGATTTC\n+CGAGACCTCAGTTTCCCCAAACTCATCATGATCACTGATTACTTGCTGCTCTTCCGGGTC\n+TACGGGCTGGAGAGCCTGAAGGACCTGTTCCCCAACCTCACCGTCATCCGGGGCTCACGC\n+CTCTTCTTTAACTATGCGCTGGTCATCTTTGAGATGGTTCACTTGAAGGAGCTTGGCCTC\n+TACAATTTGATGAACATCACCAGGGGTGCTGTCCGCATCGAGAAGAACAATGAGCTCTGC\n+TACCTGGCGACCATTGACTGGTCGCGCATCCTGGACTCTGTGGAGGATAATTACATTGTG\n+CTGAACAAAGACGACAACGAGGAGTGTGGGGACATTTGCCCAGGCACTGCGAAGGGCAAG\n+ACCAATTGCCCTGCCACCGTCATCAATGGGCAATTTGTCGAGCGGTGTTGGACGCACAGT\n+CACTGCCAGAAAGTGTGCCCGACCATCTGTAAGTCGCACGGCTGCACTGCTGAGGGCCTC\n+TGCTGTCACAGCGAGTGTTTGGGCAACTGCTCTGAGCCAGACGACCCCACCAAGTGCGTG\n+GCCTGCCGCAACTTCTACCTGGACGGCAGATGCGTGGAGACCTGCCCGCCCCCCTACTAC\n+CACTTCCAAGACTGGCGCTGCGTGAACTTCAGCTTCTGCCAGGACCTGCACAACAAATGC\n+AAGAACTCAAGGAGGCAGGGCTGCCACCAGTACGTCATTCACAACAACAAGTGTATCCCT\n+GAGTGCCCCTCAGGGTACACGATGAATTCCAGCAACTTGATGTGCACTCCGTGCCTAGGC\n+CCCTGTCCCAAAGTGTGTCACCTCCTGGAAGGCGAGAAGACCATCGACTCAGTGACATCC\n+GCCCAGGAGCTCCGAGGCTGCACCATTATCAACGGGAGCCTAATCATCAACATTCGAGGA\n+GGCAACAACCTGGCAGCCGAACTAGAGGCCAACCTTGGACTCATTGAGGAGATTTCAGGG\n+TACCTGAAAATCCGCCGATCCTATGCCCTCGTGTCACTTTCCTTCTTCCGGAAGTTGCGT\n+CTGATCCGAGGGGAGACGTTGGAAATTGGGAACTATTCTTTCTATGCCTTGGACAACCAG\n+AACCTAAGGCAACTGTGGGACTGGAGCAAACACAACCTCACCATCACTCAGGGGAAACTC\n+TTCTTCCATTATAATCCCAAACTCTGCTTGTCGGAAATTCACAAGATGGAGGAAGTTTCT\n+GGAACCAAGGGGCGCCAGGAGAGAAATGATATTGCCCTGAAGACCAATGGGGACCAGGCG\n+TCCTGTGAAAATGAGTTACTTAAATTTTCTTACATTCGGACATCTTATGACAAGATCTTG\n+CTGAAGTGGGAGCCGTATTGGCCCCCCGACTTCCGAGACCTCCTGGGGTTCATGCTCTTC\n+TACAAAGAGGCCCCTTATCAGAACGTGACGGAGTTTGACGGGCAGGATGCGTGTGGCTCC\n+AACAGCTGGACGGTGGTGGACATTGACCCGCCTACGAGGTCCAATGACCCCAAGTCCCAG\n+AACCATCCTGGGTGGCTGATGCGTGGTCTCAAGCCCTGGACCCAGTATGCCATCTTTGTC\n+AAGACTTTGGTCACCTTTTCTGATGAACGACGCACCTATGGAGCCAAGAGTGACATCATC\n+TACGTCCAGACAGATGCCACAAGTAAGCATGTC\n+\n'