[med-svn] [Git][med-team/resfinder][master] 7 commits: Build-Depends: s/python3-pdm/python3-pdm-backend/g
Andreas Tille (@tille)
gitlab at salsa.debian.org
Tue Jan 16 07:06:29 GMT 2024
Andreas Tille pushed to branch master at Debian Med / resfinder
Commits:
63e23486 by Andreas Tille at 2024-01-16T07:09:27+01:00
Build-Depends: s/python3-pdm/python3-pdm-backend/g
- - - - -
8b5b3cf9 by Andreas Tille at 2024-01-16T07:10:34+01:00
routine-update: New upstream version
- - - - -
014c078b by Andreas Tille at 2024-01-16T07:24:05+01:00
New upstream version 4.4.2
- - - - -
e7450de4 by Andreas Tille at 2024-01-16T07:24:07+01:00
Update upstream source from tag 'upstream/4.4.2'
Update to upstream version '4.4.2'
with Debian dir 9e4eb8139454bb8c1786d18306759d79665a5d77
- - - - -
03da7d8c by Andreas Tille at 2024-01-16T07:24:46+01:00
routine-update: Build-Depends: s/dh-python/dh-sequence-python3/
- - - - -
c0aba9b8 by Andreas Tille at 2024-01-16T08:02:03+01:00
Build-Depends: python3-pdm-pep517
- - - - -
4674a866 by Andreas Tille at 2024-01-16T08:04:26+01:00
Upload to unstable
- - - - -
23 changed files:
- + .pdm-python
- CHANGELOG.md
- Dockerfile_no_dbs
- bitbucket-pipelines.yml
- debian/changelog
- debian/control
- debian/rules
- dockerfile
- + dockerfile_EFSA_2022
- src/resfinder/__init__.py
- src/resfinder/cge/output/gene_result.py
- src/resfinder/cge/output/std_results.py
- src/resfinder/cge/phenotype2genotype/isolate.py
- src/resfinder/cge/phenotype2genotype/res_profile.py
- src/resfinder/cge/pointfinder.py
- src/resfinder/cge/resfinder.py
- src/resfinder/run_resfinder.py
- + tests/data/test_isolate_12.fa
- + tests/data/test_isolate_14.fa
- + tests/data/test_isolate_14_1.fq
- + tests/data/test_isolate_14_2.fq
- tests/resfinder/cge/output/test_gene_result.md
- tests/resfinder/cge/output/test_phenotype_result.md
Changes:
=====================================
.pdm-python
=====================================
@@ -0,0 +1 @@
+/Users/rolf/workcloud/Scripts-CGE/resfinder/.venv/bin/python
\ No newline at end of file
=====================================
CHANGELOG.md
=====================================
@@ -5,6 +5,50 @@ All notable changes to this project will be documented in this file.
The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
+## [4.4.2] - 2023-11-27
+
+### Fixed
+
+- Issue where mutations in promoters that are not indels (i.e. substitutions) could be missed as they were uncorrectly registered as unknown mutations.
+
+## [4.4.1] - 2023-11-02
+
+### Fixed
+
+- Issue where mutations in promoters could be missed if using fastq files as input. Issue continues to persist for PointFinder database versions < 4.0.1
+
+## [4.4.0] - 2023-10-27
+
+### Fixed
+
+- Issue where application crashed with a "KeyError: 'var_aa'", when finding mutation in promoters.
+- Issue where PointFinder database reference sequence names containing "promoter-size-\d+bp" would be missed. Issue continues to persist for PointFinder database versions < 4.0.0
+- Spelling error in json output "ref_seq_lenght". It is now "ref_seq_length"
+
+## [4.3.3] - 2023-08-22
+
+### Fixed
+
+- Issue where ResFinder crashed with the error: "UnboundLocalError: local variable \'hit_class\' referenced before assignment'" if the flag '--output_aln' was enabled along with FASTQ input
+
+## [4.3.2] - 2023-05-31
+
+### Fixed
+
+- Issue where PointFinder part crashed with the error: "NameError: name 'combined_id' is not defined"
+
+## [4.3.1] - 2023-03-03
+
+### Removed
+
+- Overlap warning from pointfinder (message was only intended for testing)
+
+### Fixed
+
+- dockerfile, so that it works with NextFlow
+- Issue with alignment output in json
+- Issue with reporting concatenated gene hits in PointFinder
+
## [4.3.0] - 2023-02-09
### Added
=====================================
Dockerfile_no_dbs
=====================================
@@ -8,19 +8,20 @@ RUN apt-get update -qq; \
wget \
ncbi-blast+ \
libz-dev \
+ procps \
; \
rm -rf /var/cache/apt/* /var/lib/apt/lists/*;
# Install KMA
RUN cd /usr/src/; \
- git clone --depth 1 -b 1.4.3 https://bitbucket.org/genomicepidemiology/kma.git; \
+ git clone --depth 1 -b 1.4.11 https://bitbucket.org/genomicepidemiology/kma.git; \
cd kma && make; \
mv kma /usr/bin/; \
cd ..; \
rm -rf kma/;
# Install ResFinder
-RUN pip install --no-cache-dir resfinder==4.3.0
+RUN pip install --no-cache-dir resfinder==4.3.1
# Setup environment
RUN cd /; \
=====================================
bitbucket-pipelines.yml
=====================================
@@ -1,10 +1,15 @@
# Testing ResFinder
-
+clone:
+ depth: full
image: python:3.10
pipelines:
custom:
- test-and-build:
+ mirror-to-github:
+ - step:
+ script:
+ - git push --mirror git at github.com:genomicepidemiology/resfinder
+ test-and-build:
- step:
name: Test
caches:
@@ -29,4 +34,4 @@ pipelines:
- export PATH=/root/.local/bin:$PATH
- pdm install
- pdm run test
- - pdm build
+ - pdm build
\ No newline at end of file
=====================================
debian/changelog
=====================================
@@ -1,3 +1,13 @@
+resfinder (4.4.2-1) unstable; urgency=medium
+
+ * New upstream version
+ * Build-Depends: s/python3-pdm/python3-pdm-backend/g
+ Closes: #1060881
+ * Build-Depends: python3-pdm-pep517
+ * Build-Depends: s/dh-python/dh-sequence-python3/ (routine-update)
+
+ -- Andreas Tille <tille at debian.org> Tue, 16 Jan 2024 08:02:14 +0100
+
resfinder (4.3.0-1) unstable; urgency=medium
* New upstream version
=====================================
debian/control
=====================================
@@ -5,11 +5,12 @@ Uploaders: Andreas Tille <tille at debian.org>,
Section: science
Priority: optional
Build-Depends: debhelper-compat (= 13),
- dh-python,
+ dh-sequence-python3,
pybuild-plugin-pyproject,
python3,
python3-cgelib,
- python3-pdm
+ python3-pdm-backend,
+ python3-pdm-pep517
Standards-Version: 4.6.2
Vcs-Browser: https://salsa.debian.org/med-team/resfinder
Vcs-Git: https://salsa.debian.org/med-team/resfinder.git
=====================================
debian/rules
=====================================
@@ -3,7 +3,7 @@
export PYBUILD_NAME=resfinder
%:
- dh $@ --with python3 --buildsystem=pybuild
+ dh $@ --buildsystem=pybuild
override_dh_install:
dh_install
=====================================
dockerfile
=====================================
@@ -9,6 +9,7 @@ RUN apt-get update -qq; \
wget \
ncbi-blast+ \
libz-dev \
+ procps \
; \
rm -rf /var/cache/apt/* /var/lib/apt/lists/*;
@@ -20,7 +21,7 @@ RUN cd /usr/src/; \
cd ..; \
rm -rf kma/;
-ENV RESFINDER_VERSION 4.3.0
+ENV RESFINDER_VERSION 4.4.2
# Install ResFinder
RUN pip install --no-cache-dir resfinder==$RESFINDER_VERSION
=====================================
dockerfile_EFSA_2022
=====================================
@@ -0,0 +1,60 @@
+FROM --platform=linux/amd64 python:3.10-slim-bullseye
+
+ENV INSIDE_DOCKER true
+ENV DEBIAN_FRONTEND noninteractive
+
+RUN apt-get update -qq; \
+ apt-get install --no-install-recommends -y -qq git \
+ build-essential \
+ wget \
+ ncbi-blast+ \
+ libz-dev \
+ procps \
+ ; \
+ rm -rf /var/cache/apt/* /var/lib/apt/lists/*;
+
+# Install python dependencies
+RUN pip3 install --upgrade pip
+RUN pip3 install -U tabulate===0.8.9 cgecore==1.5.6 cgelib==0.7.3 pandas==1.4.2 \
+ biopython==1.79 gitpython==3.1.27 python-dateutil==2.8.2;
+
+# Install KMA
+# Note: In future versions remove kma_index as it is only needed in old versions of pointfinder
+RUN cd /usr/src/; \
+ git clone --depth 1 -b 1.4.11 https://bitbucket.org/genomicepidemiology/kma.git; \
+ cd kma && make; \
+ mv kma /usr/bin/; \
+ mv kma_index /usr/bin/; \
+ cd ..; \
+ rm -rf kma/;
+
+# Install BLAST
+RUN cd /usr/src/; \
+ wget "https://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/2.13.0/ncbi-blast-2.13.0+-x64-linux.tar.gz"; \
+ tar -xf ncbi-blast-2.13.0+-x64-linux.tar.gz; \
+ mv ncbi-blast-2.13.0+/bin/* /usr/bin/;
+
+ENV RESFINDER_VERSION EFSA_2022
+
+# Install ResFinder
+RUN cd /usr/src/; \
+ git clone -b $RESFINDER_VERSION --depth 1 https://bitbucket.org/genomicepidemiology/resfinder.git;
+
+# Install databases
+RUN cd /usr/src/resfinder; \
+ git clone -b $RESFINDER_VERSION --depth 1 https://git@bitbucket.org/genomicepidemiology/resfinder_db.git db_resfinder; \
+ git clone -b $RESFINDER_VERSION --depth 1 https://git@bitbucket.org/genomicepidemiology/pointfinder_db.git db_pointfinder; \
+ cd /usr/src/resfinder/db_resfinder; \
+ python INSTALL.py; \
+ rm -rf .git; \
+ cd /usr/src/resfinder/db_pointfinder; \
+ python INSTALL.py; \
+ rm -rf .git;
+
+# Setup environment
+RUN cd /; \
+ mkdir app;
+WORKDIR /app
+
+# Execute program when running the container
+ENTRYPOINT ["python", "/usr/src/resfinder/run_resfinder.py"]
=====================================
src/resfinder/__init__.py
=====================================
@@ -1 +1 @@
-__version__ = '4.3.0'
+__version__ = "4.4.2"
=====================================
src/resfinder/cge/output/gene_result.py
=====================================
@@ -30,7 +30,7 @@ class GeneResult(dict):
self["ref_database"] = [res_collection.get_db_key(db_name)[0]]
self["identity"] = res["perc_ident"]
self["alignment_length"] = res["HSP_length"]
- self["ref_seq_lenght"] = res["sbjct_length"]
+ self["ref_seq_length"] = res["sbjct_length"]
self["depth"] = res.get("depth", None)
self["ref_start_pos"] = res["sbjct_start"]
self["ref_end_pos"] = res["sbjct_end"]
@@ -41,12 +41,7 @@ class GeneResult(dict):
self["notes"] = []
if conf and conf.output_aln:
- for ab_class, hits in alignments.gene_align_query.items():
- ab_class_keys = list(hits.keys())
- hit_key = next(key for key in ab_class_keys
- if key.startswith(self["ref_id"]))
- hit_class = ab_class
- break
+ hit_class, hit_key = self.get_hit_class(alignments.gene_align_query)
self["query_string"] = alignments.gene_align_query[hit_class][
hit_key]
@@ -56,19 +51,53 @@ class GeneResult(dict):
# BLAST coverage formatted results
coverage = res.get("coverage", None)
- if(coverage is None):
+ if (coverage is None):
# KMA coverage formatted results
coverage = res["perc_coverage"]
else:
coverage = float(coverage) * 100
self["coverage"] = coverage
- self["grade"] = GeneResult.calc_gene_grade(coverage=self["coverage"], identity=self["identity"])
+ self["grade"] = GeneResult.calc_gene_grade(
+ coverage=self["coverage"], identity=self["identity"])
self.remove_NAs()
uniqueness = self._get_unique_gene_key(res_collection)
self["key"] = uniqueness
+ def get_hit_class(self, gene_align_query):
+ """
+ Input:
+ gene_align_query: Dict from either FinderResult (KMA) or
+ Blaster (BLAST)
+ Output:
+ hit_class: Name of AB class
+ hit_key: Key used for query (hit)
+
+ Method looks through all the queries that align to a reference and
+ finds the query (hit) that matches the GeneResult object.
+ For BLAST hits the query is named after the coresponding contig
+ header (query_id), but for KMA the query key is named after the
+ reference it aligns to (ref_id).
+ """
+ lookup_key = self["query_id"]
+ if self["depth"] is not None:
+ lookup_key = self["ref_id"]
+
+ for ab_class, hits in gene_align_query.items():
+ if hits:
+ hits_keys = list(hits.keys())
+ hit_key = next((key for key in hits_keys
+ if key.startswith(lookup_key)), None)
+ if hit_key is None:
+ continue
+ hit_class = ab_class
+ break
+ else:
+ continue
+
+ return hit_class, hit_key
+
@staticmethod
def calc_gene_grade(coverage: float, identity: float) -> int:
if coverage == 100.0 and identity == 100.0:
@@ -87,7 +116,7 @@ class GeneResult(dict):
"""
na_keys = []
for key, val in self.items():
- if(val == "NA" or val is None):
+ if (val == "NA" or val is None):
na_keys.append(key)
for key in na_keys:
del self[key]
@@ -107,7 +136,7 @@ class GeneResult(dict):
If gene_key is found in res_collection. Creates a unique key by
appending a random string ton minimum_gene_key.
"""
- while(gene_key in res_collection["seq_regions"]):
+ while (gene_key in res_collection["seq_regions"]):
rnd_str = GeneResult.random_string(str_len=4)
gene_key = ("{key}{deli}{rnd}"
.format(key=minimum_gene_key, deli=delimiter,
@@ -143,7 +172,7 @@ class GeneResult(dict):
sbjct = header.split("_")
template = sbjct[0]
- if(len(sbjct) > 1):
+ if (len(sbjct) > 1):
variant = sbjct[1]
acc = "_".join(sbjct[2:])
else:
@@ -180,7 +209,7 @@ class GeneResult(dict):
query_id = self.get("query_id", "NA")
# Query id == "NA" when FASTQ
- if(query_id == "NA"):
+ if (query_id == "NA"):
res_collection["seq_regions"][gene_key]["ref_database"].extend(
self["ref_database"])
gene_key = None
=====================================
src/resfinder/cge/output/std_results.py
=====================================
@@ -2,6 +2,7 @@
import random
import string
import json
+import re
from ..phenotype2genotype.res_profile import PhenoDB
from ..phenotype2genotype.feature import ResGene, ResMutation
@@ -11,6 +12,7 @@ from .gene_result import GeneResult
from .seq_variation_result import SeqVariationResult
from .phenotype_result import PhenotypeResult
+from ..pointfinder import PointFinder
def add_gene_result_if_key_not_None(gene_result, res_collection):
'''
@@ -190,3 +192,27 @@ class PointFinderResultHandler():
res_collection, mismatch, gene_results, ref_db_name)
res_collection.add_class(cl="seq_variations",
**seq_var_result)
+
+ @staticmethod
+ def rename_promoter_results(results, method):
+ if method == PointFinder.TYPE_BLAST:
+ PointFinderResultHandler._rename_blast_promoter_results(results)
+ else:
+ PointFinderResultHandler._rename_kma_promoter_results(results)
+
+ def _rename_blast_promoter_results(results):
+ promoter_keys = [k for k in results if re.search(r"promoter-size", k)]
+ for promoter_key in promoter_keys:
+ val_dict = results.pop(promoter_key)
+ new_promoter_key = PhenoDB.if_promoter_gene_rename(promoter_key)
+ results[new_promoter_key] = val_dict
+
+ def _rename_kma_promoter_results(results):
+ promoter_keys = set()
+
+ for inner_dict in results.values():
+ promoter_keys = [k for k in inner_dict if re.search(r"promoter-size", k)]
+ for promoter_key in promoter_keys:
+ val_dict = inner_dict.pop(promoter_key)
+ new_promoter_key = PhenoDB.if_promoter_gene_rename(promoter_key)
+ inner_dict[new_promoter_key] = val_dict
=====================================
src/resfinder/cge/phenotype2genotype/isolate.py
=====================================
@@ -234,7 +234,7 @@ class Isolate(dict):
f"_{feat_res_dict['nuc_change']}")
aa_format = (f"{feat_res_dict['seq_regions'][0]}"
f"_{feat_res_dict['ref_start_pos']}"
- f"_{feat_res_dict['var_aa']}")
+ f"_{var_aa}") # var_aa can be None
return nuc_format, aa_format
# Amino acid mutation
else:
@@ -395,7 +395,7 @@ class Isolate(dict):
hit = DBHit(name=gene_info["name"],
identity=gene_info["identity"],
match_length=gene_info["alignment_length"],
- ref_length=gene_info["ref_seq_lenght"],
+ ref_length=gene_info["ref_seq_length"],
start_ref=gene_info["ref_start_pos"],
end_ref=gene_info["ref_end_pos"],
acc=gene_info["ref_acc"],
@@ -481,7 +481,7 @@ class Isolate(dict):
elif(feature.hit.match_category > best_match):
best_match = feature.hit.match_category
- output_str += "\t" + str(best_match)
+ output_str += "\t" + str(best_match)
gene_list = ab.get_gene_namewacc(tostring=True)
mut_list = ab.get_mut_namewannot(tostring=True)
=====================================
src/resfinder/cge/phenotype2genotype/res_profile.py
=====================================
@@ -610,6 +610,22 @@ class PhenoDB(dict):
return out_string
+ @staticmethod
+ def if_promoter_gene_rename(s):
+ """
+ Needed instead of if_promoter_rename() when only gene name is
+ present without variant number and acc no
+ """
+ out_string = s
+
+ regex = r"^(.+-promoter)-size-\d+bp$"
+ promoter_match = re.search(regex, s)
+
+ if promoter_match:
+ out_string = f"{promoter_match.group(1)}"
+
+ return out_string
+
@staticmethod
def get_csv_tuple(csv_string, sep=",", lower=True):
""" Takes a string containing a comma seperated list, makes it all
=====================================
src/resfinder/cge/pointfinder.py
=====================================
@@ -43,6 +43,7 @@ class PointFinder(CGEFinder):
PHENOTYPE_FILE = "phenotypes.txt"
# Deprecated variable, use phenotype (kept for legacy dbs)
RES_OVERVIEW_FILE = "resistens-overview.txt"
+ PROMOTER_REGEX = r"(.+promoter)-size-(\d+)(?:bp)"
def __init__(self, db_path, species, gene_list=None, ignore_indels=False,
ignore_stop_codons=False):
@@ -57,6 +58,8 @@ class PointFinder(CGEFinder):
self.gene_list = PointFinder.get_file_content(
f"{self.specie_path}/{self.GENE_FILE}")
+ self.promoters = PointFinder.load_promoters(self.gene_list)
+
if os.path.isfile(f"{self.specie_path}/{self.RNA_GENE_FILE}"):
self.RNA_gene_list = PointFinder.get_file_content(
f"{self.specie_path}/{self.RNA_GENE_FILE}")
@@ -82,6 +85,19 @@ class PointFinder(CGEFinder):
"the '{self.RES_OVERVIEW_FILE}' file "
"(legacy database)")
+ @staticmethod
+ def load_promoters(gene_file_list):
+ promoters = {}
+
+ for gene_filename in gene_file_list:
+ promoter_match = re.search(PointFinder.PROMOTER_REGEX, gene_filename)
+ if promoter_match:
+ promoter_key = promoter_match.group(1)
+ promoter_length = promoter_match.group(2)
+ promoters[promoter_key] = int(promoter_length)
+
+ return promoters
+
def get_user_defined_gene_list(self, gene_list):
genes_specified = []
for gene in gene_list:
@@ -255,12 +271,7 @@ class PointFinder(CGEFinder):
for gene, hit in GENES.items():
# Start writing output string (to HTML tab file)
gene_name = gene.split("_")[0] # Not perfeft can differ from specific mutation
- regex = r"promoter-size-(\d+)(?:bp)"
- promtr_gene_objt = re.search(regex, gene_name)
-
- if promtr_gene_objt:
- gene_name = gene.split("_")[0]
-
+ gene_name = PhenoDB.if_promoter_gene_rename(gene_name)
output_strings[1] += "\n%s\n" % (gene_name)
# Ignore excluded genes
@@ -421,7 +432,11 @@ class PointFinder(CGEFinder):
f"_{mis_match[3].lower()}")
mis_match_key_aa = (f"{gene_db_id}_{mis_match[1]}"
f"_{mis_match[-1].lower()}")
- elif mis_match[4].startswith('r.'):#RNA substitution
+ elif mis_match[4].startswith('r.'): # RNA substitution
+ mis_match_key_aa = ""
+ mis_match_key_nuc = (f"{gene_db_id}_{mis_match[2]}"
+ f"_{mis_match[3].lower()}")
+ elif mis_match[4].startswith('n.'): # non-coding DNA substitution
mis_match_key_aa = ""
mis_match_key_nuc = (f"{gene_db_id}_{mis_match[2]}"
f"_{mis_match[3].lower()}")
@@ -621,6 +636,7 @@ class PointFinder(CGEFinder):
gene_name = str(mutation[0].split("_")[0])
# Only consider mutations in genes found in the gene list
if gene_name in gene_list:
+ gene_name = PhenoDB.if_promoter_gene_rename(gene_name)
mut_pos = int(mutation[3])
ref_codon = mutation[4] # Ref_nuc (1 or 3?)
ref_aa = mutation[5] # Ref_codon
@@ -849,10 +865,6 @@ class PointFinder(CGEFinder):
# If alternative query overlap excist save it
if pre_qry_overlap != next_qry_overlap:
- eprint("OVERLAP WARNING:")
- eprint("{}\n{}"
- .format(pre_qry_overlap, next_qry_overlap))
-
# Save alternative overlaps
alternative_overlaps += [(next_block_start,
overlap_end_pos,
@@ -934,13 +946,8 @@ class PointFinder(CGEFinder):
qry_seq)
else:
# Check if the gene sequence is with a promoter
- regex = r"promoter-size-(\d+)(?:bp)"
- promtr_gene_objt = re.search(regex, gene)
-
- # Check for promoter sequences
- if promtr_gene_objt:
- # Get promoter length
- promtr_len = int(promtr_gene_objt.group(1))
+ promtr_len = self.promoters.get(gene_name, None)
+ if promtr_len is not None:
# Extract promoter sequence, while considering gaps
# --------agt-->----
=====================================
src/resfinder/cge/resfinder.py
=====================================
@@ -73,6 +73,7 @@ class ResFinder(CGEFinder):
min_cov=min_cov, threshold=threshold, blast=blast,
allowed_overlap=allowed_overlap)
self.blast_results = blast_run.results
+
return blast_run
def load_phenos(self, pheno_file):
=====================================
src/resfinder/run_resfinder.py
=====================================
@@ -282,37 +282,25 @@ def main():
conf)
else:
- if (conf.nanopore):
- kma_run = acquired_finder.kma(
- inputfile_1=conf.inputfastq_1,
- inputfile_2=conf.inputfastq_2,
- out_path=conf.outPath_res_kma,
- db_path_kma=conf.db_path_res_kma,
- min_cov=conf.rf_gene_cov,
- threshold=conf.rf_gene_id,
- kma_path=conf.kma,
- databases=acquired_finder.databases,
- sample_name="",
- kma_cge=True,
- kma_apm="p",
- kma_1t1=True,
- kma_add_args='-ont -md 5'
- )
- else:
- kma_run = acquired_finder.kma(
- inputfile_1=conf.inputfastq_1,
- inputfile_2=conf.inputfastq_2,
- out_path=conf.outPath_res_kma,
- db_path_kma=conf.db_path_res_kma,
- min_cov=conf.rf_gene_cov,
- threshold=conf.rf_gene_id,
- kma_path=conf.kma,
- databases=acquired_finder.databases,
- sample_name="",
- kma_cge=True,
- kma_apm="p",
- kma_1t1=True
- )
+ kma_nano_args = None
+ if conf.nanopore:
+ kma_nano_args = '-ont -md 5'
+
+ kma_run = acquired_finder.kma(
+ inputfile_1=conf.inputfastq_1,
+ inputfile_2=conf.inputfastq_2,
+ out_path=conf.outPath_res_kma,
+ db_path_kma=conf.db_path_res_kma,
+ min_cov=conf.rf_gene_cov,
+ threshold=conf.rf_gene_id,
+ kma_path=conf.kma,
+ databases=acquired_finder.databases,
+ sample_name="",
+ kma_cge=True,
+ kma_apm="p",
+ kma_1t1=True,
+ kma_add_args=kma_nano_args
+ )
# DEPRECATED
# TODO: make a write method that depends on the json output
@@ -363,37 +351,25 @@ def main():
conf)
else:
- if (conf.nanopore):
- kma_run = disinf_finder.kma(
- inputfile_1=conf.inputfastq_1,
- inputfile_2=conf.inputfastq_2,
- out_path=conf.outPath_disinf_kma,
- db_path_kma=conf.db_path_disinf_kma,
- min_cov=conf.dis_gene_cov,
- threshold=conf.dis_gene_id,
- kma_path=conf.kma,
- databases=disinf_finder.databases,
- sample_name="",
- kma_cge=True,
- kma_apm="p",
- kma_1t1=True,
- kma_add_args='-ont -md 5'
- )
- else:
- kma_run = disinf_finder.kma(
- inputfile_1=conf.inputfastq_1,
- inputfile_2=conf.inputfastq_2,
- out_path=conf.outPath_disinf_kma,
- db_path_kma=conf.db_path_disinf_kma,
- min_cov=conf.dis_gene_cov,
- threshold=conf.dis_gene_id,
- kma_path=conf.kma,
- databases=disinf_finder.databases,
- sample_name="",
- kma_cge=True,
- kma_apm="p",
- kma_1t1=True
- )
+ kma_nano_args = None
+ if conf.nanopore:
+ kma_nano_args = '-ont -md 5'
+
+ kma_run = disinf_finder.kma(
+ inputfile_1=conf.inputfastq_1,
+ inputfile_2=conf.inputfastq_2,
+ out_path=conf.outPath_disinf_kma,
+ db_path_kma=conf.db_path_disinf_kma,
+ min_cov=conf.dis_gene_cov,
+ threshold=conf.dis_gene_id,
+ kma_path=conf.kma,
+ databases=disinf_finder.databases,
+ sample_name="",
+ kma_cge=True,
+ kma_apm="p",
+ kma_1t1=True,
+ kma_add_args=kma_nano_args
+ )
# DEPRECATED
# TODO: make a write method that depends on the json output
@@ -434,38 +410,30 @@ def main():
results = blast_run.results
else:
-
+ kma_nano_args = None
+ if conf.nanopore:
+ kma_nano_args = '-ont -md 5'
+
method = PointFinder.TYPE_KMA
- if (conf.nanopore):
- kma_run = finder.kma(inputfile_1=conf.inputfastq_1,
- inputfile_2=conf.inputfastq_2,
- out_path=conf.outPath_point_kma,
- db_path_kma=conf.db_path_point,
- databases=[conf.species_dir],
- min_cov=0.01, # Sorts on coverage later
- threshold=conf.pf_gene_id,
- kma_path=conf.kma,
- sample_name=conf.sample_name,
- kma_cge=True,
- kma_apm="p",
- kma_1t1=True,
- kma_add_args='-ont -md 5')
- else:
- kma_run = finder.kma(inputfile_1=conf.inputfastq_1,
- inputfile_2=conf.inputfastq_2,
- out_path=conf.outPath_point_kma,
- db_path_kma=conf.db_path_point,
- databases=[conf.species_dir],
- min_cov=0.01, # Sorts on coverage later
- threshold=conf.pf_gene_id,
- kma_path=conf.kma,
- sample_name=conf.sample_name,
- kma_cge=True,
- kma_apm="p",
- kma_1t1=True)
+
+ kma_run = finder.kma(inputfile_1=conf.inputfastq_1,
+ inputfile_2=conf.inputfastq_2,
+ out_path=conf.outPath_point_kma,
+ db_path_kma=conf.db_path_point,
+ databases=[conf.species_dir],
+ min_cov=0.01, # Sorts on coverage later
+ threshold=conf.pf_gene_id,
+ kma_path=conf.kma,
+ sample_name=conf.sample_name,
+ kma_cge=True,
+ kma_apm="p",
+ kma_1t1=True,
+ kma_add_args=kma_nano_args)
results = kma_run.results
-
+
+ PointFinderResultHandler.rename_promoter_results(results, method)
+
if (conf.specific_gene):
results = PointFinder.discard_from_dict(
in_dict=results, wanted_list=conf.specific_gene)
=====================================
tests/data/test_isolate_12.fa
=====================================
@@ -0,0 +1,25 @@
+>erm(38)_first1-321
+GTGTCCACACCACATCACGGCCGGCACGAGCTCGGCCAGAACTTCCTGTCCGATCGGCGC
+GTCATCGCCGATATCGTCGAAATCGTCTCGCGCACAAACGGTCCGATCATCGAGATCGGG
+GCGGGCGACGGCGCGCTGACCATACCCTTGCAACGACTCGCCCGCCCGCTCACCGCCGTC
+GAGGTCGACGCGCGGCGCGCGCGGCGGTTGGCGCAGCGCACCGCGAGATCCGCCCCGGGG
+CCTGCCTCGCGGCCCACCGAGGTCGTCGCCGCCGACTTCCTGCGCTACCCACTGCCCCGC
+TCACCCCACGTGGTCGTGGGC
+>erm(38)_middel201-578
+CGGCGGTTGGCGCAGCGCACCGCGAGATCCGCCCCGGGG
+CCTGCCTCGCGGCCCACCGAGGTCGTCGCCGCCGACTTCCTGCGCTACCCACTGCCCCGC
+TCACCCCACGTGGTCGTGGGCAACCTGCCGTTCCACCTCACCACCGCGATCCTGCGGCGA
+CTGCTGCACGGTCCGGGCTGGACCACGGCCGTGCTGCTCATGCAGTGGGAGGTGGCCCGC
+CGACGCGCCGCGGTGGGCGGCGCCACCATGATGACCGCCCAGTGGTGGCCGTGGTTCGAA
+TTCGGCCTTGCCCGAAAGGTTTCCGCGGCGAGCTTCACGCCGCGGCCCGCGGTCGACGCC
+GGACTGCTCACCATCACGCGCCGCAGCCGGCCGCTGGTCGACGTCGCGGACCGGGCGCGT
+>erm(38)_last630-1161
+GGACGCGGACACGGCATGGCGCAGATCCTG
+CAACGGTTGCCCACGCCGGTGCCCCGCACTTGGTTGCGGGCCAACGGGATAGCACCGAAC
+TCCCTGCCCCGCCAGTTGTCCGCGGCGCAGTGGGCGGCGCTGTTCGAGCAGACGCGTCTA
+ACTGGTGCCCAACGGGTCGATCGTCCACGCGATGTACAGCACGGCCGCGCTCACCGTCGC
+CGTGGTGGCGAAGTCGATCGCCCGGCTACGCACCACAAGCAGACCGGCCCGGTCGTCGGT
+CAGCGCCAACCGCAGCGCGGCCGCGACGCCGACGCCGATCCCGATGACCAGCGCACCGCG
+CCGCCAGTAACCCGCCACCACCAGGGCGAACGCCGCGATGAAGATCAGGCCGACCACCAG
+GATCGGCCATTGACCGGCGAACACCTTGCGGGCGAATTCCTTTGGCGTCACGCCAGTTTC
+GACTCTTCGGCTTCGACGACGTTGGTCAGCAGGAAGGCGCGGGTCAACGGGCCCACGCCA
=====================================
tests/data/test_isolate_14.fa
=====================================
@@ -0,0 +1,60 @@
+>test_isolate_14_contig_01
+GGTGGTAACTATAACCATCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGC
+ACGAATGGCGTAACGACTTCTCAACTGTCTCAACCATAGACTCGGCGAAATTGCACTACG
+AGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTT
+GGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCC
+AGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAA
+CCCTGAATCGGGTTTAGGGACAGTGCCTGGCGGGTAGTTTAACTGGGGCGGTTGCCTCCT
+AAAATGTAACGGAGGCGCCCA
+>test_isolate_14_contig_02
+GCTTCCTGCTCTGGCATGTCATCGGCGCGAATTCGTCGGACGACGGCTACATCCTGGGCA
+TGGCCCGAGTCGCCGACCACGCCGGCTACATGTCCAACTATTTCCGCTGGTTCGGCAGCC
+CGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACG
+CCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGT
+CGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGG
+CGGCGGCCATGGTCTTGCTGACCGCGTGGATGCCGTTCAACAACGGCCTGCGGCCGGAGG
+GCATCATCGCGCTCGGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACA
+GCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGC
+CCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCT
+TGGTGCGCCGTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCG
+GCACCGTCATCCTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCA
+CCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGTATACCGAGAACCTGCGTTACT
+ACTACCTCATCCTGCCCACCGTCGACGGTTCGCTGTCGCGGCGCTTCGGCTTTTTGATCA
+CCGCGCTATGCCTGTTCACCGCGGTGTTCATCATGTTGCGGCG
+>test_isolate_14_contig_03
+GCACGCCGCGGCGGCGCGAGCGCAGCCACAGCGCTTTTGCACGCAATTGCGCGGTCAGTT
+CCACACCCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGT
+AGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCG
+AAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGC
+CGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGA
+CGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAGGGGCCAAACCCCCATTCGTATCC
+CGTTCAGTCCTGGTTACCGGAGGAAACCGGGGGATCGGGCTGGCGATCGCACAGCGGCTG
+GCTGCCGACGGCCACAAGGTGGCCGTCACCCACCGTGGATCCGGAGCGCCAAAGGGG
+>test_isolate_14_contig_04
+CGGCGCTGATCGTCGGCGGTCACACTTTCGGTAAGACCCATGGCGCCGGCCCGGCCGATC
+TGGTCGGCCCCGAACCCGAGGCTGCTCCGCTGGAGCAGATGGGCTTGGGCTGGAAGAGCT
+CGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGA
+ACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGC
+TGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCA
+CCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCT
+CGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGG
+AATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCG
+TTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCC
+CTGCGGTCAGCCACGACCTCGTCGGCGAAGCCGAGATTGCCAGCCTTAAGAGCCAGATCC
+TGGCATCGGGATTGACTGTCTCACAGCTAGTTTCGACCGCATGGGCG
+>test_isolate_14_contig_05
+GACGTATGCGGGCGTTGATCATCGTCGACGTGCAGAACGACTTCTGCGAGGGTGGCTCGC
+TGGCGGTAACCGGTGGCGCCGCGCTGGCCCGCGCCATCAGCGACTACCTGGCCGAAGCGG
+CGGACTACCATCACGTCGTGGCAACCAAGGACTTCCACATCGACCCGGGTGACGACTTCT
+CCGGCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGGCG
+CGGACTTCCATCCCAGTCTGGACACGTCGGCAATCGAGGCGGTGTTCTACAAGGGTGCCT
+ACACCGGAGCGTACAGCGGCTTCGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATT
+GGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTG
+TGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACC
+TGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCGCTGGAGGAGATGCGCACCGCCA
+GCGTCGAGTTGGTTTGCAGCTCCTGATGGCACCGCC
+>test_isolate_14_contig_06
+ACGTTGATCAACATCCGGCCGGTGGTCGCCGCGATCAAGGAGTTCTTCGGCACCAGCCAG
+CTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTG
+TCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGCGACGTG
+CACCCGTCGCACTACGGCCGGATGTGCCCGATCGAAACCCCTGAGGGGC
=====================================
tests/data/test_isolate_14_1.fq
=====================================
@@ -0,0 +1,360 @@
+ at test_isolate_14_contig_01-20
+GGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGGGTTTAGGGACAGTGCCTGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-18
+TAACTATAACCATCCTAAGGTAGCGAAATTCCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAACGACTTCTCAACTGTCTCAACCATAGACTCGGCGAAATTGCACTACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-16
+AAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-14
+GGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-12
+ACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGGGTTTAGGGACAGTGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-10
+ACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-8
+ACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-6
+ACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGGGTTTAGGGACAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-4
+CCTTGTCGGGTAAGTTCCGACCTGCACGAATGGCGTAACGACTTCTCAACTGTCTCAACCATAGACTCGGCGAAATTGCACTACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-2
+GGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-50
+GCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-48
+GCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-46
+GTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCGGCACCGTCATCCTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCACCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-44
+GCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-42
+GTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGCGCCGTCATCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-40
+TGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGTCTTGCTGACCGCGTGGATGCCGTTCAACAACGGCCTGCGGCCGGAGGGCATCATCGCGCTCGGCTCGCTGGTCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-38
+CTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-36
+CTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGTCTTGCTGACCGCGTGGATGCCGTTCAACAACGGCCTGCGGCCGGAGGGCATCATCGCGCTCGGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-34
+CATCCTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCACCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGTATACCGAGAACCTGCGTTACTACTACCTCATCCTGCCCACCGTCGACGGTTCGCTGTC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-32
+CGTTCAACAACGGCCTGCGGCCGGAGGGCATCATCGCGCTCGGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-30
+GACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCACCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGTATACCGAGAACCTGCGTTACTACTACCTCATCCTGCCCACCGTCGACGGTTCGCTGTCGCGGCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-28
+ATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-26
+CTACATGTCCAACTATTTCCGCTGGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-24
+CGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGCGCCGTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-22
+GCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGTCTTGCTGACCGCGTGGATGCCGTTCAACAACGGCCTGCGGCCGGAGGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-20
+GGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAAC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-18
+TGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGCGCCGTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCGGCACCGTCATCCTGACCGTGGTGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-16
+GAAGCCACCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGTATACCGAGAACCTGCGTTACTACTACCTCATCCTGCCCACCGTCGACGGTTCGCTGTCGCGGCGCTTCGGCTTTTTGATCACCGCGCTATGCCTGTTCACCGCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-14
+TCGGCGCGAATTCGTCGGACGACGGCTACATCCTGGGCATGGCCCGAGTCGCCGACCACGCCGGCTACATGTCCAACTATTTCCGCTGGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-12
+TGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-10
+CGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-8
+AGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-6
+GGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-4
+CACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-2
+TGGATGCCGTTCAACAACGGCCTGCGGCCGGAGGGCATCATCGCGCTCGGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-30
+CACCCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-28
+TAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-26
+GGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-24
+CACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAGGGGCCAAA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-22
+ATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-20
+CCGCGGCGGCGCGAGCGCAGCCACAGCGCTTTTGCACGCAATTGCGCGGTCAGTTCCACACCCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-18
+GCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAGGGGCCAAACCCCCATTCGTATCCCGTTCAGTCCTGGTTACCGGAGGAAACCGGGGGATCGGGCTGGCGATCGCACAGCGGCTGGCTGCCGACGGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-16
+TAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-14
+TTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-12
+AGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-10
+CACACCCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCAC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-8
+CTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-6
+TGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-4
+GGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-2
+AATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-40
+GTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-38
+GATCTGGTCGGCCCCGAACCCGAGGCTGCTCCGCTGGAGCAGATGGGCTTGGGCTGGAAGAGCTCGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-36
+GGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCGATCTATGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-34
+CGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-32
+CTCGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-30
+ACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-28
+CCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCTCGTCGGCGAAGCCGAGATTGCCAGCCTTAAGAGCCAGATCCTGGCATCGGGATTGACTGTCTCACAGCTAGTTTCGACCGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-26
+CAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-24
+TGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-22
+GCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-20
+ACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-18
+CGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-16
+GACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-14
+TTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCTCGTCGGCGAAGCCGAGATTGCCAGCCTTAAGAGCCAGATCCTGGCATCGGGATTGACTGTCTCACAGCTAGTTTCGACCGCAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-12
+CAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-10
+CACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCTC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-8
+AGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-6
+CGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-4
+GCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-2
+GGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-30
+GGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-28
+AGCGTACAGCGGCTTCGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-26
+TGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCGCTGGAGGAGATGCGCACCGCCAGCGTCGAGTTGGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-24
+AACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-22
+CAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-20
+CGCGCCATCAGCGACTACCTGGCCGAAGCGGCGGACTACCATCACGTCGTGGCAACCAAGGACTTCCACATCGACCCGGGTGACGACTTCTCCGGCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-18
+TGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-16
+GGCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGGCGCGGACTTCCATCCCAGTCTGGACACGTCGGCAATCGAGGCGGTGTTCTACAAGGGTGCCTACACCGGAGCGTACAGCGGCTTCGAAGGAGTC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-14
+CGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGAC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-12
+GTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCGCTGGAGGAGATGCGCACCGCCAGCGTCGAGTTGGTTTGCAGCTCCTGATGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-10
+GTGTTCTACAAGGGTGCCTACACCGGAGCGTACAGCGGCTTCGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-8
+CTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-6
+CGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-4
+CCGCGCCATCAGCGACTACCTGGCCGAAGCGGCGGACTACCATCACGTCGTGGCAACCAAGGACTTCCACATCGACCCGGGTGACGACTTCTCCGGCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-2
+CATCGACCCGGGTGACGACTTCTCCGGCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGGCGCGGACTTCCATCCCAGTCTGGACACGTCGGCAATCGAGGCGGTGTTCTACAAGGGTGCCTACACCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-10
+ACGTTGATCAACATCCGGCCGGTGGTCGCCGCGATCAAGGAGTTCTTCGGCACCAGCCAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-8
+AGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGCGACGTGCACCCGTCGCACTACGGCCGGATGTGCCCGATC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-6
+AGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGCGACGTGCACCCGTCGCACTACGGCCGGATGTGCCCGATC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-4
+GCACCAGCCAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGCGACGTGCACCCGTCGCACTACGGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-2
+CAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGCGACGTGCACCCGTCGCACTACGGCCGGATGTGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
=====================================
tests/data/test_isolate_14_2.fq
=====================================
@@ -0,0 +1,360 @@
+ at test_isolate_14_contig_01-20
+GCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-18
+ACCTGCACGAATGGCGTAACGACTTCTCAACTGTCTCAACCATAGACTCGGCGAAATTGCACTACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-16
+GTCTCAACCATAGACTCGGCGAAATTGCACTACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-14
+AGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGGGTTTAGGGACAGTGCCTGGCGGGTAGTTTAACTGGGGCGGTTGCCTCCTAAAATGTAA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-12
+ACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-10
+GGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGGGTTTAGGGACAGTGCCTGGCGGGTAGTTTAAC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-8
+AACCATAGACTCGGCGAAATTGCACTACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-6
+GGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAGTCGTTGTTGAAATACCACTCTGATCGTATTGGGCATCTAACCTCGAACCCTGAATCGGGTTTAGGGACAGTGCCTGGCGGGTAGTTTAACTGGGGCGGTTGCCTCCTAAAATG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-4
+CTGTCTCAACCATAGACTCGGCGAAATTGCACTACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_01-2
+GGCGAAATTGCACTACGAGTAAAGATGCTCGTTACGCGCGGCAGGACGAAAAGACCCCGGGACCTTCACTACAACTTGGTATTGATGTTCGGTACGGTTTGTGTAGGATAGGTGGGAGACTGTGAAACCTCGACGCCAGTTGGGGCGGAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-50
+GGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-48
+GCCGGCTACATGTCCAACTATTTCCGCTGGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-46
+CTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGCGCCGTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCGGCACCGTCATCCTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-44
+CGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGCGCCGTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-42
+TGCGGCCGGAGGGCATCATCGCGCTCGGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-40
+GCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGTCTTGCTGACCGCGTGGATGCCGTTCAACAA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-38
+CAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-36
+GGTCTTGCTGACCGCGTGGATGCCGTTCAACAACGGCCTGCGGCCGGAGGGCATCATCGCGCTCGGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-34
+GCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCGGCACCGTCATCCTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCACCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGTATACCGAGAACCTGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-32
+GGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-30
+CTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCGGCACCGTCATCCTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCACCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGTATACCGAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-28
+CAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-26
+GACGACGGCTACATCCTGGGCATGGCCCGAGTCGCCGACCACGCCGGCTACATGTCCAACTATTTCCGCTGGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-24
+GGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGCGCCGTCATC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-22
+ATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGTCTTGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-20
+TGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGTCTTGCTGACCGCGTGGATGCCGTTC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-18
+CCGATGCTGCGGATCTTGGTGCGCCGTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCGGCACCGTCATCCTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCACCAGGGTTCGCGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-16
+CTGACCGTGGTGTTCGCCGACCAGACCCTGTCAACGGTGTTGGAAGCCACCAGGGTTCGCGCCAAAATCGGGCCGAGCCAGGCGTGGTATACCGAGAACCTGCGTTACTACTACCTCATCCTGCCCACCGTCGACGGTTCGCTGTCGCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-14
+TCGCCGACCACGCCGGCTACATGTCCAACTATTTCCGCTGGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-12
+ATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGCTGCGGATCTTGGTGCGCCGTCATCGCCTGGTCGGCACGTTGCCGTTGGTGTCGCCGATGCTGGCCGCCGGCACCGTCAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-10
+CCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGTCTTGCTGACCGCGTGGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-8
+GACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGGCCGCCGGGCTAGTGTGCTGGCTGCTGCTGTCGCGTGAGGTGCTGCCCCGCCTCGGGCCGGCGGTGGCGGCCAGCAAACCCGCCTACTGGGCGGCGGCCATGGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-6
+TGGCCCGAGTCGCCGACCACGCCGGCTACATGTCCAACTATTTCCGCTGGTTCGGCAGCCCGGAGGATCCCTTCGGCTGGTATTACAACCTGCTGGCGCTGATGACCCATGTCAGCGACGCCAGTCTGTGGATGCGCCTGCCAGACCTGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-4
+CATCATCGCGCTCGGCTCGCTGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_02-2
+TGGTCACCTATGTGCTGATCGAGCGGTCCATGCGGTACAGCCGGCTCACACCGGCGGCGCTGGCCGTCGTTACCGCCGCATTCACACTGGGTGTGCAGCCCACCGGCCTGATCGCGGTGGCCGCGCTGGTGGCCGGCGGCCGCCCGATGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-30
+CAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-28
+AATTGCGCGGTCAGTTCCACACCCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-26
+GCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCAC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-24
+TACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-22
+GCAATTGCGCGGTCAGTTCCACACCCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-20
+CCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-18
+TAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAGGGGCCAAACCCCCATTCGTATCCCGTTCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-16
+AAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-14
+GCAATTGCGCGGTCAGTTCCACACCCTGCGGCACGTACACGTCTTTATGTAGCGCGACATACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-12
+AGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAGGGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-10
+TACCTGCTGCGCAATTCGTAGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-8
+AGGGCGTCAATACACCCGCAGCCAGGGCCTCGCTGCCCAGAAAGGGATCCGTCATGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-6
+TAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAGGGGCCAAACCCCCATTCGTATCCCGTTCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-4
+GTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCCACAGCCACTGAAGGGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_03-2
+TGGTCGAAGTGTGCTGAGTCACACCGACAAACGTCACGAGCGTAACCCCAGTGCGAAAGTTCCCGCCGGAAATCGCAGCCACGTTACGCTCGTGGACATACCGATTTCGGCCCGGCCGCGGCGAGACGATAGGTTGTCGGGGTGACTGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-40
+CGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-38
+TCGTCGGCGGTCACACTTTCGGTAAGACCCATGGCGCCGGCCCGGCCGATCTGGTCGGCCCCGAACCCGAGGCTGCTCCGCTGGAGCAGATGGGCTTGGGCTGGAAGAGCTCGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-36
+CGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-34
+CGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-32
+TCGGCCCCGAACCCGAGGCTGCTCCGCTGGAGCAGATGGGCTTGGGCTGGAAGAGCTCGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-30
+GAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-28
+CGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCTCGTCGGCGAAGCCGAGAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-26
+TACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-24
+CCACTGACCTCTCGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-22
+AAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-20
+CTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-18
+GGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCTCGTCGGCGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-16
+AATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAGCCCTGCTGGCGCTTGGCAATACACCGCCAAGGACGGCGCCGGTGCCGGCACCATCCCGGACCCGTTCGGCGGGCCAGGGCGCTCCCCGACGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-14
+GTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCTCGTCGGCGAAGCCGAGATTGCCAGCCTTAAGAGCCAGAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-12
+CGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGATCCACCGAGACATGGGTCCCGTTGCGAGATACCTTGGGCCGCTGGTCCCCA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-10
+CCGTTGCGAGATACCTTGGGCCGCTGGTCCCCAAGCAGACCCTGCTGTGGCAGGATCCGGTCCCTGCGGTCAGCCACGACCTCGTCGGCGAAGCCGAGATTGCCAGCCTTAAGAGCCAGATCCTGGCATCGGGATTGACTGTCTCACAGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-8
+CGCTGGAGCAGATGGGCTTGGGCTGGAAGAGCTCGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-6
+CCCGTTCGGCGGGCCAGGGCGCTCCCCGACGATGCTGGCCACTGACCTCTCGCTGCGGGTGGATCCGATCTATGAGCGGATCACGCGTCGCTGGCTGGAACACCCCGAGGAATTGGCCGACGAGTTCGCCAAGGCCTGGTACAAGCTGAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-4
+GGGCTTGGGCTGGAAGAGCTCGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGGAGCTGACGAAGAG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_04-2
+CGCTGGAGCAGATGGGCTTGGGCTGGAAGAGCTCGTATGGCACCGGAACCGGTAAGGACGCGATCACCAGCGGCATCGAGGTCGTATGGACGAACACCCCGACGAAATGGGACAACAGTTTCCTCGAGATCCTGTACGGCTACGAGTGGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-30
+GCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCGCTGGAGGAGATGCGCACCGCCAGCGTCGAGTTGGTT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-28
+CACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-26
+GCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-24
+AGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-22
+GTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCGCTGGAGGAGATGCGCACCGCCAGCGTCGAGTTGGTTTGCAGCTCCTGAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-20
+CGACGTGCAGAACGACTTCTGCGAGGGTGGCTCGCTGGCGGTAACCGGTGGCGCCGCGCTGGCCCGCGCCATCAGCGACTACCTGGCCGAAGCGGCGGACTACCATCACGTCGTGGCAACCAAGGACTTCCACATCGACCCGGGTGACGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-18
+GAGCGTACAGCGGCTTCGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-16
+GCGACTACCTGGCCGAAGCGGCGGACTACCATCACGTCGTGGCAACCAAGGACTTCCACATCGACCCGGGTGACGACTTCTCCGGCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGGCGCGGACTTCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-14
+TCGGCAATCGAGGCGGTGTTCTACAAGGGTGCCTACACCGGAGCGTACAGCGGCTTCGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCAT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-12
+GCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-10
+GCGGTACTCCCGGCGCGGACTTCCATCCCAGTCTGGACACGTCGGCAATCGAGGCGGTGTTCTACAAGGGTGCCTACACCGGAGCGTACAGCGGCTTCGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-8
+TACAGCGGCTTCGAAGGAGTCGACGAGAACGGCACGCCACTGCTGAATTGGCTGCGGCAACGCGGCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-6
+GCGTCGATGAGGTCGATGTGGTCGGTATTGCCACCGATCATTGTGTGCGCCAGACGGCCGAGGACGCGGTACGCAATGGCTTGGCCACCAGGGTGCTGGTGGACCTGACAGCGGGTGTGTCGGCCGATACCACCGTCGCCGCGCTGGAGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-4
+AACCAAGGACTTCCACATCGACCCGGGTGACGACTTCTCCGGCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGGCGCGGACTTCCATCCCAGTCTGGACACGTCGGCAATCGAGGCGGTGTTCTACAA
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_05-2
+GCACACCGGACTATTCCTCGTCGTGGCCACCGCATTGCGTCAGCGGTACTCCCGGCGCGGACTTCCATCCCAGTCTGGACACGTCGGCAATCGAGGCGGTGTTCTACAAGGGTGCCTACACCGGAGCGTACAGCGGCTTCGAAGGAGTCG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-10
+TCTTCGGCACCAGCCAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGCGACGTGCACCCGTCGCACT
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-8
+CAAGGAGTTCTTCGGCACCAGCCAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGCGACGTGCACCC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-6
+GTCGCCGCGATCAAGGAGTTCTTCGGCACCAGCCAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGGAGGTCCGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-4
+TGATCAACATCCGGCCGGTGGTCGCCGCGATCAAGGAGTTCTTCGGCACCAGCCAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGC
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+ at test_isolate_14_contig_06-2
+GGCCGGTGGTCGCCGCGATCAAGGAGTTCTTCGGCACCAGCCAGCTGAGCCAATTCATGGACCAGAACAACCCGCTGTCGGGGTTGACCCACAAGCGCCGACTGTCGGCGCTGGGGCCCGGCGGTCTGTCACGTGAGCGTGCCGGGCTGG
++
+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
=====================================
tests/resfinder/cge/output/test_gene_result.md
=====================================
@@ -52,7 +52,7 @@ results from a ResFinder hit.
>>> assert(gene_result_blast["ref_acc"] == "KF986263")
>>> assert(gene_result_blast["identity"] == rf_dat_blast["perc_ident"])
>>> assert(gene_result_blast["alignment_length"] == rf_dat_blast["HSP_length"])
->>> assert(gene_result_blast["ref_seq_lenght"] == rf_dat_blast["sbjct_length"])
+>>> assert(gene_result_blast["ref_seq_length"] == rf_dat_blast["sbjct_length"])
>>> assert(gene_result_blast["ref_start_pos"] == rf_dat_blast["sbjct_start"])
>>> assert(gene_result_blast["ref_end_pos"] == rf_dat_blast["sbjct_end"])
>>> assert(gene_result_blast["query_id"] == rf_dat_blast["contig_name"])
=====================================
tests/resfinder/cge/output/test_phenotype_result.md
=====================================
@@ -64,7 +64,7 @@
... "name": "blaOXA-384",
... "identity": 98.2,
... "alignment_length": 90,
-... "ref_seq_lenght": 100,
+... "ref_seq_length": 100,
... "coverage": 90.0,
... "depth": 24.8,
... "ref_id": "blaOXA-384_1_KF986263",
@@ -84,7 +84,7 @@
... 'key': 'gyrA;;1;;CP073768.1',
... 'identity': 99.92,
... 'alignment_length': 2628,
-... 'ref_seq_lenght': 2628,
+... 'ref_seq_length': 2628,
... 'depth': 21,
... 'ref_start_pos': 1,
... 'ref_end_pos': 2628,
View it on GitLab: https://salsa.debian.org/med-team/resfinder/-/compare/f08fa09bf2ec23dbde4ffd1291925bea61b723c1...4674a866569cc31dcc5069ab6f4c9c6e9afa805c
--
View it on GitLab: https://salsa.debian.org/med-team/resfinder/-/compare/f08fa09bf2ec23dbde4ffd1291925bea61b723c1...4674a866569cc31dcc5069ab6f4c9c6e9afa805c
You're receiving this email because of your account on salsa.debian.org.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://alioth-lists.debian.net/pipermail/debian-med-commit/attachments/20240116/70cda2f0/attachment-0001.htm>
More information about the debian-med-commit
mailing list