diff --git a/.github/workflows/create_release.yml b/.github/workflows/create_release.yml
new file mode 100644
index 0000000..ecb96d5
--- /dev/null
+++ b/.github/workflows/create_release.yml
@@ -0,0 +1,30 @@
+name: Create Release
+
+on:
+ push:
+ branches:
+ - master
+ paths:
+ - 'version.txt'
+
+jobs:
+ release:
+ runs-on: ubuntu-latest
+
+ steps:
+ - name: Checkout code
+ uses: actions/checkout@v2
+
+ - name: Read version
+ id: get_version
+ run: |
+ VERSION=$(cat version.txt)
+ echo "VERSION=${VERSION}" >> $GITHUB_ENV
+
+ - name: Create bundle zip
+ run: zip -r bundle.zip *.wdl
+
+ - name: Create Release
+ run: gh release create ${{ env.VERSION }} jgi_assembly.wdl bundle.zip
+ env:
+ GITHUB_TOKEN: ${{ secrets.GITHUB_TOKEN }}
diff --git a/.gitmodules b/.gitmodules
new file mode 100644
index 0000000..e69de29
diff --git a/README.md b/README.md
index 8c605a4..c04f276 100644
--- a/README.md
+++ b/README.md
@@ -1,51 +1,40 @@
# The Metagenome Assembly Pipeline
## Summary
-This workflow is developed by Brian Foster at JGI and original from his [repo](https://gitlab.com/bfoster1/wf_templates/tree/master/templates). It take paired-end reads runs error correction by bbcms (BBTools). The clean reads are assembled by MetaSpades. After assembly, the reads are mapped back to contigs by bbmap (BBTools) for coverage information.
+This workflow is developed by Brian Foster at JGI and original from his [repo](https://gitlab.com/bfoster1/wf_templates/tree/master/templates). It takes in paired-end Illumina short reads or PacBio long reads.
-## Running Workflow in Cromwell
+In short reads, the workflow reformats the interleaved file into two FASTQ files for downstream tasks using bbcms (BBTools). The corrected reads are assembled using metaSPAdes. After assembly, the reads are mapped back to contigs by bbmap (BBTools) for coverage information. The `.wdl` (Workflow Description Language) file includes five tasks: *bbcms*, *assy*, *create_agp*, *read_mapping_pairs*, and *make_output*.
+
+In long reads, the workflow uses Flye for assembly, pbmm2 for alignment, Racon for polishing, and minimap2 for read mapping and coverage analysis. The :literal:`.wdl` (Workflow Description Language) file includes six tasks: *combine_fastq*, *assy*, *racon*, *format_assembly*, *map*, and *make_info_file*.
-Description of the files:
- - `.wdl` file: the WDL file for workflow definition
- - `.json` file: the example input for the workflow
- - `.conf` file: the conf file for running Cromwell.
- - `.sh` file: the shell script for running the example workflow
## The Docker image and Dockerfile can be found here
-[microbiomedata/bbtools:38.94](https://hub.docker.com/r/microbiomedata/bbtools)
+[microbiomedata/bbtools:39.03](https://hub.docker.com/r/microbiomedata/bbtools)
-[microbiomedata/spades:3.15.0](https://hub.docker.com/r/microbiomedata/spades)
+[microbiomedata/spades:4.0.0](https://hub.docker.com/r/microbiomedata/spades)
## Input files
-1. fastq (illumina paired-end interleaved fastq)
+1. The path to the input FASTQ file (Illumina paired-end interleaved FASTQ or PacBio paired-end interleaved FASTQ) (recommended: output of the Reads QC workflow).
-2. contig prefix for fasta header
+2. Project name: nmdc:XXXXXX
-3. output path
-
-4. input_interleaved (boolean)
-
-5. forwards reads fastq file (required value when input_interleaved is false, otherwise use [] )
+3. Memory (optional) e.g., `"jgi_metaAssembly.memory": "105G"`
-6. reverse reads fastq file (required value when input_interleaved is false, otherwise use [] )
+4. Threads (optional) e.g., `"jgi_metaAssembly.threads": "16"`
-7. memory (optional) ex: "jgi_metaASM.memory": "105G"
+5. Whether the input is short reads (boolean)
-8. threads (optional) ex: "jgi_metaASM.threads": "16"
```
{
- "jgi_metaASM.input_file":["/global/cfs/projectdirs/m3408/ficus/11809.7.220839.TCCTGAG-ACTGCAT.fastq.gz"],
- "jgi_metaASM.rename_contig_prefix":"503125_160870",
- "jgi_metaASM.outdir":"/global/cfs/projectdirs/m3408/aim2/metagenome/assembly/ficus/503125_160870",
- "jgi_metaASM.input_interleaved":true,
- "jgi_metaASM.input_fq1":[],
- "jgi_metaASM.input_fq2":[],
- "jgi_metaASM.memory": "105G",
- "jgi_metaASM.threads": "16"
+ "jgi_metaAssembly.input_files": ["https://portal.nersc.gov/project/m3408/test_data/smalltest.int.fastq.gz"],
+ "jgi_metaAssembly.proj": "nmdc:XXXXXX",
+ "jgi_metaAssembly.memory": "105G",
+ "jgi_metaAssembly.threads": "16",
+ "jgi_metaAssembly.shortRead": true
}
```
@@ -54,31 +43,33 @@ Description of the files:
Below is a part list of all output files. The main assembly contigs output is in final_assembly/assembly.contigs.fasta.
```
- ├── bbcms
- │ ├── berkeleylab-jgi-meta-60ade422cd4e
- │ ├── counts.metadata.json
- │ ├── input.corr.fastq.gz
- │ ├── input.corr.left.fastq.gz
- │ ├── input.corr.right.fastq.gz
- │ ├── readlen.txt
- │ └── unique31mer.txt
- ├── final_assembly
- │ ├── assembly.agp
- │ ├── assembly_contigs.fna
- │ ├── assembly_scaffolds.fna
- │ └── assembly_scaffolds.legend
- ├── mapping
- │ ├── covstats.txt (mapping_stats.txt)
- │ ├── pairedMapped.bam
- │ ├── pairedMapped.sam.gz
- │ ├── pairedMapped_sorted.bam
- │ └── pairedMapped_sorted.bam.bai
- └── spades3
- ├── assembly_graph.fastg
- ├── assembly_graph_with_scaffolds.gfa
- ├── contigs.fasta
- ├── contigs.paths
- ├── scaffolds.fasta
- └── scaffolds.paths
+# Short Reads
+ output/
+ ├── nmdc_XXXXXX_metaAsm.info
+ ├── nmdc_XXXXXX_covstats.txt
+ ├── nmdc_XXXXXX_contigs.fna
+ ├── nmdc_XXXXXX_bbcms.fastq.gz
+ ├── nmdc_XXXXXX_scaffolds.fna
+ ├── nmdc_XXXXXX_assembly.agp
+ ├── stats.json
+ ├── nmdc_XXXXXX_pairedMapped.sam.gz
+ └── nmdc_XXXXXX_pairedMapped_sorted.bam
+# Long Reads
+ output/
+ ├── nmdc_XXXXXX_assembly.legend
+ ├── nmdc_XXXXXX_contigs.fna
+ ├── nmdc_XXXXXX_pairedMapped_sorted.bam
+ ├── nmdc_XXXXXX_read_count_report.txt
+ ├── nmdc_XXXXXX_metaAsm.info
+ ├── nmdc_XXXXXX_summary.stats
+ ├── nmdc_XXXXXX_scaffolds.fna
+ ├── nmdc_XXXXXX_pairedMapped.sam.gz
+ ├── stats.json
+ ├── nmdc_XXXXXX_contigs.sam.stats
+ ├── nmdc_XXXXXX_contigs.sorted.bam.pileup.basecov
+ ├── nmdc_XXXXXX_assembly.agp
+ └── nmdc_XXXXXX_contigs.sorted.bam.pileup.out
```
+## Link to Doc Site
+Please refer [here](https://nmdc-workflow-documentation.readthedocs.io/en/latest/chapters/3_MetaGAssemly_index.html) for more information.
diff --git a/docs/index.rst b/docs/index.rst
index be63c31..afa01ae 100644
--- a/docs/index.rst
+++ b/docs/index.rst
@@ -1,31 +1,64 @@
-Metagenome Assembly Workflow (v1.0.2)
-========================================
+:github_url: https://github.com/microbiomedata/metaAssembly/blob/master/docs/index.rst
-.. image:: workflow_assembly.png
- :scale: 60%
+..
+ Note: The above `github_url` field is used to force the target of the "Edit on GitHub" link
+ to be the specified URL. That makes it so the link will work, regardless of the Sphinx
+ site the file is incorporated into. You can learn more about the `github_url` field at:
+ https://sphinx-rtd-theme.readthedocs.io/en/stable/configuring.html#confval-github_url
+
+Metagenome Assembly Workflow (v1.0.7)
+=====================================
+
+.. image:: lrassy_workflow2024.svg
:alt: Metagenome assembly workflow dependencies
-
+
Workflow Overview
-----------------
-This workflow takes in paired-end Illumina reads in interleaved format and performs error correction, then reformats the interleaved file into two FASTQ files for downstream tasks using bbcms (BBTools). The corrected reads are assembled using metaSPAdes. After assembly, the reads are mapped back to contigs by bbmap (BBTools) for coverage information. The .wdl (Workflow Description Language) file includes five tasks, *bbcms*, *assy*, *create_agp*, *read_mapping_pairs*, and *make_output*.
+This workflow takes in paired-end Illumina short reads or PacBio long reads.
+
+**Short Reads**:
-1. The *bbcms* task takes in interleaved FASTQ inputs and performs error correction and reformats the interleaved fastq into two output FASTQ files for paired-end reads for the next tasks.
-2. The *assy* task performs metaSPAdes assembly
-3. Contigs and Scaffolds (output of metaSPAdes) are consumed by the *create_agp* task to rename the FASTA header and generate an `AGP format `_ which describes the assembly
+In short reads, the workflow reformats the interleaved file into two FASTQ files for downstream tasks using bbcms (BBTools). The corrected reads are assembled using metaSPAdes. After assembly, the reads are mapped back to contigs by bbmap (BBTools) for coverage information. The `.wdl` (Workflow Description Language) file includes five tasks: *bbcms*, *assy*, *create_agp*, *read_mapping_pairs*, and *make_output*.
+
+1. The *bbcms* task takes in interleaved FASTQ inputs, performs error correction, and reformats the interleaved FASTQ into two output FASTQ files for paired-end reads for the next tasks.
+2. The *assy* task performs metaSPAdes assembly.
+3. Contigs and Scaffolds (output of metaSPAdes) are processed by the *create_agp* task to rename the FASTA header and generate an `AGP format `_ which describes the assembly.
4. The *read_mapping_pairs* task maps reads back to the final assembly to generate coverage information.
-5. The final *make_output* task adds all output files into the specified directory.
+5. The final *make_output* task collects all output files into the specified directory.
+
+**Long Reads**:
+
+In long reads, the workflow uses Flye for assembly, pbmm2 for alignment, Racon for polishing, and minimap2 for read mapping and coverage analysis. The :literal:`.wdl` (Workflow Description Language) file includes six tasks: *combine_fastq*, *assy*, *racon*, *format_assembly*, *map*, and *make_info_file*.
+
+1. The *combine_fastq* task combines the input FASTQ files into a single FASTQ file, which is used as input for polishing and mapping tasks.
+2. The *assy* task takes in the input FASTQ files and performs assembly using Flye.
+3. The *racon* task cleans up the assembled contigs through two rounds of error correction using :literal:`pbmm2` and :literal:`Racon`.
+4. The *format_assembly* task formats the polished assembly using BBTools' :literal:`fungalrelease.sh`, creating release-ready scaffolds and contigs, along with an `AGP format `_ file and a legend file that describes the assembly.
+5. The *map* task maps the input reads back to the final assembly using minimap2 to generate coverage data.
+6. The final *make_info_file* task produces a summary file documenting tool versions, parameters, memory usage, and Docker containers used throughout the workflow.
+
Workflow Availability
---------------------
-The workflow from GitHub uses all the listed docker images to run all third-party tools.
-The workflow is available in GitHub: https://github.com/microbiomedata/metaAssembly; the corresponding Docker images are available in DockerHub: https://hub.docker.com/r/microbiomedata/spades and https://hub.docker.com/r/microbiomedata/bbtools
+The workflow from GitHub uses all the listed Docker images to run all third-party tools.
+
+The workflow is available on GitHub: `https://github.com/microbiomedata/metaAssembly`
+
+The corresponding Docker images are available on DockerHub:
+
+- `https://hub.docker.com/r/microbiomedata/spades`
+- `https://hub.docker.com/r/microbiomedata/bbtools`
+- `https://hub.docker.com/r/staphb/flye`
+- `https://hub.docker.com/r/staphb/racon`
+- `https://hub.docker.com/r/staphb/minimap2`
+- `https://hub.docker.com/r/staphb/samtools`
Requirements for Execution
--------------------------
-(recommendations are in **bold**)
+(Recommendations are in **bold**)
- WDL-capable Workflow Execution Tool (**Cromwell**)
- Container Runtime that can load Docker images (**Docker v2.1.0.3 or higher**)
@@ -33,198 +66,184 @@ Requirements for Execution
Hardware Requirements
---------------------
-- Memory: >40 GB RAM
+**Memory: >40 GB RAM**
The memory requirement depends on the input complexity. Here is a simple estimation equation for the memory required based on kmers in the input file::
predicted_mem = (kmers * 2.962e-08 + 1.630e+01) * 1.1 (in GB)
.. note::
-
- The kmers variable for the equation above can be obtained using the kmercountmulti.sh script from BBTools.
- kmercountmulti.sh -k=31 in=your.read.fq.gz
+ The kmers variable for the equation above can be obtained using the `kmercountmulti.sh` script from BBTools.
+
+ Example command:
+
+ ::
+ kmercountmulti.sh -k=31 in=your.read.fq.gz
Workflow Dependencies
---------------------
-Third party software: (This is included in the Docker image.)
+Third-party software: (This is included in the Docker image.)
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-- `metaSPades v3.15.0 `_ (License: `GPLv2 `_)
-- `BBTools:38.94 `_ (License: `BSD-3-Clause-LBNL `_)
+- `metaSPAdes v4.0.0 `_ (License: `GPLv2 `_)
+- `BBTools v39.03 `_ (License: `BSD-3-Clause-LBNL `_)
Sample dataset(s)
-----------------
-- small dataset: `Ecoli 10x (287M) `_ . You can find input/output in the downloaded tar gz file.
+**Short Reads:**
-- large dataset: `Zymobiomics mock-community DNA control (22G) `_ . You can find input/output in the downloaded tar gz file.
+- Small dataset: `Ecoli 10x (287M) `_ (Input/output included in tar.gz file)
-Zymobiomics mock-community DNA control (`SRR7877884 `_); this original dataset is ~4 GB.
+- Large dataset: `Zymobiomics mock-community DNA control (22G) `_ (Input/output included in tar.gz file)
-For testing purposes and for the following examples, we used a 10% sub-sampling of the above dataset: (`SRR7877884-int-0.1.fastq.gz `_). This dataset is already interleaved.
+- Zymobiomics mock-community DNA control (`SRR7877884 `_); this `dataset `_ is has 6.7G bases.
+
+ - The non-interleaved raw fastq files are available as `R1 `_ and `R2 `
+ - The interleaved file is `here `_
+ - A 10% subset of the interleaved file is available as a quick dataset `here `
+
+**Long Reads:**
+
+Zymobiomics synthetic metagenome (`SRR13128014 `_) For testing we have subsampled the dataset, the original dataset is ~18GB.
Input
-----
-A JSON file containing the following information:
+A `JSON file `_ containing the following information:
+
+1. The path to the input FASTQ file (Illumina paired-end interleaved FASTQ or PacBio paired-end interleaved FASTQ) (recommended: output of the Reads QC workflow).
+2. Project name: nmdc:XXXXXX
+3. Memory (optional) e.g., `"jgi_metaAssembly.memory": "105G"`
+4. Threads (optional) e.g., `"jgi_metaAssembly.threads": "16"`
+5. Whether the input is short reads (boolean)
-1. the path to the input FASTQ file (Illumina paired-end interleaved FASTQ) (recommended the output of the Reads QC workflow.)
-2. the contig prefix for the FASTA header
-3. the output path
-4. input_interleaved (boolean)
-5. forwards reads fastq file (required value when input_interleaved is false, otherwise use [] )
-6. reverse reads fastq file (required value when input_interleaved is false, otherwise use [] )
-7. memory (optional) ex: “jgi_metaASM.memory”: “105G”
-8. threads (optional) ex: “jgi_metaASM.threads”: “16”
+Example input JSON for short reads::
+
+ {
+ "jgi_metaAssembly.input_files": ["https://portal.nersc.gov/project/m3408/test_data/smalltest.int.fastq.gz"],
+ "jgi_metaAssembly.proj": "nmdc:XXXXXX",
+ "jgi_metaAssembly.memory": "105G",
+ "jgi_metaAssembly.threads": "16",
+ "jgi_metaAssembly.shortRead": true
+ }
-An example input JSON file is shown below::
+Example input JSON for long reads::
{
- "jgi_metaASM.input_file":["/path/to/SRR7877884-int-0.1.fastq.gz "],
- "jgi_metaASM.rename_contig_prefix":"projectID",
- "jgi_metaASM.outdir":"/path/to/ SRR7877884-int-0.1_assembly",
- "jgi_metaASM.input_interleaved":true,
- "jgi_metaASM.input_fq1":[],
- "jgi_metaASM.input_fq2":[],
- "jgi_metaASM.memory": "105G",
- "jgi_metaASM.threads": "16"
+ "jgi_metaAssembly.input_files": ["/global/cfs/cdirs/m3408/www/test_data/SRR13128014.pacbio.subsample.ccs.fastq.gz"],
+ "jgi_metaAssembly.proj": "nmdc:XXXXXX",
+ "jgi_metaAssembly.memory": "105G",
+ "jgi_metaAssembly.threads": "16",
+ "jgi_metaAssembly.shortRead": false
}
Output
------
-The output directory will contain following files::
+The output directory will contain the following files for short reads::
+ output/
+ ├── nmdc_XXXXXX_metaAsm.info
+ ├── nmdc_XXXXXX_covstats.txt
+ ├── nmdc_XXXXXX_contigs.fna
+ ├── nmdc_XXXXXX_bbcms.fastq.gz
+ ├── nmdc_XXXXXX_scaffolds.fna
+ ├── nmdc_XXXXXX_assembly.agp
+ ├── stats.json
+ ├── nmdc_XXXXXX_pairedMapped.sam.gz
+ └── nmdc_XXXXXX_pairedMapped_sorted.bam
+
+The output directory will contain the following files for long reads::
output/
- ├── assembly.agp
- ├── assembly_contigs.fna
- ├── assembly_scaffolds.fna
- ├── covstats.txt
- ├── pairedMapped.sam.gz
- ├── pairedMapped_sorted.bam
- └── stats.json
-
-Part of an example output stats JSON file is shown below:
-
-```
-{
- "scaffolds": 58,
- "contigs": 58,
- "scaf_bp": 28406,
- "contig_bp": 28406,
- "gap_pct": 0.00000,
- "scaf_N50": 21,
- "scaf_L50": 536,
- "ctg_N50": 21,
- "ctg_L50": 536,
- "scaf_N90": 49,
- "scaf_L90": 317,
- "ctg_N90": 49,
- "ctg_L90": 317,
- "scaf_logsum": 22.158,
- "scaf_powsum": 2.245,
- "ctg_logsum": 22.158,
- "ctg_powsum": 2.245,
- "asm_score": 0.000,
- "scaf_max": 1117,
- "ctg_max": 1117,
- "scaf_n_gt50K": 0,
- "scaf_l_gt50K": 0,
- "scaf_pct_gt50K": 0.0,
- "gc_avg": 0.39129,
- "gc_std": 0.03033,
- "filename": "/global/cfs/cdirs/m3408/aim2/metagenome/assembly/cromwell-executions/jgi_metaASM/3342a6e8-7f78-40e6-a831-364dd2a47baa/call-create_agp/execution/assembly_scaffolds.fna"
-}
-```
+ ├── nmdc_XXXXXX_assembly.legend
+ ├── nmdc_XXXXXX_contigs.fna
+ ├── nmdc_XXXXXX_pairedMapped_sorted.bam
+ ├── nmdc_XXXXXX_read_count_report.txt
+ ├── nmdc_XXXXXX_metaAsm.info
+ ├── nmdc_XXXXXX_summary.stats
+ ├── nmdc_XXXXXX_scaffolds.fna
+ ├── nmdc_XXXXXX_pairedMapped.sam.gz
+ ├── stats.json
+ ├── nmdc_XXXXXX_contigs.sam.stats
+ ├── nmdc_XXXXXX_contigs.sorted.bam.pileup.basecov
+ ├── nmdc_XXXXXX_assembly.agp
+ └── nmdc_XXXXXX_contigs.sorted.bam.pileup.out
+
+Example output stats JSON file::
+
+ {
+ "scaffolds": 58,
+ "contigs": 58,
+ "scaf_bp": 28406,
+ "contig_bp": 28406,
+ "gap_pct": 0.00000,
+ "scaf_N50": 21,
+ "scaf_L50": 536,
+ "ctg_N50": 21,
+ "ctg_L50": 536,
+ "scaf_N90": 49,
+ "scaf_L90": 317,
+ "ctg_N90": 49,
+ "ctg_L90": 317,
+ "scaf_logsum": 22.158,
+ "scaf_powsum": 2.245,
+ "ctg_logsum": 22.158,
+ "ctg_powsum": 2.245,
+ "asm_score": 0.000,
+ "scaf_max": 1117,
+ "ctg_max": 1117,
+ "scaf_n_gt50K": 0,
+ "scaf_l_gt50K": 0,
+ "scaf_pct_gt50K": 0.0,
+ "gc_avg": 0.39129,
+ "gc_std": 0.03033
+ }
The table provides all of the output directories, files, and their descriptions.
-=================================================== ================================= ===============================================================
-Directory File Name Description
-=================================================== ================================= ===============================================================
-**bbcms** Error correction result directory
-bbcms/berkeleylab-jgi-meta-60ade422cd4e directory containing checking resource script
-bbcms/ counts.metadata.json bbcms commands and summary statistics in JSON format
-bbcms/ input.corr.fastq.gz error corrected reads in interleaved format.
-bbcms/ input.corr.left.fastq.gz error corrected forward reads
-bbcms/ input.corr.right.fastq.gz error corrected reverse reads
-bbcms/ rc cromwell script sbumit return code
-bbcms/ readlen.txt error corrected reads statistics
-bbcms/ resources.log resource checking log
-bbcms/ script Task run commands
-bbcms/ script.background Bash script to run script.submit
-bbcms/ script.submit cromwell submit commands
-bbcms/ stderr standard error where task writes error message to
-bbcms/ stderr.background standard error where bash script writes error message to
-bbcms/ stderr.log standard error from bbcms command
-bbcms/ stdout standard output where task writes error message to
-bbcms/ stdout.background standard output where bash script writes error message(s)
-bbcms/ stdout.log standard output from bbcms command
-bbcms/ unique31mer.txt the count of unique kmer, K=31
-**spades3** metaSPAdes assembly result directory
-spades3/K33 directory containing intermediate files from the run with K=33
-spades3/K55 directory containing intermediate files from the run with K=55
-spades3/K77 directory containing intermediate files from the run with K=77
-spades3/K99 directory containing intermediate files from the run with K=99
-spades3/K127 directory containing intermediate files from the run with K=127
-spades3/misc directory containing miscellaneous files
-spades3/tmp directory for temp files
-spades3/ assembly_graph.fastg metaSPAdes assembly graph in FASTG format
-spades3/ assembly_graph_with_scaffolds.gfa metaSPAdes assembly graph and scaffolds paths in GFA 1.0 format
-spades3/ before_rr.fasta contigs before repeat resolution
-spades3/ contigs.fasta metaSPAdes resulting contigs
-spades3/ contigs.paths paths in the assembly graph corresponding to contigs.fasta
-spades3/ dataset.info internal configuration file
-spades3/ first_pe_contigs.fasta preliminary contigs of iterative kmers assembly
-spades3/ input_dataset.yaml internal YAML data set file
-spades3/ params.txt information about SPAdes parameters in this run
-spades3/ scaffolds.fasta metaSPAdes resulting scaffolds
-spades3/ scaffolds.paths paths in the assembly graph corresponding to scaffolds.fasta
-spades3/ spades.log metaSPAdes log
-**final_assembly** create_agp task result directory
-final_assembly/berkeleylab-jgi-meta-60ade422cd4e directory containing checking resource script
-final_assembly/ assembly.agp an AGP format file describes the assembly
-final_assembly/ assembly_contigs.fna Final assembly contig fasta
-final_assembly/ assembly_scaffolds.fna Final assembly scaffolds fasta
-final_assembly/ assembly_scaffolds.legend name mapping file from spades node name to new name
-final_assembly/ rc cromwell script sbumit return code
-final_assembly/ resources.log resource checking log
-final_assembly/ script Task run commands
-final_assembly/ script.background Bash script to run script.submit
-final_assembly/ script.submit cromwell submit commands
-final_assembly/ stats.json assembly statistics in json format
-final_assembly/ stderr standard error where task writes error message to
-final_assembly/ stderr.background standard error where bash script writes error message to
-final_assembly/ stdout standard output where task writes error message to
-final_assembly/ stdout.background standard output where bash script writes error message to
-**mapping** maps reads back to the final assembly result directory
-mapping/ covstats.txt contigs coverage informaiton
-mapping/ mapping_stats.txt contigs coverage informaiton (same as covstats.txt)
-mapping/ pairedMapped.bam reads mapping back to the final assembly bam file
-mapping/ pairedMapped.sam.gz reads mapping back to the final assembly sam.gz file
-mapping/ pairedMapped_sorted.bam reads mapping back to the final assembly sorted bam file
-mapping/ pairedMapped_sorted.bam.bai reads mapping back to the final assembly sorted bam index file
-mapping/ rc cromwell script sbumit return code
-mapping/ resources.log resource checking log
-mapping/ script Task run commands
-mapping/ script.background Bash script to run script.submit
-mapping/ script.submit cromwell submit commands
-mapping/ stderr standard error where task writes error message to
-mapping/ stderr.background standard error where bash script writes error message to
-mapping/ stdout standard output where task writes error message to
-mapping/ stdout.background standard output where bash script writes error message to
-=================================================== ================================= ===============================================================
+
+=================================================== ===================================================== ===============================================================
+Directory File Name Description
+=================================================== ===================================================== ===============================================================
+**Short Reads** Short reads assembly output directory
+/make_info_file nmdc_XXXXXX_metaAsm.info Summary information about the short reads assembly process
+/finish_asm nmdc_XXXXXX_covstats.txt Coverage statistics for assembled contigs
+/finish_asm nmdc_XXXXXX_contigs.fna Final contig sequences in FASTA format
+/finish_asm nmdc_XXXXXX_bbcms.fastq.gz Error-corrected FASTQ file from bbcms
+/finish_asm nmdc_XXXXXX_scaffolds.fna Final scaffold sequences in FASTA format
+/finish_asm nmdc_XXXXXX_assembly.agp Assembly information in AGP format
+/finish_asm stats.json Assembly statistics in JSON format
+/finish_asm nmdc_XXXXXX_pairedMapped.sam.gz SAM file with reads mapped back to assembly
+/finish_asm nmdc_XXXXXX_pairedMapped_sorted.bam Sorted BAM file with reads mapped back to assembly
+
+**Long Reads** Long reads assembly output directory
+/finish_lrasm nmdc_XXXXXX_assembly.legend Mapping file from contig to scaffold names
+/finish_lrasm nmdc_XXXXXX_contigs.fna Final contig sequences in FASTA format
+/finish_lrasm nmdc_XXXXXX_pairedMapped_sorted.bam Sorted BAM file with reads mapped back to assembly
+/finish_lrasm nmdc_XXXXXX_read_count_report.txt Read count report for validation
+/make_info_file nmdc_XXXXXX_metaAsm.info Summary information about the long reads assembly process
+/finish_lrasm nmdc_XXXXXX_summary.stats Summary statistics for assembly
+/finish_lrasm nmdc_XXXXXX_scaffolds.fna Final scaffold sequences in FASTA format
+/finish_lrasm nmdc_XXXXXX_pairedMapped.sam.gz SAM file with reads mapped back to assembly
+/finish_lrasm stats.json Assembly statistics in JSON format
+/finish_lrasm nmdc_XXXXXX_contigs.sam.stats SAM file statistics for contigs
+/finish_lrasm nmdc_XXXXXX_contigs.sorted.bam.pileup.basecov Base coverage information for contigs
+/finish_lrasm nmdc_XXXXXX_assembly.agp Assembly information in AGP format
+/finish_lrasm nmdc_XXXXXX_contigs.sorted.bam.pileup.out BAM file pileup output for contigs
+=================================================== ===================================================== ===============================================================
+
Version History
---------------
-- 1.0.2 (release date **03/12/2021**; previous versions: 1.0.1)
+- 1.0.7 (release date **11/14/24**; previous versions: 1.0.6)
Point of contact
----------------
diff --git a/docs/lrassy_workflow2024.svg b/docs/lrassy_workflow2024.svg
new file mode 100644
index 0000000..e2bb5be
--- /dev/null
+++ b/docs/lrassy_workflow2024.svg
@@ -0,0 +1,234 @@
+
+
\ No newline at end of file
diff --git a/docs/workflow_assembly.png b/docs/workflow_assembly.png
index c555347..8a6022c 100644
Binary files a/docs/workflow_assembly.png and b/docs/workflow_assembly.png differ
diff --git a/imports.zip b/imports.zip
new file mode 100644
index 0000000..088cd40
Binary files /dev/null and b/imports.zip differ
diff --git a/input.json b/input.json
index 8ea0673..2b3ee69 100755
--- a/input.json
+++ b/input.json
@@ -1,8 +1,7 @@
{
- "jgi_metaASM.input_file":["/global/cfs/projectdirs/m3408/ficus/11809.7.220839.TCCTGAG-ACTGCAT.fastq.gz"],
- "jgi_metaASM.rename_contig_prefix":"503125_160870",
- "jgi_metaASM.input_interleaved":true,
- "jgi_metaASM.input_fq1":[],
- "jgi_metaASM.input_fq2":[],
- "jgi_metaASM.outdir":"/global/cfs/projectdirs/m3408/aim2/metagenome/assembly/ficus/503125_160870"
+ "jgi_metaAssembly.input_files":["/global/cfs/cdirs/m3408/www/test_data/SRR13128014.pacbio.subsample.ccs.fastq.gz"],
+ "jgi_metaAssembly.proj":"nmdc:503125_160870",
+ "jgi_metaAssembly.memory": "105G",
+ "jgi_metaAssembly.threads": "16",
+ "jgi_metaAssembly.shortRead": false
}
diff --git a/jgi_assembly.wdl b/jgi_assembly.wdl
old mode 100755
new mode 100644
index 94d80fa..2b9a163
--- a/jgi_assembly.wdl
+++ b/jgi_assembly.wdl
@@ -1,187 +1,174 @@
-workflow jgi_metaASM {
- String? outdir
- String? memory
- String? threads
- String input_file
- String proj
- String resource
- String informed_by
- String rename_contig_prefix="scaffold"
- Float uniquekmer=1000
- String? git_url="https://github.com/microbiomedata/mg_annotation/releases/tag/0.1"
- String? url_root="https://data.microbiomedata.org/data/"
- String bbtools_container="microbiomedata/bbtools:38.96"
- String spades_container="microbiomedata/spades:3.15.0"
- Boolean paired = true
-
- call stage {
- input:
- container=bbtools_container,
- input_file=input_file
+version 1.0
+import "shortReads_assembly.wdl" as srma
+import "make_interleaved_WDL/make_interleaved_reads.wdl" as int
+import "https://code.jgi.doe.gov/BFoster/jgi_meta_wdl/-/raw/bc7c4371ea0fa83355bada341ec353b9feb3eff2/metagenome_improved/metaflye.wdl" as lrma
+
+workflow jgi_metaAssembly {
+ input {
+ Boolean shortRead
+ String proj
+ String prefix=sub(proj, ":", "_")
+ # shortReads parameters
+ String? memory
+ String? threads
+ # longReads parameters
+ Array[String] input_files
+ String flye_container = "staphb/flye:2.9.2"
+ String flye_parameters = "--meta -o flye -t 32 --pacbio-hifi"
+ String smrtlink_container = "bryce911/smrtlink:12.0.0.177059"
+ String racon_container = "staphb/racon:1.4.20"
+ String minimap2_container = "staphb/minimap2:2.25"
+ String minimap2_parameters = "-a -x map-hifi -t 32"
+ String samtools_container = "staphb/samtools:1.18"
+ String bbtools_container = "microbiomedata/bbtools:39.03"
+ String spades_container="staphb/spades:4.0.0"
}
- call bbcms {
- input: input_files=stage.assembly_input, container=bbtools_container, memory=memory, paired = paired
- }
- call assy {
- input: infile1=bbcms.out1, infile2=bbcms.out2, container=spades_container, threads=threads, paired = paired
- }
- call create_agp {
- input: scaffolds_in=assy.out, container=bbtools_container, rename_contig_prefix = rename_contig_prefix, memory=memory
+ if (shortRead) {
+ if (length(input_files) > 1) {
+ call int.make_interleaved_reads {
+ input:
+ input_files = input_files,
+ container = bbtools_container
+ }
+ }
+ call srma.jgi_metaASM {
+ input:
+ memory = memory,
+ threads = threads,
+ input_file = if length(input_files) > 1 then make_interleaved_reads.interleaved_fastq else input_files[0],
+ proj = proj,
+ bbtools_container = bbtools_container,
+ spades_container = spades_container
+ }
+
}
- call read_mapping_pairs {
- input: reads=stage.assembly_input, ref=create_agp.outcontigs, container=bbtools_container, memory=memory, threads=threads, paired = paired
+ if (!shortRead) {
+ call lrma.metaflye {
+ input:
+ proj = proj,
+ input_fastq = input_files,
+ flye_container = flye_container,
+ flye_parameters = flye_parameters,
+ smrtlink_container = smrtlink_container,
+ racon_container = racon_container,
+ minimap2_container = minimap2_container,
+ minimap2_parameters = minimap2_parameters,
+ samtools_container = samtools_container,
+ bbtools_container = bbtools_container
+ }
+ call finish_lrasm {
+ input:
+ proj = proj,
+ prefix = prefix,
+ container = bbtools_container,
+ contigs = metaflye.final_contigs,
+ bam = metaflye.final_bam,
+ scaffolds = metaflye.final_scaffolds,
+ agp = metaflye.final_agp,
+ legend = metaflye.final_legend,
+ basecov = metaflye.final_basecov,
+ sam = metaflye.final_sam,
+ output_file = metaflye.final_output_file,
+ stats = metaflye.final_stats,
+ summary_stats = metaflye.final_summary_stats,
+ pileup_out = metaflye.final_pileup_out
+ }
}
-
- call finish_asm {
- input:
- proj=proj,
- start=stage.start,
- git_url=git_url,
- url_root=url_root,
- container="microbiomedata/workflowmeta:1.1.1",
- informed_by=informed_by,
- resource=resource,
- input_file=stage.assembly_input,
- fasta=create_agp.outcontigs,
- scaffold=create_agp.outscaffolds,
- agp=create_agp.outagp,
- bam=read_mapping_pairs.outbamfile,
- samgz=read_mapping_pairs.outsamfile,
- covstats=read_mapping_pairs.outcovfile,
- asmstats=create_agp.outstats
- }
-
output {
- File contig=finish_asm.outcontigs
- File scaffold=finish_asm.outscaffolds
- File agp=finish_asm.outagp
- File bam=finish_asm.outbam
- File samgz=finish_asm.outsamgz
- File covstats=finish_asm.outcovstats
- File asmstats=finish_asm.outasmstats
- File objects=finish_asm.objects
- }
-
- meta {
- author: "Chienchi Lo, B10, LANL"
- email: "chienchi@lanl.gov"
- version: "1.0.0"
+ # long reads output
+ File? lr_contigs = finish_lrasm.final_contigs
+ File? lr_bam = finish_lrasm.final_bam
+ File? lr_scaffolds = finish_lrasm.final_scaffolds
+ File? lr_agp = finish_lrasm.final_agp
+ File? lr_legend = finish_lrasm.final_legend
+ File? lr_basecov = finish_lrasm.final_basecov
+ File? lr_sam = finish_lrasm.final_sam
+ File? lr_output_file = finish_lrasm.final_output_file
+ File? lr_stats = finish_lrasm.final_stats
+ File? lr_summary_stats = finish_lrasm.final_summary_stats
+ File? lr_pileup_out = finish_lrasm.final_pileup_out
+ File? lr_asminfo = metaflye.asminfo
+
+ # short reads output
+ File? sr_contig=jgi_metaASM.contig
+ File? sr_scaffold=jgi_metaASM.scaffold
+ File? sr_agp=jgi_metaASM.agp
+ File? sr_bam=jgi_metaASM.bam
+ File? sr_samgz=jgi_metaASM.samgz
+ File? sr_covstats=jgi_metaASM.covstats
+ File? sr_asminfo=jgi_metaASM.asminfo
+ File? sr_bbcms_fq = jgi_metaASM.bbcms_fastq
+
+ #Both
+ File? stats = if (shortRead) then jgi_metaASM.asmstats else finish_lrasm.asm_stats
}
-
}
-task stage {
- String container
- String input_file
- String? memory = "4G"
- String target = "staged.fastq.gz"
- String output1 = "input.left.fastq.gz"
- String output2 = "input.right.fastq.gz"
-
- command <<<
- set -e
- if [ $( echo ${input_file}|egrep -c "https*:") -gt 0 ] ; then
- wget ${input_file} -O ${target}
- else
- ln ${input_file} ${target} || cp ${input_file} ${target}
- fi
-
- reformat.sh -Xmx${default="10G" memory} in=${target} out1=${output1} out2=${output2}
- # Capture the start time
- date --iso-8601=seconds > start.txt
-
- >>>
-
- output{
- Array[File] assembly_input = [output1, output2]
- String start = read_string("start.txt")
- }
- runtime {
- cpu: 2
- maxRetries: 1
- docker: container
- }
-}
-task finish_asm {
- Array[File] input_file
- File fasta
- File scaffold
- File? agp
- File bam
- File? samgz
- File? covstats
- File asmstats
- String container
- String git_url
- String informed_by
- String proj
- String prefix=sub(proj, ":", "_")
- String orig_prefix="scaffold"
- String sed="s/${orig_prefix}_/${proj}_/g"
- String resource
- String url_root
- String start
-
+task finish_lrasm {
+ input {
+ File contigs
+ File bam
+ File scaffolds
+ File agp
+ File legend
+ File basecov
+ File sam
+ File output_file
+ File stats
+ File summary_stats
+ File pileup_out
+ String container
+ String proj
+ String prefix
+ String orig_prefix="scaffold"
+ String sed="s/~{orig_prefix}_/~{proj}_/g"
+ # String start
+ }
command<<<
- set -e
+ set -oeu pipefail
end=`date --iso-8601=seconds`
- # ln ${fasta} ${prefix}_contigs.fna
- # ln ${scaffold} ${prefix}_scaffolds.fna
- # ln ${covstats} ${prefix}_covstats.txt
- # ln ${agp} ${prefix}_assembly.agp
+ ln -s ~{output_file} ~{prefix}_read_count_report.txt
+ ln -s ~{stats} ~{prefix}_contigs.sam.stats
+ ln -s ~{summary_stats} ~{prefix}_summary.stats
##RE-ID
- cat ${fasta} | sed ${sed} > ${prefix}_contigs.fna
- cat ${scaffold} | sed ${sed} > ${prefix}_scaffolds.fna
- cat ${covstats} | sed ${sed} > ${prefix}_covstats.txt
- cat ${agp} | sed ${sed} > ${prefix}_assembly.agp
-
- ## Bam file
- samtools view -h ${bam} | sed ${sed} | \
- samtools view -hb -o ${prefix}_pairedMapped_sorted.bam
- ## Sam.gz file
- samtools view -h ${samgz} | sed ${sed} | \
- gzip -c - > ${prefix}_pairedMapped.sam.gz
-
- # Remove an extra field from the stats
- cat ${asmstats} |jq 'del(.filename)' > stats.json
-
-
- /scripts/generate_object_json.py \
- --type "nmdc:MetagenomeAssembly" \
- --set metagenome_assembly_set \
- --part ${proj} \
- -p "name=Metagenome Assembly Activity for ${proj}" \
- was_informed_by=${informed_by} \
- started_at_time=${start} \
- ended_at_time=$end \
- execution_resource=${resource} \
- git_url=${git_url} \
- version="v1.0.3-beta" \
- --url ${url_root}${proj}/assembly/ \
- --extra stats.json \
- --inputs ${input_file[0]} ${input_file[1]} \
- --outputs \
- ${prefix}_contigs.fna "Final assembly contigs fasta" "Assembly Contigs" "Assembly contigs for ${proj}" \
- ${prefix}_scaffolds.fna "Final assembly scaffolds fasta" "Assembly Scaffolds" "Assembly scaffolds for ${proj}" \
- ${prefix}_covstats.txt "Assembled contigs coverage information" "Assembly Coverage Stats" "Coverage Stats for ${proj}" \
- ${prefix}_assembly.agp "An AGP format file that describes the assembly" "Assembly AGP" "AGP for ${proj}" \
- ${prefix}_pairedMapped_sorted.bam "Sorted bam file of reads mapping back to the final assembly" "Assembly Coverage BAM" "Sorted Bam for ${proj}"
+ cat ~{contigs} | sed ~{sed} > ~{prefix}_contigs.fna
+ cat ~{scaffolds} | sed ~{sed} > ~{prefix}_scaffolds.fna
+ cat ~{agp} | sed ~{sed} > ~{prefix}_assembly.agp
+ cat ~{legend} | sed ~{sed} > ~{prefix}_assembly.legend
+ cat ~{basecov} | sed ~{sed} > ~{prefix}_contigs.sorted.bam.pileup.basecov
+ cat ~{pileup_out} | sed ~{sed} > ~{prefix}_contigs.sorted.bam.pileup.out
+
+ ## Bam file
+ samtools view -h ~{bam} | sed ~{sed} | \
+ samtools view -hb -o ~{prefix}_pairedMapped_sorted.bam
+ ## Sam.gz file
+ samtools view -h ~{sam} | sed ~{sed} | \
+ gzip -c - > ~{prefix}_pairedMapped.sam.gz
+
+ # stats file
+ bbstats.sh format=8 in=~{scaffolds} out=stats.json
+ sed -i 's/l_gt50k/l_gt50K/g' stats.json
+ cat stats.json |jq 'del(.filename)' > stats.json
>>>
output {
- File outcontigs = "${prefix}_contigs.fna"
- File outscaffolds = "${prefix}_scaffolds.fna"
- File outagp = "${prefix}_assembly.agp"
- File outbam = "${prefix}_pairedMapped_sorted.bam"
- File outsamgz = "${prefix}_pairedMapped.sam.gz"
- File outcovstats = "${prefix}_covstats.txt"
- File outasmstats = "stats.json"
- File objects = "objects.json"
+ File final_contigs = "~{prefix}_contigs.fna"
+ File final_bam = "~{prefix}_pairedMapped_sorted.bam"
+ File final_scaffolds = "~{prefix}_scaffolds.fna"
+ File final_agp = "~{prefix}_assembly.agp"
+ File final_legend = "~{prefix}_assembly.legend"
+ File final_basecov = "~{prefix}_contigs.sorted.bam.pileup.basecov"
+ File final_sam = "~{prefix}_pairedMapped.sam.gz"
+ File final_output_file = "~{prefix}_read_count_report.txt"
+ File final_stats = "~{prefix}_contigs.sam.stats"
+ File final_summary_stats = "~{prefix}_summary.stats"
+ File final_pileup_out = "~{prefix}_contigs.sorted.bam.pileup.out"
+ File asm_stats = "stats.json"
}
runtime {
@@ -190,205 +177,3 @@ task finish_asm {
cpu: 1
}
}
-
-task make_output{
- String outdir
- File contigs
- File scaffolds
- File agp
- File bam
- File samgz
- File covstats
- File asmstats
- String contigs_name=basename(contigs)
- String scaffolds_name=basename(contigs)
- String agp_name=basename(contigs)
- String bam_name=basename(contigs)
- String samgz_name=basename(contigs)
- String covstats_name=basename(contigs)
- String asmstats_name=basename(contigs)
- String container
-
- command{
- if [ ! -z ${outdir} ]; then
- mkdir -p ${outdir}
- cp ${contigs} ${scaffolds} ${agp} ${bam} \
- ${samgz} ${covstats} ${asmstats} ${outdir}
- chmod 764 -R ${outdir}
- fi
- }
- runtime {
- docker: container
- memory: "1 GiB"
- cpu: 1
- }
- output{
- File? outcontigs = "${outdir}/${contigs_name}"
- File? outscaffolds = "${outdir}/${scaffolds_name}"
- File? outagp = "${outdir}/${agp_name}"
- File? outbam = "${outdir}/${bam_name}"
- File? outsamgz = "${outdir}/${samgz_name}"
- File? outcovstats = "${outdir}/${covstats_name}"
- File? outasmstats = "${outdir}/${asmstats_name}"
- }
-}
-
-task read_mapping_pairs{
- Array[File] reads
- File ref
- String container
- String? memory
- String? threads
- Boolean paired = true
- String bbmap_interleaved_flag = if paired then 'interleaved=true' else 'interleaved=false'
-
- String filename_unsorted="pairedMapped.bam"
- String filename_outsam="pairedMapped.sam.gz"
- String filename_sorted="pairedMapped_sorted.bam"
- String filename_sorted_idx="pairedMapped_sorted.bam.bai"
- String filename_bamscript="to_bam.sh"
- String filename_cov="covstats.txt"
- String system_cpu="$(grep \"model name\" /proc/cpuinfo | wc -l)"
- String jvm_threads=select_first([threads,system_cpu])
- runtime {
- docker: container
- memory: "120 GiB"
- cpu: 16
- maxRetries: 1
- }
- command{
- set -eo pipefail
- if [[ ${reads[0]} == *.gz ]] ; then
- cat ${sep=" " reads} > infile.fastq.gz
- export mapping_input="infile.fastq.gz"
- fi
- if [[ ${reads[0]} == *.fastq ]] ; then
- cat ${sep=" " reads} > infile.fastq
- export mapping_input="infile.fastq"
- fi
- bbmap.sh -Xmx${default="105G" memory} threads=${jvm_threads} nodisk=true ${bbmap_interleaved_flag} ambiguous=random in=$mapping_input ref=${ref} out=${filename_unsorted} covstats=${filename_cov} bamscript=${filename_bamscript}
- samtools sort -m100M -@ ${jvm_threads} ${filename_unsorted} -o ${filename_sorted}
- samtools index ${filename_sorted}
- reformat.sh -Xmx${default="105G" memory} in=${filename_unsorted} out=${filename_outsam} overwrite=true
- ln ${filename_cov} mapping_stats.txt
- rm $mapping_input
- }
- output{
- File outbamfile = filename_sorted
- File outbamfileidx = filename_sorted_idx
- File outcovfile = filename_cov
- File outsamfile = filename_outsam
- }
-}
-
-task create_agp {
- File scaffolds_in
- String? memory
- String container
- String rename_contig_prefix
- String prefix="assembly"
- String filename_contigs="${prefix}_contigs.fna"
- String filename_scaffolds="${prefix}_scaffolds.fna"
- String filename_agp="${prefix}.agp"
- String filename_legend="${prefix}_scaffolds.legend"
- runtime {
- docker: container
- memory: "120 GiB"
- cpu: 16
- }
- command{
- fungalrelease.sh -Xmx${default="105G" memory} in=${scaffolds_in} out=${filename_scaffolds} outc=${filename_contigs} agp=${filename_agp} legend=${filename_legend} mincontig=200 minscaf=200 sortscaffolds=t sortcontigs=t overwrite=t
- if [ "${rename_contig_prefix}" != "scaffold" ]; then
- sed -i 's/scaffold/${rename_contig_prefix}_scf/g' ${filename_contigs} ${filename_scaffolds} ${filename_agp} ${filename_legend}
- fi
- bbstats.sh format=8 in=${filename_scaffolds} out=stats.json
- sed -i 's/l_gt50k/l_gt50K/g' stats.json
- }
- output{
- File outcontigs = filename_contigs
- File outscaffolds = filename_scaffolds
- File outagp = filename_agp
- File outstats = "stats.json"
- File outlegend = filename_legend
- }
-}
-
-task assy {
- File infile1
- File infile2
- String container
- String? threads
- String outprefix="spades3"
- String filename_outfile="${outprefix}/scaffolds.fasta"
- String filename_spadeslog ="${outprefix}/spades.log"
- String system_cpu="$(grep \"model name\" /proc/cpuinfo | wc -l)"
- String spades_cpu=select_first([threads,system_cpu])
- Boolean paired = true
- runtime {
- docker: container
- memory: "120 GiB"
- cpu: 16
- }
- command{
- set -eo pipefail
- if ${paired}; then
- spades.py -m 2000 -o ${outprefix} --only-assembler -k 33,55,77,99,127 --meta -t ${spades_cpu} -1 ${infile1} -2 ${infile2}
- else
- spades.py -m 2000 -o ${outprefix} --only-assembler -k 33,55,77,99,127 -t ${spades_cpu} -s ${infile1}
- fi
- }
- output {
- File out = filename_outfile
- File outlog = filename_spadeslog
- }
-}
-
-task bbcms {
- Array[File] input_files
- String container
- String? memory
- Boolean paired = true
-
- String filename_outfile="input.corr.fastq.gz"
- String filename_outfile1="input.corr.left.fastq.gz"
- String filename_outfile2="input.corr.right.fastq.gz"
- String filename_readlen="readlen.txt"
- String filename_outlog="stdout.log"
- String filename_errlog="stderr.log"
- String filename_kmerfile="unique31mer.txt"
- String filename_counts="counts.metadata.json"
- runtime {
- docker: container
- memory: "120 GiB"
- cpu: 16
- }
-
- command {
- set -eo pipefail
- if file --mime -b ${input_files[0]} | grep gzip; then
- cat ${sep=" " input_files} > infile.fastq.gz
- export bbcms_input="infile.fastq.gz"
- fi
- if file --mime -b ${input_files[0]} | grep plain; then
- cat ${sep=" " input_files} > infile.fastq
- export bbcms_input="infile.fastq"
- fi
- bbcms.sh -Xmx${default="105G" memory} metadatafile=${filename_counts} mincount=2 highcountfraction=0.6 in=$bbcms_input out=${filename_outfile} > >(tee -a ${filename_outlog}) 2> >(tee -a ${filename_errlog} >&2) && grep Unique ${filename_errlog} | rev | cut -f 1 | rev > ${filename_kmerfile}
- if ${paired}; then
- reformat.sh -Xmx${default="105G" memory} in=${filename_outfile} out1=${filename_outfile1} out2=${filename_outfile2}
- fi
- readlength.sh -Xmx${default="105G" memory} in=${filename_outfile} out=${filename_readlen}
- rm $bbcms_input
- }
- output {
- File out = filename_outfile
- File out1 = if paired then filename_outfile1 else filename_outfile
- File out2 = if paired then filename_outfile2 else filename_outfile
- File outreadlen = filename_readlen
- File stdout = filename_outlog
- File stderr = filename_errlog
- File outcounts = filename_counts
- File outkmer = filename_kmerfile
- }
-}
-
diff --git a/labels.json b/labels.json
new file mode 100644
index 0000000..68140e3
--- /dev/null
+++ b/labels.json
@@ -0,0 +1,6 @@
+{
+ "sample_id": "503125_159845",
+ "project_id": "503125",
+ "pipeline": "metaAssembly",
+ "submitter": "vli"
+}
diff --git a/make_interleaved_WDL/make_interleaved_reads.wdl b/make_interleaved_WDL/make_interleaved_reads.wdl
index affdb7e..6c1fd5b 100644
--- a/make_interleaved_WDL/make_interleaved_reads.wdl
+++ b/make_interleaved_WDL/make_interleaved_reads.wdl
@@ -1,9 +1,16 @@
+# Make interleaved workflows for QC, etc.
+version 1.0
workflow make_interleaved_reads {
+ input{
Array[File] input_files
String output_file = "interleaved.fastq.gz"
-
+ String container="microbiomedata/bbtools:38.96"
+ }
call interleave_reads {
- input: input_files = input_files, output_file = output_file
+ input:
+ input_files = input_files,
+ output_file = output_file,
+ container = container
}
output {
File interleaved_fastq = interleave_reads.out_fastq
@@ -21,24 +28,28 @@ workflow make_interleaved_reads {
}
task interleave_reads{
-
+ input{
Array[File] input_files
String output_file = "interleaved.fastq.gz"
-
+ String container
+ }
command <<<
- if file --mime -b ${input_files[0]} | grep gzip; then
- paste <(gunzip -c ${input_files[0]} | paste - - - -) <(gunzip -c ${input_files[1]} | paste - - - -) | tr '\t' '\n' | gzip -c > ${output_file}
+ set -euo pipefail
+ if file --mime -b ~{input_files[0]} | grep gzip; then
+ paste <(gunzip -c ~{input_files[0]} | paste - - - -) <(gunzip -c ~{input_files[1]} | paste - - - -) | tr '\t' '\n' | gzip -c > ~{output_file}
else
- paste <(cat ${input_files[0]} | paste - - - -) <(cat ${input_files[1]} | paste - - - -) | tr '\t' '\n' | gzip -c > ${output_file}
+ paste <(cat ~{input_files[0]} | paste - - - -) <(cat ~{input_files[1]} | paste - - - -) | tr '\t' '\n' | gzip -c > ~{output_file}
fi
>>>
runtime {
- #docker: "my_image"
+ docker: container
+ memory: "1 GiB"
+ cpu: 1
}
output {
File out_fastq = output_file
}
-}
\ No newline at end of file
+}
diff --git a/shortReads_assembly.wdl b/shortReads_assembly.wdl
new file mode 100755
index 0000000..22c92bc
--- /dev/null
+++ b/shortReads_assembly.wdl
@@ -0,0 +1,543 @@
+version 1.0
+workflow jgi_metaASM {
+ input {
+ # String? outdir
+ String? memory
+ String? threads
+ String? input_file
+ String proj
+ String prefix=sub(proj, ":", "_")
+ String rename_contig_prefix="scaffold"
+ # Float uniquekmer=1000
+ String bbtools_container="microbiomedata/bbtools:39.03"
+ String spades_container="staphb/spades:4.0.0"
+ String workflowmeta_container="microbiomedata/workflowmeta:1.1.1"
+ Boolean paired = true
+ }
+
+ call stage {
+ input:
+ container=bbtools_container,
+ input_file=input_file
+ }
+
+ call bbcms {
+ input:
+ input_files=stage.assembly_input,
+ container=bbtools_container,
+ memory=memory,
+ paired = paired
+ }
+ call assy {
+ input:
+ infile1=bbcms.out1,
+ infile2=bbcms.out2,
+ container=spades_container,
+ threads=threads,
+ paired = paired
+ }
+ call create_agp {
+ input:
+ scaffolds_in=assy.out,
+ container=bbtools_container,
+ rename_contig_prefix = rename_contig_prefix,
+ memory=memory
+ }
+ call read_mapping_pairs {
+ input:
+ reads=stage.assembly_input,
+ ref=create_agp.outcontigs,
+ container=bbtools_container,
+ memory=memory,
+ threads=threads,
+ paired = paired
+ }
+
+ call make_info_file {
+ input:
+ bbcms_info= bbcms.outcounts,
+ assy_info = assy.outlog,
+ container=bbtools_container,
+ prefix=prefix
+ }
+ call finish_asm {
+ input:
+ proj=proj,
+ prefix=prefix,
+ # start=stage.start,
+ container=workflowmeta_container,
+ fasta=create_agp.outcontigs,
+ scaffold=create_agp.outscaffolds,
+ agp=create_agp.outagp,
+ bam=read_mapping_pairs.outbamfile,
+ samgz=read_mapping_pairs.outsamfile,
+ covstats=read_mapping_pairs.outcovfile,
+ asmstats=create_agp.outstats,
+ bbcms_fastq = bbcms.out
+ }
+
+ # call make_output{
+ # input:
+ # outdir= outdir,
+ # contigs = finish_asm.outcontigs,
+ # scaffolds = finish_asm.outscaffolds,
+ # agp=finish_asm.outagp,
+ # bam=finish_asm.outbam,
+ # samgz=finish_asm.outsamgz,
+ # covstats=finish_asm.outcovstats,
+ # asmstats=finish_asm.outasmstats,
+ # contigs_name=basename(finish_asm.outcontigs),
+ # scaffolds_name=basename(finish_asm.outcontigs),
+ # agp_name=basename(finish_asm.outcontigs),
+ # bam_name=basename(finish_asm.outcontigs),
+ # samgz_name=basename(finish_asm.outcontigs),
+ # covstats_name=basename(finish_asm.outcontigs),
+ # asmstats_name=basename(finish_asm.outcontigs),
+ # container = workflowmeta_container
+ # }
+
+ output {
+ File contig=finish_asm.outcontigs
+ File scaffold=finish_asm.outscaffolds
+ File agp=finish_asm.outagp
+ File bam=finish_asm.outbam
+ File samgz=finish_asm.outsamgz
+ File covstats=finish_asm.outcovstats
+ File asmstats=finish_asm.outasmstats
+ File asminfo=make_info_file.asminfo
+ File bbcms_fastq = finish_asm.outbbcms
+ }
+
+ meta {
+ author: "Chienchi Lo, B10, LANL"
+ email: "chienchi@lanl.gov"
+ version: "1.0.0"
+ }
+
+}
+
+task stage {
+ input {
+ String container
+ String? input_file
+ String memory = "4G"
+ String target = "staged.fastq.gz"
+ String output1 = "input.left.fastq.gz"
+ String output2 = "input.right.fastq.gz"
+ }
+
+ command <<<
+ set -euo pipefail
+ if [ $( echo ~{input_file}|egrep -c "https*:") -gt 0 ] ; then
+ wget ~{input_file} -O ~{target}
+ else
+ ln -s ~{input_file} ~{target} || cp ~{input_file} ~{target}
+ fi
+
+ reformat.sh \
+ ~{if (defined(memory)) then "-Xmx" + memory else "-Xmx10G" }\
+ in=~{target} \
+ out1=~{output1} \
+ out2=~{output2}
+ # Capture the start time
+ date --iso-8601=seconds > start.txt
+
+ >>>
+
+ output {
+ Array[File] assembly_input = [output1, output2]
+ String start = read_string("start.txt")
+ }
+ runtime {
+ cpu: 2
+ maxRetries: 1
+ docker: container
+ }
+}
+
+task make_info_file {
+ input{
+ File assy_info
+ File bbcms_info
+ String prefix
+ String container
+ }
+
+ command<<<
+ set -euo pipefail
+ bbtools_version=`grep BBToolsVer ~{bbcms_info}| awk '{print $2}' | sed -e 's/"//g' -e 's/,//' `
+ spades_version=`grep 'SPAdes version' ~{assy_info} | awk '{print $3}'`
+ echo -e "The workflow takes paired-end reads runs error correction by bbcms.sh (BBTools(1) version $bbtools_version)." > ~{prefix}_metaAsm.info
+ echo -e "The clean reads are assembled by metaSpades(2) version $spades_version with parameters, --only-assembler -k 33,55,77,99,127 --meta" >> ~{prefix}_metaAsm.info
+ echo -e "After assembly, Contigs and Scaffolds are consumed by the *create_agp* task to rename the FASTA header and generate an AGP format (https://www.ncbi.nlm.nih.gov/assembly/agp/AGP_Specification/) file which describes the assembly" >> ~{prefix}_metaAsm.info
+ echo -e "In the end, the reads are mapped back to contigs by bbmap (BBTools(1) version $bbtools_version) for coverage information." >> ~{prefix}_metaAsm.info
+
+ echo -e "\n(1) B. Bushnell: BBTools software package, http://bbtools.jgi.doe.gov/" >> ~{prefix}_metaAsm.info
+ echo -e "(2) Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017 May;27(5):824-834." >> ~{prefix}_metaAsm.info
+ >>>
+
+ output {
+ File asminfo = "~{prefix}_metaAsm.info"
+ }
+ runtime {
+ memory: "1 GiB"
+ cpu: 1
+ maxRetries: 1
+ docker: container
+ }
+}
+
+task finish_asm {
+ input {
+ File fasta
+ File scaffold
+ File? agp
+ File bam
+ File? samgz
+ File? covstats
+ File asmstats
+ File bbcms_fastq
+ String container
+ String proj
+ String prefix
+ String orig_prefix="scaffold"
+ String sed="s/~{orig_prefix}_/~{proj}_/g"
+ # String start
+ }
+
+ command<<<
+
+ set -euo pipefail
+ end=`date --iso-8601=seconds`
+ # ln -s ~{fasta} ~{prefix}_contigs.fna
+ # ln -s ~{scaffold} ~{prefix}_scaffolds.fna
+ # ln -s ~{covstats} ~{prefix}_covstats.txt
+ # ln -s ~{agp} ~{prefix}_assembly.agp
+
+ ##RE-ID
+ cat ~{fasta} | sed ~{sed} > ~{prefix}_contigs.fna
+ cat ~{scaffold} | sed ~{sed} > ~{prefix}_scaffolds.fna
+ cat ~{covstats} | sed ~{sed} > ~{prefix}_covstats.txt
+ cat ~{agp} | sed ~{sed} > ~{prefix}_assembly.agp
+ ln ~{bbcms_fastq} ~{prefix}_bbcms.fastq.gz || ln -s ~{bbcms_fastq} ~{prefix}_bbcms.fastq.gz
+
+ ## Bam file
+ samtools view -h ~{bam} | sed ~{sed} | \
+ samtools view -hb -o ~{prefix}_pairedMapped_sorted.bam
+ ## Sam.gz file
+ samtools view -h ~{samgz} | sed ~{sed} | \
+ gzip -c - > ~{prefix}_pairedMapped.sam.gz
+
+ # Remove an extra field from the stats
+ cat ~{asmstats} |jq 'del(.filename)' > stats.json
+
+ >>>
+ output {
+ File outcontigs = "~{prefix}_contigs.fna"
+ File outscaffolds = "~{prefix}_scaffolds.fna"
+ File outagp = "~{prefix}_assembly.agp"
+ File outbam = "~{prefix}_pairedMapped_sorted.bam"
+ File outsamgz = "~{prefix}_pairedMapped.sam.gz"
+ File outcovstats = "~{prefix}_covstats.txt"
+ File outasmstats = "stats.json"
+ File outbbcms = "~{prefix}_bbcms.fastq.gz"
+ }
+
+ runtime {
+ docker: container
+ memory: "1 GiB"
+ cpu: 1
+ }
+}
+
+
+task read_mapping_pairs{
+ input {
+ Array[File] reads
+ File ref
+ String container
+ String? memory
+ String? threads
+ Boolean paired = true
+ String bbmap_interleaved_flag = if paired then 'interleaved=true' else 'interleaved=false'
+
+ String filename_unsorted="pairedMapped.bam"
+ String filename_outsam="pairedMapped.sam.gz"
+ String filename_sorted="pairedMapped_sorted.bam"
+ String filename_sorted_idx="pairedMapped_sorted.bam.bai"
+ String filename_bamscript="to_bam.sh"
+ String filename_cov="covstats.txt"
+ String system_cpu="$(grep \"model name\" /proc/cpuinfo | wc -l)"
+ String jvm_threads=select_first([threads,system_cpu])
+ }
+ runtime {
+ docker: container
+ memory: "120 GiB"
+ cpu: 16
+ maxRetries: 1
+ }
+ command<<<
+ set -euo pipefail
+ if [[ ~{reads[0]} == *.gz ]] ; then
+ cat ~{sep=" " reads} > infile.fastq.gz
+ export mapping_input="infile.fastq.gz"
+ fi
+ if [[ ~{reads[0]} == *.fastq ]] ; then
+ cat ~{sep=" " reads} > infile.fastq
+ export mapping_input="infile.fastq"
+ fi
+
+ bbmap.sh \
+ ~{if (defined(memory)) then "-Xmx" + memory else "-Xmx105G" } \
+ threads=~{jvm_threads} \
+ nodisk=true \
+ ~{bbmap_interleaved_flag} \
+ ambiguous=random \
+ in=$mapping_input \
+ ref=~{ref} \
+ out=~{filename_unsorted} \
+ covstats=~{filename_cov} \
+ bamscript=~{filename_bamscript}
+
+ samtools sort \
+ -m100M \
+ -@ \
+ ~{jvm_threads} \
+ ~{filename_unsorted} \
+ -o ~{filename_sorted}
+
+ samtools index ~{filename_sorted}
+
+ reformat.sh \
+ ~{if (defined(memory)) then "-Xmx" + memory else "-Xmx105G" } \
+ in=~{filename_unsorted} \
+ out=~{filename_outsam} \
+ overwrite=true
+
+ ln -s ~{filename_cov} mapping_stats.txt
+ rm $mapping_input
+
+ >>>
+ output {
+ File outbamfile = filename_sorted
+ File outbamfileidx = filename_sorted_idx
+ File outcovfile = filename_cov
+ File outsamfile = filename_outsam
+ }
+}
+
+task create_agp {
+ input {
+ File scaffolds_in
+ String? memory
+ String container
+ String rename_contig_prefix
+ String prefix="assembly"
+ String filename_contigs="~{prefix}_contigs.fna"
+ String filename_scaffolds="~{prefix}_scaffolds.fna"
+ String filename_agp="~{prefix}.agp"
+ String filename_legend="~{prefix}_scaffolds.legend"
+ }
+ runtime {
+ docker: container
+ memory: "120 GiB"
+ cpu: 16
+ }
+ command<<<
+ set -euo pipefail
+ fungalrelease.sh \
+ ~{if (defined(memory)) then "-Xmx" + memory else "-Xmx105G" } \
+ in=~{scaffolds_in} \
+ out=~{filename_scaffolds} \
+ outc=~{filename_contigs} \
+ agp=~{filename_agp} \
+ legend=~{filename_legend} \
+ mincontig=200 \
+ minscaf=200 \
+ sortscaffolds=t \
+ sortcontigs=t \
+ overwrite=t
+
+ if [ "~{rename_contig_prefix}" != "scaffold" ]; then
+ sed -i 's/scaffold/~{rename_contig_prefix}_scf/g' \
+ ~{filename_contigs} ~{filename_scaffolds} ~{filename_agp} ~{filename_legend}
+ fi
+ bbstats.sh format=8 in=~{filename_scaffolds} out=stats.json
+ sed -i 's/l_gt50k/l_gt50K/g' stats.json
+
+ >>>
+
+ output {
+ File outcontigs = filename_contigs
+ File outscaffolds = filename_scaffolds
+ File outagp = filename_agp
+ File outstats = "stats.json"
+ File outlegend = filename_legend
+ }
+}
+
+task assy {
+ input{
+ File infile1
+ File infile2
+ String container
+ String? threads
+ String outprefix="spades3"
+ String filename_outfile="~{outprefix}/scaffolds.fasta"
+ String filename_spadeslog ="~{outprefix}/spades.log"
+ String system_cpu="$(grep \"model name\" /proc/cpuinfo | wc -l)"
+ String spades_cpu=select_first([threads,system_cpu])
+ Boolean paired = true
+ }
+ runtime {
+ docker: container
+ memory: "120 GiB"
+ cpu: 16
+ }
+ command <<<
+ set -euo pipefail
+ if ~{paired}; then
+ spades.py \
+ -m 2000 \
+ -o ~{outprefix} \
+ --only-assembler \
+ -k 33,55,77,99,127 \
+ --meta \
+ -t ~{spades_cpu} \
+ -1 ~{infile1} \
+ -2 ~{infile2}
+ else
+ spades.py \
+ -m 2000 \
+ -o ~{outprefix} \
+ --only-assembler \
+ -k 33,55,77,99,127 \
+ -t ~{spades_cpu} \
+ -s ~{infile1}
+ fi
+ >>>
+
+ output {
+ File out = filename_outfile
+ File outlog = filename_spadeslog
+ }
+}
+
+task bbcms {
+ input{
+ Array[File] input_files
+ String container
+ String? memory
+ Boolean paired = true
+ String filename_outfile="input.corr.fastq.gz"
+ String filename_outfile1="input.corr.left.fastq.gz"
+ String filename_outfile2="input.corr.right.fastq.gz"
+ String filename_readlen="readlen.txt"
+ String filename_outlog="stdout.log"
+ String filename_errlog="stderr.log"
+ String filename_kmerfile="unique31mer.txt"
+ String filename_counts="counts.metadata.json"
+ }
+ runtime {
+ docker: container
+ memory: "120 GiB"
+ cpu: 16
+ }
+
+ command<<<
+ set -euo pipefail
+ if file --mime -b ~{input_files[0]} | grep gzip; then
+ cat ~{sep=" " input_files} > infile.fastq.gz
+ export bbcms_input="infile.fastq.gz"
+ fi
+
+ if file --mime -b ~{input_files[0]} | grep plain; then
+ cat ~{sep=" " input_files} > infile.fastq
+ export bbcms_input="infile.fastq"
+ fi
+
+ bbcms.sh \
+ ~{if (defined(memory)) then "-Xmx" + memory else "-Xmx105G" } \
+ metadatafile=~{filename_counts} \
+ mincount=2 \
+ highcountfraction=0.6 \
+ in=$bbcms_input \
+ out=~{filename_outfile} \
+ > >(tee -a ~{filename_outlog}) \
+ 2> >(tee -a ~{filename_errlog} >&2) \
+ && grep Unique ~{filename_errlog} \
+ | rev | cut -f 1 | rev \
+ > ~{filename_kmerfile}
+
+ if ~{paired}; then
+ reformat.sh \
+ ~{if (defined(memory)) then "-Xmx" + memory else "-Xmx105G" } \
+ in=~{filename_outfile} \
+ out1=~{filename_outfile1} \
+ out2=~{filename_outfile2}
+ fi
+
+ readlength.sh \
+ ~{if (defined(memory)) then "-Xmx" + memory else "-Xmx105G" } \
+ in=~{filename_outfile} \
+ out=~{filename_readlen}
+
+ rm $bbcms_input
+
+ >>>
+
+ output {
+ File out = filename_outfile
+ File out1 = if paired then filename_outfile1 else filename_outfile
+ File out2 = if paired then filename_outfile2 else filename_outfile
+ File outreadlen = filename_readlen
+ File stdout = filename_outlog
+ File stderr = filename_errlog
+ File outcounts = filename_counts
+ File outkmer = filename_kmerfile
+
+ }
+}
+
+# task make_output{
+# input{
+# String? outdir
+# File contigs
+# File scaffolds
+# File agp
+# File bam
+# File samgz
+# File covstats
+# File asmstats
+# String contigs_name=basename(contigs)
+# String scaffolds_name=basename(contigs)
+# String agp_name=basename(contigs)
+# String bam_name=basename(contigs)
+# String samgz_name=basename(contigs)
+# String covstats_name=basename(contigs)
+# String asmstats_name=basename(contigs)
+# String container
+# }
+# command{
+# if [ ! -z ~{outdir} ]; then
+# mkdir -p ~{outdir}
+# cp ~{contigs} ~{scaffolds} ~{agp} ~{bam} \
+# ~{samgz} ~{covstats} ~{asmstats} ~{outdir}
+# chmod 764 -R ~{outdir}
+# fi
+# }
+# runtime {
+# docker: container
+# memory: "1 GiB"
+# cpu: 1
+# }
+# output{
+# File? outcontigs = "~{outdir}/~{contigs_name}"
+# File? outscaffolds = "~{outdir}/~{scaffolds_name}"
+# File? outagp = "~{outdir}/~{agp_name}"
+# File? outbam = "~{outdir}/~{bam_name}"
+# File? outsamgz = "~{outdir}/~{samgz_name}"
+# File? outcovstats = "~{outdir}/~{covstats_name}"
+# File? outasmstats = "~{outdir}/~{asmstats_name}"
+# }
+# }
\ No newline at end of file
diff --git a/submit_curl.sh b/submit_curl.sh
new file mode 100755
index 0000000..5676e73
--- /dev/null
+++ b/submit_curl.sh
@@ -0,0 +1,15 @@
+#!/bin/bash
+#
+
+
+# zip bundle.zip *wdl
+WD=/global/cfs/cdirs/m3408/aim2/dev/kli_training/mg_annotation/
+
+curl --netrc -X POST "https://nmdc-cromwell.freeddns.org:8443/api/workflows/v1" \
+ -H "accept: application/json" \
+ -H "Content-Type: multipart/form-data" \
+ -F "workflowSource=@jgi_assembly.wdl" \
+ -F "workflowInputs=@input.json;type=application/json" \
+ -F "workflowDependencies=@imports.zip" \
+ -F "labels=@labels.json;type=application/json"
+
diff --git a/test_output/choose_wdl.json b/test_output/choose_wdl.json
new file mode 100644
index 0000000..65b61ee
--- /dev/null
+++ b/test_output/choose_wdl.json
@@ -0,0 +1,5 @@
+{
+ "jgi_metaAssembly.shortRead": false,
+ "jgi_metaAssembly.input_fastq":["/expanse/projects/nmdc/edge_app/test/kli/metaAssembly/test_output/SRR13128014.pacbio.subsample.ccs.fastq.gz"],
+ "jgi_metaAssembly.proj":"nmdc:xxxxxxx"
+}
\ No newline at end of file
diff --git a/test_output/long_reads.json b/test_output/long_reads.json
new file mode 100644
index 0000000..1abbbca
--- /dev/null
+++ b/test_output/long_reads.json
@@ -0,0 +1,14 @@
+{
+ "metaflye.input_fastq":["SRR13128014.pacbio.subsample.ccs.fastq.gz"],
+ "metaflye.flye_container":"staphb/flye:2.9.2",
+ "metaflye.flye_parameters":"--meta -o flye -t 32 --pacbio-hifi",
+ "metaflye.smrtlink_container":"bryce911/smrtlink:12.0.0.177059",
+ "metaflye.racon_container":"staphb/racon:1.4.20",
+ "metaflye.minimap2_container":"staphb/minimap2:2.25",
+ "metaflye.minimap2_parameters":"-a -x map-hifi -t 32",
+ "metaflye.samtools_container":"staphb/samtools:1.18",
+ "metaflye.bbtools_container":"microbiomedata/bbtools:38.96",
+ "metaflye.bc_container":"doejgi/bc:1.07.1",
+ "metaflye.combined_fastq_filename":"reads.fastq.gz",
+ "metaflye.flye_memory":"115G"
+}
\ No newline at end of file
diff --git a/test_output/test.1.fastq b/test_output/test.1.fastq
new file mode 100755
index 0000000..916c510
--- /dev/null
+++ b/test_output/test.1.fastq
@@ -0,0 +1,41472 @@
+@:114:A016U:1:1:16853:1322 1:N:0:1
+AGTGCCTGCAAAGAGAGGGTGCAATCCTGCTGCAGTAAGGCCAGCAGCTTACGGTCAATTTTATCTAACATAGCCCTTCCACAGAGAATTTTTTTCTCTTTTAATTTTATTTTAAAGGTGAAATCACGACAATTTATTCTGTGCGGTTTAA
++
+CCCFFFFFHGGHHGIIJJJHIIIIJJIIIIIJIJJGIIJJJJIJJJIIJJJJIJIJIGJJJJJJJGHHHHHHFFCDFEEEDEDBDDCCDD@CCBDDCDDDDCDDCECDEDDEDCC>AC@CCCCCDCCDDBBBDCCDEDDEEDC>38(6>(6>>ADAA?=AAA>:38><><<>55>8:>>AA>A:>>9>;;?>9AAA:028<>5&25343+4(5&&09>A((4:><9
+@:114:A016U:1:1:17297:1520 1:N:0:1
+TTTCCATCGCTTTCTTTACTGGGGTGTGATTATTCGCGACGTGCCCTTCCTGCCCCACCCGATCTCTCACTGTTAGACCCGCTTTGTTATGCCGCAACGACCTGGCGAATGACCATCATTGCCGGCCTTCCCGTTGAATATAACGATCGCT
++
+CCCFFFFFHHGHHJJJJJJJJJJGDFFHHGIJJJIIJJJJGIIIJHHHHHFFFFFDDDDDDBBDDDEDCDDD>AACDDDDDDBDBCCCDEDDD@>@DDBDBBBABDDBBB?CDDACA@:>>>C@7BDDCCCCBB8?8:@CCCA:(282&+2
+@:114:A016U:1:1:13758:1594 1:N:0:1
+CCGCAATGGCTGACCACAGTAGCGCCATTGCGTCGCGTTCCAGCTCTTTTTCACCATCCTGGCGGATATGCTCATGAATTGCCATCGCCCGGGAGGGGAGTCGGTCTTCATCTATTTCTATTTTTTTGCCGCGCTCTTTTTCTTCGCTCTC
++
+CC@FFDFFFDHHHGGGIJJGIFDFHIDEHIII6DFIIHHEHGIJJGCEEEE?>BDEDECECCB@@8>;CCCCA:@A:>AC@@>CC?BBB7<@@-55<<)0>@B@AAD@A@B><@AB@DB@@DCBBC
+@:114:A016U:1:1:9947:1618 1:N:0:1
+AGGCTAACCGCCGTTTGCATTCCCGCGCCCCGCTGTGTGAACGCCTGTCACCTGTACTCAACGGCTTACAGAATTTAACCCTGCTACGTGATATCGAATTGCGGGTGTATGACACTGATGATGAAGAGAATCATCAGGAGTTTACCTGCTG
++
+<@@B?D?BHF?DFCGHFHGHGGCHC?@A@GBFAABAEHECD;?/9;=A553;@?>@55@@>ABB>8>>@A>:CA:A:@::8<<<:@@BDC@?C@@0<)05:>C@@CC>AADD@>:>:@:+((:++>@::A@BAABB(55>ACCDDDDDCCB9:>A@ACDD<>4>A@@9&5
+@:114:A016U:1:1:11855:1707 1:N:0:1
+TAAGTATAAGTATATCGATATAGATCAGTGTGATTCGCTATGAGCGACTTCGGCATTGCCATGATTCTGGCAGCCAACTTCTTCACAACGTTGCATCATGTAACCCAGCGCGACCAGCACTACTGTTGATAACATCACGCTGGTGGTGGTG
++
+BBCFDEFFHHHHHJFIGGGIIIJJIJJIHIIJJJJIJIIJIJJJJJJJJJJIJJIIJJJJJHHHHHHHFFFFFDDEEDDDCDDDECDDDD?ADDDECCCADDCCAA@D@BBD>BBDDBDBCCDDDDCDCDDCCD@CC@BDBDD?AB5?B9<
+@:114:A016U:1:1:7770:1708 1:N:0:1
+CAATCATATCCTTGATGAAGAAGCGATTAAATTTCTGGAGAAATTAGATGCCGAGAGAATTATTGTTTTGGCACCTTATCACATCAGTAAACTAAAAGCTAAAGCGCCTATTTATTTTGTTAGCCGCAAAGAAAGTATCAAAAATCTTCTT
++
+B@@FFFFFHGFHHJIJIIIIJIIJIIJIGIJIJJJIJJIJJIJJJJIJJJJJIGEHIJGJJIIIIHHHHGEFDFCCEECECCDDDDDADCDDCDCDCCDCDDCCDDBBBDDC@DEDEEDACCDCD099A<<&+4(:@>C755592(:>BCCC34>C0B?CC@:>@CCD:@C:AC::>>@BB589A8)50)5884@CDCC?<@CCC:>22))88((38<<&+:4
+@:114:A016U:1:1:15305:1814 1:N:0:1
+AAAGATATCGACATACTCAAGCCCCATACGCTTCAGGCTTTGGTCGAGGCTGGCGAGCAGGTATTTACGTGAACCGCCAGAGCCGTAAGGGCCGGGCCACATATCGTAGCCAGCCTTGGTAGAGATAATCAGTTCATCGCGATAAGCGGCA
++
+@@CFDDFFHDFHHJJGIIJIEIJIJGIJDHHIIIIJJA@DHIGJGJGAFBHIGEH<>BD8:8@?09<5<@?ACAACD94::AC:@@@CA@@4:3>@AC@9>592+:999>
+@:114:A016U:1:1:15782:1885 1:N:0:1
+GAATTAATGCAGGGTTTTGACGCTGAAGTGCTCTTACGTAATGACGAAGGCACCGAGGCTGAAGCCAACAGCGTTATTGCGCTGCTGATGTTGGATTCTGCCAAAGGACGGCAGATTGAAGTTGAAGCGACCGGTCCACAGGAAGAGGAAG
++
+@@@FFFDBFHHHHJHIIIJIJJJJGIJIIJIGJJIGIJHGHJJJGJJIIIIJJJJIEHFDDDEEDDDDDB@DD8?BDCCABBDBBCCCDDEDDDCCCCDCDDDDDDDBBDB<9>CDD?B<<<(22812
+@:114:A016U:1:1:19171:1926 1:N:0:1
+TGCGCTGGTGCAAAAACTGACAGCAAAAGATGTGCAAGTCGAAATGACTCTGCGCTTCGCCACGCCGCGCACGTCACTACTGGAAACCAAAATCACCAGCAATAAACCGCTGGATCTGGTGTGGGACGGCGAACTGCTGGAAAAACTGGAA
++
+CCCFFFFFDFHHHJJIIJJJJJJJJJJJIJJJIIJJIIIJJJJJJIJJJJJIIIJJJIHHHFFFDDDDDBDDDBBBDDDDDDDCCCDDDDDDACDCCDDDDCDDDCCCBDD@B@CCDDDCCBBD@@DDDBB95<(4:@CC1@@AA@DEDEDCCCB7?BDDECDD@BBDDDD=:?CCACA@@4(4>?9?A3&
+@:114:A016U:1:1:11647:1941 1:N:0:1
+AAGACTCATTCTTCAACTCTCGCCTGGAACGTGATCTCGCCACCTTCAAACAACAAGCCGACGTCATTATCTCTAACCGAATGGCAGAAGAGCTTAAGGATGTGGCAGATAAGGTATACACCCGCGATCTCTTTGGCAGCGACTAACATCT
++
+CC@FFFFFHHFBFHGIJGHEHIJJIJJJIJIGHGHIIHIGHGIGIIIIGIJJJIJJIGHHHFF@C?BDDDDDDDCCCBB=@DDDDC@@CCA?C@A@4>A:@>A>@ADBBCD3>A4@DDDDCD9@@B5@@ACD>CACA3?<<943>AC
+@:114:A016U:1:1:6643:1965 1:N:0:1
+TGGGCATACGCCTGCTGTAACACTTGCGCGACTTGCGCCTGGGTCACACCCGATTTCAGGCGGCAGGTAATGGTTTCGAGAATGCCGCGCGGGAAATTGCCCAGATGTGGGGTGAAGATAACGTCAGCACCGAGGTGTGTGGCGATCTCTG
++
+@@CFDDFFHHHHFIJJJGIIGGEHHHIGIG6FGG@FHIIGIJIIJIJGHFFEDDDCCDACCB@DB@BCCAC<@DB7@<>3?CC@C>CCDA>CACA>9(4+(202:>?21()05<@BCDEACACCAC@CD
+@:114:A016U:1:1:21204:2084 1:N:0:1
+CACTTAAATTGATACAACTTTTTAGAGTAGTCATTAGTGAACAGATAATTGATAACTCAGAACCAGTAAAATGCAAAAGACGCACCACGCTGGGAAAATCGCTGCTATTCAGCAGAGTTTTGTATCGAGAATGGTGGCCAGCCACTTATTG
++
+@@CFFFFDHFDHHJJJGIJJJJIIIJFCHHHGJJJJJGIJJGIHIIIIIIGGGJJJJIIGGIIHIEFGHIJEGEHJJGIJIEHFFDEAB?BBDD@CDCD@B??<5>ADCCCCA:<CDDCDDDCDBDB>?CACA>?@ACDA9595;;6;;@C?82=?>CCC>:@9&&&0)0508(85338+9
+@:114:A016U:1:1:19838:2164 1:N:0:1
+TACCAATTCACCTGCGGTATTAACTGCTACGGTTTTCCACTGCGGATTGTTGCTTTCGCCCAGTCCATCCACCAGGTCAGATGCGCGGATCATGCGCCCGGGAACGTAGCTACCATCGTCGCGAGGCTCCAGCATTACCAGCATCGGCATG
++
+@@CFFFFFHHGHGJJIJ@D9B@@BB>BDD?AC:?:>CBABBBB>5@B1?83>CCAA4:@C@&.595?
+@:114:A016U:1:1:17629:2217 1:N:0:1
+GTCCTGCATGATATGGACAAACGGACGTGCACGACTGTGTTTCCAGGCACCAGCCAATTCGTTAATCAGCCCTTCCGCCGGATTGAAGCCGTTGCCTTTATAGCCCTGCTCCGCCCAACGTTCGAGGATCAAGCCCGCTTCGGTATGACAG
++
+CCCFFFFFHHHHHJJJJJJJJJJIJJJJJJJJJJJJJIIHJIJIJJJJJJJJHJJHGHHHFFFFDEEEEDDDDDDDDDDDD>BDDDDCDDDDDDDDDDDDDDDDDDDDDCCC>
+@:114:A016U:1:1:23298:2302 1:N:0:1
+AAAATACGAATGGGAATAAAACGCATTACCGCTATTTCTTTAATCAGCGAAGATCCCCTCGCTGATTAAAATATTTAGAGAAAGACAGGGAAGGCGGAATGATGGCACTGATATGACAAAGTGAATTCCGTATTCAAACATAAGGAACCCG
++
+@CCFFFFD?=ACFGCGAHIIIHHIHIGGHGHGIIIGIIIIIGGH5<@CCCCCCCCC4:@>>ACCCCC:>C@@3:C>CBA9BC:::C(:>?(5&
+@:114:A016U:1:1:20776:2348 1:N:0:1
+AATGCCGCCACCAGATCGTGCGGCAGGTGGCGAATAAACCCTACGGGATCTGCAAGTACGGTTTCACCGACATCTGCAACGTCAATACGCCGCAACGTCGGGTCGAGGGGGGCAAACGACAGGTGTGCCGCGTAGACCCGGGCTTCGGCGG
++
+@@@FFDADDFABFGHIGGEEAFAA?@G/*?BFEA';C=3=CE?B><'33=CC5@@C:@(58555>>::<>>>->CC::@AB@?&5<&&)&0(+5@>&&+((+&&))01+&)&)&
+@:114:A016U:1:1:18138:2369 1:N:0:1
+TAATTAACGCCGCGATTTCTGGCGGGATTTGTTGCGGATTAATGATTTTAATTTGTTTGTGGCGGGCAAGTTCGCCTTTTTTAATCACCACCTGGCTTTTGGCAACCCGGAATTGCTTAACGAGAAACTTCACCAGATGACTGTTGGCCTG
++
+@CCFFFFFHHHGHJIJJJGJJJIIIH',;@>;B9(9@'2=:CA@3@BDD3>>@D@:?C>A1(5&000+:((++((42(5>9CC(::::@(+(+4344(+4((+((+
+@:114:A016U:1:1:15808:2404 1:N:0:1
+CGGTACCGGACCCTGGTCGGTTTCAGCCAGGAAACGCATACCTACCTGCAGTTCATCAACGCCCATAAATACGTCTTTAGGAACACGTTGCACCAGGTTTTCGTCGTACTGACCGTAAGCGTCGTTCGCGCCAACAGCGACATCAAATTTG
++
+CCCFFFFFHHHHHIIJGIIJFHIGGIIJJIJEIGIGJJJJJJJJJJJIHHHHGHFFFFFFDDDDDDDDDEEDBABDDDDDCCDDDDDABBBDDDDBD>ACDDDDDDDDDDDDDDD@DDDDDBDB)5?>BB>A@@>@:CCCCC;@@>5C@CB3<4>44:(:@:48A(:::34::(495.0<3:.59
+@:114:A016U:1:1:13999:2432 1:N:0:1
+CGATCTTCCTTTCGCGAAAAGTGCCGCACTGGGAGTGAAGAATCCCTACCTGACTTACGCGCGCATGGCGCAAATTTTAGATACCACGCCGCAGCCCGCGCAGAACATTGCACCCAGTGCGGTGATCGACGCGACGGCGAAGCTGGGGAAC
++
+@@@DDDDDFFD?FBFFFIDI31?@FF<0DF@F(BFFDFDFF9=EEEFIFFEEEDDDDD;>938=>BBBBBB;5;BB@ABBABEBABBB<<>B@9B?850<<3<&9?3:@BBB:>?8?:A>9<)0(4328>055&0&09&&)9.>@(2&+)(
+@:114:A016U:1:1:22470:2440 1:N:0:1
+ATCACAAAGCCAAAAAAACCGGATACGAAGGTCATTGCCAGCGTTTCCCATACGCCACGAACCAGCAGCCACATCATCGGCTCAGACATAACCCAGTACCTCTACTTTTACATGGTGTTCCTGCAGCCAGGCAATGGCGGCTTGCGTATCT
++
+CCCFFFFFHHHHHIIJJJJIJJJFGIJJJJJIGJJIIIIIIJJIIHHHHHHFFFFDDDDDDDDDBDDDDDDDDDCDDDDDD<>>CDD@CCCCC?A<<@@<B.588(::(+4<<<9<<<509.&&&+(:
+@:114:A016U:1:1:10900:2542 1:N:0:1
+ACCGGCGGCGAAGCGGTAGTAGAAGCGGCGCGTAAACACACCAATAAACGTCTGATTGCCGCAGGCGCTGGCAACCCGCCGGTAGTGGTGGATGAAACCGCCGACCTCGCCCGTGCCGCTCAGTCCATCGTCCAAAGCGCCTCTTTCGATA
++
+@CCFFFFFHDFHFHGI@EBCCCEH@@DDBDDB@@DDCCBD7B?ACDDD>?CBB7)5>B>7;B?AA<8@@@@0905<>CA@BD(:@:>>:2559>9.50>?@D<&)05<59<<594+>A:(28<(8((3<9<&004+>(2(:?
+@:114:A016U:1:1:18393:2548 1:N:0:1
+AACGTAAACAACGTGCCTGGCAGGCACTGCTGCGTTATCTGCCTGCGGAAGTGAAAATTCATCATAATGACAGTGGTTACTTTCTCTGGTTGGAGCTCCCCGAGCCGTTAGATGCCGGCGAATTAAGCCTGGCGGCACTGACGCATCATAT
++
+@@@DDFFFBDH=DFEGIIII@HGGB3EGG9?CA@@DAHCFDDD?FGIIHGIEEHHFHFB?@DFCD@ACC@CDD@CA>@CCCCD@ACCC@:AB>CCCBCCAA>><<<<5>@B@C>C3:1:@:
+@:114:A016U:1:1:13470:2581 1:N:0:1
+CATAGCTGGTTGCACAGTTAACGTCGTTCAGGCCACGTTCCAGACCGTAGTTTGCCAGCGTCAGATCATAAACGGTGGTCACCAGGGCGGTGCTGCCATCAGCCAGTTGCAGGCGTTTCACCGGCAGTTTGTGCAGCAGCACGTTTTCCAG
++
+@@CDBDDFGHFFHEGIIIIIIIGFGDEGHIIIIIIIIIIIIIGCHIIIIIGHGHIIIIGH=BEDECCCCDCDDDD>@@CCD@CCCCDDDDDC:ACD>ABBB39>>CCBBDB59(:>89>CCCD?<1?52832A(+4
+@:114:A016U:1:1:16005:2593 1:N:0:1
+TTCAGCCATGTCACCTGTTTTGCTGTCAGACGCATTTAGAAATTCTCCCGCAAATCCAGGCCACGATACAGCGATGCCACCTGGGCGGCGTAACCATTAAACAGTAGCGTTGGCTGGCGCTGTACATCGAGGATGCCGCCTGAACCAATAA
++
+CCCFFFFFHAHFFHJIIGHHJJJEGHGIIIJIGGIIIIIEIIIJIGIIIJJGGIIIJGGGIFHHHFFDAECEDBB?BDCDDDDDDD>9>>5>CCDDDDBDDCDBBBBDDDDCABC?C<@3>@ACB@9>BC
+@:114:A016U:1:1:13755:2611 1:N:0:1
+TGACCGAAGTACAAGTAGCTCGCTTTGCCGTCAATCAGTACCGTTTTCCGGGTGTCGAAGTTAAAGGCTATAAACGTCGTTACTATCCTTACGGTTCGGCGTTGACCCACGTCATCGGCTATGTGTCGAAAATCAACGATAAAGACGTCGA
++
+CCCFFFFFHFHHHIJHIGJHIIIJJJJJJJJJJJJJJJJJJJHHIJJJJIJIHFFFDDDDDDEEDDDDDDDEDEDDBDBBDDDDEEDDEDDCDBCB@BBDD>BDDCDDDDBDDDDCCBDD@DDAC@ABAB?BCDCD?BAC925<9
+@:114:A016U:1:1:26525:2624 1:N:0:1
+AAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCGGGAG
++
+@@CFFFFFHHFFBIJIJBGHIJIJJIGJJIFIJJIIIIDHIJJJIDIIFHIIJJJIEGID@BHHEFFFFFE9(8>?CDCCC9>B@>B?BD<5<0:>A:4:9&)&&&
+@:114:A016U:1:1:22560:2658 1:N:0:1
+CTGCTGAATAGTGCGTACGTGGTAGCAACGCCCATGCTATAGCCGAATCCGCGTACAAAGCTGTCAACAAACGGTTCTCCGGCTGGTGGTGGCAGATAAACCGGCGGATACGCAGTATTGGCCCAGTTCCCGTAAACCACGGGTGGGGTGG
++
+@CCFFFDFFHHHHJJHIIJIIIGIIJJJJJJJJIJJJJJJIJJJJIIJJJJJHFFFDEEEEDDCDDDDDBBDBDDDDDDDDDDDDDCACC5?@0?:CCCCBB9>59>CD:3>@7>C3::8
+@:114:A016U:1:1:24598:2755 1:N:0:1
+CAGCATATCCATCATCGCGTCAGTACGCTGATCGTCTTCATCATCGATAGTCAACTGCACCAGACGGCGAGTGTTCGGATCAAGCGTGGTTTCGCGCAATTGCATCGGGTTCATTTCCCCCAGACCTTTAAAACGCTGGACGTTCGGCTTG
++
+CCCFFFFFHHHHHJJJJJJJJJJIJJJJJJJJJJHJJJJJJJJJJJJJJJJJJJJJJJJJJEEEHFFDDDD@BCDDDB@DDDDDDDDDDB@BCDDBDBDDDCDCDDBDD?BBDCDDECDDDDDBDDDDDDCCCDBB?@DD(098<@>9@<&
+@:114:A016U:1:1:21792:2814 1:N:0:1
+CAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGAGGCACTGGCG
++
+CCCFFFFFHHHHHJJJJJIJJJIJIJJIJJJIJIJJJJJJHGIJJJJIIHGHHHFCEEEECDDDDDDDDDDDDDDDDEDCDDCCCCCCCCBCDDCCCCCACADDEDDCDCDACDDD>@>BB@BCCDCC>CCCCC@?ACBB@DDDB?9
+@:114:A016U:1:1:22055:2831 1:N:0:1
+GATGTTGCCACATCTGATCATACTCAGCGCGACGTTGTTCATCACTTAACACTTCCCAGGCTTCAGCGACCTCTTTGAAGCGGGCTTCGGCATCCGGTTCTTTGCTGACATCAGGATGGTATTTGCGGGCAAGTCGACGATAGGCGGTCTT
++
+C@CFFFFFHGHHHJJJJJIJIJIIJJIIGIIIJGHEHIJIJIJJJJJJJJJJIIIHEHHFFFFFCEEDDDDDDDDDDACDDDDDDDDDDDDD@CDDDDDDDCC@CCDCDC>A@CCCAABCDCDD@CDDD@B@BA@?@ABB?AABBBA;>BB?AAD:@>B<9?B>:4>A??BBBBAB:@B::>>@3AABD@>@:ABB3@>3
+@:114:A016U:1:1:7776:3009 1:N:0:1
+TTTTGAACGGCTGGGCGTAGGAGAACACCAACGGCCCCAATGGGGACATCCATTGTAATGCGATACCCGCAGACATACGGATATTGCTTGGATCACTATAGTCCGGATATCCAGAATATTGGCTGGAATCCCAGTTTGTATACCAAACGGT
++
+@C@FDFFDDHHAFEGIGGEGGGGIGEEDBHHI?HDHFBGIBCCHEDDD@;A@ECECD@CDCBBDD8>CDDD99?BCAC?@<90@CDCCDCB@>@8A@ACCAC>C@BBBBC:ADD>ACDCCCCD@85+:4((+:++4@CACC(4>>:2(02+
+@:114:A016U:1:1:16110:3027 1:N:0:1
+AGTCGTGATGTGAACTAACGCAGGAACCGATCAGATTGAGAGTTACCGTTCCAGAGAGGGGGACCGAATCCTTATATAAACACTGAGGTAACTCTCATGCTTCATATCCAGTATTCATGATGCGGGCTTTTTGCTGATATTTAACCCGAAC
++
+BCCFFFFFGHFHHJJJIJIGIGIJIJJJIHFHIIJJIIIIIJIIGIJJJJJJJJHFHHFFDDDDDDBDBDDDDDDDEEEECDDDDDDDCDDDEDDDDEDDDEDDEEEDDACCCDDDEEDDDDCDBBB.>>>ABB50
+@:114:A016U:1:1:20207:3073 1:N:0:1
+AAAAGTAGAACTCAGGGTCGCAGAGCGAAAATGGTTGTCTTCACTCCAGAGAATATCGCGCGTCCAGCGGCTGCCATCTTTTATTGCCAGCGCGTTACGCAGAGGGTCCTGCCCGGAATTTGTCACTTCCAGCAGATTATTATTAAGTTTG
++
+CCCFFFFFHHHGHJJJIHGIJJGIJIJIJJJIJIHHIJJJIJJJIGIJIJJHIIJIHHHFFDDDDDDDDDDDDDDCDDDDACCDD>@CDAD<@B9<@BDDBD@?BB?CDDCD?A>@BBDCCCDDDD@AAC:@ACDC@C?BBDCDCACCCDDC@CDDDCACC44(5?B@CC3:>4>@AC>:@+
+@:114:A016U:1:1:8205:3212 1:N:0:1
+AATGGCTCAAAGCCCATATCCAGACCGGTCACCGTCGCGCCTTCGCGCGCCATACTCTCGGCCAGAATGCCGCCGCCACAACCGACATCGAGCACCTTTTTGCCAAATAAACCGCCAGCACGCTCGGCAATATAGCACAGAAGCAGCGGGT
++
+CCCFFFFFHHHGHJJJIIGIIJIIGHGIDFHGJJHIJJJFGIJFJHHFCBDDBDDDDDDDDBDDDBCCCACBBDB@5>BBA??B@BBD(8??BB@?9@CCA@:AAA?A>A>(99>>BDDC<>909CDD>3@(4@C:@>5BACDDD>CBD?AABDDAB@CCCCCC??@8?BDAC9?C:A3:A3>>CC5
+@:114:A016U:1:1:26079:3256 1:N:0:1
+CTGAGCTTCTGAACATTCACGGGTACAGCCGGAGACACCGAACTTCATTTTGTGCGGCGTACGGATGCCTTTGTAGCGGTTTTCCAGTTCCACGCCGAGGCCGACGCTGTCGCCAACGCCGTAGCGGCACCAGGTGCTCCCCACGCAGGTT
++
+@C@DFFFFHHGHHIJJJIIJJIIHGIIIIJJJJJIJJJIIIIJGHIGHIIIAEEEEECDDDDDDDBBDCCCDADEAC@BBDDDACA:AC>CDDBBBDD@9B(9@5)50?39<&5<@D@?ACBACACC?A9C?CBBDDDDCCCD((4:?C9>8ABDFF>BDEEECCECDDDDD:@AA>@:AC
+@:114:A016U:1:1:21329:3381 1:N:0:1
+TATCAACGCTTATTTAAATCAGACTGAAGACTTTATCTCTCTGTCATAAAACTGTCATATTCCTTACATATAACTGTCACCTGTTTGTCCTATTTTGCTTCTCGTAGCCAACAAACAATGCTTTATGAATCCTCCCAGGAGACATTATGAA
++
+@CCFFFFFGHHFHIJJJIJIIJJJJIJIJIJJJIGGHIGIIIJJIJJJHJJJJJIJJJGGIIIIJJJJIJJIJGIGJIEHHHEEHFFFBDEEECEECCCDDDDBDDDBDDCD@BDDD@CACCEDDDCC>>CCDDCBDDBB@ABCDDDCDCD
+@:114:A016U:1:1:22146:3388 1:N:0:1
+GACGGCCTTTCCCCACGCATCGTACTGTGTCCCGTTTCTGCTGCCGTCCGGGGCAATGACGGCGGTCAGGTCACCGGCGATGTTGTATTCATACCGCGTTTCATGGCCCTGCGTGTCTTTCACGGCCATTAACTGTCCACGGCTGTCGTAT
++
+CCCFFFFFHHHHHJIJIJIIJJJJJIJIGHJJJIFGJJIEIJJJJIIIJJGGFDDDDDDDDBDDDBDDDDDCDDDBDDB@@BDDCCCCECCD@CD@BB@3(+4>@?9>598>0++(
+@:114:A016U:1:1:8098:3460 1:N:0:1
+GACTCTGCGCAGCGGCTTAGCGAGGTATGTCAGTGGCTGGCGGAAGCTGCACCGACGTAAAAAGACGGTAAGTATCGCTTTCAGTCTTATGAATATCGCAATCGGCGAATACCTCTGGTCGTAGAGTTTCAGGATAAAGAGGGAGGTCTCC
++
+@@@FFFDDHDFHHIDHGIIIJIFDHDGGIJJJJAC:@>@AAACDDC>B8A>(((44@<.<<8>(:@>BCCC@>CCCDBBDDD?9ACDEDBBBBDCBB<5>ACDCDDDCDCADA@CDCACCDDACDCCC@AC>3:C>:A<98A1>@
+@:114:A016U:1:1:13486:3499 1:N:0:1
+GGGTTATCTGTATTTTTATGCTCGCCGTAGATCTGCGTGAAGTAGTAGTTTTTATCTGCAAACCACCGATTAAGGTAGTAATCTGTTAACCATTCGCGCCAGCGAATAGTCAGTAGTTTAATTAACCATGTTTTGTTTACGGATATTAATA
++
+CCCFFFFFHHFHHJJJJIIJJIJIJJJIIJJJJJGJJIIJJGHHIIJJHIJJJJJJJJJJJJJHHHHFFDDDCDD@CDCFFEEDEEDDEDDDDDEEDDDDDDBDBDDDDDEDDDEDEEDECDDDDDCDDDDDDCBDDDAA?>095@CDDEC
+@:114:A016U:1:1:13475:3523 1:N:0:1
+TGTAATTAGATTTGACTGAAATCGTACAGTAAAAAGCGTACAATAAAGGCTCCACGAAAGTGGGGCCTTTTTTAGCGCGAGAGCCTTTTTTGTCAGCTATCTATACTTTAAGGTGACTGCCAAAACAGACTCGACGTAGCCTTCGAGCTGC
++
+CCCFFFFFHHHHHJJJJJJJJJJJJJJJJIIJJJJIJIIJJJIJJJJGHIHIJJJJJJIJGIJGHFFDEEEEDDDDDDDDBBDBDDDDDDDBDDEDDDCDDDEEEDEDDDCDCCCDDDDDCDDDDDDDDDDCBDBD@BDDDCCC@BDBDC9
+@:114:A016U:1:1:23479:3534 1:N:0:1
+ATTACCGGGTGGAGATAAAAGGGAAATCAAAGGTAATTATATCAGGAAGATATAACGCATTACATTTATTGTGTGTAGAAGAAAATCCATTCCTGTTCTAATTATCCAATTTAAACATCTTAGCATAAAACAAATGATGAATAAGGAATTT
++
+@CCFFFFFHDFHFGGIJJJJJJJIJJJJJJIIJFHIJJJJJIJIJJHIIJJJJJJJJJJHHHHFFFFFFFEEEEEEEDEDDDDDDDDDDDFEEDDDDDEEDDEEEDCCCDDDDECCDDDDDEDBDEDCDDDBDDDDCDCDCCDDC:>>A@>
+@:114:A016U:1:1:24447:3535 1:N:0:1
+GTGTGGACTGGCTTCTTGATGTTTTTGCTACCTGGCTCTACGGCTTAAAAGTAATCGCGATAACGTTAGCGGTCATCATGTTCATCAGCGGGCTGGACGATTTTTTTATTGATGTCGTCTACTGGGTACGCCGCATTAAAAGCAAGTTGAG
++
+@@@FDFFDFHFHFHJJJIIJHHIJJIGFIBEHGGIFGIFEIJJIGIIJIIIFEGIJJJJFHEBEFCCECEBB/==A>@@A@A>C@CCCCBBB>99<>C3:A@BB??((4@<<@><<>BBD.5:::+>239<4>@>4
+@:114:A016U:1:1:18521:3587 1:N:0:1
+AACCTTCATCAATGTTTTTACCCGAGAGGGCTGGGGTGAGGCTTTTCCCTCACCCAAACCCTCTCCGTGAGGAGAGGGGCAGCGCCGGATGCGGCGTAAGCGCCTTATCCGGCCTACAAGAACGTGCAAATTCAACAAATTGTGATGCATC
++
+CCCFFFFFHHHHHJIJJJJJJJJJJEHECHDGIJJJHHJJJJJJJJJJJJJJJHHHHFFFDEDEEDDDDDDDDDDDDDDDDDDBDDBBBBBB@DB<09?AC@BB@@DCDC<<>9CB?C:(+(4@@@CCC
+@:114:A016U:1:1:11141:3587 1:N:0:1
+AGCTTTCGGGCGGTGAACAACAGCGTGTTGGCATTGCCCGCGCGGTGGTGAACAAGCCCGCGGTACTGCTGGCGGACGAACCGACTGGTAACCTGGACGACGCGCTGTCGGAAGGCATTTTACGTCTGTTTGAAGAGTTTAACCGCGTTGG
++
+CCCFFFFFHHGHHHIIJJJJJJJJJGHIIGJJJJJIJGHHFDDDD;?DBDDDDDDDDBDDDDDBDDDDDDDCDDDB95).&5&42<&<@?18>C@A(48289(4:@9(:>><@>>::+493<&09?
+@:114:A016U:1:1:25369:3595 1:N:0:1
+GCTTTGTCAGAAACACTGAGCAGAGCTCGGCGGACTGGACGACGTTGTTGCATGGTAAATCCCCTGGATTTGACTATTACAGAGAGCGTTAGCTGAATTTTTCGCGAAAAACTCAGCTAACGCCCCCAACGGGGCATCCTTATTTTTCCCC
++
+@@CFFFFDHFAFHGGGJIIJIJGIEIIJGIIJJAFAHGHGG:BDA?CBDDDCEDBCDDDD5@CBC@B8@CDDE@CACDDDDAC>A0+9>?BBAACDA>ACD@3&&599&0<5>(4:@:(:?@D5>5&0<88905>)3>>::3>AA:@9((+
+@:114:A016U:1:1:16824:3599 1:N:0:1
+TGAAGTCATACTGGGAAAGCCATCAGCATGAGATCCAGCCGTTAACCCGCAATACCATTACCGAGCTGCCCGCGATGTTCGACGAACTGGCGCACATTCGTGAAAGCGGAGCGGCGATGGACAGAGAAGAAAACGAACTCGGCGTCTCCTG
++
+CCCFFFFFHHHHHJIIJJJJJJJJJJJJJJJJJJJJJJJJJIIIIIJJIIIIIJJJJJJJHHHHFFDEEEDDDDBDDDDEDDDDDDDDDDDDDDDBDDCDD@CBDCCDDBBDDDDDDDDBDDCBDCDCCCACCABBB7B@?BBB9@BCC>>
+@:114:A016U:1:1:23567:3665 1:N:0:1
+CCTTCACCCTGCGTGTGCAGGACTTCACGGCGGTTCAGGCTGTCGGTGATGGTGATGCGGTCCTTCTCATACTGATACGTGTAACTTAAGCCTGCCGGGTTTTGCTGTTCTGTCACCCGCCCGTCGCTGTCGCAACGGGAACAGATTTCCG
++
+CCCFFFFFHHHHHJJJJJJJJGJJJJJJHGJJJ;F;=)@76667;9'5(6;A->-;;52?23;:A>44+:34:(4+:A<(+3(:@3:(>:4?3+8?59<&8?(83(+4>(:44(((+8&5<>50855095:0(&&&050&.FHDHFHGGGEEFFFFFDDCDDDDDBBCCDDDDDDDACDDDDBDCCDDDDACDD@>BBBCCCCC??BCDDCB?B@BBBBD9@ACB>BB<<8??:???9?B@8>C@(:>:::@C.
+@:114:A016U:1:1:15768:3772 1:N:0:1
+GTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGCTTCTCTGACTGCGGTTAG
++
+?@@DBDDD?AFHDIGIEACHBBEGEF99:?D:C??@?@FFFA.=;=FFAAE;6.=DA=CH3?;C?C>BB>AC;>CDCC@C;;;@5;>>C;5:AC:A@>:2(4:>(+>493<4+25?3.9C:8((+0>>ACC@A(+&099>8
+@:114:A016U:1:1:14829:3811 1:N:0:1
+TTTCGCCTTTACACACCAGATCGCAAACCGGGAGTGAACTGCGTTCCTGATGGATCGACAGGCCACTAAGACAGAAGGTGCCAGTGAGTGCGAAAAATGGCACCTTTTCCATGCCTTAGCGTGAATTGACGCTCCCGTCGCTGTGCCTGTC
++
+@7?D?D6DFDDFD4?EEF)).,,389;5(9(5>3>(((5(+(8+(((44+43>++()093<(2(((++(9:4(4(4(+(+4:4((0&&(++:@3(&)&&)()285)09)((((+++
+@:114:A016U:1:1:13671:3857 1:N:0:1
+GTGCCATGGGCGACGTCGAGCATTGCGGTGATGTTTGTACTGTTGTTCCTCTGCGGTTGGTTCCAGGGGATGGGGTGGCCGCCGTGTGGTCGCACTATGGTGCACTGGTGGTAGCAGAAAGAACGTGGCGGGATTGTGTCAGTGTGGAACG
++
+;@@D1ADD8;AD<6@CCAAFBC>=CDE3>@B;@@:;>@>889>9<<>8?@+:9&5994@::4::3:>?@@??&(4>@@<9::383<8??B@B&)&+94:<+4:(4442((3&
+@:114:A016U:1:1:6763:3907 1:N:0:1
+TTCACCTAATTTGCCCATGATGATATCTCCGTATTACCCCGCGCCGGGAATGCGCGGCCGCCAATTTTAGTTTAGGTTCCCGGTAAAGATGACATTAAGGAAAACGTGCTGAATCCTCAATAATCATCCGCTGCTGTTTATTATCTTTTTT
++
+CCCFFFFFHHHHHJJJJJJIJJIJIJJIJJIGIJJJIIJJJIIJJIJIHHFFFFDDDDDDDDDDDDDDACACEDDDCC@C@BDBDBCDDDDDDDDCCCDDDDDDCB<@>ADCDDCDDCDDDD>BDDBDDDDA@CCCACCCCCCD
+@:114:A016U:1:1:26769:4002 1:N:0:1
+TTTTACTGGTGCTGACGCTGAATACGATGACCATTCTGTTGTCGATTGCCGCGCTAGCGCTGGCGTGGGTGTACCCGTTTATGAAGCGGTATACCCATCTACCGCAAGGGGTGCTGGGCGCGGCGGTTGGCTGGTCGATTCCAACGGCTTT
++
+?@@FFFFFHFDHHJJJJJJFIJGJIJIGDGIJGHIJJJJIIHGHGGJJJJJJJJHBEFFDDDBBDDDB=050>>C@<5AC?CACD(:@@99<@&0050+48>((9.8>@DDC((&&)&3
+@:114:A016U:1:1:15506:4030 1:N:0:1
+ACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATGTTCGTGACGGTTATGGCATTCGTTGTGCGCTCACCTC
++
+CCCFFFFFHHHHHJJJJJIIJJJJJJJJJJJJJJJJJJJJJJJJJIJIJJJIJJIJJJJJJJJHHHGHFFFFFEEEFEFEEDDCDDDDEDDDDBDDDDDDDDDDDDD@BDDEDC@DDDBBCDBB@DFHH@FFHGGACFHBAHEIEEBDFB@@ACC@CDEDEDCDDCCCCCDCCDDDDDDDDD@@CC4@AC@ADCAC@9CDC@:@CCCAC>CCA3>>A8?CA::<<3:::>>>(4>9>(+0
+@:114:A016U:1:1:12925:4046 1:N:0:1
+GAGGCGATTGCCGCGTCCGCCGTCGTGAATCAGGTGGCACCGAAGTGGTGGTCACCTTTATTCCCGAAAAAACTTTCACAGACGTCCAAGGAGATACCCATGAGTAATCAGGAACCGGCTACTATCCTGCTGATTGACGATCACCCGATGC
++
+@?CCD>@>C>+4:@C@7AC9@ACDDC>CCB>BB@@@DBB
+@:114:A016U:1:1:14887:4061 1:N:0:1
+ACATCAATTTAACATTGAGTTAACCAAGACAAGGTCACAGAGCTGGAAAAAAAATGGTCTGACCGGTAGCTAAAGAGATAGACGAAAACGAAAAGCCCGCTTAATAACTGTTCACAGAAGCAGCGCGCAAAAATCAGCTGCCACACAACAC
++
+CCCFFFFFHHHHHJJJJJJJJIJJJJJJJJJJJJGIJJIJJJIJJJIIJIIJJJJJJHJHHHHHFDDDDDDEDDDDDDDDDDDDDDDDDDDDDDCBDBDBDDBDDECDDDCDDEDDDDDDDDDDDDAC
+@:114:A016U:1:1:23499:4147 1:N:0:1
+GATCTCGGGCAAAAACTTGGCGTGGATACCAGTACAGCTTCCAGTTTACTGGCAGAACAATTGCCGAAGATTATTGATGCGCTCTCACCGCAAGGTGAAGTGTCACCACAAGCCAATAACGATCTGCTTTCCGCAGGCATGGAACTGCTGA
++
+CCCFFFFFHHHHHJJJJJJJJIHHGCHIHIIGHIGIIJJJJJIJJJJHGHGIJJJDIGIHFHEHHFFCCCDCC>B3;A;@@B@BDDCDD?BBD@D:@CDD>CDECDDDDDD8<@?CCCDD??8??@>:@@CDDD@BB>B>AAC>?>>C:(+
+@:114:A016U:1:1:11367:4212 1:N:0:1
+AAGCGTTGGTTTATGATAAGCCGCAGCGATATGAACAAGAGCTGTATCAACAGAAATGACAAAATCACTATACTTGGTCAACGCAACGGTATAAATAAACTCATCAAATGGTAGTGTTTCGATTTCAAGTATTGGTATTGTCAATAAATCT
++
+CCCFFFFFHHHHGJIJJJJJJJJJJJJIJIJJIJIJJJJJJJJJHIIJJJJJJJJIHHHHHHFFFFFEEEEEEACCDDDCDDDDDDDDD?A?DDEEEDDDACDDDDCACCC@CC@CCCCC??DDDDEDCDCDDCACDDCCDDDDDD<@C3@
+@:114:A016U:1:1:24206:4262 1:N:0:1
+AAATGTAGACCAGACGGGGGGCGAAAGGGAAGACTCCAGGTACGACAAATCAGGCGTTAAATCACGTTTTCTGGTCTTATTGCTGAAAAATTATACCAATATTGCTCTCACCCGCCTCACGCAGGTCTGCCCTAAGGCCTTTTATCAGCTC
++
+CCCFFFFFHHHHHJIJJJJJDBDDDDDDDDDBDDDDDDDDCCDDDDDBCDDDDCCDDDDDDEDDDD2C@CCDDBDDBBBBDDCCDCCD>CC@@BCCDACB<@@B>95@095>C@CBBD:>@CB<<@>A@@?>32()&0999(:4(4@CBCA.&0)+:<&
+@:114:A016U:1:1:26275:4289 1:N:0:1
+GTAAATACGTTTCCATAGATACGCTGGGTTCGCTTGATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGCTGTACACCATCG
++
+@@CFFFFFDFHHHIHHIIEDIDFHGIIIFE?FHD=FFCCGEIIIHFDHDGIEGIIIGGGGCEHHAB?ABDD99@C>:<8A??>?B>B55:@BDBBB?B@BB9.9BB&&5((+:((2:>8
+@:114:A016U:1:1:7806:4317 1:N:0:1
+TTTAAGTAATAGCGAACGTCCTGATGCGCCATAGAAAACTGGGTCGCTTCAATCAGACGCTTCATCGTTGCCTGCGGCAGGGTAAATTCGACTTCACTCTGCCAGTCATCGAGGTTTGGGACATCCCCCGCTGGCAGGGGAGACAGCCACA
++
+@@CFDBDBHHAFFHIJJGDDBHIJIJIJFDHBHC@HHHGIIHFGHI:EHIJHHHHHGFFB>?BCDDDDD<;CDDBBDB99-5+3(:(:(+((++(:44+4((4::::@C@3((5&(+(()((2(((+(0050&&((+559>.58<<(+((+
+@:114:A016U:1:1:26946:4318 1:N:0:1
+CACGCCGCCCGCGAGGGCAAACGCCAGCACATGTTGGCGTTCAGCCGCGACACGATCGTTGACCTTCATATCCGGCGCGCTAATGCTGCGTAAACGCGGCAACAGCCACAGTAGCAGCGGCGCACCGATCAACGCAGTGACCGAACCGGTG
++
+@@CFFFFFHHHHDGIIIBHIIIIIBHHIHEHEHDEB=C>ABDDDD?BDDDDDDDDBD@??CBDDDDDDDDEDCBD900<@B59@CACDDB08?B:<>55<<@8??::>C@<.059950&5&55C35&50(4:(:0.5<&&)&)
+@:114:A016U:1:1:17473:4377 1:N:0:1
+CATACCGGAGCCGTTATCACCGAACATCGGTTTTGGCATAAAGGTCGCGGTTTTACCGAAGCGGTGCGCTACGTTGTGCACAACATATTTGTAGATCTGAATTTCGTCAGCTTTTTTGGTCATGGTATTGAAGCGGGTAGCCACTTCGTTC
++
+@@@FFFDDH)@?FFAFBGH>C@GGEHDHBH?@DG@@FHIG;CECGGEHFC9?BB@CCB?B'7;B;>;5>;;;0A<(+2>::>?B8A:4:CC:A>>>CD:CC@@A@>:CDCD>C>3:&&&09834(<>@CB?C8
+@:114:A016U:1:1:12649:4389 1:N:0:1
+ACTCAAGCTGGTGGATCAGGTAATTCAGCGCCTTCAGCGGAGAAATCACGCCGCAGGTAGAGACTTCAATATCGGCGCGGAAGGTACATAAACCGCCTTCAGGATGACTTTCCGGGTAGGTATGTACGCAAATATGACTTTTATCAAGATG
++
+CCCFFFFFHHHDHIGIJJJJHHJJJIJIJJIJIJJIJJJJJIJIJJJJJJJJHFFDDADDDDDDDDDDDEEDEDDDDDDDDD@@94:ACCDC>A@9B9B>@@CCCCD@CCCC@?>B)5??+:>CDD@BD?.5>C>33:>CCCCC:::>A:A
+@:114:A016U:1:1:5752:4415 1:N:0:1
+TTCAAACGATTCGGCAATCACCACACGAATACCAAGCAGACGCGGACCTTTTGCCGCCCAGTCACGACTGGAGCCTGATCCATACTCTTTCCCGGCAATCACCGCCAGCGGCGTTTGCTCCTGCTTATAGCGCCTCGCAGCATCATAAATA
++
+CCCFFFFFGGHHHJIHJJIJJIJIJJJJIFIHIJIGHGIJJJJIIFFCEEEDCDDDDDDBDDDDDBBDD@BCDDDDB?C@CD:@CCC@CACCC>DB@3AAC4989>7<<<@BB99@(::ACACD:ACACCDD<&)9>99.<CDD
+@:114:A016U:1:1:15250:4420 1:N:0:1
+TAATATTCGCTTGCCGTGTGGCGCAAAGCGGTCGGCAGCGTAGCGGATATTATCGATAAATACTGCCCGGTACCGTTCAGCATCTTCGCCAGCGGGCACGACGCCTGCCATCACATGGACTTGTTCACAATTGAGCGCCCATGCATATTCC
++
+@@@DDDDDDDDDFHGEHH?F@AEFHBBHGEGIG8FGE@HH;?9A<<<>CDC>::CBBB9;@CB22?CCC@C::A@A<9?>@-05>9<9059<>88(+44?>>(4:3@A34@4+4+++(((+(+(&0&&)33>:(+++:(
+@:114:A016U:1:1:14952:4449 1:N:0:1
+GCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGG
++
+CCCFFFFFHHGHHJJJJJJJJJJJJJJJJJJJJJJJIJJJJJIJIJJJHHHHHHHFFDDDDDDDDDDDDDDDDDDDDCCDDDDDDDDDDDDDD@BBD<?<5(9C??B@B9(+)0&50(
+@:114:A016U:1:1:28353:4561 1:N:0:1
+GCAAATGCGCATGCATAATGTCGCTGGCCGGGCGACCGTGCGCCAGCAAATCAGCCATTGCTTTAAGATATGGCGGTAGATGGCGGAAATAACGCTGATACCCGGCACACAAATAATTCAGTCCCGGTTTGCCGCGGGCATCGAGCATGAA
++
+@@@FFFFFG>FHGIBHGGGGHHEGGIJJJIGBGIIIIFEEDB>CDBDDDDDCDCCDBDCCCCCCACCAC>CCCDBB)55>:>C9559-9>CC(832<>@A>::053093:(+:>>@7&99<5&44::
+@:114:A016U:1:1:26837:4634 1:N:0:1
+CATGAAGCATGAAAGCAACTACCTGCGCGGCACCATTGCGGAAGATTTAAACGACGGTCTGACCGGCGGCTTTAAGGGCGACAACTTCCTGCTGATTCGCTTCCACGGCATGTATCAGCAGGATGACCGCGACATCCGCGCCGAACGCGCG
++
+=@@DADDD>FDFHCEFFGE3AFEAFE?:@DF>?DAC)/)-@554=73??;7?DA>:88@',::+085>B>>&(4:@>?<&))09(:>:(:(>3@>CCBB0925A>>B<99&2(+(+4>C321((4>C:0&&&))5(()(&&&)&&&)09<&
+@:114:A016U:1:1:14705:4651 1:N:0:1
+CCCAGCGCCAGCAGCGAATCAATCGCAACAGCGGTACGTTTCGGGCAACTCAGAGCATGAACAAAGTGCAGGAGTGGCGAGGCGAAATATCCTTGCGCGGCATAACGTAAATAACTGACGCTCACCGCTGTGGTAACGAGTTGCAGATTGT
++
+CCCFFFFFHHHHHJJJJJJJJJJJJJJJJJIJJJHIJJHIJJHFFFDDDEEDDDDDDDDDDDDDDDCDDDDDDDDDDDDDDDDDBD@DDACDDCDDDDDDB>BDCCDDDDDEEDDCC@A>@BBDAB@BAB<>??CCD:3:@::
+@:114:A016U:1:1:14619:4683 1:N:0:1
+TTTCACCTAAAATACCCTTATCACCTGATATACAGTAAATACCTTACCTGGCGCGAATTTCTTTGAATAACGGCCTACCGCCCCGCTTATTTGCGCATCTTTTGTCTCGATATACATTTTCTGACATCCTCATCCGCTATGCTACAGGCAC
++
+CCCFFFFFHHHHHJJIIJJJJJJJJJJJJJJJJJJIIJJJJJJJJJJJJJJIJJJJIJJJJJHHHHHHHFFFBDDDCDABBBDDDDBDDDEEEDBBDDDDDDCCBDDDDDCDD>CDEEEDEDC@CACDDDDDDDCBBDBBDDCAC@CCDD?
+@:114:A016U:1:1:9879:4689 1:N:0:1
+AACTATAAGGTTTATATGACCAGCGACGTCAACGCATGGGCGATGGCGAACGGCTGTGTTCGTGTCTACAGTGGCCTGATGGACATGATGAATGACAACGAAATTGAAGGCGTTCTGGGCCATGAACTGGGCCACGTCGCGTTGGGTCACT
++
+B@@FFFFFHHFHGJIJJFIIJJJJJIJJJJJJJJJJIJIGIJHJJJHHFFDDDDDDDCCCDDDDDDDEDDDEDDCDDDDDDDCDDDDDCDCCDDCCACBDDDDDDDCCCC?@BDDDDDCDD3<9>9A>BCA>:
+@:114:A016U:1:1:7176:4699 1:N:0:1
+TTTTATTCTTGGCCTGTTTATTTCTGTGGGTGCCGTGCTGGCGCGTGCGATGTTGCGTCGTGGTGTAGAAGCCCCGGAACAACTGGAAGAGCACGGCATCAGCGTTTATGCCACTATCCCAATGTCCGAGTGGCTGGATAAAAGCACCCGT
++
+@@CFFFFEFFFDHJIGIJIIJJJIGIHJJJEHIIIGIIIJIIIJJGHHIGEFCC>B<8=;>6??@ACC>@59ABBBD<@BB?BCCC?CCB?>9(89ACD:A>(49>>?@:?<>>4@B.&(2<(+8<(((+4>((33(+&)
+@:114:A016U:1:1:25566:4740 1:N:0:1
+TACATTCAGGCGCACCGCGGTGGTATCACCAATGACCTGATTGACGTCCAGCGTGCCGCGGCGGAACCAGGCGCTGCCAATACTGGCGGAGGCGTCAATGCCGGAATCATTGCGCGGCTGCTTGCTGATCATATTGATCGAGCCTGTCGGG
++
+CCCFFFFFHHHHHJJJJJJJGHIHEHIJJJJJJJJJJJJJJJJJIIHJJHHHDDFDDDDDDDDDDDDDDDDBBBBDDDDCDDDDEDDDDB9<00)5>@AA?CDDD>?@CCC3(&)9<>3>BCC<1(4::>CC@>C>:3>2+805><:((0&
+@:114:A016U:1:1:29039:4740 1:N:0:1
+AGTTCTGGAAGGAGATCGAAGCCAGCGTTTGGTTTTCGTTCTGGATCTGCACACCTTCTTTCGCTTCCACAGCCTGGTGCAGACCATCGGACCAGCGACGGCCCTGCATGGTACGACCGGTGTGTTCGTCAACGATGATAACTTCACCATC
++
+@?@DDDDDHFFB+A>&25.&0)&)82?:(+8<9:&+++:BC(4(4>>@>:>(
+@:114:A016U:1:1:14396:4753 1:N:0:1
+TTTACAGCCCGATGTGATCTTTCCCCATGATCCACTATCACAGAATGCTCTATCATTGGCGCTTAAACTGAAACGTGGATAATTTCGTTGTGAATGTTACTTGCTTGCGAAGTTATCATTTTGAAACCTAAATCAGGTAATCACGCCCATG
++
+CCCFFFFFHHHHHJJJJJJJIJJJJJJJJJGJJJJJJJJJJJJJJJIJJJIJJJJJIDHIEIJJHHHHHHFFFFFDDDDDDDDDEDDDDDDDDEEDDEEDCDCDDDCDB<@DCAACCCDDDDD@DCEDDDDBDBBBDDB@DD<5?CCC>5@BB
+@:114:A016U:1:1:23080:4812 1:N:0:1
+GCCATTTGCGCTGACGCCCGTCTACGACATGCTGCCGATGGTCTATGCACCAAACAGCGCTGGAATGCTGCGTGATGCTGCCATTGAGGTGAAGTTTGATCTTAACGTCAGTAAAAGCGCTTGGTTAACGGCGATCCCGCTGGCGCAGCAG
++
+CCCFFFFFHHHHHIJGIJJJIJJIJIJIJJIJJIJJJJJJJJJJJJJHHHHHFFFFEEDDDDDDDDDDDDDD>BDDDDDDDDCDDDDDD>ADDDC@CDDDDEDEEDBCDDDCACDDDBDDDDDDDDDDC@BDBBBD@?BB@@DDDDBDB?A
+@:114:A016U:1:1:18875:4821 1:N:0:1
+AATCGCTTCGTACACCGCTTTTAGCGTACGGGAAGCGTGCGCGCTGGCTTTTTCCGCGCCATCTTTCATCACCTGTTGCAGGAAGGCTTCATCGTTGCGGAAACGGTGATAGCGTTCCTGCAATTCAGTCAGCATACCGGAAACGGCATCA
++
+@@@FFDD?DAFHHIGIGEEFHAHIGIFHIIIIIEGGE@@F@GABA6>BCBCCC>C?@B@C?CCCCCCC@<<<<@ACCCCED?AB?@<359>CBBBB8>@@CBCC??ACAC3@:@A()))05988<>5+4
+@:114:A016U:1:1:26541:4842 1:N:0:1
+CTCTTTCATTAATAACGCTGGCGGGGTGTACGGTTACTCGCCAAGCCCATGTGAGCGAAGTTGACGCGGCAACAGGAATTGTACGGTTAGTTTATGACCAGGCTTTTTTTCAGCATGCCCATACCGATCGTTTTTGTAGTCGCGGGGACGC
++
+@CCFFFFFHHHHHJJJIJIJJIIIJI-=9;?BB6?>;;>;@B???9@975@;2>33+555;5;B<82A(+:>>(4(2<(8514>>4:@>@CA??AACDDDB&+++::@@CC?(39398(+(&(+44(&&))0&&)55
+@:114:A016U:1:1:12356:4854 1:N:0:1
+TTGCACGATTCACCCAACGTTTTGGTCGTTTTGTTAGCACGAAAAAGGGGGCGTGATGAAAGTTCTTAACTCTCTGCGTACCGCAAAAGAACGCCATCCAGACTGTCAGATTGTGAAGCGAAAAGGACGGCTATATGTGATTTGTAAATCT
++
+CC@FFFFFHHHHHJJJIJJJJJJIIIJJJJJJJGIIIJJJJJJIJIJJJJHDD?BDCCDDDDCDEEDDEDDDDDDDDDBBDDDBDDDDDDDDBDDDDCDDDDDDCCCDDCCCDDDD@CCD?BBDDB@?BDB@8::>B>>C@>A>@C>BCC3
+@:114:A016U:1:1:14590:4861 1:N:0:1
+CCAAAATTATATCACATTTCAGCAGGGTTATTATTTCTATAAACCTGTTACATACATCGACCTGGTCAAAATTATCCTTTCTAAACCGAAGGTGAAGGTTGTGGTTGAGTGAAAATTGATCAGTAAGGCCATAGTGCGGTGTAATTATAGA
++
+CCCFFFFFHHHHHIJJJJJJJJJJJJJHIJJJGIJJJJJJJJJJJJJJJJJJJJIJJJIIJJJIJJJIIIJJJJHHHHHHHFFFFFFDCBDD?CCDDDACDBDDDDDDDCCDDDDDDDDEDEDEEDCDDDDCDDCDDDD>B@DDEEEEECA
+@:114:A016U:1:1:22487:4868 1:N:0:1
+GGTTCGCCGTCTAACATGATTTTAAGACCAGCACGAAAATCGTTGCTATAGTACGTTGCCATAAGGCCCTCTGAAATTTGTTAATTGGTAGCTAAGCCACAAAATGGCGCATATTGTAACCCTAAATACCCCATCCAGAGAAGATTGGTTA
++
+CCCFFFFFHHHHHIIJJJJIJJJJJGIIIJGIIJJJJJJIIIIIJJJJJJJIIJI?EEHFFFFFFEEDDDDDCDDDDDDDEDDEEEDDCCDDDEDDCDCDBBBBCACBDDDDBDEEDDDCDBDDDDDEDCCDBDDDDCCBC>CDDDD@>?A
+@:114:A016U:1:1:13246:4918 1:N:0:1
+TATCCTTGCCTGATTAACGTTTTCCGCTGTGCAACCAGAACGGTTCCTGATTGTTTGCCTCATTACGGAAATGTTCAATTTCGATGCGTTGCCCGGTCGCAATAGCGTGATGCAGCCCCTGCCAGTTCTCCAGCCAGGCTTGCGCGACGAC
++
+CCCFFFFFHHHHHJJJJJGHIIIIIFJJJHGGIIIJIJIJJJIIGIIJIIIGIHIJJJJJJJJIGHHHDDDEEEEDEEEEDDDDDEDDDDDDDCDDDDB@DD@CDDDD<@DDDCCCBBBBBDBDDDDCD>ACCDDD8AA8@ACD:?B?>ADBBD7?@B080@:(
+@:114:A016U:1:1:16830:4926 1:N:0:1
+ACCAGACGCTCCAGTAATGACATCTGGCGATATTTTGGATTTTGTTCAAGAAAAGGCCAGAGTGTGGTCATTACTTCACGAACGGCTTGCGCGAACTCGGTTTGATTCGGGTCGCGCTTTTGGACATGGTTGAGGAATGACTCCAGAGAAT
++
+BCCFFFFFHDHHHJHCGGFHJJJIJJJJJFJJIJJJJJGIJJJIIJJJJJJJJJIJIJJJHHHHHFFFFFEEEEEEEDDDDDDDDDBDDB@B@5@BBCBB?<9A>C?<>C:A3::::4>>@C<99(
+@:114:A016U:1:1:16201:4934 1:N:0:1
+AAACGCTGGACACCGCCGACGCGCTCGAACAAAAAGAGATGCCGGAAGAGCTGCCGCTCGATGCCAGTTGGGACACCATTTACACCGCTGGTACACCATCCGGCACCAGCGGTGACTACATTGACGACGAGCTGCCGGTCTACCAGGGCGA
++
+CCCFFFFFHGHHHIAHIIIIFIAHGHIGGFHHFFFCCCCECCCBBBCBCCCCCCCCC5@?BBBB<@ACCA>8.<>>
+@:114:A016U:1:1:8650:4936 1:N:0:1
+TTCAGCCAACAGGCCAGATTGAGCGTAGCCACCAAAGCGGTCGTGTGCCCAGAAGATAATGTCAGGGCCATCGCCAGTTGCCGCAACCTGTGGGAATTTCTCTTCCAGTTTATCCGGATGCTCAACGGTGACTTTAATTCCGGTATCTTTC
++
+CCCFFFFFHHHHHJJJJJJJJIJIJHHIJJJIJJJJJIIIIJIIJIIJJJHHFHHFFFFFFFEEEEDDDDDDD8=B=@CDCCDDBBBBCDCDDCBA?AC>CCDDDCCCAC@D>CDDBBBDBDDDDDBD8<@CCDCCDCE>:9>@>BDE>A@
+@:114:A016U:1:1:8186:4950 1:N:0:1
+ATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGTTACCTAATTTTTCCAGCACAGAACGCGCAACAATCACGTTCAGTTCAATGTCTTCCACCAGCACCACATTCAGCGCCGGTAAAGGCATATCGTCTTCATCAAACGCAT
++
+8=@+===3:?@DF6E?8?91?@?6?9*08('<558=@E@4(5,;3@>6;@CEE;A;(;;3?'0-5&00(8:9@:>4:(:44:4:(+>(2?@>@BA.&00)&)(4448?(:(+23?(8@@ADA3:>.502
+@:114:A016U:1:1:12744:4961 1:N:0:1
+AAGAAGAGGTTGCCGCCATTACGCGTGTTTCGTAAGCCACGCTGGTTGATTTTGTTGGCAGTGGCGACCGCCGGGCTGTTTGGGAACTTCATCCTGTTCAGCTCATCCTTGCAATACCTGAGTCCGACCGCTTCGCAGGTGATTGGGCAAC
++
+CCCFFFFFHHHFHJJJJJJJJJIIIGIGGIJJJJJIJJJJJJJJIHHHHHFFFFFFBDDEDCDCDDDDDBDDDDDDDD?CCCDDBDDDCDDEEDDDDDDEACDDDDDDDCCCCDDDDDDCDCDDCBBDBBBBD@BDBB@+4>AC@A<ABDDBBB@?BC?@CACCA?<<@55<9<5<@?@C>?282?>+85?>3
+@:114:A016U:1:1:16122:5114 1:N:0:1
+TTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCG
++
+CCCFFFFFHHHHHJJIJIIJJIIIJJJIJJJJJJJJJJJJJHFDDDDDDDDDDDDBDBDBBD3<>BDDDD&08?C@AC>AA>CC@<
+@:114:A016U:1:1:8649:5115 1:N:0:1
+CTGCGCCAGACCGTTCGGGGTATCAAACTTCTCTTTAAACAGATAAGAAGTCCCCTGATCGGCATCTGAAGTATCGCGGCAAGGCCACTGAATGAAGCCCAGTTCGCCCATTTTCTCGTAAGTCGCACCGTAGAAATCCGGGCACAGATGA
++
+@@CFFDFFHHHHHFGGGIIG0@FHEHIDDEGGIIIGGGHGIIJJIGIJJGA>@GIGHEECB/3=ACDDCCDCCCCA88?A@CCBD055?:A>ACCDBB22:>B>BBBD?A?<3434>.59&0093944>
+@:114:A016U:1:1:24795:5123 1:N:0:1
+TAACCGATCAGGATCGGCGTTATAACGTGCGGGCGTGTCGGCTTCAATATCACGAACGCGGTGGACATAGTCATAATCGTGGGCGTGACGACCATCGTGACGCACATCGTCATCTTTGATGCGATCATCATCCAGCTTGCTAAGGATCTGC
++
+CCCFFFFFHHHHHJJJJJJJJGIGIIHHIIJJJHFCCBBBD@BDCDDCDCCDDBDDDDDDB>CC?BBDC<@C>34445(+<:A@<>CCC>4
+@:114:A016U:1:1:24365:5219 1:N:0:1
+ACGCGTGAGAGTTGAAATCGCAGTACGGGCATTTCTGCACGCACCACGGGATGTGAATGTAGAGACTCAGCGGAGGTAATTTAACCATTACGTAAAGCGTCCAGCAGCAGTTTCAACGCCTGACCACGGTGGGAAATGGCGCTCTTTTCTT
++
+CCCFFFDDHGHHHJIJJJGEGIJEGGHIJIJJJJHIJJJJJJJJGIJIJHEBE;BDEEEC@CECCA>CCCCBBDDD8@ADDEDE@CDDDDDDACBDACB??<@A><@8(29@B?<:>((5005DDCCDEDA@@CD
+@:114:A016U:1:1:17824:5259 1:N:0:1
+TTGTCGCGTTCCAGAGTGAAGGTACTGTCTTTCAGTTGGGCAACACCGGCGAATTCAGTCCCTGTTGCATGGGTAAAGCCAAACATCTTGTCGGAGGATGACAGCTGCACTCGCATCAAGCCATCGCCTTTGAGCGCATTGGTCAGCGTGT
++
+@CCFFDFFHHHGHIIJCBHGGIEHHIJIJIJJIJJGIIGGIIJIIII>AHHFFDDEEDDDDDDCCCDDDEDDDCCCDCDDDDDDBDDCDCDDDBCCACDDBDDDACCCD>B?@CC::?CC?8?B05@ACBCBBB.0::(::<@>?D9
+@:114:A016U:1:1:18035:5345 1:N:0:1
+GTGTCGATGACGCGACTCACTATCAGTTGAACATCGGCACCAGTGACAAACACACCAGCGTTGATGGCTACTACAGCCATGATGGTTCGCTGGCGCAGGTTGACCTCAGTGCGAACTACCATGAAGGGCAATACACCTCTGCGGGCTGGTC
++
+@@?DDDDDDDHHH@BBBBBAC?9A:4>@9AC@4:><>B9&5<@
+@:114:A016U:1:1:8237:5373 1:N:0:1
+CAGCACCGCGATAGGGAGCTGTAACGCCCTACCGAGTCGCTGGAAAAAACCTAAAATATTCATCTTATTCCCCCTACGAGAACCCTATTTGGCTCGTTTCAAGCCGTATTTTTATTTTGCTGCAAATTGTACTGCCGATGTTCTGTAATCA
++
+CCCFFFFFHHHHHJJJJJJJJJJJJIIJJJJJJGJGHIIIIIHHHHFFCDDDDDDDEDEEEFEEDDEDE@DDDBDDDDDDDDCBBBCCCDEDDDD@B@BCCCDCDDB>ACEEDDDDEDCCDDCCACCCCCDDDACD>@@B?ADDDEDDCCA
+@:114:A016U:1:1:6511:5398 1:N:0:1
+GAGAATTTAGCAACCAAAAAAAGAATTCAGACCATTCCTCTTCGCGTAACTGAAGCATTGCCGAACGGGTATTTTGGCTTAATCTGCCAGAAGGCTTGATCAATTCTAAATGTTCAATAAGCGCAGCCAAAACCGCCTCACTATCCGCCCC
++
+@@@D1DDD??FHDIGDEEHD?CDGIGGIGAFGEBFHHIIIGIGEBBG@5@GGE<;AABE@C>A<;@B8?=8@CCECBCCB>ACCCDCC>:(+><((35CAD@((4::3@::44(4>9B>BBB?3(2<>B<>AC:ACC<@?B@
+@:114:A016U:1:1:12033:5403 1:N:0:1
+GGTGACCAGCGCCGCATTATGACGCCAGAACAGGCGTTGTCGGCTGGTGTTGATTATATGGTGATTGGTCGCCCGGTAACGCAATCGGTAGATCCAGCGCAGACGCTGAAAGCGATCAACGCCTCTTTACAGCGGAGGGCATGATGAGTGA
++
+@CCFFFFFGHAFHIIEGFHIHIEEEGEG>GEBBGFGFGGHGHF@C>;>A,:(5;:9;BB>@-0?>?@5@>?BDDD7A9@@&5A@@>5<829C3&))29:A>>A:::55&0&)9>:4:(:(
+@:114:A016U:1:1:12003:5459 1:N:0:1
+GGGGTGGCCGGGGCGCTACAGGCTCCTACATTGAGCTTATTTCTGAGTCGTGAGGTTGGCGCGCAACCTTTCTGGATCGGCCTCTTTTATACGGTGAATGCTATTGCTGGGATCGGCGTAAGCCTCTGGTTGGCAAAACGTTCTGACAGTA
++
+CCCFFFFFHHGHHJIJJJJJJJJIIJFHGHGHFFFFFFEEEEEEEEDDDDDDDDDABB>BDDDDCADDCDDDBBBD>CCC@CDCDDAA?C@@@CCDDDDDCC@8:?B>BB@BDBCCDC>CC@:<:8??<8DFFCECCB95>?<<;AC?>55?B:4<0<<5&5?(:A+:4+4>A>C?B&+8?:3>3:@(&0)000&5:(4:4(43+3(((:>3&&0
+@:114:A016U:1:1:29020:5604 1:N:0:1
+GTTCGATGCAGATGCTGCGCGGCATGGGCGTGAAAGTGATTGAAATCCCAACCGATCCAGAAACTGGCATCAGCGTTGAAGCGCTGGAACTGGCGCTGGAACAGTGGCCGATTAAAGGCATCATTCTGGTGCCAAACTGTAATAATCCGAG
++
+CCCFFFFFHFHHHJJJJJJJJJIJJJJJIIDHEHIG=CFHHIGIIJIFHHHFFFDDDDDDDDDDCDDDCCDDDCDBDDDDDDDB@BDDDDCCCDDDBDDBCDBDDDCD25550<:4>A@?ACCCC>CDD49A9@C<1>A@B@CD>>A:&(&
+@:114:A016U:1:1:29647:5671 1:N:0:1
+TGTCGCCTGCTGGACGTTAGCCCCCCACTTTCTCCTTTCTCTCCCCTTTCTTTCCCCCTCCCTCCTCTCCCCTCTCTTCTCTTCTCCTTTTCTTTCTCTCCCTTTTTCTTCCCTTTCCCTCTCTCTCTCCTCCCCCTTCCCCTCCTTCTTC
++
+CCCFFFFFHGHGHJJJJJJJJ4<11B0)09B00?B999?=BFHIIIGI=CHCCGIHHH?DDC6;?CC5;ACD5=9A::>>@CCDCDC4:A:>@:>>:A@CCD(:@ACC@:::<@44>CC3<3>(::A@CC(:ACCD@DDBDDDDCCBD?CCDDDB<@<>3<@:::>((++
+@:114:A016U:1:1:12580:5721 1:N:0:1
+TGGGCTTACGGCAACGCCAATGGTAGAAGCTGCGTTGAAGAAGCCCATCGGGCGTCCGCGCTCTTTGATGTTGAACCAGTCGTTAATGACTTTTACACCGCATGGGTTCATCGGCGCTTCACCAATCCCCATACCGATACGCACCAACACG
++
+C@CFFFFFHHHHHJJIJJJJJJJIIIJJIJJJJJJJIJJIJJJJJJJJJHEHFDDDDDDDD@BDDDDDDDEDCDDDDDDCDDDDDDEDDCDDDDDDDDC@B@BBDD?@BCDDCDB<@>BDDDDDDDCDCDDDDDCB3>@@?@DDB8?BDD5
+@:114:A016U:1:1:23349:5747 1:N:0:1
+GATTATAAAATGACGCCAGGTGAGAGGCGCGCGACCTGGGGTTTAGGGACCGTATTCTCGTTGCGCATGCTGGGCATGTTTATGGTTCTGCCGGTTCTGACCACGTACGGCATGGCTCTGCAAGGTGCCAGCGAAGCATTAATCGGTATTG
++
+@@CFFFFFHHHHGJJJJIIJEHIIIIJIJIIIIIGJJHGHF4=BECEDDDDDBDDCDDC<@BDDDBDDBDCDDDBCDDCDDCCDDCC@:>ACBD59.(++(49>ACDD?(4:85)5.&+((+(4:00(+(((
+@:114:A016U:1:1:24391:5777 1:N:0:1
+TGACTCATCAAGCTCTCCATATTTGTCCTACAGATAATTGCGAATAAATACAGAGAACGCCATTTATCCTTTTAGCATACCTGCCGGATATTTTTATTTGTATTAATATTAATTGCATCTCCAACTCCGAGCATATCCATCATTAAAAACA
++
+CC@FFFFFHHHHGIJJJIJJJJJJIHIJJJJJJJJJIIIGJJIJIIIEIIIJIJIJJJJJJJIJIJGCGHGHHHFCFFFFFECECDDBBBCDDDDBDDDEDDEEEDEEEDEEEDCCA@C@C4>?CCDCBDB>3>C:>>CCC>2
+@:114:A016U:1:1:27066:5785 1:N:0:1
+TTTTTAAACGGCGTAAACAGCGCGGAGCGGGTGCCGATCACAATCGCCGCCTCACCATTTTTCGCTTTCAGCCACGCCGAAAGACGCTCGCTGTCGTTCAGGCCGGAATGCAGAACTTCCACGGGGGCATTAAAACGTTCACGAAAACGGG
++
+?@@FADDFFFFHHEHIIIIBFG8GHBGGHFH4=AC@>;;@>AC9><<@BBB5<@8ACCCEDC95;7BB>@@:?AC90<59@55?C>&)5-9255(88<89(43855<>?CC>::><>:>>(:<5>955&9((+4>A><3+(4(<<5<3.8.
+@:114:A016U:1:1:16294:5786 1:N:0:1
+AGAAGATCAGCGATAAAAAGCCAAACACGGCATCGCGTTCAACGCCAAAACCAAACTGGCCGGACAAACATTCACGAAGTGTATATAACGGGCTGGTACCAATATCGCCGTAGACAACTCCAATCGCCGCGAGGGTAATCGCGGGCAATGA
++
+CCCFFFFFHHHGGJJJJIIGIJJJGGJIIIJIIIGGGHHIJJIIJJHHHFFFDDDDCCDDDDBBDDDDDDDDDEDDDDD<@:@DEDDDC@BDDB@BACCCDCDDDEDDDB>B@BDCCBCACCCC9?BDDB>>BBB+9>@C?<0<9@9<>C@
+@:114:A016U:1:1:23064:5834 1:N:0:1
+GGTGGCCCGATAACCCAGGCAAGCGAGATCTGAGCACGTAAAAATGTACTGAACATGACCGTCTCACGGCCTGTTCTGTCGGCGTGTTCACGGGCGAGGGCGAACATTTGCGGGTTTTCGGGGGGAGGAAAAATACTCAGAAGAACGCCAG
++
+@CCFFFFFGHHHHGIJIIJGJIEHDGGGIIIIIGGGIIHGIJIJJICDHIJIJIGHGHEEFFDACDED?BBBDAC@@ACCB@DD>B9A:A@(5@0559@B.<.09<:AD:49@B&(++(5&&)5&)0&(299((:444>(+:3(3(&&5&&
+@:114:A016U:1:1:7283:5858 1:N:0:1
+GTGACCAGTGATATCGCTTACGCCTTTGGCACGCCGCCAAGCGACGCCGAAGCGATTAAAGTTCGCCACGGTTGTGCGCTGGGTTCCATCGTTGGAAAAGATGAGAGCGTGGAAGTGCCGAGCGTAGGTGGTCGTCCGCCACGGAGTCTGC
++
+@@CFFFFFFFHGHIIEI@GIIEGGBCHIIIIIIFDF@FGHG=GIEEDCA?BBBBBBBBCCC@CDCBBCCCC<@BBB@<<>>B@?@?@CCCBA?35<<>@>0&+(++0:<9(8<&29>:?++EB::ACDEFICEA)?:)7?BDD*)=@@AADDDD@<5>AAA3('',858>-((9<<85-5&&)&)&0(2+&))55(((9>>3810+(((0&)(+(:(4:+<:+(&&&&((+&2&&&)
+@:114:A016U:1:1:26063:5926 1:N:0:1
+AACCCATAACCAGGGCGTATTTGATGTTTACTCACCGGATATGCTGCGCTGCCGTAAATCTGGCGTGCTGACCGGTTTACCAGATGGCTATGGCCGTGGGCGCATTATCGGTGACTATCGCCGCGTAGCGCTGTATGGCATCAGTTATCTG
++
+CCCFFFFFHHHHHJJJJIIJJJIIJJJJJJIJJJJJJJIJIJJJJJJJJJIHHHFFFDFDEEEEDDDDDDDDDDDBDDDDDDDBDDDDCDDDCDDD@DDDDDDDDDDEDDB88?CDDDC?B>BB<993:
+@:114:A016U:1:1:18051:6047 1:N:0:1
+TACCCTGCTGAATATGAAATTGGAGATATCGTCTTTACATGTATAGGTGCTGCCTTATTTGGTCAAATATCAGCTGCATCAAATTGCTGGAGTAATCACGTCGGGATCATTATCGGTCATAACGGTGAAGACTTCCTGGTTGCAGAAAGCC
++
+CCCFFFFFHHGHGJJJJJIJJJIGIGJIIIJJJJIJJJJJJIIJJJIIIIJJJJJIJJJIIHHIJJJJJJJIJIJJJJJHGHHHHFFFFFFDCCEFEDDDDDDDBBBDBDDC@CB<@DDCEDDD@CC?CC>8
+@:114:A016U:1:1:17032:6091 1:N:0:1
+CTCCTTACCGCTGGTGTAGCGTTTGAAAAGATCGTTGTAAGCCAGCGGGAATCGCTCTTCCAGATGATCCCACAGCTTCCCGGCGATCGTTGCCTGAACGCTAAGATCCCACCATTGGTTGTTGAAAGCATACCGCGCCAGTTGAGCCTGC
++
+C@CFFFFFHHHHHJEGHJJJIJJJIJJJJJJIIJIJJJJJIJJJJJHGIHHFFFFFDDEEEDDDDDDDDDDDDDDDDDDDDBBBBDDDDCDCDDDCDDDDD?BBDDCCDDDDDDCA>@>ADDDDACCDDC:>A<<>9>D<<+:>A@CCA?<
+@:114:A016U:1:1:4965:6094 1:N:0:1
+CCCCTGGGCGATAGCGTCAGCCATCATGCGGGTGTTATTCGACATGGTGTCGTAGAAAATGGTGATTCTGTCTTCCTGATAATCAGCCGCCCATTTCAGGTACAGCTCGACAATTTGCGTCGGGTTATCGCGCCATACCACGCCGTGGGAA
++
+CCCFFFFFHHHHHJJJJJJJJJJIIGIIJJJJ;@FEHJJJJJJHHHHFFFFDDEDDDDDDDC@CDCCFDEEDDDACCDDDCDDDDDCDDB>@BDCDEAACCCDD@CACBDDDA?.99@8@CCCCBBD<<.5<23
+@:114:A016U:1:1:22538:6094 1:N:0:1
+AATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACCAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCA
++
+@CCFFFFFFHHHHJJJJJJJJJJJJIJJJJJJJIJJIJJJJIJJIJJJJJJHIIJJJJJJJJIJHHHHFFFFDCCDEEDDDDDDDDDDDDCCDDDDDDEDDDDDDCDBDDDB>ADDDDDDDDDEEDDDDDDDDDAD>@@@ACDAAD3:ACC
+@:114:A016U:1:1:16397:6188 1:N:0:1
+TGCGCATGTTGATGGCGATGACGAAGAGTAGTCAGCAGCATAAAAAAGTGCCAGTATGAAGACTCCGTAAACGTTTCCCCCGCGAGTCAAATGTATGTTGATGTTATTTCTCACCGTCGCCATGGTGCACATTGTGGCGCTTATGAGCCCC
++
+CCCFFFFFHHHHHJJIJIIFIJIIJJIJGGIJJIJJIJJJJJJJIIJJHHIJIJIIJJHHHHHHFFFFEACEDDCCDDDDDDDDDDDDCA@AC>CDEDEDCDDEDCDDEEEDECDDB@BDB@DDDDDDCC:@CCDD@BDDBB.@ACCD?85
+@:114:A016U:1:1:9477:6239 1:N:0:1
+GTTGATTGATCCTTGCGCGCGGCAAATTGACGGGGAGTTTAAAGATAACCCGCTGCTGCACGCCGGTCATAATGAACCTGACTATCGCGACAACGCCGCCATTGCGCCGAAATGCGTAGTCGTGGTTGATCACTATGACTGGGAAGCTGAT
++
+@@@D=?D>?8DAHG>9C8?G@D:@@G@DA>FAFAG/=5?@@;@AACC>;>?8=8=?CA@C@25?;308&>:@@:@:@291>AA@A>@B&)059<>&0)000?4:40&)50&&+4352(02(&(&)(28+((4(+:3(++((28238(+:((
+@:114:A016U:1:1:18919:6242 1:N:0:1
+ATGGCGTCAAATCAAAATCGCTTTTTTTGTCTTGCGCAAAATCGGCAGCTATTTCCCCTAAAACTGACGCAAATTTAAAACCGTGCCCACTCAGGCCGGTAATGAGCAGTGTATTATCGTGGCCGGGTAGGGTATCGATAATAAAATCTTC
++
+CCCFFFFFGHHHHJJJJJJIJJJJJJJJHIJJIEGGIJIJJJJHHHFFDEDEEEFEDDDB?CDDDDDDDDDBBDDDDDDDDDBDBDDD@BDCDACDB5@BDD(>>ABBBD?@ACC@AA@AA
+@:114:A016U:1:1:6779:6250 1:N:0:1
+CGGTAACACCCGCAATACGGCGTATTTCGGCCATCAGACTGCTGAAACTCTCAAACTCCAGTTCAGCAAAATTGAGGTAGATTCGCCCAATGGGATCAATCTCAATACCGCGTAAATCAATGCCTCTTAGCACGAGTAGATCGAGTAATTC
++
+BBCDDFEFDHHHHJIGIIJGIIGGIIJJGGGHEHHGJIJIJIGGIHHGHHFFFFFFECECCCDDEDCCCDDD@CCDDCCC:ADCBB>BD@?CCD?CBBDCCDDDAADCCBBD>982::>CC:>ACCCCCCACA9<>?8AC@3<C>>
+@:114:A016U:1:1:14329:6297 1:N:0:1
+GTTTTCAGCCCTTCTTCTTCACCTTCAATGGATTCGTGGATTTCCTTGCTGGCTTTCCAGATAAGGAACAAGCCACCCAGAAGCAGAATCAAATCACGGGCAGAAATTTCCTGACTGAATATTGTAAAAAGCGGATTTGTCAGGCGCGTAA
++
+CCCFFFFFHHHHHJJJJJJJIJJJIIJIJJJJJIJIJJJJJJJJJJJCHIGHHJJJJJIIJGIIIJJJIIIJHHHHFFDD?CEEDDCCDDDDDDDDDB?DDDDDDCDDDDEDDDDDDDACDEEEDEDEDDDBB9BBBDDEDE@ACD@BBB2
+@:114:A016U:1:1:23266:6355 1:N:0:1
+AATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAAC
++
+CCCFFFFFHHGHHJJJIJIJJJJIJJJJIJJGHIIJJJJIIIJIHHHHHFFFFFEEEEEDDDDBDDDDDDDDDDCDDEDDDDDDDDDCDDBBB>BDDDDBBDDDDDDDDDB@BDEECCBDBCACBDDDDDBBDB@>3>ACCCCCCCCEDE:
+@:114:A016U:1:1:28071:6363 1:N:0:1
+AAGAATTGCCGCCGCCGTTGCGACCACAGTATTCACAAATGCCAGTGCCGCGATTTCATTCGCCGTGCCCGCTGACCCGGCGTTAAAGCCAAACCAACCGATATAGAGAATGGCAGTCCCGGTGAAGACCATCGGCAGGTTGGGGGGGTTA
++
+CCCFFFFFHHHHHGIIJIJJJGIIIAGHJGIIJIIIJJJJJHHHGHFFFFDDDDDDDDEFEABDBBDDDDDDBD@DDD@@BBDDBDCCCDCDDDDDDDDDD@BDECC@CCCCAB9093>@>BBD?CDDDCD>9@>3BDDBDDDDCDDDBD@>CDCCCD@B>C>B>>@DDCCDDDDD@BBBBCC@AC:@A?8>508(
+@:114:A016U:1:1:26795:6374 1:N:0:1
+TGGAAGCGGCACAATCCAGTTTACCTATGACGCCGATGACCCTGCGGGGCATGTATTTTTATATGGCTGATGCGGCGACCTTCACTGATTGCGCGACCGGGAAACGTTTAATGGTAGCGAATAACGCAGAGCTGGGGCGTAGCTACCTGGC
++
+@@@DDDDDHF?@FGB>GDCBFGH>:@@CC0>3::@1(+(:>>@B-0559@>3:>>A4:@@>505<<<@B&052982?A((4>CC>3@>&.90:((0&&&(((3(())5<05(:@(<441(
+@:114:A016U:1:1:6650:6381 1:N:0:1
+ATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGTGGCCCCCCCCGTCCACGAGCTTTTTGCCCTTTTACATGCTGTTTCGGGGGACGGCGGGTTTTTTCTTT
++
+CCCFFFFFHHHHHJJJJJHIJJIJJJJJJJJJIIJJJJJGHFHCHFDDBDDDDDCDD?B((&((,(+((+((+((((4((((&&&&&0&&&&)&&234(+
+@:114:A016U:1:1:16007:6397 1:N:0:1
+TAGATTTCAGGGGGGAGATCAAAAAATTTCGCTCTGTGCCAGAGCGGACATTCACGGAGCTGGTTCATTACCAATGAGGTTGGGCTTTTGAGGATAAATCAATGATCAGACGCCAACGGAAATCAAAAGCAACCATGGAACGGAATTGCTA
++
+?=?D?DDB?==?FF;06-?CCCCCBB9(>:@CDACC@(4+3:?(8?59?(<9:>889(9AC>>3:>@C>>@C(4>0)559<(&((0@@:@C:(@DCDCDDDCD>CBDCDDDDD9CCD@@CCB>589>9<>B&3A@(++:AAC
+@:114:A016U:1:1:17301:6405 1:N:0:1
+TATTCTTTGCTCAGCTTGCTGTGCGGAACCGCAACCAGTTTCGCATCTTCACCGGCTTCGTCGGTCATTTTCAGAACGCCAACCGGACGGCAACGGATCACAGAACCCGGCTGCAGCGGGTACGGAGTTGGGACCAGTACGTCAACCGGGT
++
+CCCFFFFFHHHHHIJJJJJJJIIHJJJJJJJJJJJJJJIJJJJJJJJJJJJJJHHFDDEDDDDDDBDDEEEEDDDDDDDDDDDDDBDDBBDDDDDD@5
+@:114:A016U:1:1:10506:6424 1:N:0:1
+AACCAAAGCGCAGATTCTGGAAAGAGAAGGTCGCCTGCCGGATGCCGTTATCGCCTGTGTTGGCGGCGGTTCGAATGCCATCGGCATGTTTTCTGATTTCATCCAAGAAACCCACCTCCGCCTGATTTGTGTGGTGGCCGGGGGGCTCTGG
++
+CCCFFFFFHHHGHIJJJJJJJJJJJJIJJJHIJJJIGJJJJIJJJIIJHHHFFFDDDEECDDDDDDDDD;BDBDDDCDDDDDDDDDD?+4@((+4((4>((4((:(+29?(((2(&(2(5)+)+(+:4:+9<(8(++((09.05&&)((+4
+@:114:A016U:1:1:9558:6433 1:N:0:1
+CATTAAGCCCTTACGACTACACACTGTATTGCAGGAAGTGGGTAAAACCAAATCTGACTTGTACGCCACAGGATGATTCCCAGTATCCAGGGATGGATACAAAGACAGCAAGTGATTACCTGCCACAGACAAATTGGCCAACAACTCAATT
++
+CCCFFFFFHHHHHJJJJJJJJJIJJJJJJJIJJJJJJJHIJJIJJJJJJJJJJJJJJIJJJIJJJHHHFFFDEEEDFEFEDDDDDEEDDDDDDDDDDDEDDDDDDDDDDDDDCCCCDEDDDDDDCCDDDDDCDCDACCCDDDDA@CDD9A@
+@:114:A016U:1:1:16441:6449 1:N:0:1
+CGGCGAGGTACTGGTGACAGACTTAAAATTATGGTGGAACGCCGATACGCTCATGTTGGCTTCTGCCGCCAGTTGCTCGACGCTCAGGTTTTCGGTGTATTTATTCTCAATCCGTTTCAGCACGCGGCTAATCAGACTGAAGTGAGTCTGG
++
+BCCFFFFFHHHHHJAGHIJJJIJJJJJJJJJJJJHIJJJJJJJJJIHHHFFFDEEEEEDDDDDDDDDDDDDDCDDDDDDDDDDDDDDD>AC@CDD??<995>A(>98
+@:114:A016U:1:1:8917:6530 1:N:0:1
+AAATGGCTTACCGGGCAGCCGGGTGAGCGCATTGTGCGTAAAGACCGCTATGGTCGCGTAATTGAAGATATTTCTTCTACTGACAGCCAGGCAGCGCACAACCTGGCGCTGAGTATTGATGAACGCCTGCAGGCGCTGGTTTATCGCGAAC
++
+CCCFFFFFHHHHHJJJIJJJJJJHIIJJJJGIJJIJJJIGGGHHHHFDCDDDDDDDDDBBBDDDCDCCDCDEEEEDDDDDEDDDDDDBBDDDDBDDDBBBDDDCBABB@>B<
+@:114:A016U:1:1:26826:6539 1:N:0:1
+TTGGCTTTAACGTACGTGCTGATGCCTCTGCACGTAAAGTGATTGAAGCGGAAAGCCTGGATCTGCGTTACTACTCCGTCATCTATAACCTGATTGACGAAGTGAAAGCGGCGATGAGCGGATTGCTGTCTCTGGGACTGCACTAACACAT
++
+@?15>;(--59?B@>5(5(:@C>@(4>4::0+0((+:3>:(0<.&&504((05&&&(4>>CCC3(4((2(+((++(+(4(((4
+@:114:A016U:1:1:5792:6551 1:N:0:1
+TTTTCGTAAACAGAAAACTGTAAGCAAATCGGGAAAATGGCTCGCCCGCCCGGCAGTAAGCTGTCTCGTTCTCCCCATTCTTTCTTTTGTTTGCAGCCGGAGCATCTCGCCCTGGCAGACGCTTAAGATTGTAGATGAACATCAGCCACAA
++
+@C@FDDDDAFDFHFHEBGIIC@HAFGH@48<3(44438&5&)0244(+++&))+(29((2).025((44>:(4(4>((4+(4(++((++2(
+@:114:A016U:1:1:9486:6605 1:N:0:1
+CAAAGGTCATTGTGGCGAAAGTCAGCTGCTTGCCGACTTTCAGGCCGCGCTTTTGCGCTTCCTCCAGCCACGCTTGCGGCACTTCGCGTGAACTGCGCAACGCCCGATCGCCCGCCATAAACTCACGGCTTTGCTTGCTTAAGCCCTTCTC
++
+CCCFFFFFHHHFHJIIIJJIJIJIIJJIJIJGJJJIIIJIIIFHIJGHEFEDDECEDBBBDDDDDDDDDDD@B@DBD>@B@BCDC@DBDDBBDC@CB@9@D@BDBD@>D>BDD@B>B<@>CCCC>9(09&&9AC3+:?3>::>3>??<<>>
+@:114:A016U:1:1:17513:6694 1:N:0:1
+CTAACCTGTTAAAGGTGCAAAAGATTGCACCTTTTTCAGATTACCTCGGGCCTGGTTACGTCGCTTTCGGCTTACTGATTTTCCTCAATTAGCAACGCTTCAAGCAAATCGAGATCGTGCAGTAATTTTTGCAGCGTTTCGTTGCTGATCT
++
+CCCFFFFFHHHHHHIHIJJJJJJJJJIJJJJJJJJJJJJJIJJJJJJJIJJJJJJIJJJJGHHHHFFFFDDDDDDDDDDCEDEDDDDDDEDCDDDDDDBDDDDDDDDDDDDDDDDDDCDDDD:ACDDDDDDDDCDDDDDCBBDCDDDDDD@
+@:114:A016U:1:1:22972:6695 1:N:0:1
+ATGCGCTACGCTTATCAGGCCTGAATTTCCATACAATATATTGAATTCGCAAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCACACTTTGCCAACAATCTAAATCTGTTTAATTGCCGGATGTCAGACAT
++
+CC@FFFFFHGHHHJJJJJJJJHIJJJJJJJIJJJJJJJJJJJGJJIJJJJJJJJJJJJJJJIGHHHFDDDDDCDDDDDDDDDDDDDDDDDBDDDDDDDDCBCCDDDDDDDCDCDDDDDDCCCDDCDCDDDDCACDCC@>@DBBCCC@CCB3
+@:114:A016U:1:1:25091:6782 1:N:0:1
+GAAAACGACCAGGTAGAAGCGGTTAAACGCTTTTTTGCTGAAAATGGCAAAGCACTGGCTGTTGGGGTGATTTTGGGCGTTGGCGCACTGATTGGCTGGCGCTACTGGAACAGCCATCAGGTTGATTCTGCACGCTCCGCTTCTCTTGCCT
++
+CCCFFFFFHHHHHHHGFIGGIIHIJJJJIGJJJJJJHGIJJJJJJJIJHHHHHFFFFFDCDEDDDDD;??CDECCBBBDDBBDDDDBDDDCDDCDDDBDDBBBDDCDCDDDBCDBDDC@CA@ACCCDD<@@CDBD5>><@G@FC>B0BGFF9@FC4=BB=CCFFG4C;@CF@EFF@EEDDDD>BB8=2/(5>>@>39?A>@A:A3@B@4:>AAABAAAA:>A>:+(+29(:20959.0(8<14(+((4>@@B:A3(+53&&
+@:114:A016U:1:1:6902:6810 1:N:0:1
+GCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCA
++
+CCCFFFFFHHHHHJJJJJJJJHIJIJIJJJJIJJIJJIIJJIJJJJJJJJHHGHFFFFFFDDDDDDBDDDDDDD>>BDDDDBDDDDDDDDDDDDDDDDDDCCCDDDDDDBBBDDEDDDDADCDDDCDDDDDBDEEDECBD@BD>>B@@BB8
+@:114:A016U:1:1:23688:6834 1:N:0:1
+GGTGGCGCAATCTATTCTTCTGGTACTAACGATACTGGTGCCGTCGATTTACGTGTCACTAACGCCATGTTTCGCAATAACATCGCTAATGATGGCAAATGTGGCGCAATTTATACCATTAATAATGACGTTTATTTAAGTGATGTTATTT
++
+@@?DFFADHHFHBGIIGGHIIGIFHIGEGIGGGIHIIIGIIIIGIIIIIIIIIHGGHHHHHFBFFDDDCCDDDDDDDDDDDDACBB<@BDDDCDCDDC?DC@DC?>BDB>CDDEDDDCCDCCAD@DCDDDDDC?ACAC3>>ACDD:>>CCC
+@:114:A016U:1:1:15962:6837 1:N:0:1
+TCGACTGGAACCGACCACCGCGCCAAATGAACAGCACCTCGTTTTTCTACAATCATTCCAGCCAATGGCGCTATGAAAAAGTCTCTGCTCAGGAGTTACTTTCACCGCTCGCCGATGCCAGTAAGTACAGCGGTCATCTGATTGATTTCAA
++
+CCCFFFFFHHHHHJJIFJJGIIIJJJJJIJJIJIIIJJGIIHIJHHHHHFFFFFFEEEEEEEDDDCDDBDDDDDDDEDDDDDCEDDDDCDDDDCBAAA@CDEDDDDDDBBDDDD<>BBDDD:@ACCDCCEC><@BBDDDAD::>>CDCCAC
+@:114:A016U:1:1:18021:6876 1:N:0:1
+TATTAATGTCGGGTAAAACCACAATAAACGGTGCAATTCCCGTTTGTTGGCGAAAAAGCTCGGCTTTATGCGCGATCTCTTCAGCATCTGGTACAGTACCGGTAACGTCCGTAAAAAGTGTCTTCATTATTTGAATTAATTGCTTTCTCGG
++
+???DDDDDDEDDD:@AEI:;5:;??-&28A:>:3:>AAA:4:44>A:>>:>5>05<:<98<&+19:>B(:(44:(:4(4+(4((4(((
+@:114:A016U:1:1:23634:6900 1:N:0:1
+TTTTCTGCATCAGTTCCGGCGTACGGTAGATACCGCAGCCTTCTTCCATAGCCAGGCCCATTTCGTCGCGGATCTTCGCCCAGTTTTCGCCGCCATCCTGGTTAACCAGATCTTTCAGACGTTGTTCAACGCCAGCTGCCTGCGCTTCAAT
++
+@@@FFFFFGHHGGJJJJIIGEFHIJJIJIJJJJJGHIJJJJJGJJIIJGJHHHGHFFFBDEEEDCDDDDDDD;BDCDBBB99<@@@C@0>>BDDBDDDCCCCCCCDCDB@CACCDDDC@CBDDDD?ACAD>B<5.>4(+
+@:114:A016U:1:1:19671:6905 1:N:0:1
+AGGTATTTGCGTAGCCGATCCGCGTCATTGACGCTGGCTTTGCCCTGGCGCGAAACGTCAAAAAGCTGGCGTCCGGCTGCGGAAAGCGACTTTGCCTGGCGGCAGATAGCGATAACGTGTTCCAGTTGCATACGGTCGAAGAGATCGATGT
++
+CCCDFFFFHHHHHIJJJJIJJIIIHHGIJJJJGHFHJIJIJJIJJJJHHFFDDDDDDADDDDDDDDDDDDDDDDDDDDDDDDBDDBDDBDDBDDDDDDDDDBBDDDDDDDDDDDDDDAACDDDDDCCCCCCAC<@DDB.<<@BC>AB@BAB
+@:114:A016U:1:1:10011:6914 1:N:0:1
+GTATAACGTGCCCAGTCAAAATCGGGGAGTCCGTGGTTTAACTCAGCAATATTCACACACGCCAGGCCGATGCCAATACAAATCAGTGCTGGGGCGAGTTTTCCCCACGCGAGTTTGCCGTCGCGCTGGCTATAAAGCAGATTGGCAAAGA
++
+CCCFFFFFHHHHHIJJJIJJJJJJJJJJFGJIJIJJFGIIJJJJJJJJJJIJJJJIGHHFEFDDDDDDDDDDDDDCCDDDDCCDDDDDDEDDDDBDDD9@DCACCCDBDDBBB@BCCCDDDB7@DDB@DBCAC@CDCCCC<3>C8>::839
+@:114:A016U:1:1:11406:6925 1:N:0:1
+CCGTCCTCAATGCCAATCATCTCGATATCTGTTGCATAACCTTTGTGATAATAATCTTTAATTTCAGCAAATGTTTTAGTTTTATCTGCATCCGCTTTCTTCTGAAAGACGGGGTCAAGTTTGCCACTTTGCAGTAAAGGATAAACGGATT
++
+CCCFFFFFHFHHHJIGIJGHIGIJFHIJIJJJJJJGIJJJJJJJIHGGHGIIJIIJJJJJGIJJJIJJJGIJICGIGIIHHHHHHFFFFFFFDCCDDDCCDDDDCCDACCB;<>BB?@CACDDCCACDCCDDDDDDD>3@BCCCCCDBD@>
+@:114:A016U:1:1:5627:6928 1:N:0:1
+GCTTTATGACGTACGCAGTGGAGTTTGACAACGTCTCGCGGTTGTACGGTGACGTGCGCGCAGTAGATGGCGTCAGTATTGCGATAAAAGATGGTGAGTTCTTCTCTATGCTGGGGCCGTCCGGCTCCGGCAAAACCACCTGCCTGCGCCT
++
+CCCFFFFFHHHHHJJIJJHHIGJJIJIJJJJJIIJGJIJJJHIIJJJIIGEHHHFFFDDDDDDBDCCECDCDDDDDCDEEDDD@BDDDCCDCCDDDD@CCCDDDDCCDDECDACB<@@@9@BDDDB>BDDBD.8CBBDDD>D;@;>6;3>ABB/8?/,58A((9<&0&0)2:C44@:>C(9@?CBB24299<<1<<((+(+(&)+(((((+((+4+((+8((+:&&
+@:114:A016U:1:1:9225:6942 1:N:0:1
+TCACCGGTACATTCTGAGCGGTAAGGATCTCCAGCATTTGCTGAATGTTGTAGCCTTCCCCTGGTTTTTGCCCGCGCAACCCGGCTGTGACCGCATCAGACATGAGGAACAGACGCAGATCCAGATTGCTCTCCTGCTCTCGTAACGCAAT
++
+CCCFFFFFHHHHHJJJJJJJJJJJJJJJJIJJJJJGJJJJJJJIHIJIIJJIIJJJJJJJJJJGCEHHFFFEEDDDDBDDDBDDDDDDDEDDDDDDDCCDDBDDCDDCBD@DCBDBDBBCCCCDDDDDDDDDDDDCDDD@@@?8AB?@<.9
+@:114:A016U:1:1:3640:6983 1:N:0:1
+ACATTCAAAAAAACTCATCCCCAGGTCATCGAGCCATTGACGCAGCGTATCCAGACCAGGAACGACCAGCGATGTCATAATTTTGTACAACCTCTTTCGATAAAAAGACCGGCACAGCTTACGCACAAAGCGCAGGCAAAACCATGATCAG
++
+@@@FFFFFGHHDFIJIAHHGDHEIHGHHGECEFHH@*?BGIHIGGIIAEFEFECDE?C;?2??>BDBD?@;;2&2@@:@CC4>:A8@CDCC@AA@>C:@>CDDB&52<9?0<@5..(8+32((+4>:A@A
+@:114:A016U:1:1:21075:6983 1:N:0:1
+TTCAGAGAACTGCGTAAGTATTACGCATGTTTTCCCTGTCATTCATCCAGATTATTCCTAATCACCAGACTAATGATTCCATCAATCCTGGCGCATTTTAGTCAAAACGGGGGAAAATTTTTTCAACAAATGCTCAACCAGCATTGGGTAT
++
+CCCFFFFFHHGHHIEHIJHIIIIJIJJIJHIJJJJJJJJJJJJJJJJJJJJJJJJJGHJJJJJIJJJIGIJIHHHHHHHFFFFFFFEEECDDB@DDDDEECDDDDDDDDBDDBDDDDDDDDD@CBBDD
+@:114:A016U:1:1:6958:6990 1:N:0:1
+TTTATCGTTACTTATGTCAGCATTCGCACTTCTGATACCTCCAGCATGCCTCACAGCACACCTTCGCAGGCTTACAGAACGCTCCCCTACCCAACAACGCATAAGCGTCGCTGCCGCAGCTTCGGTGCATGGTTTAGCCCCGTTACATCTT
++
+CCCFFFFFHHHHHIJIJJJJJJJJJIJGIIJJJJJJJJJJJJIJJJJJJJJJJJJJIJJIJJJJJJHHHFFFCDCECECDBDDDDDDDDBDDCBBDDD>BBBD@C@B@?>BBDBBDDD<@CCDD@8?@CDDD>CBC>C@ABB)<@?CDDC>
+@:114:A016U:1:1:29514:7001 1:N:0:1
+CTGCTGCCCGTCTGTTTACAGGTTTCGCCGGAGATCCTCCGCCAGCGCCTGGAAAACCGTGGCCGTGAAAACGCCAGTGAAATTAACGCCCGCCTGGCGCGCGCCGCCCGCTATACTCCACAGGACTGCCATACGCTCAACAATGACGGCA
++
+@@@DFFFFHHFGBEFHIEEBDGFHIJIIGGIJGIHH9;9>>4>:3@>>@@4<-9905<075<<5@@<5>>30::ACA:<428<(++(9ACC9+985@@<(2+(4><3&&
+@:114:A016U:1:1:26852:7013 1:N:0:1
+CAGTTTGTCATCTGTTTTACGGGAAGCGAAGCTGCCTTTCTTGAGGAAAGTGGCATGCTCCGCATCAGAAGCAATAACAGAGATGCGTGTGTCATGCTCCCCCAGTTTTCCTTGATTGGATATCCACCATGTGCTGTATCCTGCTTTTGCT
++
+CCCFFFFFHHHHGJHJJJJIJJJJJIJIIIHIIJJJJJJJIJJGIJIJIJBFHIIJIC4BGBABDFE6>@>;CD;@C3>5555::5;50)22:A>>@>:4<>@&++:CD:>:>(+:::A3+:>>(:>?<:@:A(44>:CDCC3:>>>CC3(
+@:114:A016U:1:1:23707:7056 1:N:0:1
+TATCTTGAAGCCAATGAGTTGTTAACTGGCAAGACAGATGTGCTGGTTTGTGACGGCTTTACAGGAAATGTCACATTAAAGACGATGGAAGGTGTTGTCAGGATGTTCCTTTCTCTGCTGAAATCTCAGGGTGAAGGGAAAAAACGGTCGG
++
+CCCFFFFFHHHHHJJJJJJJJJJJJIJJJJJJJJIJJIJJIJJJJJGHIJJIJIJJJJJJIJJJIGFHHHHFFFFFFFCEEEEDDDDDACDD>@CDBDDEDDDDDDDEDDDDDDDDDDDDDDCDDCDDCDD9@ABBABDBDB<CBDDDDBDDDDDDCCCDDADDDDB5<@CCCCCDDD>.9?B>4>A@<)
+@:114:A016U:1:1:13916:7106 1:N:0:1
+CAAGGAAGATCTTCCGTTCAATCTGCCTTAGCTGAGCCAGACTGATTCCCAGAATACGCTGGTCAACATCAGCCACGACGGGCATCCCTTCCTTGTCATAGAAGCGACCACAAATGACACCTACTGCGCCTAAATCCCGATACGTCTGCAT
++
+CC@FFFFFGHHGHJJJJJJJJJJJJJJJJIHIJIJJJJJJJJJIJJJJJJJJJIJJJJIIIJIJJJJJJHHHHHFFFBCDBDDDDDDDDDDDDDDDDEEDDDDDDDDDDDDDDBDDDDDCDDDDCCCBD>BBDDDD>ADD&5>BC>B@B(
+@:114:A016U:1:1:25882:7129 1:N:0:1
+CAGGCTCGCGGGTAGCCAGGAAAAACCAAGTTTCTGCTGGTGGCTCTTCAAGCGTTTTCAGCAATGCGTTAGCCGCGGCGTCGGTTAGTAAGGCAGCATCGGTTACCCAAACGACTTTCGCACCACCTAAGCGTGCGTGCTCATTCAGCTT
++
+CCCFFFFFGHHHCGGIJGIJIJJJIGIIGIEHIGIIIGJJGGHG@CDFFFFEEDBDDDCDDDDDDDDDDDDDDDDDDBBD>BBD9>@?C<@BBB79@0:C@CD>4(+4
+@:114:A016U:1:1:6958:7183 1:N:0:1
+AGCTCAGGTAACGAATAGCCGCTTTCTCTTCAGCGCTTAAATCCAGCGCGCGCAGCATACCTGCATTTTCCAACTGGGGCAGGCATTTTTCCAGGACCGCCAGCTCAGCTTTCGCGTCTTTATCGCCACCTTTGGCTTTCTTCTGTACGCG
++
+<@@DFDFFCFDHHE;EEFFGIIGHIGIIGGHCDEHGGIEIAGEIC>AHHHHDBA???CA@CDCB<>C:CCDEDDCC?A>BB<@@B@@C>C>C>33?A?@><<9BDDDDCCDDC05@59@<>CC3(8)(0++4443>@9<33(:>34@CA>CCBD7<3ACBD8AD34(::>ACA35&&)(4:@3
+@:114:A016U:1:1:9644:7202 1:N:0:1
+CGCAAGTTTGAGACTCACAAATCCACGATGACAGGCGTTCATTAACGTTAACAACAGTGGTCCGCACAGACCTGCCAGCACACCTGTACTGATAATCAGCGCATAGTCACGAGCCTGAACCGTCACTGAGAGCTGTACGTTGTAGAGTAAC
++
+CCCFFFFFHHDHHIIJIJJJJJJJJIIIJJJJIGGGIIJJIJJJJJJJJJJJJJIJJGHJJJHHFFDDEDDDDDDDDDDDDDDDDDDDDEDCCDDDDCDCBDBDDDDDDDDDDDBCDCCCBDDDDBCC>@CC(:ACC@A?CDDCCEDAA>:
+@:114:A016U:1:1:9728:7211 1:N:0:1
+CTTCTCAGCTGTGTTGATTTTCGTGGTGGTATGGCTGACGCTCTCTTACATTCCGATTGCGCATATGGTGTGGGGCGGTGGTTTGCTGGCTTCTCACGGTGCGCTGGATTTCGCGGGTGGCACCGTGGTGCACATTAACGCCGCAATCGCC
++
+CCCFFFFFHHHFHIJJJJJJJJJJJJHIJEHIJJJJJJJJJJJJJJJJJIJJJJJIJJJJJGFFDEFFCEDDDDDDDDDDD@DDDDDDDBDBCCADDDDDDBDDBDB?CCCDDDDDDDDDBDDDDDD@A?CDD@C@>CC@.9>BBDBDCC:@@8A???CCACC?@(283
+@:114:A016U:1:1:8975:7259 1:N:0:1
+ATTGATATCAACCATGATATTGCCAATGTCACCGAACAATGTCGGCAGGCTGTGCTGGCGATACGACAAAACAGAATATGTGCGAAAGAGGGCAGCGCCTCAGATCAGCGCTGCGAATGATTTTCAAAAATCGGCTTTCAACACCACGCGG
++
+@BCFFFFFHHGHHGIJIIFIJHHIIGJHJJJIIGIGIJCHIIJJJJGIIJJJJIJGIIHHFFDDDDD885?B5>5@CABBBDDDDDDDDC?CCDCDCB@<@DD5<9?:@>CCDAA:A8?CBDDBBCC@@CCB?8?BD<<
+@:114:A016U:1:1:18574:7292 1:N:0:1
+TTATGCATAGGGATGGTGCCAACACATTTCGCCAAACCGTGGCTTAATGAGGGGAAGTGGGTAGCATTAGAACTGGAAACCCCCTTTCCGGATTCGGCATGTTGTCTGACATGGCAGCAAAATGATATGTCGCCTGCGTTAACCTGGGTGC
++
+?1?DDD?DFFFFF:EDCBEGDGEB?>AA(5<=BB8AA:<9-5<>??B>>B>>:9>ABB@:3>??>B94>:3>:(::?B<@8195(08::>CACCCC>C<>BBB(943>:8<@A
+@:114:A016U:1:1:17779:7305 1:N:0:1
+CTCCTGATTCTGGCACTTTATTCTATGTCTCTTTCGCATCTGGCGAAAAGTCGTGTACCGGCAAAGGTGCAGTCGTTATATACATGGAGATTTTGATGGCAAGCAGTAACTTGATTAAACAATTGCAAGAGCGGGGGCTGGTAGCCCAGGG
++
+CCCFFFFFHHHFGJIIJJJIJIIIJJJJIIJJJJIIJIJJJJJJJJJIJJIIJIJIGIGIHHFFEDECCDDDDDDBDBBDDEEDDDDDDCCDDDDDDDDDDDBDDDCDDDEEDDDDECCDDDDDDDDDDDDDBBBBBBBBB9@CC>>?B?<
+@:114:A016U:1:1:13112:7307 1:N:0:1
+TCACGCGTGGAAGCCACACCGCAGGACATTATCGCCTGGACCGAAGGTAACGCGCTGGTCGCCACGGGCAGCCCGTTTAATCCAGTGGTATGGAAAGATAAAATCTACCCTATCGCCCAGTGTAACAACGCCTTTATTTTCCCGGGCATCG
++
+CCCFFFFFHHHHGJJIJIIJJJIJJJJGJJJGIIJJJJJJJJJHFFFCCEEDDDDDDDDDDDDDD@DDDD@BDDD@BCBDDDDDCCCDCDDDDDDDDDDDDDAACDDCCBDCDDDBDDDBACD@ACCA?>@>8@CDC@DDD:@BBB><<9<
+@:114:A016U:1:1:21363:7316 1:N:0:1
+ATTCTGCCGTCTGGCATTGCACCTAGCTACTTGGAGCAATATTCCAGATTTGAATGACGCCTACGTAAACATTTATCACTGACGGAAGAGTCTCTCTGGCAAAAATGTGCACAAAAACTTATTGCCGCAATTCCACATATTCCAGAATGGC
++
+CCCFFFFFHHHHHJIJJJIJIJIJJJJJJJJJJJJJJJJJJJJJJIJJJJJJJJJIJJIJJJIJJIHHHHHFFFFFFFEEDCEDDDDBBDDDCDDCCDCDDCDDDDD@CCEDDDDDCDDDDDECDDDDDBDDCDDDDCCDCDDDCCCCC:(
+@:114:A016U:1:1:15732:7342 1:N:0:1
+GATATCCGCCCATCCTTCACGGTCCAGTACGATACCAATCGCCTCCGGTTTGTGTCGGAGGACAAAACTTAAAAACTTGCTCGTATCCGCCAGCTCTTTTTCGTTGTACTTTGCCATCTCGTTTTCGTCCCTGGATTTTCTGTTTTTATTT
++
+CCCFFFFFHGHHGIJJJJJJJJHJJJJHIJIJGGIIIJIJFGIGIGGEEHIGHHEFFFCDBDDCDDDDDDDDDDDDDCDDDCDDDDBDDDDDDDCDCDDDDCDBD@DDDABCCAC>CCCDD<@C?8@ACCDCDC>CCB2>:@
+@:114:A016U:1:1:8807:7344 1:N:0:1
+CCTTATCCGCTGGAAACCATGCTACGCATTCACTGCATGCAGCATTGGTACAACCTGAGCGATGGCGCGATGGAAGATGCTCTGTACGAAATCGCCTCCATGCGTCTGTTTGCCCGGTTATCCCTGGATAGCGCCTTGCCGGACCGCACCA
++
+CCCFFFFFHHHHHJJJJJJJJJJJIGIJIJJJJJEIJJJJJGJJIJJJIGIIGIJIJJHIIJHHGFFDDDDDDDDDDDDCDDDDDEDBDDDDDDBDBDDDDDDD<CDDDCCDD@9:B:>@ABA@ACCCAC:CDA>C::::CA<<>C@C::::@C@@C8?CACDDDDACCCCA@CCD:@CACCDC:@>(2+8@BDEEE:@A:>>CC
+@:114:A016U:1:1:14078:7399 1:N:0:1
+TACTACAACCGGGCTGACACCCGCGGCCCGAGCCAGTTGCAGTTCAGAAAACGTTTCCCCTGGGCGTAACTGGCGCTGGTTTATCATCGTGAAAAAAGCCTGCTCAAACACTTTGTTTTGTTCACTCATTGAAGCGGTAGTACAGGCGAAA
++
+CCCFFFFFHHHHHHGIIIJJJJIJGIIJJBEEFF@BEEEEDDDDDDCDDDDDDDDDDDDDD?BDDBDDDDEDCDDDDBD9@BDDDDEDDDDCDCDDD59@@DD@>9?@BB&
+@:114:A016U:1:1:5559:7408 1:N:0:1
+GACCCAGGCACGGCAAAGCCTGACACTGGAAGAGCTTTATCAGATGTTCCGCGATGGCGCGGACGGCGTAATCAGTCGCGGGAAAGCCGAACCTGGCGATAAAACCATGTGTGATGTGTGGGTGCCGGTGGTGGAATCGTTACGTCAGTCC
++
+CCCFFFFFHHHHGJJJJIIJJJIIJJJJJEIIIGIJJJJJJJJJJIIGHJJJJIHHFFFDDDDDDDDDDDDDEDDEDDDDDD<>B?CDB<95:A@B>>@+4:8?95<8ADB
+@:114:A016U:1:1:7922:7411 1:N:0:1
+ACGTAGCGAAACATTTTTCACGCAGAGCATTAGCGGAGTCCTTGTCTGACATAACCTACCCATTTTGCGACTAACGCGGTCAGGGCAAAAATAATAAGCGGTAATAGCAGTTGCCATAAAGCCTGGGCGGCGAGAATACCGTTGCTGCCGC
++
+CCCFFFFFHHHHGJJJJJJJJJJJJJJJIJJJIJIJIJIJJJJJJJJJIJJJJHHHHHHFFFFEEEDDDDDDDDDBDDDDDDDDDDDDDDDADCEEDDCDDDDDDCDDDD@CCDDDDDDDDC@CB<5@CA@B@B@@CDCB3<
+@:114:A016U:1:1:25950:7418 1:N:0:1
+GCAAAGAGAAACTGCCGAAACCGCGGATTTCAATACGCTCGCCCTGCGCAAGAGTCGAGGCCATATGCTCCAGCATCTCTTTTACTGCATCTTCAACCGTCTTGGCGGGAAAGTGGGGTTGGTGGGGGGGAAGGTTTTCTTTCAATTTTTC
++
+CCCFFFFFHHHHHJJJIJJJJJIFIJHIIJJJJJJJIJIIIHHFFFDDDDDDCDCDDDDDBDDDDDEDCDDDDDCDD@CDDDDDCDDDCDDCDEECABDBDDDDCCBB<&0(((+(())00&((2955>&&&((+++8(:(++(((((+(+
+@:114:A016U:1:1:26470:7461 1:N:0:1
+GCGACGGCGGCTTTAGCGCCCGGGGGCAAGCGGGCCAGATGCTGAAAAAAGAAGGCGTCACGCTGATGGTTGAAGCCATTCACGCCAGCCCGGGCAGTGGGTTTACCGTCACCAAATACCCCACCCTGGGGGTGATCAACCAACTGCAAAA
++
+??AA>>AA<<=(55<7<7=;>7>A:(489AAA?<:>:>A8>?>>>>?950;5>+(++0283>(52980((+2443(:((5&0(9<2.&098(4>>9<>80(((+:A9
+@:114:A016U:1:1:5847:7548 1:N:0:1
+CGGTGAATATAATTCCATCCAGGCGACGTAATGAAGCTGCGTGTCCGGCAATATGACGGGCAATTCGGTGAACAAAGGTTTTAATTGCCAGTTGCGCGCGTTCGTGACCTTCATGCCAGGCTTTTTCCAGAACACGTAAATCCGAAGAAAG
++
+@CCDFFFFHHHHHJJIJJJJJJJJJJJJIJJJJIJIJJJJJJJJJJJIIJHGHHHFFFCDDDDDDDDDDDBDDDDDDD4@CDDDDED>ACDCCDDDDDDDDDDDBDDBCDDCCCCCCCC?CDDADBD@CCDDBDCDEEEEDACDDD@CCDD:@CCD>>C:@CC@@3?>@B>9B99@>ACA>:::A>C>>>CD?CAC>C:>>@A8<@>39&++4422<<&9?>(+:CACA(3(((4:(.0&)28<5(+(+(
+@:114:A016U:1:1:19742:7606 1:N:0:1
+CGCCGCCTGGCTGCCAGGCTTTAAATTTCAGACCATCACGCATACCACCGGCGTAATCTTCGAGGATCTCGCGTGCGGTGGTGCCGAACGGCAGTTCCCACAGTCCCGGATTTTTCACCCGACCGGAGAAGCCCATCAGCTTGGTGCCAGC
++
+CC@FFFFFHHHHHJJJJJJJJJJJJIJJIJIJJJJJJJIIIIJJJJJJJHHIFECDDEDDEDDDDDDDDDDDDBDDDDC:@C+49<48:<>C?
+@:114:A016U:1:1:13676:7613 1:N:0:1
+CTTCAATTATCTACCGTTTATTGCTGATTGGCCTCGGATCACCAATGGTGATTTACGGCCTGGTTCGCCCGCTCACCATCGAAACGCGAGATTAAACTCCTGACAAAACACGCCCCAGAAACGTACTTCTGTTGGCGTAGTCTGGGTTATT
++
+CCCFFFFFHHHHHJJJIIJJJJJJJJIIJIJJJJIJJGIJJJJJIJJJFHIIJIJIGGIIGHHEHHFFDDDDDDDDDDDD?CDBDBDDDDBBDCCDBDDDBBCDA??BDACDCDED<<@B@DD@DEDCA(2<>>
+@:114:A016U:1:1:23337:7628 1:N:0:1
+CTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGACTACATCCC
++
+@CCFFFFFHHHHHJJJJJIIEHGIGHJJJJJGJJJJIJIJIGIJJIJJJIJJJGHHHHHFFFFFFFEDEEEDDDDDDDDDDDDDDDDDDDBDDCCDCCDDDCCC<@B9CBB>@BDCB2<@B<58?B(((+4:(+(
+@:114:A016U:1:1:3424:7641 1:N:0:1
+GCATTTGCCGCACAGTCGCTTTGTCGCCAACCAACTGCACCAGTTGATCGATACGATCCTGCTCCAGCAGGTCTTCCGCCCGCTGTTGCTGAAGCGCCAGCCCGGTATAGAGGAATGTCAGACGCGAGCGTTTACCGCGCCCGACTTACGC
++
+@@@DFDFFHHHGHJIIGGIIIJFDHGIJIIAHIJJEGGHIEEGFHEIGIIIIJHDHFFEECEEDCEACAB?A@CDDCDDBBDDB99CCD@@8?ACCDEDDCDBBB5<9@98>?::>BDB@<9@@&((+85
+@:114:A016U:1:1:21156:7682 1:N:0:1
+AATCTGTGGCTGGTTAATCCACTGGTGAAACCCGACTCCCAGGCGATTGCCCTGGCAGCAGGGGCCAGCGATGGGATTGGTAATGTTCACATCGTCGCTACTCGCGACGAAGCGTTCGCCGGTTGTAGCCTGGTGGTGGGCACCTGCGCAC
++
+@@@DDAD4CCB?8CE:AEC?3AH<@++9(+(:3(4(44+8(&&029(:())0599.+0&(2(&)9@953(+4(+8:<()87.(((8(++&)5.
+@:114:A016U:1:1:26391:7690 1:N:0:1
+GGTCATTGAGCCGAAAACCACGGTACCGCTGAATGTGAAAATGGATAGCGTACCGGTGCTGACTTACGTCAACGACTTTGGCGCACGGATGCCGCTGTTCTTCCAGTGCAATGGCAACAGCTGCCAGGGGGATGAAGAGCAGAGCCGCAAG
++
+@CCFFFFFHGHHHIIHIIJGIIGHHJJEGGHIJJJHJIJJ@GIIFCGIIJHJIHHFFFDEECCDDDDBBDDDCB?@@DCCCCDD>BDB<:3+:AC:@>ABDDDD(24:C<@>?&05..2:::((+433<>1&090&
+@:114:A016U:1:1:13571:7695 1:N:0:1
+AAGAAAACGTACTTATCGGTATTACCCAGGCTCCTTATATTATGGGGCTGGCTGATGGTTTAATGAAAAACGATATACCATACAAATTAATATCAGAAGGAAATTGTACGGGATATCATCTACCAGCCAAACAAACCATTACGCTTATTGA
++
+CCCFFFFFHHHHHJIJJJJIJJJJJJJJIJJJJJIJIJJJJJJJJJJJJJJJJJIJJJHHIJGHHHHHHFFFDDDDEEDDDDDDDDDDDDEEFEEECCCDDDDDDDCDDDBBBDDDDEEDDDDDDCDDDCDDDDDDBBDDDD@@B8?BCC??CDDD8>B9>B
+@:114:A016U:1:1:14873:7780 1:N:0:1
+CTCGCAATCGTGGCTGGCTGCTGCTACCCATTGTCTATGCCTTTGCGCTTATCAACCTGGCGCTGGCAACCTTCATGCCCAATGAATACATCACCCATCTGGAAGCTACGCCTGGGATGCTGGTGCACATTGGCTTATCGCTCTTTTCCTA
++
+CCCFFFFFHHHHHJJJJJGJIIIIHIIJJEHIGIIJJIIJJFIJIJJJJIGJJIJJJIHHHFFDCDDCCDDDDDDECDDDDDDDDCEEDCCDDDCDDDDDCCCDDBCCDBBDBBBD@DDDDCCCDDDDDDDD>CCDADDDDBB<@CD>CC(
+@:114:A016U:1:1:21523:7846 1:N:0:1
+TTTCGTCAGCATCGCTGGGACCTGGAAGAAGAAGCGGAAAACTTAATTCAGGAACAAAGTGAAGATGATCAAGGCTGGGTCTGGTTACCCTGATCCAGATATTCGTCCTTCCATTTCACGTAATTATTCGCGGAATAGCGTAACCCAGCCT
++
+@@@FADD;FDDHHE?GIIIGII@BFFHEHIG@EFGHG@G;CHIIADHGGHGIIHHHFHC>>CBCFCCCACEDDDBDD8@CCCC?:@@@CACCCDCC>CCCDDDEDD@B?C>?CACCDECD?<@DDDCDDDEDDD>@BAC>@@DDDDBBCC<9<8@>CBB?8C>::@AC3:
+@:114:A016U:1:1:28887:7858 1:N:0:1
+CCTCGAAAATCGCGCCGATATCGCCGTACACTCAATGAAAGATGTGCCGGTTGAATTCCCGCAAGGTCTGGGACTGGTCACTATTTGTGAGCGTGAAGATCCTCGCGATGCCTTTGTGTCCAATAACTATGACAGTCTGGATGCGTTACCG
++
+CCCFFFFFHHHHGDGJIIJGJJGIIIIJJIJGHGGGIJJJIHHHHFHFFDDDDDDDEDDDDDD;BDCDDDDDAABDB@CCCDDDEEDCD@@>55@CCCAC:@>>>CCDCD:A85??A0+@>BBBDD5<@DCBD?>BBDCCAA8?ACCCCCDDDD58A3@3:88?5>.>8ACCACC22>CCDABBC
+@:114:A016U:1:1:25479:8033 1:N:0:1
+AGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCG
++
+CCCFFFBFFHHHHJJJJJGJJJJIJIGJHIJIJJJHEGGIJIJEHIJGIIJJJEHFEEHFFFFDDEEEED?CDDACCCDC,>@?BDBDDBBDDBDDCDEDAD4>>A:(3(4(4:<9AC@?@>39.99
+@:114:A016U:1:1:21544:8052 1:N:0:1
+GATTATCTGGCGGACCTGTTTAAGATTGATAACTCCCACGGCGGGCGGTTCAAAACCCTGCTGTTAACCCTCCTGCCCCCCGCGTTTTTTTTTCTGATCTTTCCGAACGGCTTTTTTTTCTGGAATGGTGGGGGCCGGCGGGGTGCCCCCC
++
+CCCFFFFFHHHGHJJJIIIJJJJJJJJJJJJIIHIJFIJJJJIF?50&09((+9A(2((+(+(+:(4(4(2+2(((8(+9&)0&&&&+2)0&&(+(+44:>(+>(&+)))&&)94:@<&3+:3((((++(2&&&)&00&)&)&&)2++2&0
+@:114:A016U:1:1:24871:8061 1:N:0:1
+CATCAAGCCGCTGTATACCGAAGCCGATCTCGATAATCTGGAGGTGACAGGTACCCTTCCTGGTTTGCCGCCCTACGTTCGTGGCCCGCGTGCCACTATGTATACCGCCCAACCGTGGACCATCCGTCAGTATGCTGGTTTTTCAACAGCA
++
+CCCFFFFFHHHHGHIIJJJJJJJJJJEGHGIJIJJJJJJJJJJJHIJJJJJEHHHHHFFFFFEEEEDDDDDDDDDDDDDDBDBDDDD>B5095A:A:(:3+4>CAC<@BDDA2902AB&+8(9@@@>:@@A9>A<8+(4::(8+
+@:114:A016U:1:1:6805:8063 1:N:0:1
+GTTATGAAAAAGCTCATCGATATCCGTAACGGTAAAGGCGAAGTCGATGATATCGACCACCTCGGCAACCGTCGTATCCGTTCCGTTGGCGAAATGGCGGAAAACCAGTTCCGCGTTGGCCTGGTACGTGTAGAGCGTGCGGTGAAAGAGC
++
+CCCFFFFFGHGHHIJJJJJJJIJJJJJJJJJJHJIIIJJJJIIIJJIJJJJJJJHHHFFFDEECBDDDDDDBDDDBDDCDDDDDDCDDBDB3>@BCDBDDBDBDCDD@CDCDB79>CC9
+@:114:A016U:1:1:26095:8095 1:N:0:1
+GCGCTTGACCGCCAGCCGGAGCTGATCACCGAAGCGTTGAAAAACAACATCATGCTGGTTAGCCCGACTACGCTGCTGGTGGCGCTGCGCACTATCGCCAACCTGTGGCGTTATGAGCATCAAAGCCGCAACGCCCAGCAAATCGCCGATC
++
+@@CFFFFFHGFHGJIJJJJHGHIIJIJIGHHGIIGGFHIIHHHGFFEFDCECEDDDCDDCDACDDB=BDDDDDBB<:ACDDDDDBBBD?B?BCBDCCABDBB<@C@CAC@ACD3??5<0&598<55&
+@:114:A016U:1:1:7393:8099 1:N:0:1
+AAGCGGTGCGCGACGGCATGGTGGTGAATCTGTTGCTGGAGGCAATTCCGGATCTGCATTTGCGTCTTGACGCAAATCGCGCCTGGACACCGCTGAAAGGTCAGCAGTTTGCCAAATACGTTAACCCGGATTATCGCGACCGCATCGCGTT
++
+CCCFFFFFHHHHHJIJJJJJJFHH@GHIJJJJGIIJJJJHHHFFFFFEEDDDDDDDCDDDDCDDDDDDDDDDBDDDDDBDDBBDBDDDD@DDBDBBCCDD@CDDDDDACDCDDDCDCDDDCCCBDDDBDBDD:>C?@<>BDD@>DD5<>5<
+@:114:A016U:1:1:27220:8127 1:N:0:1
+ACTTGGGATTTGATATCTTTTCCTGAAATACCTTTTTGACCAGAAGTTGTTTTCACGCAGTTCATCATTGCATTTCGTGCTGAGGGGGATGAAAAAAATATTTCGATATATTCTGGTAAAGCATCTTTGGTTAATCGAGCTCGAATAAGTT
++
+@CCFFFFDHHHHHJFIGIJJJJJJJJJJIJJJJJJJJIIJGIJJIJIIJJJJJIIJJJJIJJJIJJIJIGGIJHHHGHFFFFFDDDDD@BBDDDDDDDDDEDEECB@BBDCCBCD<:>CACCDDDDDCD(8?CCDDBD?@BDDDCCDC<(4
+@:114:A016U:1:1:24859:8129 1:N:0:1
+ATACAAGCCGATTGGCCCGGCGGTGCACAGCACCAGATCTACCCAGTCCAGCGCCTTGTCAGATGCCAACAACGGATCGCTTTCTCCGGTCAACGCTTCGGCTTCATCTTCGTTCATCGCAAGGATAGAGACGTGATCTTTGAGGAATTGC
++
+CCCFFFFFHHGHHIJIJJIGHGJAHIGHIIJIIHGGEFFFFFFECDCEDDDDDDDDBDCDECCCCDDDDB?BDBBB5>@BBBDDDCA>BDDDDD>BBDBDDDD@CCDEDCD1<9A?::9@@BBDAACC>:>?B28BDDDDBD<@CDD@BBCCDDDCDBD5<@B@BC:ACAA:ACC::>88
+@:114:A016U:1:1:7250:8289 1:N:0:1
+CTTAAGGAGACAGCTGATGGCAATGTTCGGTTTTCCTCACTGGCAGTTAAAATCGACCTCTACAGAGAGCGGCGTGGACGCGCCGGGGGGACGAGAGCCGTTGGGTCAGAAGGCAGTGGTGGGGGCGTAACACCAGGAGGGGCGGTGTGTG
++
+B@@FFFFFHHHHFJIJIGIJJJJJIHJIGGHGIJGIJJIJEBEEGGBGHGGJIIGJIGGHGHIFEEAHEDDDDDD&0((&&0))0))&0&)&&)&&&((&&)((&(+(((++(++(+((((()&0&0&)22:>(((+(()))&)&)&)((+
+@:114:A016U:1:1:18589:8298 1:N:0:1
+GCCGCGCTGCACCACTTTAGTATCGCCAGTAACGATGGCAATGCCTGCCGCGCGGGCGGTTTCTGCCATGCTGGTCACTACGGCTTTCAGTGTCTCCATCGGCAATCCTTCTTCGAGGATAAAGCCACAGGAGAGATAGCGCGGAAGAGCG
++
+@?@BDDD@FF?FHAFEHGBFCGIC=FIGIBFDFDG?B;CGGCHIIIEHGGEBC>;=BBB7@BCCCC@44:>>>>CCCCC@C<@ACCCCDC@CCCC<@99CCCB<9@09&8A<9
+@:114:A016U:1:1:28618:8377 1:N:0:1
+GGGTGAGGCGTTAAAAGATGCGGTTGCACGGGCCGATGTGGTGCTCGACTGTACCGACAATATGGCGACTCGCCAGGAGATTAATGCCGCCTGCGTGGCACTCAACACGCCGCTTATCACCGCCAGCGCGGCCGGATTTGGCGGTCAGTTG
++
+@@CBDDFFHFHFHJJJJIIJIIJJJJJJIJIIIJIGHHAEH>BDEFDCDDDDEDDB@DDDDDDEDDDDD9@?BB@BDB@B?C3>(9<4<5<@09><<44:AC09B99BD&0&)))0:>A@<<&)03+4(
+@:114:A016U:1:1:8252:8389 1:N:0:1
+TTCGCCTTCACATCTTCCGGCAGAGTAAATGTTACCTCAGTGTTGGCTATCGCATTGCCCTCTGTATCAGCGACTGTTGCTGTTAGTGTCGTGAGATCGTTATTGTCGGCAATCACCGGATCCTTCGAGGCTGCCAGCGTGATTTCAGCAT
++
+CCCFFFFFHHHHHJJJJJJJJJJJJFHIJIJJJIJJIJJJJJIJIJJJIJJJJIJJJJJJJHHHHGHFFFFFDDDDDCDDDDDDCDCDEEDBBCBCCDD?BBDDCDD<@BBBDDDDDDBDBDDDDD@@99C
+@:114:A016U:1:1:14873:8438 1:N:0:1
+ACCGGCGCGGGCAAAGCGCAGCAGTTTTTCTTCAACATCGGCAGGAATCGATGTGTCATCGGCATCCTGAACGTCATGACCGTAAATGGAGAATGCTGGGATGCCTTTCTGGCTGTGAGCTGCCCGAGCCGCCGCCCGGGAAACAGCGCCG
++
+CCCFFFFFHHHGDIJGIIJEEHAHEHHHHFHFDBEDDCEADDBDD2:3:3:++55059?4(4@:22<(2(+4>&)+0?:4:(23?3(+4:@@((+(::@C(4+:4<7++(43>9(<()&))5&)&&&50&5&&&2>A9<9050
+@:114:A016U:1:1:4129:8474 1:N:0:1
+ATCTGCGGAACACCGACGCCAGTCACGATACGAGTTGTACAGATAGAGCCAGGGCCAATGCCGACTTTAACCGCACTGCAACCAGCTTCTGCCAGAGCGCGTGCACCTGCAGCTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGAT
++
+CCCFFFFFHHHHHJIJJIJJJIIJJJJJJJIJJIIIIJJJHHGHHHFFFFDEDDDDDDDDDDDDDDDDCDDDBDBDDDDDDDDB?BCCDDACCDCC4@A>(+4@?325&9
+@:114:A016U:1:1:16373:8500 1:N:0:1
+GATCGTAATCATCGACATCTTCGCGAACCAGCGCGATATCAATTTTACCCGGACCGACTTTCCAGTTTTCCAGACCTACACCCGCTGCGGCATCAGTACGCTGCGTTTTCCAGTCAAGCATCTGGATTTCAATTTTCGGCGCACCGTGTTT
++
+CCCFFFFFHHHHHJJJJJIJJJGJJJJJJIGIGIJHIJHHGHHHFFFFFFDDDDDDDDDDDDDDDDDEDDDEDDDDDDDDDDDBDDBDDDDDDDCDDDDDBDDDD9B@DDCCDDDDACDDDCCCDDCDDDDCCDDDE9>BB9@BDBA?A@C
+@:114:A016U:1:1:21659:8560 1:N:0:1
+ACCACCGGTGAAGTCGCTGACGATGCCGCCCGCTTCACGAACCAGCAGCTCGCCTGCGGCGAAGTCCCACGGGCGCAGACCGATTTCAAAGAAACCGTCAACACGACCCGCAGCGACGTAAGCCAGCTCCAGCGCCGCACACCCGGGACGC
++
+@CCFFFFFHHHHHGHJIJIJIGGI?FHGIIIJI@FFCGHIIGFFFFFAE(>@BBBDDD@B;;B@BCDDCD>BDBBD@>BBDB<5<@DDDCCCBCCDB89BBC4?&595>5@9&&&590)88<>>98(+(:4+3)59.0.&)(()0&)&0&&
+@:114:A016U:1:1:29672:8572 1:N:0:1
+ACACAGACGTTAAGCCAGCTTGAAAACAGCGGCGCTTTTATTGAACGCCATATCGGACCGGACGCCGCGCAACAGCAAGAAATGCTGAATGCCGTTGGGGCACAATCGTTAAACGCGCTGACCGGCCAGATTGTGCCGAAAGATATTCAAT
++
+@@@DFDDDHHHFHGBGIJIIHIHJIJIJHIEGA@FGIIJIIJHHG>=BBBDDBBDC?CCC@CCCCC>>::(59@B@&59<<92++8BC99<9.5+2>((+4(())55&3(+(4(>:(
+@:114:A016U:1:1:23895:8600 1:N:0:1
+AAAGCCCTGAGTGATTTGTTCCGGCTCCAGGTGCCAGAGGCCGAACTCACGCACTCCTGGGGAGTCAATCACATCACCGCCGTGCGGGAAGTGATACAGCCGAGCGGCGGTCGTGGTGTGCTGACCGAGCCCCGAGTTGTCCGAGATATCG
++
+CCCFFFFFHHGFFEGIIDIGIJIFG8EEFHG@BGHHGIJIJJIJFGDCCH?EAEAHFFFFFCABDDDCC>CCCA@ADDDDBB@<@DDCD<@(+4>>ACCD<@09<(8+8::>:?A>3.95&)(2@4(((&)&+:>:(
+@:114:A016U:1:1:19835:8615 1:N:0:1
+TTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCG
++
+CCCFFFFFHHHHHJJJJJJJIJJIJJIIJGIJIEIIIIIJDHIJJHIJJHGJJJJJJJJIHHFDDDDDDDEDDDDDDDDDDDDDDDDDDDDDDD>CBDDDDDDEDDDCDDDDDDDD02>CCDDDEEEEDCDDBDD@BD<>ACDD?>?CCC5
+@:114:A016U:1:1:13445:8638 1:N:0:1
+TGACTTCCAGCACTTTGAGCGAACCGAGGCTGGTCATCATCTGGAAAGCGCGGGTATAAATAAAATCGCCTACCAGCACGCTGGCGGCATTGCCAAATGCGGCGTTGGCGGTAGCTTTACCCCTGCGCATATCTGATTCATCCACAACGTC
++
+CCCFFFFFHHHHHJJJJJJJJJJIJGIIIJJJJGGIIJJJJJJJJIJIGIHFDD?BDEEEEEDDDCC?BDDDDDDDDDBBDDBDDB9B@CCCCCCCB??BDBD@?CC@@CDCDDDDCDD>889A<
+@:114:A016U:1:1:29588:8669 1:N:0:1
+TTCGTATGTGCATATTAACGTGGTGACATCGTTTGCCATCGACAATACCCACAGCGTTGCGCTGCAGCAGGATTTCTTCGATAAAGTCTTCCCGAACATTTTACCGATGGCCTACACCCTGCTGATGTATTACTTCCTGCGGGCGAAAAAA
++
+;;@DBDDDHHHFGIBEGGEFGHIBEHGIGHHICG?ADFHAHGHFBDGIBBGFGGIBEGIHEDDACCE@?=;;>@CCEACBBBABC>>>CD>CCBBBB@CCE:::>BB.>BCCC?CCCC.2(234:@CDDDCD:>@3::>CB@B&599@B@.
+@:114:A016U:1:1:5487:8688 1:N:0:1
+CCAGAATGACGCGGTAACCCTCTTCACGCAGGGCTTTACACGCTTGCGCGCCAGAGTAGTCAAACTCACACGCCTGACCGATAACAATCGGGCCCGCACCCAGAATCAGGATACTTTTTATATCTGTACGTTTTGGCGTGGCTCTTTTACT
++
+@@@FFFFFHHFHHIJJJCHGGIBHIIFGGGGEIIGG;FHIJGI@GGIGHF@B?BBB>CD(;AACCDDDCCC9;:CCB<@95>B>BBDDDD8>>AC3>AC@@:@:>>?CD>C@C@@C??CC?<&:5>5532:+:C:>A:
+@:114:A016U:1:1:19677:8712 1:N:0:1
+ACGTTGTAGAATCTTTCTTGCGTTATTTATTGTCACCATAAAACACCTACAGATTAATGATGAGGTCTACGGCCTCAAATCGCCAATATTCAATATCTATTTGGACAATATTGCCTTGCTGGTCGGCGCGATGCTTTTCAATTAAGATCGC
++
+CCCFFFFFHHHHHJJJIJJJJJIJIJIJGIJIIJJJJJJJJJJJIJIJIJIJIIJJJJIIJIIIJJIIIIGIHHFFFFEEEBDDBDDDDEEEEEEDEEEFEDEDDCDBDCCDDCDDDDDDDCCD>@BDDDDDDBCDCDDCCCDC>ACCD=8
+@:114:A016U:1:1:26176:8740 1:N:0:1
+CCGCATCATACTGGCTCTGAATAAAGCGTTTATCAATTAACTGTGCCAGCGTTTCGCTCTCTATCACGAAGCCCATTTCAGTATTGAGCAATGTTGAACGCGGATCGAAATTGAAAGAGCCGATAAACACCGTTTTACCATAGATGCTAAG
++
+@?@DDDDDBHFG@CEGCHEEABB?CCABCCDDAC-5@C(5:@>>;>:C?B>>509@BCBBC4@>A:>959?@:@B<0+29@C:ACA((4(+:A@3(
+@:114:A016U:1:1:18969:8753 1:N:0:1
+CAAACACCGCCGTCGATATGAACTCTTGGGCGGTATCAGCCTGTTATCCCCGGAGTACCTTTTATCCGTTGAGCGATGGCCCTTCCATTCAGAACCACCGGATCACTATGACCTGCTTTCGCACCTGCTCGCGCCGTCACGCTCGCAGTCA
++
+CCCFFFFFHHHHHJJIJJJIJJJIJIJJIIHIJHHGHHHHFFFFFEEEDEDDDDDDDDDDDDDDDDDDDDDDDDDDDDBDDDDDDDDEEEEDDDDDDD@BDBDDDDDDDCDDDDDDDDDDDDDDDDCBCDDDBB@9<@B@9B5>DDB<BBBBD>B9CCC9>BB99<>C5<@0<>@59AC?B55@838?9@C7>?<
+@:114:A016U:1:1:14416:8772 1:N:0:1
+CCGTCGCTGTTGTCCTGCCAGCTATAGCCAGCCCAGAAACGAACACTACTGACCGCTTCATTCGGCGCTTTGACTTCCAGACGATAAGCGCCCCACTCTACCGGAAAACTGACCTTGCCGGTTTCGTCCGCTTTCAGATCCAGAGTTTGTT
++
+@@@DDDA?FHHHFEHGDFECEHGGGGGCEGBBGGF@??BDDEHH@FHGIEHCEHBHEDF<;C>;6>B=2=BCCCD>ACC>@?BBD8?>:<9<<@BBDDC@CDBB9BDDDA:A>9344>909099:88(8.95&++>34:33:>(8:::@89
+@:114:A016U:1:1:14193:8777 1:N:0:1
+AACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATC
++
+@@CFDDDADB>AFGHGFEBF>FHCHEFGGEHFGDHGGIIGFEHE@BGHICEHA@CCDDFCEB=@DDD@CCDCCA?AACAA?BB>BBDA@:@@CD@C?A?@D<@>92<+4>38>>CCCC3>CA
+@:114:A016U:1:1:10661:8787 1:N:0:1
+TACCTGCAAAATCTGGTGCCGTTTAAAGAGGGCGATGAGTACGAATCGAAAGATCTGGCAGAACTGAACCGCCGACTTTCTGCTACCGGCTGGGTTAACTCGGGGGGGGTGGCTCCACAATTTGATAAAGCGCGCGAAACGAAAGAATTTC
++
+CCCFFFFFHGHHHJJJHIJIJJIJJJJJIJJJJJJIJJJJIJJJJJJJJJJJJIJHHGHHFFFFEEEEEDDDDDDDDDDDDDDDDDCDDDBD@29?8>::@08&0<&&5&091(4<:9(+4>A(4>AAC3(5500&0&5&()(5.(+4(((
+@:114:A016U:1:1:22888:8790 1:N:0:1
+CGGCCATCAAGCACGGTGTCATTCATCACTACCAGCACGCCACGGTTGGCGGAGGCTTTATCAGCTGCGGTCACTACCGCGTTATACAGGTTGAATGGACCGTCTGCGCTCATAGACGTGGACGGACGCATTGCGCCGACCATCACCACCG
++
+CCCFFFFFHHHGGJJJGGIJJIJIIIJJJJIJJJIJJJJJJJJIIHIJIJJGFDDDDDDDDDEEDDDDDDDDDDDDDDDDB@>BBDDEDD@@CCCACDCCCBBBDDDD<>BBDCDDDBBBBBDD<9>>DD?<>@9095&5<54>>>9@?89
+@:114:A016U:1:1:22833:8799 1:N:0:1
+TGGTATTAATGCTGATTATTTAAATCTTTTAAAAAGAGCGTTAAATATCAAATTAACACTCCGGGAATACGCAGATCATCAAAAAGCAATGGACGCGCTTGCAGAAGGTGAAGTCGATATAGTGTTATCACATTTAGTTACTTCGCCGCCT
++
+B@@FFFFFHHFHDGIGHGIJJGIJICEHJJGIJJIDFEHGIIIJJIGIIIJJJJJJJJIGIJIJJIJHHHHFFFDEDEDDCDDCDDDDDDDDDDDD@D@BDDDDCDDDACCDC:>??&
+@:114:A016U:1:1:7280:8815 1:N:0:1
+ACTTCGTTCATTTCACCACCTCTGAGTAGAGGCGCTATTCTACCACTGCTGGAGAGGAAGAAAATCTAGTGCTGAAAAAATGATATCACCACGGGATAAACTGGTTAATGGCACAAATTCCGCTGAAAAATTGTACATCCTGATCGCACAT
++
+CCCFFFFFHHHHHJJJJJJJIIJJJJIJJJIJJJJJIJJJJJJJJJJJJJJIHIIJJJJHHHHHFFFFFDEEEEDEEDDDDDDDEEDECDDDD@:
+@:114:A016U:1:1:29919:8851 1:N:0:1
+AAGGTTTCTTTCGCGCCGAGCATATTACCGTGACGCACTGCGCGACGAATGATACGACGCAGTACATAACCACGGTTTTCATTGGACGGCATTACGCCATCCGCGATCAGGAACGCACAAGAACGAATGTGGTCAGCGCTTACGCGCAGCG
++
+1??DDDDDDDDD1EDEAAEDDDEEIIIIDDDDDIIIIIIEEADAD?????A>>>A?=<57798>>>AAABAAAA>99&5:>><;9><>AA>>?5;;>><<<99;.8<>4>>>>>>>;&+5&:230&&55)
+@:114:A016U:1:1:20227:8900 1:N:0:1
+TTAACGAAAGCGTCCCTTCATCGGGACGAGTGGAGATGAGCGTGGAGTCCATAATGATCCTGTTTGTACTGTCTTTTTATAATTAGTATGGAAAAGGTAACTCACATTGCCGTCGATCTATAAACGAACCTTATGAGACTACCTCGACAGA
++
+CCCFFFFFHHHGHJJJJIJJIJJJJJJJAGHGIHIGJIJJJJHIGHHFHFFFFFFFECEEEEDDDDDEEFCEDDDDDDBDDADEDC@DDACCCDDDBCCDDDCDCDDDCDDDDDDDBBDDEEDED<9@BBCCDDCDACDCDC?>>B>
+@:114:A016U:1:1:23534:8913 1:N:0:1
+CTCAGGGTAATTGGTGTAAAAAGCCACTTTAGGTCGAGTGAATGGTACCGGGTGGCGGACACATAACGAATAGTAATCGTACTCCAGCTGCTGAGCCTGAAGCTCAATTTCATGGTAGACCCCTTCTGCGGCCTCCATCCCCTGAACACGC
++
+???;BDB+2CBD:?:AE>DEFEIII@FEDED491?C?0?BAEDEI9?DDDDA;5<@'=@??59>5::BA<2128:>:(9A>9:>>8AAA<((+::AA>A>AA>A:(:A2+2<>A::;3&0&+83(++(8+(34(:3)&
+@:114:A016U:1:1:12930:8967 1:N:0:1
+AACTACTAAAAGAGCTTCATACATAATCAACCTGTTCCTTGCGGAGTTGCCGCGTACCAATGCTTCAACCAATAAAGCGGGAGTTTTTAATCTTTCCCACTGAAGCGGGTGTGAATACTAACCAAAGCGAACGACCTTCGCAAGCGCAATT
++
+8??DDB;?CDBDD?A?EDEE?A>@6.5;?BCC:5@>;/0&0<9::88A>CA:@4(++2?3:C::&)&0&)8(4((+(+4+8C(<+(020&)5&054(299><(85&05(
+@:114:A016U:1:1:17767:8969 1:N:0:1
+TCAGTTTCTGGCGTGGGGCGAAGAGAAAATCGGCATTGCTATTGATACGGGAAATAATAAAGTCAGCGAATCACTATTAGAGCAGCGCAATATTATTTTTTATCACGACAGTGAACATGAATCATTTATAGAAATGATACCGGGCAGCTAC
++
+CCCFFFFFGHHHHIIIJJJJIIIIIGJJJJIJFIGHHHFFFFFFFEEEDDDBDDDDDEEDDCDCEDDDDDDDDDDDDEEEDDDDDDDDDDDCDEDEEEEDDDDDEABBDDDCCDDDDDDDDEDDDDDDDECC(:@CCCED:@3<@BDBDDC
+@:114:A016U:1:1:22003:8995 1:N:0:1
+TTCACCTGTTCCCTTTGATATCTCGCTTCAGAATTGTATTCGTGTGCGTAATATTGAAACAAAACTCGTCACTGGAAAAGTAGGTACTCAAAACACACAACTGCTTGGTAATACGCTTACTGGAAGCACTGCCGCCAAAGGCGTTGGCGTA
++
+CCCFFFFFHHHHHJJJJIJIJJEHGIJJJJJJJJJIHIJJJJJJJJJJJJIHJJJJJIJJJJJJJJJIIJHHHHFFFFFE;AEEECCDDDDDDDBDDDDDDDDCDDDDDDDEEDDBDDBDDCCDDCADCDDDC@BDBDCC??B9@BBCB59
+@:114:A016U:1:1:4359:8995 1:N:0:1
+CCAGTTGGCGATGTGCTGGCTGGAGATACGCAGCGTTGCGCGGTCTTCCATCAACGCCACATTGTGAATATCCGGCACTTTTGAACAACCAATCCCCTGCTCCACCCAGCGCACCACGTACCCCAGAATCCCCTGCACGTTGTTATCCAGC
++
+CCCFFFFFHHHHHIJJJJJIJJJGGIJJJJJJIIJIJJJJJJJHGFFFFFFEEEDDDDDDDDDDDDDEEEEEDDDDDDDDDDDDDDDCBDDDDCCDDDDDDDDCCDD@ADBDBBDDBBBCDDDDDDACCD>C@??BCC?2CA
+@:114:A016U:1:1:10888:9001 1:N:0:1
+ATCCTGAGCACTTTCCCGGAGTGGGGTAGCTGGTTGAATGAAGAAATTGAACAAGAACAGGTCGCTCCTGGCACATTTGCGATGTGGTGGCTTGGCTGCACCGGGATCTGGTTGAAATCGGAAGGTGGCACCAACGATTGCGTTGATTTCT
++
+@@CFFDDDBDFBDGHHIG@EH?FEHEFGCGIIEBDF@FFHIEGCBHDGGGGIEHH>@GCEH=?BDDD<>CCDCCCDCCDCB@AB??CD?C>A@.>9.8>>C:@BCCDDDB>??>@(43
+@:114:A016U:1:1:5706:9066 1:N:0:1
+ATGCGTTGGGGAAAATGGATAAATGGTTGTACGCCGATATCACGCATTTTTCCCAGTTCTGGCATTACCTGAACGAACAGGATGAAACGCCCGGATTTGCCGATGACATGACCTGGGATTTTATCAGCAACGTCAATAGCATAACCCGCAA
++
+CCCFFFFFHHHHHJIJJJIIJJJIJJHHJJJJJIJJJIJJJJJJJJJHHHHHFFFFFEEEEEEEDDDDDDDDDDDDDDDDDDDDDDDDD@A3>B9
+@:114:A016U:1:1:26993:9067 1:N:0:1
+TGGCTGCCCAGCAGATCACAGGTTTTGAACCTTAGTGTCCTTCGTGGCTGTAGTAAAAGACCGGTAATCCCAGCTTAAGGCGTAGCGCCAGTAAGCGGGCAAAAAAGCCGAATACCAGGGTGGAGATAATTACCACATCATGGTTAGAAAC
++
+CCCFFFFFHHHHHJJJJJJJJJHHIJGIJIJJJIJGJJJJJJJJJJJJJJIIJIGIJJIIJJIGHDFDFEEEEEDDCCCCDDDDBDDDDDBDCCDCB?DDDCCCCD@>ACBC:ACCCACCBCC@(
+@:114:A016U:1:1:24070:9090 1:N:0:1
+AATTGTGGACATGGGAAACCCTGTCTGAAAGTTATAGTTAACCTACCCCATCCGTAGATGGGGGGATGTATGGGTACGTTGTAATTAGGGATTTAACGAATTAGCGCCAGGCGTCAAACCAGCCAAGCCCTTCTTCGGTGAGGCCACGAGG
++
+@@@FDFBDDFFHHGHIJGIIIIIGIIIIJGHCGGCGHIHIJJJGGJIJJJIJJJGFH7FGIIJID;>BACCDDDCCDB?CB?::@DDCDD>9?CCCCDDBBBC:CBDBBDBB<9>BBCCBABD7A<@?9BB
+@:114:A016U:1:1:15197:9103 1:N:0:1
+ACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACG
++
+CCCFFFFFHHHDDGGIJJJJJJJIIGIJJIJJJJIJJIJJJIIJJJJJJJJJJIJIHHHHFFDEDDDDBDDDDDDDDDDDDDDDDDDDDDDDDDDCDDDCDDDBDDBCDDDDDDEDCDDDDDDEECDDDDADDCDCCCDA>@C:@ACCA>2
+@:114:A016U:1:1:7868:9182 1:N:0:1
+TCCATGACGAGCGTCTTAGCACTGTGGAAGCCCGTTCCGGTCTGTTTGAACAGGGCGGCTATCGGGCGCTCAACAAAGGCAAAGTTGACTCTGCCTCTGCGGTTATTATTCTCGAAAACTATTTCGAGCAGGGATATTAAGGCGATTTAAA
++
+CCCFFFFFHHHHHIJJJJJJIJJIHIJJJIJJJIIIJJJJJJJJIJIJJJJJHHHHFDDBDDDDBDDDDDDDDDDDDDDDDDDD@CDDDDCDDDDDDCDDDD>>DDDEDDDDD?CDB&9ACC@DDDB<>BDDBBC:?955??CCC
+@:114:A016U:1:1:25771:9231 1:N:0:1
+TTCGTTACCTGCTCAAATCACTTTCCAGAAGCAATGTGCGACATCGTGAGCCATGTCAGTAAGTGGTTAAATTCGCAGTGGTAGCATGCCTCCCGTGATTTCTTACTAACCGGCGAAGAAATGGAACTGCTTTTATTGAGTAACTCGACGC
++
+@@@DDDDDHHBHDGEHDGHFCDHGGHCAHHEHIJGE:DHIBHDHGIFHHGAF>GI8=@F)>@@7=D@DE>;AC>D5;;>AC;AC>CB<=A08>4>@CCDDC@3>195&00.0&+3+++(:>@CC3ACC:>>+4:@4>((48(559
+@:114:A016U:1:1:27258:9263 1:N:0:1
+TATGTTTAAAGTGAGTGTCATATCACGGTGATAACGCGAAACAGCTCACAAAAATACACGCTTACGGGGTCTTTGCAACACAAAATAGCGCCGTTTGGCAAACTGAAGGGTTTATTGCTGAATGCCTGCTCCCCTCTCGTTCGAGGGGAAC
++
+@CCFFFFFHHHDFIJIIJJJJJIIGIJJFGIJJJIIJJJJJJJJJGIJJJJJJJJJJHHHFFFDEDDDD=BDDDDDEDDDDDDDDDDDDDDBD@DDDDDDDDDDDDDCDDDDDDDCDDDCDCDDDDDDD>CCDD?B?C<<@DDB<>BB@@0
+@:114:A016U:1:1:9093:9280 1:N:0:1
+CGTGCGATTATTTGCTCAATTAAGCTGGTATTTCCGTCGGGAATGGCGTCGCTATCTCGGGGCTGTCGCCTTGCTTGTCATTATCGCGATGCTGCAACTGGTTCCGCCAAAAGTGGTTGGTATTGTTGTCGATGGCGTGACAGAACAACAC
++
+@@@BDDADHAFHH@FBF:EEHCHEHGGHDGGIIBGBCFGHIAGEGHHAFFGEBBBEBC6=BBB?BA@ACCCAAC4@CCA@@>:>9@@9?A:@:>@@B5@<B(5<@A@A>CAA>>88:9?>95(<>A9:?3<1(2
+@:114:A016U:1:1:27685:9292 1:N:0:1
+AAACATGCCGTTAATCCAGCCGCCCAGCGAACCGCCAATCGCCACGCCAAGAAACTGGCTGGTGGAGTAAACACCCATCGCCGTACCTTTGTAACCTGCTGGCGACTCTTTACTGATAAGTGAGGGCAGGAGGGCTTCCATCAAATTAAAC
++
+CCCFFFDFGGHHHJJIJJJJJIGHIIHIIIGGJGJIHGEHEHFFFACDBDDBDDDDDCDCDD9C?:::>C@>>2<34:?..9554@C@@>34:::3(:>>A@380.(()&&+++:3(++:A:::@((
+@:114:A016U:1:1:20132:9298 1:N:0:1
+CGTGGATGAGCAGGGTAATGCCATTGACGTGGTCGACCCGATGCTGGCGGAGTTCCAGAAGATCAACGCGCAGTATCAGGGCGCAGACCGCGTGAAAGCGCTGCTGGGCCTGAGCGGTATTTTTGCCGATGATCTGCCGCAGAATGCCGAC
++
+@@CFFFDDA8DHHIIGHGEDBFFHEHFHGG?F8CFFGGGDA;FA=FIGEH<9=C>CDCC?CCCCDDC8?05;@(:@A:A99))8>@C:)3:5@>>+4>:>CCBBBDB9::>:950
+@:114:A016U:1:1:14853:9323 1:N:0:1
+AAGGCAGTCACTGCAACAGTTGGCGGTTGGGTACTGACGAACCGCTATCACTGGAAGGCGCTCAGGTGGCATCACCGGCGCTCACCGATTTACGGGCAAACCCAACGGCACGCGCCGCGTTATGGCAACAAATTTGCACATATGAACACGA
++
+CC@FFFFFHGHHHJJJJJJJJJIJIJHHJJJGIGIJIJJJJJJJHEFFFFFEEDDEDDDDDDDDDCCDDDDDDDDDDDDDD@BDDBBBBDDDDBDDDBDBABDDDD?BBBDDB><@BB@9>>?::<9@AABB@AACAACAC::CDEDCC59
+@:114:A016U:1:1:12207:9363 1:N:0:1
+GATGCAGTAATAAATCAGCGGCTGCCATTAATTCTGACACATCATTGCGACCGGAGAAGAAATGCACATTGCTCCGCACGCCGAGTTTTTCTGCCAGCGCTTCAAATTTTCGCGGCTTATCCTGACCAACAACAAATAAAAGCGTATTGTG
++
+CCCFFFFFGHHHGJJJJIIGIJIJIIIJIIJJJJJJJJJJIJJJJJJJJJJJJHHFFEEEEDDDDDDDDDDDCDDDDDDDDDDBD@BDDDDDDCDDCDDDDDDDDDDDDDD@BD>CCCCACCD7)82>ACC((:4:+4((+20(
+@:114:A016U:1:1:7663:9402 1:N:0:1
+GAATTGTTCAGGGCGATGTGCCGGAAGTGATGGCTGACTGTACGATTTACTCTCTCGATATCGGTTCTCTGTTAGCGGGCACAAAATATCGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAG
++
+@CCFFFFFHHHHHJJJJIIIJJJJGGIDHGIJIIGIHJJJJJJJFIIIJJJJIJIGGGHFFFFDCDBCDDCDDDDCDDDD@@CBB8ADEEBDB>BDDD>@C@@CDDDCCC18AC9:::
+@:114:A016U:1:1:6640:9403 1:N:0:1
+GAGCAGCCGGGTATCGTCGGCCTGTGGGCGATATTTTTTCAGCGCCATCGTCACCTGAGAGCGTAAATTCACAATCCGCTTCTGCTCCGGGAACCACGTTAAGAACTGCTCCTCCGCCTGGGTGCGGCTTTGCGCCACCTGTTCGCTGACG
++
+?@@D=DDDDHH@FHGEHGBHGIIG@DDHGIGFHGCGHHE;@CEEDECCC5=A=@CACC<(5<37???@>CCC3(4>3029B@C:::>>-0-5-<(2<@+2(24((+4@(:(:9?309?58A?<5?9@@<99(40<5<0889>4>(5<&)()
+@:114:A016U:1:1:15094:9494 1:N:0:1
+CCAGCGGGATATCCTGCGGTTCAAAACGCGCTTCTTCTTCAATCTCAGCGTTGATGGCAACCTGCTCGCCACCCAATACTTTTTTTTTCGGCTTATCACAAACTTTGCCGTTAACCAGCACTCGCTGGTCGAGGATCCATTCTTTTATTCG
++
+CCCFFFFFHHHHHJJIIJJIIIIJJIJJJIIHIJJIHIJIIHHHHHDHFFFDDCEDCCCDDDDDDDDDBDDDDDBDCDDDDEDDBDDDD>>@BBDDDDCACDDDDDCCDDDBDDDDDDDDD@CBDBBBCCD9@DBDDCC@@CDCDCCDCA?
+@:114:A016U:1:1:22788:9497 1:N:0:1
+TCGCTGAAATTGTCTTTGTGCCGACCGGAGAGAATTTAGCGATTATCGAACATACCGATGTCGATGAAAGCCTGAAAGGGCAAGGGATTGGTAAACAGCTGGTTGCGAAAGTCGTGGAAAAAATGCGTCGGGAAAAACGAAAAATTATCCC
++
+CCCFFFFFHHHHHJJJJIJFIIIJJJJJGHGHIIJJJJIJJJJJJJJJHHHFFFFFDCDDDDDDDEDDDDDDDDDDCCDDDDDDDBBDCDDDCCCCCDDDCDDCBBB>B9@::?<@BBCDDDB?>(525BBD?(4+(:(
+@:114:A016U:1:1:27941:9499 1:N:0:1
+GCTTTTACCCGCGACGAGCGTTTGAAATAGGGGACGATCAATGCGGTTGCCAACAAAGTAAACAGCAACTGGCGAACGGCTATGCCATCTTCCACGGGTGAAAAAGCGAATGCTTTTACTGCCAGCAATACGGAAATAAGCAGCCAGATAA
++
+CCCFFFFFHHHHHJIJGIJJJIJIJHIIIIJJJGIHFFFDDEEEDBDDDDDDDDDDDDCCDDDDDDDCBDDDDDDDD@BDBDACD@CCDDDDCCC?@75>BDDDDB<ACDDDDDDDBDCC?C?3>C?BDDB95?CDCCCBBAC?BBBDB3>BDCCA?<1>C3>CDD@>@9
+@:114:A016U:1:1:20499:9538 1:N:0:1
+CTGGCGGTTTCTGGCGTAGAAGCGGACGACGTTATCGGTACTCTGGCGCGCGAAGCCGAAAAAGCCGGGCGTCCGGTGCTGATCAGCACTGGCGATAAAGATATGGCGCAGCTGGTGACGCCAAATATTACGCTTATCAATACCATGACGA
++
+@@@FFFFAHHHHGIIFIIJJJIJJJGJIGGIHHFFFFDDDDDDDDDDDDDBBBDBDDDDDDDD((55>>CAA35:@CCCB<@CC>A?C9>>(:@@>ACA@<>BC?CC:4::>B@.8>(&28>(2<23:A:3(+&&&())00BBDDDBDDDCCCBCDDCCC?B@C3::AA
+@:114:A016U:1:1:8977:9633 1:N:0:1
+TTGCTGCACCAGAAGTGATTGACGCTGCACTTCGCGATGGCGCTGGTTTCCCGATGGGGCCGCTGGAATTAACCGATCTGATTGGTCAGGACGTCAATTTTGCTGTCACCTGTTCGGTGTTTAACGCTTTCTGGCAGGAGCGTCGTTTTTT
++
+CCCFFFFDHHDHFIGEFGJJJJJIJIJJIJJJJJIJDHIJJJJIJICHIGHHHEFFCDCDDDDDDD@CCCCCCDBBDBBD@ACDDCCCCCDDDBD@>D@CD@CCCACCCCDD@ACDBBDBBCCCCDDDBDCCD:>?BBCCDACDDDDDDDBB9>@DBD?
+@:114:A016U:1:1:28945:9898 1:N:0:1
+TCTTTTAACGACGCGCCCTCAGGGGAAACATCCACATCGTGGAAGCAATAAAATGGCACATGTAACTTGTGGAAAAACTCAAATGCGACATCTGCTTTACGCTTCGCCAACGCCAGTGCCTCACCAGGCTGCTGCCACGGACGATTAAACG
++
+<@@FFFFFHHHHHFGIFJEHIJJJJJIJJJJJJJGJIJGEHJFHHFGFFFFFED@CEEDCDDCCDCEDDDCBDCDDDBDDDDC<>A@BD@>CCCCDDDCDDD?@B4(
+@:114:A016U:1:1:23743:9915 1:N:0:1
+GACCGTTGCCAGTCGGCGTCATCCCCATCAAACGCGCGTTGAGTTTATCCTGCATGTAGCCTTTCAGAATGCCGTTCTCAATCAGCACGTTGTACTGGCCTGGCGTACCTTCGTCATCAATCGCCACCGAACCTCGGCGATCGACCATCGT
++
+CCCFFFFFHHHHHJJJJJJJIIJJJJIIJIIJJIJIJJHHFFFDDEEEEEEDDDDDDEEDDDDDDCDDDDDDDDDDDBDEDDDDDDDDDBDDDDEECDDDDDDDDDBDBDDD@>35?@<@>D:?@8
+@:114:A016U:1:1:16264:9965 1:N:0:1
+TGGAGCAGGCACGCAAAACAGAATTGAATGTTGCGGGGAAGAAAAATGCCGCTCTGTTTGAGTGGCGGCATTTGTGAAGGCGGAAAAGTAAGAATTGGCGTTACGAGTTATTCGATGTTCTGAATCTGCTCGCGCATCTGCTCAATCAACA
++
+CCCFFFFFGHHHGJJJJJJJJJJJJIGJJJHIIJJJJIEHFFFFFEEEDEDDDDDDCDDCCDCDBDBDB9>CDDCCCDDCDDB99@CC>CCCCC:A>>CB09<8?<52<::C>?:3>?C>CCCC>CA@>(42595055@(4(:(44@(>9+
+@:114:A016U:1:1:25716:10045 1:N:0:1
+GTTTGATGTTCTAACGTTGACCCGTAATCCGGGTTGCGGACAGTGTCTGGTGGGTAGTTTGACTGGGGCGGTCTCCTCCTAAAGAGTAACGGAGGAGCACGAAGGTTGGCTAATCCTGGTCGGACATCAGGAGGTTAGTGCAATGGCATAA
++
+CCCFFFFFHHHHHJJJJJIIJJJJJJJJJJJJJHHJJJJJJIJGHHIJJIAEHF?DDEEEDEDDDDDDDDDBDDDDDDCDCDCDDDCCDDDDDDDDDDDDDBDBDACDBDDDCDDDDD<>CBB>BB>ACCA@@>21:<(:A>C:3:8((::
+@:114:A016U:1:1:24242:10045 1:N:0:1
+ACGTGGGGAAGTTGCCCGATATTCATTACGCAAATTACCAGGCTGGTCAGTACCCAGATTATCGCCATCAACGGGACAAAGCCCTGCATGAGCCGGGCGACGCCATGAAGACCGCGAGTGATTGCCAGCAGAGTAAAGACAGCGAGAATAA
++
+CCCFFFFFHHHHHJJJJJJIJIJJJJJJJJJJJJJJJJJJJJJIJJGHIJIJJHHHHFFFFFFDDDDDDDDDDDBDDDDDDCDDDDDDDDDDDDBBBBDBDDDD@BCC:>CDCBDDDD9?C@ACDDDCD9@>:ACCBDB@A:
+@:114:A016U:1:1:21333:10091 1:N:0:1
+GGTTTGCTGGGGCCAAACGGCGCGGGTAAATCGACCACCTTTAAGATGATGTGCGGTTTGCTGGTGCCGACTTCCGGCCAGGCGCTGGTGCTGGGGATGGATCTGAAAGAGAGTTCCGGTAAAGCGCGCCAGCATCTCGGCTATATGGCGC
++
+CC@FFFFFHHHGHJIJJJJJIGIGIJBDDDBBBD<<@ACDDB?@BD9::CC>ACD??A@4>ACBD>BDDDDDDDBDBBD?ADCD???@?CCDDD<99
+@:114:A016U:1:1:27049:10124 1:N:0:1
+TTCCGGAGCAGTATCATTGCTGGCGTAATATTCCGCGCGCTGTTTTGCAATCCAGGCCAGTATTTCCGCCCGCTCTTCATCACTCACCTTACCCGCCCAGTGTAGCGCGGTATAACGCGTTGGTGGCATAGTTTCATACTGCATCACCCAT
++
+CCCFFFFFHHHHHIIJIJJJJJIJJHGIIIJJEIIJJJIJJIGHHHHHFFFFFFEEDDDDCDDEEFDBDDDDD@BBCDEEEDDCCDDDDCCCDCBDDD>>CCDDDD@09BDDCDDDB>>A8?BDDDCC>BCA:>A@>CCCDCACC<<9
+@:114:A016U:1:1:22335:10142 1:N:0:1
+CCACGACGGCAGAATGGTCCTTAATGATGGCGATAACGGAGAAGATATTCCCTATAGCTATCAACGGGAAGGTTTTGCTGACGGGCAACTGGTCGGTGAAGACGATCCTTGGCGTTTTGTCTGGATGACCTCACCTGATGGCAAATATCGC
++
+CCCFFFFFHHHHHJJIJHIJJJJJIIJJIJJJIGIJIJJJJIJJHHHHHHHFFFFFFEEEEEDEDDDDDDDD?CCDDCDDDDDDDDDDDDDDDDDD@BDCDDDDDDDDDDDDCD@<@DDBDDDDDCCCDDDDDDDD?CC@C?9>ACA3:?0
+@:114:A016U:1:1:4418:10160 1:N:0:1
+GGGGCTACATCCACAGGTGGTCATTGAGGCGAACTCAATTAGCGCGGTTCTGGAGCTGATTCGCCGCACTTCCCTTTCCACATTGTTACCAGCAGCGATTGCCACACAACATGACGGGCTTAAAGCTATTTCTCTTGCCCCGCCACTACTG
++
+CCCFFFFFHHHHHJJJJCGIHIJJJJIIJJJIJIJJJJJJJJJJJJHEFFFFFEEEEEDDDFDDDDDDBDDDDDDDDCDEDDDDDDCDDDDDDCBDDDBBBDDDDBDDDDDDDDDDBD>BDD@ACCDDC@CDDDDCC>CA@@BDDB:AC@>
+@:114:A016U:1:1:7108:10163 1:N:0:1
+CAGATTTACCCATTTCCTGCTGTTCGTAAGTCGAACGCGGAAGTTCTAAGCCAAATTGCTCACGGAATGTACGCTGTACGAAACCAGAACAATCGATACCTTTTTTAGTGCTGCCGCCCAGACGATAACGTACGCCTTTCCAGTCAGCATA
++
+CCCFFFFFHHHHHJJJJJJJJJJJJJIJIJIJJJJJJJJJIFIHHIJJIJJJHHHHHHFFFFFFDDDDDDFFDDDDDDDDDDDDDDDDDDDDDDDDDDDDCCDDDDDDDDDEDDDDDDDDDBDDDDBDCDCDD@DDBBB<>.@CCC@C>CDBD@BDDA
+@:114:A016U:1:1:21329:10193 1:N:0:1
+GGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCC
++
+CC@FFFFFHHHGHIJJIIJJJJIIIJJJHIIIJIJIHIJIJJIIJIJJJJJJJJJJJJFGIJIJJJJJJJJHIJIJGFHFFEFDEEEDABACEEDCDEEDDDDDDDDDBDDDDBB><@59AB@:>@CDD:9
+@:114:A016U:1:1:16672:10199 1:N:0:1
+GGAATGGCTATTGCAGGTTAACAAATTGGTGCGGGGCAACCCGCACATTTCAGGATGTTAGGAATGACCCATTTTTCACAGCAAGATAATTTTTCTGTCGCGGCGCGCGTGCTGGGCGCGCTGTTTTATTACGCTCCAGAGAGCGCGGAAG
++
+DD<3>ADDDBAAAA;ABA<::5:>AA??AADA>34>>A?>(3AA>?>?.9>;>505&+4>>?>0;9099AA>A>AAA29<9>AA(2<>(&59>&&0
+@:114:A016U:1:1:12276:10218 1:N:0:1
+AACGACCATTTATGACTACATCGCCAGGCGGCATCCACAACGCGCACAATGTGTGACTGATTTTATGAGTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGG
++
+CCCFFFFFHHHHHIJJJJJJJJGIJGIJJIJJIJJJJJJJJJHFFDDDDDDDEFDDDDDDDDEDEDDCCCDDDDBDBCDEDECDDDDDDDDDDDBD>BDDDDDDDDBDBDBDD?CCCDACDDDDDCCAC@CCACDBBBDDDC@CCCAB>&0
+@:114:A016U:1:1:11035:10223 1:N:0:1
+GCGTCTGTGCGGCGTTATTTTTGACTGGGATAAATATAATCTCGAGCTGTGCGAAGAAATTAGCAAAATGAACGAGAACCTGCCGTTGTACGCGTTCGCTAATACGTATTCCACTCTCGATGTAAGCCTGAATGACCTGCGTTTACAGATT
++
+CCCFFFFFHHHGHJJJJJJJJJJJJJJJJHJJJJJJJJJJJJJJJJJJHGHHFFFDAEDDDDDDDDDDDDDDDDBDDDDDDDDDBBBDDDEDBDDDDBDDDDDDDBCCDDDDDDDDDCDCDDBDDDDDDDDDDDDDCA?BB&8>?CCCC3:
+@:114:A016U:1:1:5645:10234 1:N:0:1
+ACAGTACGCCGCTAACGCAATTGCTACCTCTGGCATAACAAGTATATCGGGTAAGGGTTTCTGTTCCGCACACGCAGACGCAGAGTATCGTTAAGATGTCCATATTGTTGTTTTAGGCCCGCTAGTAATGCGCTACGGGTATTTAATATTG
++
+CCCFFFFFHHHHHJJJJJJJJJIJJJJJHIIJJJJJJJJJJJIHIJIIJIJGHHHHF@DEEEEEEEDDDDDDDDDDDDDDDDDDDCDDEBC?BDDDDDCDEDEDDEECCBCBCCDDDDDDBBDDDCDDDCDBDBBBDB@BCC@C@>
+@:114:A016U:1:1:21643:10269 1:N:0:1
+AGACCCAGTTTTACCCGTTTCGGTACGCTGCGTTCGCTCAGAATCGGCGCGGTGGAGATCAGCGCCAGCACGCGCAGTAACGGCCAGAAGTACAGGTTTAACCAGGATAGCCATTGTTCGCTTGTCACCTGCAACATAGTACGGCTACCCG
++
+CCCFFFFFHHHHHJJJJJJJJJJJJJJJJGIJGIJIIIJIJJJJJIIIHFDD;@DBDDDDDDDDCDDDDDDDDBDDDDDDDDDDDBDDDDCCCCCC+>@CDDDDDDDDDDDDDDDDEDDD@@BBADDDDCC?D>CBC@::@DBB7?@CAB?BCDB><.9@CA
+@:114:A016U:1:1:14413:10321 1:N:0:1
+ATGGTGATGTTTGACCTGCCATTCAACCCGGATCTACTGGAGCAGCGTATTGGTCGTCTGGATCGTATCGGCCAGGCGCACGATATTCAGATCCATGTGCCTTATCTGGAGAAAACCGCTCAGTCGGTGCTGGTGCGCTGGTATCACGAAG
++
+CCCFFFFFHHHHGJJIFIIJIJJJIIGJJJJGHIJJJJIJIJJJJJIHEIJHIHIJJHGHHHFFFFEDDDDDDDDDDDDDDDDDDBDDDDDDDDDDDDDDDDDDDDDEDDDDD?CDDDDB@>>@CDD@BDDDCCCCBB>@D3:@CCC7?9?
+@:114:A016U:1:1:18465:10327 1:N:0:1
+TCTAACGAGAAATTCCGCGGCGTACGTCTTTATGTCAGCGAAAACCAGCTGAAAATCACCGCCAACAACCCGGAACAGGAAGAAGCGGAAGAGATCCTCGACGTTACCTATAGCGGTGCGGAGATGGAAATCGGCTTCAACGTCAGTTATG
++
+CCCFFFFFGHHHHJJJJIJJIJIJJJJJJJJJJJJJIIHHHFFDEDDDDDDDDDDDDDDDDDDDDDDDDDDDBDDDDDBBDDDDDDBDBDDBDDDDDCDBDDDBBCBDCDDDECDD)5@DDBBDCDDC>CCDBDBDDDCCCDDDBDCDDDC
+@:114:A016U:1:1:28095:10338 1:N:0:1
+GATGATATCGGTCGTTGGGCCGTTAAGCTGAACATTACGCGGCTGCTGCTGTTCCTGGTAATGACTGATAGCGCGCTCAACCAGCGCCACTGCTTCGTCGATATCAAACGGTTTGGGCAGATAATCAAACGCCCCTTGTTGATAGGCGCTG
++
+CCCFFFFFHHHHHJIJJJIJJJIJJJJJJJIJJJJJIIGGGHHFDFEDEEEDDDDDDCCCCDDDDDDCDDCDDDDDBBDDDDDBBDDD@BA>CDCC@?BB9?ABCDC@C9?<@BBD@BABDDDDDDD>(4&0.323:@C3>5053
+@:114:A016U:1:1:8743:10398 1:N:0:1
+TGGCCGCCAAAGTTCAGTGCATCGTCTTTGTTCAGCCCCAGGTTACGGCGTTTACGCTCACTACCTGCCCGTGCCGCAGCGGCGGAATCGTGGAATGGCGGGGTACACAAACACGCGGCCGATTTCTCTGTTTTTTTGGTGTTTGCGTTTA
++
+B@CFFFFFHHHHHIJJIIJJJIJJIJIJJJIJJJJJJJJJJJIJFHBFHIEFFFFDDCBDDDDDDDDDDDDBDDDDDDDDDDDDDBBDDBDDDDDDCCBBB0&0.+44(8(((())&)0&&&)+(((4((+(4<(05&&+((++(()&+58
+@:114:A016U:1:1:15687:10414 1:N:0:1
+AGCATTTTTAGTCATGCTCTAATAAAATTAGGGATTGGCATATCCAATAACCCAGACAAAACCGATGTATAAGTCAACATATCCTGAATCAGACATACAATATCGCAATGAAAATCAATAATATTTTAAGGAATATCTTCATGAAATCAAA
++
+@??DDA=DDHBA:AEA@?@DHHIGEEHAGBH>?A9CE9CGHAHHHCHIBG@?B;D@GFH@F9FG8;B?>>?>CDEDBEECA>@>CCACCCCC<>@CACCCCCCD@?B39AC:>ACC>@C((3+4>>CD(:(++@A@CC@344>3>CCCC:@
+@:114:A016U:1:1:11062:10467 1:N:0:1
+GCAGAATATGGCAAAGCGCTTGCCAGCTATGCAGGTGAGCGTGAAACTCCTTTTCATCCGGTCTGCCGTTAGATTCACACAGCCACAGCGTGCGATCGATAGTCGCGCATAACGTCGTTTGTAATATCTGCGCGCCTTGCATACCGCTAAT
++
+CCCFFFFFHHHHHJJJJJJJJJIJIJJIJJJJJJIFHIJJJJJJIJJJJJJJJJJJJJJIIHHHFFFFDDDDDCEDDDDDDDDDDDDDDDBDDDDDDDDBDDCDDBDDBB?CDDCDDB>?BBEEEDEDEDDC>C>CC3AA?CB<<C>CB@B(2>CCCECCCACB<8??AC?B?0<98(:?>0)05
+@:114:A016U:1:1:20571:10559 1:N:0:1
+GCTGCCGAAGCAGATTGAACGTAATGCACTGAAAAAAGGTGAGCTGACCGTCGACGATAGCGCCATTATCGGCATTATTCGTTACTACACCCGTGAGGCGGGCGGGCGTGGTGTGGGGCGTGGAATCTCCAAACTGTGTCGCAAAACGGGT
++
+CCCFFFFFHHHHHJJIJJJIJHIJJJIIIIJIJJJJIIJFHIJJJJJJJJJJHHFFDDDDDDDDDDDDEDDDDDDDDEEEDDDBDDEDDDC>B>@BDBDDD<05&)&059B8(+80&.&009&2:444>@33++((++++&(&&)+(0290
+@:114:A016U:1:1:13177:10565 1:N:0:1
+GGTGCTGGCTGGCTGTATTGCGGAAGGGACGACGGTGGTTGATCGTATTTATCACATCGATCGTGGCTACGAACGCATTGAAGACAAACTGCGCGCTTTAGGTGCAAATATTGAGCGTGTGAAAGGCGAATAATCGTCTGAGAGCTGGCAG
++
+CC@FFDFFHHHHGJIDHIJJEGIGIIJJIGGIIIJ5CEBGHHFFEBDFFDEDCDEEDCBCBDDDDDDDCCBDDDDDDDDC@CCDDCCCBDDCDDDBDDBCDDCCACCCCCDEEDDCB>>>BDC@ABBDDBDDDE<@CCCDDDDDDDDDDDC:
+@:114:A016U:1:1:15254:10731 1:N:0:1
+CTTCCGGGTGAACGGTAGACGCGTCCAGCGGGTTATCACCGTGGTTAATGCGCAAGAAGCCCGCGCACTGCTCGAAGGCTTTCGGCCCCAGACGGCTCACTTTTAACAGTTGCTGACGGTTCTGGAACTGACCGTTCTCATCGCGCCAGGC
++
+CCCFFFFFFHHHHJJHIJJIJIJJJJJJJJJJEHHFHHGHFEECACEEEDDDDDDDDDDDDDBDBDDDBDDCDDDDDDDDDDDDDDDDDBDDDBDB>BDDCCCCCCDCD>CCDDDCCB<5?BD:@::ACDCDDD08A?CDCDDB9>@59?8
+@:114:A016U:1:1:8387:10765 1:N:0:1
+GCTGCCCTTCGCGCTCTTCCCGCGTCGCCCATTTGGTCACCGGTGCACCGGAAGATTGCGCCAGCGAAGCCAGCCGCCAGAGTTTGCTTTGCGCCTCGAACATCCCCAACGCGCGGGTCATCTGCAAAATCAACCGATGCAGACGTACAAT
++
+@@CDBDFBHHHGHI8DCFGIG@FC;@BCBBBBCCCCBBBBBB@@BB;AC1BB5A@4:@>CACACC99<>B<>CB@@CCBCCBB8>&055&)5<3>ACCCCCC3:ACC(5059