diff --git a/README.md b/README.md index 1d6b4ab..c99a503 100644 --- a/README.md +++ b/README.md @@ -3,6 +3,7 @@ # `decombinator` [![License](https://img.shields.io/badge/license-MIT-blue)](https://raw.githubusercontent.com/innate2adaptive/Decombinator/master/LICENCE) + @@ -29,6 +30,13 @@ Running `decombinator` is easy. All that is required is HTS read data and a few decombinator pipeline -in XXXX.fq -c b -br R2 -bl 42 -ol M13 ``` +## System Requirements + +`decombinator` runs on Linux (Ubuntu and CentOS tested). +Windows and MacOS are not officially supported but are expected to function correctly. +For Python and package requirements see `./pyproject.toml`. +Currently Python versions: 3.9, 3.10, 3.11, and 3.12 are supported. + ## Installation To install `decombinator` and all required packages, simply create a fresh virtual environment, activate it, and run: @@ -37,19 +45,30 @@ To install `decombinator` and all required packages, simply create a fresh virtu pip install decombinator ``` -`decombinator` also requires a number of additional files, which contain information regarding the V and J gene sequences, their tags, and the locations and motifs which define their CDR3 regions. By default `decombinator` downloads these files from [the git repository where they are maintained](https://github.com/innate2adaptive/Decombinator-Tags-FASTAs), which requires a working internet connection. In order to run `decombinator` offline, these files must be downloaded to a local location, and either stored within the directory where you wish to run `decombinator` or specified using the appropriate command line flag (`-tfdir`). These additional files can also be downloaded via `git clone`: +This process should take around 2 minutes, depending on internet connection. `decombinator` also requires a number of additional files, which contain information regarding the V and J gene sequences, their tags, and the locations and motifs which define their CDR3 regions. By default `decombinator` downloads these files from [the git repository where they are maintained](https://github.com/innate2adaptive/Decombinator-Tags-FASTAs), which requires a working internet connection. In order to run `decombinator` offline, these files must be downloaded to a local location, and either stored within the directory where you wish to run `decombinator` or specified using the appropriate command line flag (`-tfdir`). These additional files can also be downloaded via `git clone`: ```bash git clone https://github.com/innate2adaptive/Decombinator-Tags-FASTAs.git ``` -The current version of `decombinator` has tag sets for the analysis of alpha/beta and gamma/delta TCR repertoires from both mouse and man. In addition to the original tag set developed for the 2013 Thomas *et al* Bioinformatics paper, an 'extended' tag set has been developed for human alpha/beta repertoires, which covers 'non-functional' V and J genes (ORFs and pseudogenes, according to IMGT nomenclature) as well as just the functional. +The current version of `decombinator` has tag sets for the analysis of alpha/beta and gamma/delta TCR repertoires from both mouse and man. In addition to the original tag set developed for the 2013 Thomas _et al_ Bioinformatics paper, an 'extended' tag set has been developed for human alpha/beta repertoires, which covers 'non-functional' V and J genes (ORFs and pseudogenes, according to IMGT nomenclature) as well as just the functional. ## Basic Usage +`decombinator` can be run in 4 different modes depending on your use case. +Each mode is outlined below. + +### Example + +For a ready to run example of `decombinator`, see the `examples/` directory. +Within this directory there is a `.sh` script which runs `decombinator pipeline` on a small set of paired-end reads. +This script aligns both alpha and beta chains from this file, and provides expected outputs for comparison. +It should take under 30s on a laptop computer. +Please use this to check your installation is configured correctly. + ### `pipeline` -`decombinator` can be run in 4 different modes depending on your use case. Firstly, `pipeline` mode can be used via: +Firstly, `pipeline` mode can be used via: ```shell decombinator pipeline -in XXXX.fq -c b -br R2 -bl 42 -ol M13 @@ -92,7 +111,6 @@ decombinator pipeline -in XXXX.fq -c b -br R2 -bl 42 -ol M13 | `-uh`, `--UMIhistogram` | Creates histogram of average UMI cluster sizes | | `-npf`, `--nonproductivefilter` | Filter out non-productive reads from the output | - In this mode, all three main components of the pipeline are applied to the data sequentially: `decombine`, `collapse`, and `translate`. @@ -109,27 +127,27 @@ decombinator decombine -in XXXX.fq -c b -br R2 -bl 42
decombinator decombine CLI arguments -| Option | Description | -|--------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------------------------------| -| `-h`, `--help` | Show this help message and exit | -| `-s`, `--suppresssummary` | Suppress the production of summary data log/file | -| `-dz`, `--dontgzip` | Stop the output FASTQ files automatically being compressed with gzip | -| `-dc`, `--dontcount` | Stop/Block printing the running count | -| `-op OUTPATH`, `--outpath OUTPATH` | Path to output directory, writes to directory script was called in by default | -| `-c CHAIN`, `--chain CHAIN` | TCR chain (a/b/g/d) | -| `-pf PREFIX`, `--prefix PREFIX` | Specify the prefix of the output DCR file. Default = "dcr_" | -| `-in INFILE`, `--infile INFILE` | Correctly demultiplexed/processed FASTQ file containing TCR reads | -| `-br BC_READ`, `--bc_read BC_READ` | Which read has bar code (R1,R2). If used, ensure read selected is present in the same directory as the file specified by -in. | -| `-dk`, `--dontcheck` | Skip the FASTQ check | -| `-ex EXTENSION`, `--extension EXTENSION` | Specify the file extension of the output DCR file. Default = "n12" | -| `-or ORIENTATION`, `--orientation ORIENTATION` | Specify the orientation to search in (forward/reverse/both). Default = reverse | -| `-tg TAGS`, `--tags TAGS` | Specify which Decombinator tag set to use (extended or original). Default = extended | -| `-sp SPECIES`, `--species SPECIES` | Specify which species TCR repertoire the data consists of (human or mouse). Default = human | -| `-N`, `--allowNs` | Whether to allow VJ rearrangements containing ambiguous base calls ('N'). Default = False | -| `-ln LENTHRESHOLD`, `--lenthreshold LENTHRESHOLD` | Acceptable threshold for inter-tag (V to J) sequence length. Default = 130 | -| `-tfdir TAGFASTADIR`, `--tagfastadir TAGFASTADIR` | Path to folder containing TCR FASTA and Decombinator tag files, for offline analysis. Default = "Decombinator-Tags-FASTAs". | -| `-nbc`, `--nobarcoding` | Option to run Decombinator without barcoding, i.e. so as to run on data produced by any protocol. | -| `-bl BCLENGTH`, `--bclength BCLENGTH` | Length of barcode sequence, if applicable. Default is set to 42 bp. | +| Option | Description | +| ------------------------------------------------- | ----------------------------------------------------------------------------------------------------------------------------- | +| `-h`, `--help` | Show this help message and exit | +| `-s`, `--suppresssummary` | Suppress the production of summary data log/file | +| `-dz`, `--dontgzip` | Stop the output FASTQ files automatically being compressed with gzip | +| `-dc`, `--dontcount` | Stop/Block printing the running count | +| `-op OUTPATH`, `--outpath OUTPATH` | Path to output directory, writes to directory script was called in by default | +| `-c CHAIN`, `--chain CHAIN` | TCR chain (a/b/g/d) | +| `-pf PREFIX`, `--prefix PREFIX` | Specify the prefix of the output DCR file. Default = "dcr\_" | +| `-in INFILE`, `--infile INFILE` | Correctly demultiplexed/processed FASTQ file containing TCR reads | +| `-br BC_READ`, `--bc_read BC_READ` | Which read has bar code (R1,R2). If used, ensure read selected is present in the same directory as the file specified by -in. | +| `-dk`, `--dontcheck` | Skip the FASTQ check | +| `-ex EXTENSION`, `--extension EXTENSION` | Specify the file extension of the output DCR file. Default = "n12" | +| `-or ORIENTATION`, `--orientation ORIENTATION` | Specify the orientation to search in (forward/reverse/both). Default = reverse | +| `-tg TAGS`, `--tags TAGS` | Specify which Decombinator tag set to use (extended or original). Default = extended | +| `-sp SPECIES`, `--species SPECIES` | Specify which species TCR repertoire the data consists of (human or mouse). Default = human | +| `-N`, `--allowNs` | Whether to allow VJ rearrangements containing ambiguous base calls ('N'). Default = False | +| `-ln LENTHRESHOLD`, `--lenthreshold LENTHRESHOLD` | Acceptable threshold for inter-tag (V to J) sequence length. Default = 130 | +| `-tfdir TAGFASTADIR`, `--tagfastadir TAGFASTADIR` | Path to folder containing TCR FASTA and Decombinator tag files, for offline analysis. Default = "Decombinator-Tags-FASTAs". | +| `-nbc`, `--nobarcoding` | Option to run Decombinator without barcoding, i.e. so as to run on data produced by any protocol. | +| `-bl BCLENGTH`, `--bclength BCLENGTH` | Length of barcode sequence, if applicable. Default is set to 42 bp. |
@@ -139,11 +157,12 @@ decombinator decombine -in XXXX.fq -c b -br R2 -bl 42 This function performs the key computation of the pipeline, as it searches through demultiplexed reads for rearranged TCR sequences. It looks for short 'tag' sequences (using Aho-Corasick string matching): the presence of a tag uniquely identifies a particular V or J gene. If it finds both a V and a J tag (and the read passes various filters), it assigns the read as recombined, and outputs a five-part Decombinator index (or 'DCR'), which uniquely represents a given TCR rearrangement, plus some addtional UMI related information. All DCR-containing output files are comma-delimited, with the fields of that five-part classifier containing, in order: -* The V index (which V gene was used) -* The J index -* Number of V deletions (relative to germline) -* Number of J deletions -* Insert sequence (the nucleotide sequence between the end of deleted V and J) + +- The V index (which V gene was used) +- The J index +- Number of V deletions (relative to germline) +- Number of J deletions +- Insert sequence (the nucleotide sequence between the end of deleted V and J) The V and J indices are arbitrary numbers based on the order of the tag sequences in the relevant tag file (using Python indexing, which starts at 0 rather than 1). Also, note that the number of V and J deletions just represents how many bases have been removed from the end of that particular germline gene (as given in the germline FASTA files in the additional file repo); it is entirely possible that more bases were deleted, and just that the same bases have been re-added. Additionally, there are low frequencies of (predominantly alpha chain) recombinations where there is no detectable insertion, and where the nucleotides at the junction between the germline V and J genes could have derived from either. In such circumstances, the nucleotides will arbitrarily be deemed to have derived from the V gene, and thus count towards deletions from the J, however it is impossible to know which gene originally contributed these residues. @@ -157,18 +176,19 @@ Various additional fields may follow the five-part classifier, but the DCR will Which corresponds to a rearrangement between TRAV1-2 (V index **1**, with **9** nucleotides deleted) and TRAJ33 (J index **22**, with **0** deletions), with an insert sequence (i.e. non-templated additions to the V and/or the J gene) of '**CTCTA**'. For beta chains, the insert sequence will contain any residual TRBD nucleotides, although as these genes are very short, homologous, and typically highly 'nibbled', they are often impossible to differentiate. Produces a list of lists with the following entries: - 1. V index - 2. J index - 3. Number of V deletions - 4. Number of J deletions - 5. insert - 6. ID - 7. TCR sequence - 8. TCR quality - 9. barcode sequence - 10. barcode quality - -** NB The TCR sequence given here is the 'inter-tag' region, i.e. the sequence between the start of the found V tag and the end of the found J tag. + +1. V index +2. J index +3. Number of V deletions +4. Number of J deletions +5. insert +6. ID +7. TCR sequence +8. TCR quality +9. barcode sequence +10. barcode quality + +\*\* NB The TCR sequence given here is the 'inter-tag' region, i.e. the sequence between the start of the found V tag and the end of the found J tag. @@ -181,43 +201,42 @@ decombinator collapse -in XXXX.n12 -c b -ol M13
decombinator collapse CLI arguments -| Option | Description | -|--------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------------------------------| -| `-h`, `--help` | Show this help message and exit | -| `-s`, `--suppresssummary` | Suppress the production of summary data log/file | -| `-dz`, `--dontgzip` | Stop the output FASTQ files automatically being compressed with gzip | -| `-dc`, `--dontcount` | Stop/Block printing the running count | -| `-op OUTPATH`, `--outpath OUTPATH` | Path to output directory, writes to directory script was called in by default | -| `-c CHAIN`, `--chain CHAIN` | TCR chain (a/b/g/d) | -| `-pf PREFIX`, `--prefix PREFIX` | Specify the prefix of the output DCR file. Default = "dcr_" | -| `-in INFILE`, `--infile INFILE` | File containing raw verbose Decombinator output, i.e. 5 part classifier plus barcode and inter-tag sequence and quality strings | -| `-mq MINBCQ`, `--minbcQ MINBCQ` | Minimum quality score that barcode nucleotides should be to for that rearrangement to be retained. Default = 20 | -| `-bm BCQBELOWMIN`, `--bcQbelowmin BCQBELOWMIN` | Number of nucleotides per barcode whose quality score are allowed to be below -mq and still be retained. Default = 1 | -| `-aq AVGQTHRESHOLD`, `--avgQthreshold AVGQTHRESHOLD` | Average quality threshold that barcode sequences must remain above for rearrangements to be retained. Default = 30 | -| `-lv PERCENTLEVDIST`, `--percentlevdist PERCENTLEVDIST` | Percentage Levenshtein distance that is allowed to estimate whether two sequences within a barcode are derived from the same originator molecule. Default = 10 | -| `-bc BCTHRESHOLD`, `--bcthreshold BCTHRESHOLD` | Number of sequence edits that are allowed to consider two barcodes to be derived from same originator during clustering. Default = 2 | -| `-ex EXTENSION`, `--extension EXTENSION` | Specify the file extension of the output DCR file. Default = 'freq' | -| `-N`, `--allowNs` | Used to allow VJ rearrangements containing ambiguous base calls ('N') | -| `-ln LENTHRESHOLD`, `--lenthreshold LENTHRESHOLD` | Acceptable threshold for inter-tag (V to J) sequence length | -| `-di`, `--dontcheckinput` | Override the input file sanity check | -| `-bd`, `--barcodeduplication` | Optionally output a file containing the final list of clustered barcodes, and their frequencies | -| `-pb`, `--positionalbarcodes` | Instead of inferring random barcode sequences from their context relative to spacer sequences, just take the sequence at the default positions. Useful to salvage runs when R2 quality is terrible. | -| `-ol OLIGO`, `--oligo OLIGO` | Choose experimental oligo for correct identification of spacers ["M13", "I8", "I8_single", "NEBIO", "TAKARA"] (default: M13) | -| `-wc`, `--writeclusters` | Write cluster data to separate cluster files | -| `-uh`, `--UMIhistogram` | Creates histogram of average UMI cluster sizes | - +| Option | Description | +| ------------------------------------------------------- | --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- | +| `-h`, `--help` | Show this help message and exit | +| `-s`, `--suppresssummary` | Suppress the production of summary data log/file | +| `-dz`, `--dontgzip` | Stop the output FASTQ files automatically being compressed with gzip | +| `-dc`, `--dontcount` | Stop/Block printing the running count | +| `-op OUTPATH`, `--outpath OUTPATH` | Path to output directory, writes to directory script was called in by default | +| `-c CHAIN`, `--chain CHAIN` | TCR chain (a/b/g/d) | +| `-pf PREFIX`, `--prefix PREFIX` | Specify the prefix of the output DCR file. Default = "dcr\_" | +| `-in INFILE`, `--infile INFILE` | File containing raw verbose Decombinator output, i.e. 5 part classifier plus barcode and inter-tag sequence and quality strings | +| `-mq MINBCQ`, `--minbcQ MINBCQ` | Minimum quality score that barcode nucleotides should be to for that rearrangement to be retained. Default = 20 | +| `-bm BCQBELOWMIN`, `--bcQbelowmin BCQBELOWMIN` | Number of nucleotides per barcode whose quality score are allowed to be below -mq and still be retained. Default = 1 | +| `-aq AVGQTHRESHOLD`, `--avgQthreshold AVGQTHRESHOLD` | Average quality threshold that barcode sequences must remain above for rearrangements to be retained. Default = 30 | +| `-lv PERCENTLEVDIST`, `--percentlevdist PERCENTLEVDIST` | Percentage Levenshtein distance that is allowed to estimate whether two sequences within a barcode are derived from the same originator molecule. Default = 10 | +| `-bc BCTHRESHOLD`, `--bcthreshold BCTHRESHOLD` | Number of sequence edits that are allowed to consider two barcodes to be derived from same originator during clustering. Default = 2 | +| `-ex EXTENSION`, `--extension EXTENSION` | Specify the file extension of the output DCR file. Default = 'freq' | +| `-N`, `--allowNs` | Used to allow VJ rearrangements containing ambiguous base calls ('N') | +| `-ln LENTHRESHOLD`, `--lenthreshold LENTHRESHOLD` | Acceptable threshold for inter-tag (V to J) sequence length | +| `-di`, `--dontcheckinput` | Override the input file sanity check | +| `-bd`, `--barcodeduplication` | Optionally output a file containing the final list of clustered barcodes, and their frequencies | +| `-pb`, `--positionalbarcodes` | Instead of inferring random barcode sequences from their context relative to spacer sequences, just take the sequence at the default positions. Useful to salvage runs when R2 quality is terrible. | +| `-ol OLIGO`, `--oligo OLIGO` | Choose experimental oligo for correct identification of spacers ["M13", "I8", "I8_single", "NEBIO", "TAKARA"] (default: M13) | +| `-wc`, `--writeclusters` | Write cluster data to separate cluster files | +| `-uh`, `--UMIhistogram` | Creates histogram of average UMI cluster sizes |
collapse details -Takes the output files of Decombinator (run using the barcoding option) and performs collapsing and error correction. This version is a modified version of KB's script collapsinator_20141126.py (That was itself an improved version of the CollapseTCRs.py script used in the Heather et al HIV TCR paper (DOI: 10.3389/fimmu.2015.00644))* -** Version 4.0.2 includes improved clustering routines measuring the similarity in both barcode and TCR sequence of TCR repertoire data - -** NOTE - from version 4.2 this optionally looks for barcode 6NI86N at the beginning of the read; instead of M13_6N_I8_6N_I8 - (i.e. only one spacer). - This makes it compatible with the multiplex protocol in which the barcode is incorporated in the RT step and is found at the beginning of R1. ** From version V4.2 there is a required additional command line parameter -ol (see below for allowed inputs) +Takes the output files of Decombinator (run using the barcoding option) and performs collapsing and error correction. This version is a modified version of KB's script collapsinator_20141126.py (That was itself an improved version of the CollapseTCRs.py script used in the Heather et al HIV TCR paper (DOI: 10.3389/fimmu.2015.00644))\* +\*\* Version 4.0.2 includes improved clustering routines measuring the similarity in both barcode and TCR sequence of TCR repertoire data + +** NOTE - from version 4.2 this optionally looks for barcode 6NI86N at the beginning of the read; instead of M13_6N_I8_6N_I8 +(i.e. only one spacer). +This makes it compatible with the multiplex protocol in which the barcode is incorporated in the RT step and is found at the beginning of R1. ** From version V4.2 there is a required additional command line parameter -ol (see below for allowed inputs) The barcode sequence is contained in one of the additional fields output by `decombine`, which contains the first 42 bases of R2 (if `-br R2 -bl 42` specified). As Illumina sequencing is particularly error-prone in the reverse read, and that reads can be phased (i.e. they do not always begin with the next nucleotide that follows the sequencing primer) our protocol uses known spacer sequences to border the random barcode bases, so that we can identify the actual random bases. The hexameric barcode locations (N6) are determined in reference to the two spacer sequences like so: ``` @@ -228,20 +247,23 @@ GTCGTGATNNNNNNGTCGTGATNNNNNNnn The collapsing script uses the spacer sequences to identify the exact position of the barcode sequences. `collapse` performs the following procedures: -* Scrolls through each line of the input object containing DCR, barcode and sequence data. -* Removes TCR reads with forbidden errors, e.g. ambiguous base calls (with user input parameters provided to modify strictness). -* Groups input reads by barcode. Reads with identical barcodes and equivalent inter-tag sequences are grouped together. Equivalence is defined as the Levenshtein distance between two sequences being lower than a given threshold, weighted by the lengths of the compared sequences. Reads with identical barcodes but non-equivalent sequences are grouped separately. -* Each group is assigned the most common inter-tag sequence/DCR combination as the 'true' TCR, as errors are likely to occur during later PCR cycles, and thus will most often be minority variants (see [Bolotin *et al.*, 2012](http://dx.doi.org/10.1002/eji.201242517)). + +- Scrolls through each line of the input object containing DCR, barcode and sequence data. +- Removes TCR reads with forbidden errors, e.g. ambiguous base calls (with user input parameters provided to modify strictness). +- Groups input reads by barcode. Reads with identical barcodes and equivalent inter-tag sequences are grouped together. Equivalence is defined as the Levenshtein distance between two sequences being lower than a given threshold, weighted by the lengths of the compared sequences. Reads with identical barcodes but non-equivalent sequences are grouped separately. +- Each group is assigned the most common inter-tag sequence/DCR combination as the 'true' TCR, as errors are likely to occur during later PCR cycles, and thus will most often be minority variants (see [Bolotin _et al._, 2012](http://dx.doi.org/10.1002/eji.201242517)). After this initial grouping, the script estimates the true cDNA frequency. UMIs that are both similar and are associated with a similar TCR are likely to be amplified from the same original DNA molecule and to differ only due to PCR or sequencing error. Consequently, groups with similar barcodes and sequences are then clustered via the following procedure: -* The barcode of each group is compared to the barcode of every other group. -* The expected distribution of distances between UMIs can be modelled as a binomial distribution. Experimentation with simulated datasets found the best threshold for allowing two barcodes to be considered equivalent is when they have a Levenshtein distance of less than 3; a value of 2 is set by default. This can be modified through the user input parameter ```-bc```. -* Groups with barcodes that meet this threshold criteria have their inter-tag sequences compared. Those with equivalent sequences are clustered together. Sequence equivalence is here taken to mean that the two sequences have a Levenshtein distance less than or equal to 10% of the length of the shorter of the two sequences. This percentage can be modified through the user input parameter ```-lv```. -* Upon this merging of groups, the most common inter-tag sequence of the cluster is reassessed and taken as the 'true' TCR. + +- The barcode of each group is compared to the barcode of every other group. +- The expected distribution of distances between UMIs can be modelled as a binomial distribution. Experimentation with simulated datasets found the best threshold for allowing two barcodes to be considered equivalent is when they have a Levenshtein distance of less than 3; a value of 2 is set by default. This can be modified through the user input parameter `-bc`. +- Groups with barcodes that meet this threshold criteria have their inter-tag sequences compared. Those with equivalent sequences are clustered together. Sequence equivalence is here taken to mean that the two sequences have a Levenshtein distance less than or equal to 10% of the length of the shorter of the two sequences. This percentage can be modified through the user input parameter `-lv`. +- Upon this merging of groups, the most common inter-tag sequence of the cluster is reassessed and taken as the 'true' TCR. Finally, the clusters are collapsed to give the abundance of each TCR in the biological sample. -* A TCR abundance count is calculated for each TCR by counting the number of clusters that have the same sequence but different barcodes (thus representing the same rearrangement originating from multiple input DNA molecules). -* An average UMI count is calculated for each TCR by summing the number of members in each cluster associated with the TCR sequence and dividing by the number of those clusters. This gives a measure that can be used to estimate the robustness of the data for that particular sequence. + +- A TCR abundance count is calculated for each TCR by counting the number of clusters that have the same sequence but different barcodes (thus representing the same rearrangement originating from multiple input DNA molecules). +- An average UMI count is calculated for each TCR by summing the number of members in each cluster associated with the TCR sequence and dividing by the number of those clusters. This gives a measure that can be used to estimate the robustness of the data for that particular sequence. Collapsinator outputs 7 fields: the 5-part DCR identifier, the corrected abundance of that TCR in the sample, and the average UMI count for that TCR @@ -256,22 +278,21 @@ decombinator translate -in XXXX.freq -c b -ol M13
decombinator translate CLI arguments -| Option | Description | -|--------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------------------------------| -| `-h`, `--help` | Show this help message and exit | -| `-s`, `--suppresssummary` | Suppress the production of summary data log/file | -| `-dz`, `--dontgzip` | Stop the output FASTQ files automatically being compressed with gzip | -| `-dc`, `--dontcount` | Stop/Block printing the running count | -| `-op OUTPATH`, `--outpath OUTPATH` | Path to output directory, writes to directory script was called in by default | -| `-c CHAIN`, `--chain CHAIN` | TCR chain (a/b/g/d) | -| `-pf PREFIX`, `--prefix PREFIX` | Specify the prefix of the output DCR file. Default = "dcr_" | -| `-in INFILE`, `--infile INFILE` | File containing 5 part classifier plus barcode and inter-tag sequence and quality strings | -| `-sp SPECIES`, `--species SPECIES` | Specify which species TCR repertoire the data consists of (human or mouse). Default = human | -| `-tg TAGS`, `--tags TAGS` | Specify which Decombinator tag set to use (extended or original). Default = extended | -| `-npf`, `--nonproductivefilter` | Filter out non-productive reads from the output | -| `-tfdir TAGFASTADIR`, `--tagfastadir TAGFASTADIR` | Path to folder containing TCR FASTA and Decombinator tag files, for offline analysis. Default = 'Decombinator-Tags-FASTAs' | -| `-nbc`, `--nobarcoding` | Option to run CD3translator without barcoding, i.e. so as to run on data produced by any protocol. | - +| Option | Description | +| ------------------------------------------------- | -------------------------------------------------------------------------------------------------------------------------- | +| `-h`, `--help` | Show this help message and exit | +| `-s`, `--suppresssummary` | Suppress the production of summary data log/file | +| `-dz`, `--dontgzip` | Stop the output FASTQ files automatically being compressed with gzip | +| `-dc`, `--dontcount` | Stop/Block printing the running count | +| `-op OUTPATH`, `--outpath OUTPATH` | Path to output directory, writes to directory script was called in by default | +| `-c CHAIN`, `--chain CHAIN` | TCR chain (a/b/g/d) | +| `-pf PREFIX`, `--prefix PREFIX` | Specify the prefix of the output DCR file. Default = "dcr\_" | +| `-in INFILE`, `--infile INFILE` | File containing 5 part classifier plus barcode and inter-tag sequence and quality strings | +| `-sp SPECIES`, `--species SPECIES` | Specify which species TCR repertoire the data consists of (human or mouse). Default = human | +| `-tg TAGS`, `--tags TAGS` | Specify which Decombinator tag set to use (extended or original). Default = extended | +| `-npf`, `--nonproductivefilter` | Filter out non-productive reads from the output | +| `-tfdir TAGFASTADIR`, `--tagfastadir TAGFASTADIR` | Path to folder containing TCR FASTA and Decombinator tag files, for offline analysis. Default = 'Decombinator-Tags-FASTAs' | +| `-nbc`, `--nobarcoding` | Option to run CD3translator without barcoding, i.e. so as to run on data produced by any protocol. |
@@ -280,44 +301,44 @@ decombinator translate -in XXXX.freq -c b -ol M13 This step outputs `.tsv` files in the form of: -| Field | Description | -|:---:|---| -| sequence_id | A unique identifier for a given rearrangement | -| v_call | V gene used (or multiple, if they cannot be distinguished, comma delimited) | -| d_call | Blank required field (mostly cannot be assigned for TCRb, and rarely useful even then) | -| j_call | J gene used | -| junction_aa | CDR3 junction amino acid sequence | -| duplicate_count | Rearrangement abundance, from Collapsinator | -| sequence | Inferred full-length variable domain nucleotide sequence | -| junction | CDR3 junction nucleotide sequence | -| decombinator_id | Five-field Decombinator identifier | -| rev_comp | Whether rearrangements are reverse complemented (T/F) - this is never the case post-Decombining | -| productive | Whether rearrangement is potentially productive (T/F) | -| sequence_aa | Inferred full-length variable domain amino acid sequence | -| cdr1_aa | Amino acid sequence of CDR1 of the used V gene | -| cdr2_aa | Amino acid sequence of CDR2 of the used V gene | -| vj_in_frame | Whether or not the rearrangement is in frame (T/F) | -| stop_codon | Whether or not the rearrangement contains a stop codon (T/F) | -| conserved_c | Whether or not the rearrangement contains a detectable conserved cysteine (T/F) | -| conserved_f | Whether or not the rearrangement contains a detectable conserved phenylalanine or equivalent (T/F) | -| legacy_v_call | What older versions of Decombinator (i.e. <= v3) referred to this V gene as: e.g. 'TRBV12-3,TRBV12-4' in v4 was previously referred to just as 'TRBV12-4' | -| legacy_j_call | What older versions of Decombinator (<= 3) referred to this J gene as | -| v_alleles | List of V gene alleles covered by this rearrangment's tag | -| j_alleles | List of J gene alleles covered by this rearrangment's tag | -| v_gene_functionality | [IMGT predicted functionality](http://www.imgt.org/IMGTScientificChart/SequenceDescription/IMGTfunctionality.html#P1-2) of V gene (or genes) used in this rearrangement (F/ORF/P, comma delimited) | -| j_gene_functionality | [IMGT predicted functionality](http://www.imgt.org/IMGTScientificChart/SequenceDescription/IMGTfunctionality.html#P1-2) of J gene (or genes) used in this rearrangement (F/ORF/P, comma delimited) | -| sequence_alignment | Format required field - left blank | -| germline_alignment | Format required field - left blank | -| v_cigar | Format required field - left blank | -| d_cigar | Format required field - left blank | -| j_cigar | Format required field - left blank | -| av_UMI_cluster_size | The average UMI count for this particular sequence | - -You can also use the 'nonproductivefilter' flag (`-npf`) to suppress the output of non-productive rearrangements. - -As the hypervariable region and the primary site of antigenic contact, the CDR3 is almost certainly going to be the region of most interest for many analyses. By convention, the [CDR3 junction is defined as](http://dx.doi.org/10.1016/S0145-305X(02)00039-3) running from the position of the second conserved cysteine encoded in the 3' of the V gene to the phenylalanine in the conserved 'FGXG' motif in the J gene. However, some genes use non-canonical residues/motifs, and the position of these motifs varies. - -In looking for CDR3s, we also find 'non-productive' reads, i.e. those that don't appear to be able to make productive, working TCRs. This is determined based on the presence of stop codons, being out of frame, or lacking appropriate CDR3 motifs. +| Field | Description | +| :------------------: | -------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- | +| sequence_id | A unique identifier for a given rearrangement | +| v_call | V gene used (or multiple, if they cannot be distinguished, comma delimited) | +| d_call | Blank required field (mostly cannot be assigned for TCRb, and rarely useful even then) | +| j_call | J gene used | +| junction_aa | CDR3 junction amino acid sequence | +| duplicate_count | Rearrangement abundance, from Collapsinator | +| sequence | Inferred full-length variable domain nucleotide sequence | +| junction | CDR3 junction nucleotide sequence | +| decombinator_id | Five-field Decombinator identifier | +| rev_comp | Whether rearrangements are reverse complemented (T/F) - this is never the case post-Decombining | +| productive | Whether rearrangement is potentially productive (T/F) | +| sequence_aa | Inferred full-length variable domain amino acid sequence | +| cdr1_aa | Amino acid sequence of CDR1 of the used V gene | +| cdr2_aa | Amino acid sequence of CDR2 of the used V gene | +| vj_in_frame | Whether or not the rearrangement is in frame (T/F) | +| stop_codon | Whether or not the rearrangement contains a stop codon (T/F) | +| conserved_c | Whether or not the rearrangement contains a detectable conserved cysteine (T/F) | +| conserved_f | Whether or not the rearrangement contains a detectable conserved phenylalanine or equivalent (T/F) | +| legacy_v_call | What older versions of Decombinator (i.e. <= v3) referred to this V gene as: e.g. 'TRBV12-3,TRBV12-4' in v4 was previously referred to just as 'TRBV12-4' | +| legacy_j_call | What older versions of Decombinator (<= 3) referred to this J gene as | +| v_alleles | List of V gene alleles covered by this rearrangment's tag | +| j_alleles | List of J gene alleles covered by this rearrangment's tag | +| v_gene_functionality | [IMGT predicted functionality](http://www.imgt.org/IMGTScientificChart/SequenceDescription/IMGTfunctionality.html#P1-2) of V gene (or genes) used in this rearrangement (F/ORF/P, comma delimited) | +| j_gene_functionality | [IMGT predicted functionality](http://www.imgt.org/IMGTScientificChart/SequenceDescription/IMGTfunctionality.html#P1-2) of J gene (or genes) used in this rearrangement (F/ORF/P, comma delimited) | +| sequence_alignment | Format required field - left blank | +| germline_alignment | Format required field - left blank | +| v_cigar | Format required field - left blank | +| d_cigar | Format required field - left blank | +| j_cigar | Format required field - left blank | +| av_UMI_cluster_size | The average UMI count for this particular sequence | + +You can also use the 'nonproductivefilter' flag (`-npf`) to suppress the output of non-productive rearrangements. + +As the hypervariable region and the primary site of antigenic contact, the CDR3 is almost certainly going to be the region of most interest for many analyses. By convention, the [CDR3 junction is defined as]() running from the position of the second conserved cysteine encoded in the 3' of the V gene to the phenylalanine in the conserved 'FGXG' motif in the J gene. However, some genes use non-canonical residues/motifs, and the position of these motifs varies. + +In looking for CDR3s, we also find 'non-productive' reads, i.e. those that don't appear to be able to make productive, working TCRs. This is determined based on the presence of stop codons, being out of frame, or lacking appropriate CDR3 motifs. The process occurs like so: @@ -336,12 +357,12 @@ CASSLEAGINSPLHF ``` In order to do so, a third kind of supplementary data file is used, .translate files, which provide the additional information required for CDR3 extraction for each gene type (a/b/g/d, V/J). They are stored in the [TCR tag repository](https://github.com/innate2adaptive/Decombinator-Tags-FASTAs) and meet the same naming conventions as the tag and FASTA files and consist of four comma-delimited fields, detailing: -* Gene name -* Conserved motif position (whether C or FGXG) -* Conserved motif sequence (to account for the non-canonical) -* IMGT-defined gene functionality (F/ORF/P) - +- Gene name +- Conserved motif position (whether C or FGXG) +- Conserved motif sequence (to account for the non-canonical) +- IMGT-defined gene functionality (F/ORF/P) +
## Advanced Usage @@ -393,8 +414,8 @@ For any support, commerical inquiries, or other requests, contact us [here](m.co The published history of the development of the pipeline is covered in the following publications: -* [Thomas et al (2013), Bioinformatics: *Decombinator: a tool for fast, efficient gene assignment in T-cell receptor sequences using a finite state machine*](http://dx.doi.org/10.1093/bioinformatics/btt004) -* [Oakes et al (2017), Frontiers in Immunology: *Quantitative Characterization of the T Cell Receptor Repertoire of Naïve and Memory Subsets Using an Integrated Experimental and Computational Pipeline Which Is Robust, Economical, and Versatile*](https://doi.org/10.3389/fimmu.2017.01267) -* [Uddin et al (2019), Cancer Immunosurveillance: *An Economical, Quantitative, and Robust Protocol for High-Throughput T Cell Receptor Sequencing from Tumor or Blood*](http:/dx.doi.org/10.1007/978-1-4939-8885-3_2) -* [Uddin et al (2019), Methods in Enzymology: *Quantitative analysis of the T cell receptor repertoire*](https://doi.org/10.1016/bs.mie.2019.05.054) -* [Peacock et al (2020), Bioinformatics: *Decombinator V4: an improved AIRR compliant-software package for T-cell receptor sequence annotation*](https://doi.org/10.1093/bioinformatics/btaa758) +- [Thomas et al (2013), Bioinformatics: _Decombinator: a tool for fast, efficient gene assignment in T-cell receptor sequences using a finite state machine_](http://dx.doi.org/10.1093/bioinformatics/btt004) +- [Oakes et al (2017), Frontiers in Immunology: _Quantitative Characterization of the T Cell Receptor Repertoire of Naïve and Memory Subsets Using an Integrated Experimental and Computational Pipeline Which Is Robust, Economical, and Versatile_](https://doi.org/10.3389/fimmu.2017.01267) +- [Uddin et al (2019), Cancer Immunosurveillance: _An Economical, Quantitative, and Robust Protocol for High-Throughput T Cell Receptor Sequencing from Tumor or Blood_](http:/dx.doi.org/10.1007/978-1-4939-8885-3_2) +- [Uddin et al (2019), Methods in Enzymology: _Quantitative analysis of the T cell receptor repertoire_](https://doi.org/10.1016/bs.mie.2019.05.054) +- [Peacock et al (2020), Bioinformatics: _Decombinator V4: an improved AIRR compliant-software package for T-cell receptor sequence annotation_](https://doi.org/10.1093/bioinformatics/btaa758) diff --git a/examples/Logs/2025_08_29_alpha__Decombinator_Summary.csv b/examples/Logs/2025_08_29_alpha__Decombinator_Summary.csv new file mode 100644 index 0000000..f05e387 --- /dev/null +++ b/examples/Logs/2025_08_29_alpha__Decombinator_Summary.csv @@ -0,0 +1,49 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,_alpha +OutputFile,_alpha +DateFinished,2025_08_29 +TimeFinished,13:50:47 +TimeTaken(Seconds),0.01 + +InputArguments:, +species,human +chain,a +extension,n12 +tags,extended +dontgzip,False +allowNs,False +orientation,reverse +lenthreshold,130 +bc_read,R2 +bclength,42 + +NumberReadsInput,106 +NumberReadsDecombined,35 +PercentReadsDecombined,0.33 + +ReadsAssignedUsingHalfTags:, +V1error,0 +V2error,1 +J1error,0 +J2error,0 + +ReadsFilteredOut:, +AmbiguousBaseCall(DCR),0 +AmbiguousBaseCall(Barcode),0 +OverlongInterTagSeq,0 +ImpossibleDeletions,0 +OverlappingTagBoundaries,0 + +ReadsFailedAssignment:, +MultipleVtagMatches,0 +VTagAtEndRead,0 +VDeletionsUndetermined,0 +FoundV1HalfTagNotV2,7 +FoundV2HalfTagNotV1,9 +NoVDetected,55 +MultipleJTagMatches,0 +JDeletionsUndermined,0 +FoundJ1HalfTagNotJ2,0 +FoundJ2HalfTagNotJ1,0 +NoJDetected,0 diff --git a/examples/Logs/2025_08_29_beta__Decombinator_Summary.csv b/examples/Logs/2025_08_29_beta__Decombinator_Summary.csv new file mode 100644 index 0000000..113a91b --- /dev/null +++ b/examples/Logs/2025_08_29_beta__Decombinator_Summary.csv @@ -0,0 +1,49 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,_beta +OutputFile,_beta +DateFinished,2025_08_29 +TimeFinished,13:49:52 +TimeTaken(Seconds),0.01 + +InputArguments:, +species,human +chain,b +extension,n12 +tags,extended +dontgzip,False +allowNs,False +orientation,reverse +lenthreshold,130 +bc_read,R2 +bclength,42 + +NumberReadsInput,106 +NumberReadsDecombined,48 +PercentReadsDecombined,0.453 + +ReadsAssignedUsingHalfTags:, +V1error,1 +V2error,1 +J1error,0 +J2error,0 + +ReadsFilteredOut:, +AmbiguousBaseCall(DCR),0 +AmbiguousBaseCall(Barcode),0 +OverlongInterTagSeq,0 +ImpossibleDeletions,0 +OverlappingTagBoundaries,0 + +ReadsFailedAssignment:, +MultipleVtagMatches,0 +VTagAtEndRead,0 +VDeletionsUndetermined,0 +FoundV1HalfTagNotV2,1 +FoundV2HalfTagNotV1,2 +NoVDetected,55 +MultipleJTagMatches,0 +JDeletionsUndermined,0 +FoundJ1HalfTagNotJ2,0 +FoundJ2HalfTagNotJ1,0 +NoJDetected,0 diff --git a/examples/Logs/2025_08_29_beta__Decombinator_Summary2.csv b/examples/Logs/2025_08_29_beta__Decombinator_Summary2.csv new file mode 100644 index 0000000..71d5c82 --- /dev/null +++ b/examples/Logs/2025_08_29_beta__Decombinator_Summary2.csv @@ -0,0 +1,49 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,_beta +OutputFile,_beta +DateFinished,2025_08_29 +TimeFinished,13:50:00 +TimeTaken(Seconds),0.01 + +InputArguments:, +species,human +chain,b +extension,n12 +tags,extended +dontgzip,False +allowNs,False +orientation,reverse +lenthreshold,130 +bc_read,R2 +bclength,42 + +NumberReadsInput,106 +NumberReadsDecombined,48 +PercentReadsDecombined,0.453 + +ReadsAssignedUsingHalfTags:, +V1error,1 +V2error,1 +J1error,0 +J2error,0 + +ReadsFilteredOut:, +AmbiguousBaseCall(DCR),0 +AmbiguousBaseCall(Barcode),0 +OverlongInterTagSeq,0 +ImpossibleDeletions,0 +OverlappingTagBoundaries,0 + +ReadsFailedAssignment:, +MultipleVtagMatches,0 +VTagAtEndRead,0 +VDeletionsUndetermined,0 +FoundV1HalfTagNotV2,1 +FoundV2HalfTagNotV1,2 +NoVDetected,55 +MultipleJTagMatches,0 +JDeletionsUndermined,0 +FoundJ1HalfTagNotJ2,0 +FoundJ2HalfTagNotJ1,0 +NoJDetected,0 diff --git a/examples/Logs/2025_08_29_beta__Decombinator_Summary3.csv b/examples/Logs/2025_08_29_beta__Decombinator_Summary3.csv new file mode 100644 index 0000000..fc4fc49 --- /dev/null +++ b/examples/Logs/2025_08_29_beta__Decombinator_Summary3.csv @@ -0,0 +1,49 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,_beta +OutputFile,_beta +DateFinished,2025_08_29 +TimeFinished,13:50:54 +TimeTaken(Seconds),0.01 + +InputArguments:, +species,human +chain,b +extension,n12 +tags,extended +dontgzip,False +allowNs,False +orientation,reverse +lenthreshold,130 +bc_read,R2 +bclength,42 + +NumberReadsInput,106 +NumberReadsDecombined,48 +PercentReadsDecombined,0.453 + +ReadsAssignedUsingHalfTags:, +V1error,1 +V2error,1 +J1error,0 +J2error,0 + +ReadsFilteredOut:, +AmbiguousBaseCall(DCR),0 +AmbiguousBaseCall(Barcode),0 +OverlongInterTagSeq,0 +ImpossibleDeletions,0 +OverlappingTagBoundaries,0 + +ReadsFailedAssignment:, +MultipleVtagMatches,0 +VTagAtEndRead,0 +VDeletionsUndetermined,0 +FoundV1HalfTagNotV2,1 +FoundV2HalfTagNotV1,2 +NoVDetected,55 +MultipleJTagMatches,0 +JDeletionsUndermined,0 +FoundJ1HalfTagNotJ2,0 +FoundJ2HalfTagNotJ1,0 +NoJDetected,0 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_alpha_CDR3_Translation_Summary.csv b/examples/Logs/2025_08_29_dcr_TINY_1_alpha_CDR3_Translation_Summary.csv new file mode 100644 index 0000000..7964111 --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_alpha_CDR3_Translation_Summary.csv @@ -0,0 +1,30 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_alpha +OutputFile,TINY_alpha +DateFinished,2025_08_29 +TimeFinished,13:50:48 + +InputArguments:, +species,human +chain,a +tags,extended +dontgzip,False + +NumberUniqueDCRsInput,26 +NumberUniqueDCRsProductive,25 +NumberUniqueDCRsNonProductive,1 + +FunctionalityOfGermlineGenesUsed, +P_V-F,0 +P_V-ORF,0 +P_V-P,0 +P_J-F,0 +P_J-ORF,0 +P_J-P,0 +NP_V-F,0 +NP_V-ORF,0 +NP_V-P,0 +NP_J-F,0 +NP_J-ORF,0 +NP_J-P,0 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_alpha_Collapsing_Summary.csv b/examples/Logs/2025_08_29_dcr_TINY_1_alpha_Collapsing_Summary.csv new file mode 100644 index 0000000..83b5619 --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_alpha_Collapsing_Summary.csv @@ -0,0 +1,38 @@ +Property,Value +Version,5.0.0.dev0 +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_alpha +OutputFile,TINY_alpha +DateFinished,2025_08_29 +TimeFinished,13:50:47 +TimeTaken(Seconds),0.03 + +extension,freq +dontgzip,False +allowNs,False +dontcheckinput,False +barcodeduplication,False +minbcQ,20 +bcQbelowmin,1 +bcthreshold,2 +lenthreshold,130 +percentlevdist,10 +avgQthreshold,30 +positionalbarcodes,False +oligo,M13 + +InputUncollapsedDCRLines,35 +UniqueDCRsPassingFilters,26 +TotalDCRsPassingFilters,30 +PercentDCRPassingFilters(withbarcode),0.857 +UniqueDCRsPostCollapsing,26 +TotalDCRsPostCollapsing,26 +PercentUniqueDCRsKept,1.0 +PercentTotalDCRsKept,0.867 +AverageInputTCRAbundance,1.154 +AverageOutputTCRAbundance,1.0 +AverageRNAduplication,1.154 + +BarcodeFail_ContainedNs,0 +BarcodeFail_SpacersNotFound,0 +BarcodeFail_LowQuality,5 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary.csv b/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary.csv new file mode 100644 index 0000000..3a1c89d --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary.csv @@ -0,0 +1,30 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_beta +OutputFile,TINY_beta +DateFinished,2025_08_29 +TimeFinished,13:49:53 + +InputArguments:, +species,human +chain,b +tags,extended +dontgzip,False + +NumberUniqueDCRsInput,36 +NumberUniqueDCRsProductive,35 +NumberUniqueDCRsNonProductive,1 + +FunctionalityOfGermlineGenesUsed, +P_V-F,0 +P_V-ORF,0 +P_V-P,0 +P_J-F,0 +P_J-ORF,0 +P_J-P,0 +NP_V-F,0 +NP_V-ORF,0 +NP_V-P,0 +NP_J-F,0 +NP_J-ORF,0 +NP_J-P,0 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary2.csv b/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary2.csv new file mode 100644 index 0000000..5cfc503 --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary2.csv @@ -0,0 +1,30 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_beta +OutputFile,TINY_beta +DateFinished,2025_08_29 +TimeFinished,13:50:00 + +InputArguments:, +species,human +chain,b +tags,extended +dontgzip,False + +NumberUniqueDCRsInput,36 +NumberUniqueDCRsProductive,35 +NumberUniqueDCRsNonProductive,1 + +FunctionalityOfGermlineGenesUsed, +P_V-F,0 +P_V-ORF,0 +P_V-P,0 +P_J-F,0 +P_J-ORF,0 +P_J-P,0 +NP_V-F,0 +NP_V-ORF,0 +NP_V-P,0 +NP_J-F,0 +NP_J-ORF,0 +NP_J-P,0 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary3.csv b/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary3.csv new file mode 100644 index 0000000..df0d3ad --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_beta_CDR3_Translation_Summary3.csv @@ -0,0 +1,30 @@ +Property,Value +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_beta +OutputFile,TINY_beta +DateFinished,2025_08_29 +TimeFinished,13:50:55 + +InputArguments:, +species,human +chain,b +tags,extended +dontgzip,False + +NumberUniqueDCRsInput,36 +NumberUniqueDCRsProductive,35 +NumberUniqueDCRsNonProductive,1 + +FunctionalityOfGermlineGenesUsed, +P_V-F,0 +P_V-ORF,0 +P_V-P,0 +P_J-F,0 +P_J-ORF,0 +P_J-P,0 +NP_V-F,0 +NP_V-ORF,0 +NP_V-P,0 +NP_J-F,0 +NP_J-ORF,0 +NP_J-P,0 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary.csv b/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary.csv new file mode 100644 index 0000000..b4d3757 --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary.csv @@ -0,0 +1,38 @@ +Property,Value +Version,5.0.0.dev0 +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_beta +OutputFile,TINY_beta +DateFinished,2025_08_29 +TimeFinished,13:49:52 +TimeTaken(Seconds),0.05 + +extension,freq +dontgzip,False +allowNs,False +dontcheckinput,False +barcodeduplication,False +minbcQ,20 +bcQbelowmin,1 +bcthreshold,2 +lenthreshold,130 +percentlevdist,10 +avgQthreshold,30 +positionalbarcodes,False +oligo,M13 + +InputUncollapsedDCRLines,48 +UniqueDCRsPassingFilters,39 +TotalDCRsPassingFilters,44 +PercentDCRPassingFilters(withbarcode),0.917 +UniqueDCRsPostCollapsing,36 +TotalDCRsPostCollapsing,37 +PercentUniqueDCRsKept,0.923 +PercentTotalDCRsKept,0.841 +AverageInputTCRAbundance,1.128 +AverageOutputTCRAbundance,1.028 +AverageRNAduplication,1.189 + +BarcodeFail_ContainedNs,0 +BarcodeFail_SpacersNotFound,1 +BarcodeFail_LowQuality,3 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary2.csv b/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary2.csv new file mode 100644 index 0000000..6d4f49f --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary2.csv @@ -0,0 +1,38 @@ +Property,Value +Version,5.0.0.dev0 +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_beta +OutputFile,TINY_beta +DateFinished,2025_08_29 +TimeFinished,13:50:00 +TimeTaken(Seconds),0.04 + +extension,freq +dontgzip,False +allowNs,False +dontcheckinput,False +barcodeduplication,False +minbcQ,20 +bcQbelowmin,1 +bcthreshold,2 +lenthreshold,130 +percentlevdist,10 +avgQthreshold,30 +positionalbarcodes,False +oligo,M13 + +InputUncollapsedDCRLines,48 +UniqueDCRsPassingFilters,39 +TotalDCRsPassingFilters,44 +PercentDCRPassingFilters(withbarcode),0.917 +UniqueDCRsPostCollapsing,36 +TotalDCRsPostCollapsing,37 +PercentUniqueDCRsKept,0.923 +PercentTotalDCRsKept,0.841 +AverageInputTCRAbundance,1.128 +AverageOutputTCRAbundance,1.028 +AverageRNAduplication,1.189 + +BarcodeFail_ContainedNs,0 +BarcodeFail_SpacersNotFound,1 +BarcodeFail_LowQuality,3 diff --git a/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary3.csv b/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary3.csv new file mode 100644 index 0000000..b4547e4 --- /dev/null +++ b/examples/Logs/2025_08_29_dcr_TINY_1_beta_Collapsing_Summary3.csv @@ -0,0 +1,38 @@ +Property,Value +Version,5.0.0.dev0 +Directory,/home/mcowley/projects/decombinator/Decombinator/examples +InputFile,TINY_beta +OutputFile,TINY_beta +DateFinished,2025_08_29 +TimeFinished,13:50:54 +TimeTaken(Seconds),0.04 + +extension,freq +dontgzip,False +allowNs,False +dontcheckinput,False +barcodeduplication,False +minbcQ,20 +bcQbelowmin,1 +bcthreshold,2 +lenthreshold,130 +percentlevdist,10 +avgQthreshold,30 +positionalbarcodes,False +oligo,M13 + +InputUncollapsedDCRLines,48 +UniqueDCRsPassingFilters,39 +TotalDCRsPassingFilters,44 +PercentDCRPassingFilters(withbarcode),0.917 +UniqueDCRsPostCollapsing,36 +TotalDCRsPostCollapsing,37 +PercentUniqueDCRsKept,0.923 +PercentTotalDCRsKept,0.841 +AverageInputTCRAbundance,1.128 +AverageOutputTCRAbundance,1.028 +AverageRNAduplication,1.189 + +BarcodeFail_ContainedNs,0 +BarcodeFail_SpacersNotFound,1 +BarcodeFail_LowQuality,3 diff --git a/examples/cli_example.sh b/examples/cli_example.sh new file mode 100755 index 0000000..15aaf58 --- /dev/null +++ b/examples/cli_example.sh @@ -0,0 +1,10 @@ +# This file shows an example of how to run decombinator on a small test input +# Ensure you have created a virtual environment and installed decombinator via "pip install decombinator" + +echo "Example pipeline aligning alpha chains..." +decombinator pipeline -in ./input/TINY_1.fq -c a -br R2 -bl 42 -ol M13 +echo "Example pipeline completed aligning alpha chains." + +echo "Example pipeline aligning beta chains..." +decombinator pipeline -in ./input/TINY_1.fq -c b -br R2 -bl 42 -ol M13 +echo "Example pipeline completed aligning beta chains." diff --git a/examples/expected_output/Logs/2023_01_31_alpha_TINY_1_Decombinator_Summary.csv b/examples/expected_output/Logs/2023_01_31_alpha_TINY_1_Decombinator_Summary.csv new file mode 100644 index 0000000..efdcfae --- /dev/null +++ b/examples/expected_output/Logs/2023_01_31_alpha_TINY_1_Decombinator_Summary.csv @@ -0,0 +1,49 @@ +Property,Value +Directory,/home/yutanagano/Projects/Decombinator-Test-Data +InputFile,TINY_1.fq.gz +OutputFile,dcr_TINY_1_alpha.n12.gz +DateFinished,2023_01_31 +TimeFinished,14:59:14 +TimeTaken(Seconds),0.01 + +InputArguments:, +species,human +chain,a +extension,n12 +tags,extended +dontgzip,False +allowNs,False +orientation,reverse +lenthreshold,130 +bc_read,R2 +bclength,42 + +NumberReadsInput,100 +NumberReadsDecombined,35 +PercentReadsDecombined,0.35 + +ReadsAssignedUsingHalfTags:, +V1error,0 +V2error,1 +J1error,0 +J2error,0 + +ReadsFilteredOut:, +AmbiguousBaseCall(DCR),0 +AmbiguousBaseCall(Barcode),0 +OverlongInterTagSeq,0 +ImpossibleDeletions,0 +OverlappingTagBoundaries,0 + +ReadsFailedAssignment:, +MultipleVtagMatches,0 +VTagAtEndRead,0 +VDeletionsUndetermined,0 +FoundV1HalfTagNotV2,7 +FoundV2HalfTagNotV1,9 +NoVDetected,49 +MultipleJTagMatches,0 +JDeletionsUndermined,0 +FoundJ1HalfTagNotJ2,0 +FoundJ2HalfTagNotJ1,0 +NoJDetected,0 diff --git a/examples/expected_output/Logs/2023_01_31_beta_TINY_1_Decombinator_Summary.csv b/examples/expected_output/Logs/2023_01_31_beta_TINY_1_Decombinator_Summary.csv new file mode 100644 index 0000000..8817067 --- /dev/null +++ b/examples/expected_output/Logs/2023_01_31_beta_TINY_1_Decombinator_Summary.csv @@ -0,0 +1,49 @@ +Property,Value +Directory,/home/yutanagano/Projects/Decombinator-Test-Data +InputFile,TINY_1.fq.gz +OutputFile,dcr_TINY_1_beta.n12.gz +DateFinished,2023_01_31 +TimeFinished,14:59:16 +TimeTaken(Seconds),0.01 + +InputArguments:, +species,human +chain,b +extension,n12 +tags,extended +dontgzip,False +allowNs,False +orientation,reverse +lenthreshold,130 +bc_read,R2 +bclength,42 + +NumberReadsInput,100 +NumberReadsDecombined,42 +PercentReadsDecombined,0.42 + +ReadsAssignedUsingHalfTags:, +V1error,1 +V2error,1 +J1error,0 +J2error,0 + +ReadsFilteredOut:, +AmbiguousBaseCall(DCR),0 +AmbiguousBaseCall(Barcode),0 +OverlongInterTagSeq,0 +ImpossibleDeletions,0 +OverlappingTagBoundaries,0 + +ReadsFailedAssignment:, +MultipleVtagMatches,0 +VTagAtEndRead,0 +VDeletionsUndetermined,0 +FoundV1HalfTagNotV2,1 +FoundV2HalfTagNotV1,2 +NoVDetected,55 +MultipleJTagMatches,0 +JDeletionsUndermined,0 +FoundJ1HalfTagNotJ2,0 +FoundJ2HalfTagNotJ1,0 +NoJDetected,0 diff --git a/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_alpha_CDR3_Translation_Summary.csv b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_alpha_CDR3_Translation_Summary.csv new file mode 100644 index 0000000..343edad --- /dev/null +++ b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_alpha_CDR3_Translation_Summary.csv @@ -0,0 +1,30 @@ +Property,Value +Directory,/home/yutanagano/Projects/Decombinator-Test-Data +InputFile,dcr_TINY_1_alpha.freq.gz +OutputFile,dcr_TINY_1_alpha.tsv.gz +DateFinished,2023_01_31 +TimeFinished,14:59:17 + +InputArguments:, +species,human +chain,a +tags,extended +dontgzip,False + +NumberUniqueDCRsInput,26 +NumberUniqueDCRsProductive,25 +NumberUniqueDCRsNonProductive,1 + +FunctionalityOfGermlineGenesUsed, +P_V-F,0 +P_V-ORF,0 +P_V-P,0 +P_J-F,0 +P_J-ORF,0 +P_J-P,0 +NP_V-F,0 +NP_V-ORF,0 +NP_V-P,0 +NP_J-F,0 +NP_J-ORF,0 +NP_J-P,0 diff --git a/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_alpha_Collapsing_Summary.csv b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_alpha_Collapsing_Summary.csv new file mode 100644 index 0000000..d8bae09 --- /dev/null +++ b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_alpha_Collapsing_Summary.csv @@ -0,0 +1,38 @@ +Property,Value +Version,4.0.4 +Directory,/home/yutanagano/Projects/Decombinator-Test-Data +InputFile,dcr_TINY_1_alpha.n12.gz +OutputFile,dcr_TINY_1_alpha.freq.gz +DateFinished,2023_01_31 +TimeFinished,14:59:16 +TimeTaken(Seconds),0.0 + +extension,freq +dontgzip,False +allowNs,False +dontcheckinput,False +barcodeduplication,False +minbcQ,20 +bcQbelowmin,1 +bcthreshold,2 +lenthreshold,130 +percentlevdist,10 +avgQthreshold,30 +positionalbarcodes,False +oligo,M13 + +InputUncollapsedDCRLines,35 +UniqueDCRsPassingFilters,26 +TotalDCRsPassingFilters,30 +PercentDCRPassingFilters(withbarcode),0.857 +UniqueDCRsPostCollapsing,26 +TotalDCRsPostCollapsing,26 +PercentUniqueDCRsKept,1.0 +PercentTotalDCRsKept,0.867 +AverageInputTCRAbundance,1.154 +AverageOutputTCRAbundance,1.0 +AverageRNAduplication,1.154 + +BarcodeFail_ContainedNs,0 +BarcodeFail_SpacersNotFound,0 +BarcodeFail_LowQuality,5 diff --git a/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_beta_CDR3_Translation_Summary.csv b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_beta_CDR3_Translation_Summary.csv new file mode 100644 index 0000000..4acd707 --- /dev/null +++ b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_beta_CDR3_Translation_Summary.csv @@ -0,0 +1,30 @@ +Property,Value +Directory,/home/yutanagano/Projects/Decombinator-Test-Data +InputFile,dcr_TINY_1_beta.freq.gz +OutputFile,dcr_TINY_1_beta.tsv.gz +DateFinished,2023_01_31 +TimeFinished,14:59:18 + +InputArguments:, +species,human +chain,b +tags,extended +dontgzip,False + +NumberUniqueDCRsInput,35 +NumberUniqueDCRsProductive,34 +NumberUniqueDCRsNonProductive,1 + +FunctionalityOfGermlineGenesUsed, +P_V-F,0 +P_V-ORF,0 +P_V-P,0 +P_J-F,0 +P_J-ORF,0 +P_J-P,0 +NP_V-F,0 +NP_V-ORF,0 +NP_V-P,0 +NP_J-F,0 +NP_J-ORF,0 +NP_J-P,0 diff --git a/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_beta_Collapsing_Summary.csv b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_beta_Collapsing_Summary.csv new file mode 100644 index 0000000..32568f5 --- /dev/null +++ b/examples/expected_output/Logs/2023_01_31_dcr_TINY_1_beta_Collapsing_Summary.csv @@ -0,0 +1,38 @@ +Property,Value +Version,4.0.4 +Directory,/home/yutanagano/Projects/Decombinator-Test-Data +InputFile,dcr_TINY_1_beta.n12.gz +OutputFile,dcr_TINY_1_beta.freq.gz +DateFinished,2023_01_31 +TimeFinished,14:59:16 +TimeTaken(Seconds),0.0 + +extension,freq +dontgzip,False +allowNs,False +dontcheckinput,False +barcodeduplication,False +minbcQ,20 +bcQbelowmin,1 +bcthreshold,2 +lenthreshold,130 +percentlevdist,10 +avgQthreshold,30 +positionalbarcodes,False +oligo,M13 + +InputUncollapsedDCRLines,42 +UniqueDCRsPassingFilters,35 +TotalDCRsPassingFilters,38 +PercentDCRPassingFilters(withbarcode),0.905 +UniqueDCRsPostCollapsing,35 +TotalDCRsPostCollapsing,35 +PercentUniqueDCRsKept,1.0 +PercentTotalDCRsKept,0.921 +AverageInputTCRAbundance,1.086 +AverageOutputTCRAbundance,1.0 +AverageRNAduplication,1.086 + +BarcodeFail_ContainedNs,0 +BarcodeFail_SpacersNotFound,1 +BarcodeFail_LowQuality,3 diff --git a/examples/expected_output/dcr_TINY_1_alpha.freq b/examples/expected_output/dcr_TINY_1_alpha.freq new file mode 100644 index 0000000..a719398 --- /dev/null +++ b/examples/expected_output/dcr_TINY_1_alpha.freq @@ -0,0 +1,26 @@ +3, 7, 4, 5, GCT, 1, 2 +20, 11, 5, 4, GGGCGCCTT, 1, 1 +28, 26, 1, 1, ACCTTA, 1, 1 +15, 37, 4, 12, TT, 1, 1 +1, 22, 4, 0, , 1, 2 +1, 17, 2, 7, CGACTCGG, 1, 1 +8, 2, 0, 7, CGGGGG, 1, 1 +17, 12, 0, 8, TG, 1, 1 +39, 14, 3, 5, CCGG, 1, 2 +13, 7, 0, 1, TCGGGCGAGGGG, 1, 2 +43, 32, 0, 7, GCT, 1, 1 +18, 16, 4, 6, CCTCCGGCCC, 1, 1 +32, 47, 2, 11, TCCTT, 1, 1 +1, 23, 3, 3, G, 1, 1 +15, 36, 5, 0, CGG, 1, 1 +4, 7, 3, 8, CAGGG, 1, 1 +43, 33, 4, 5, CCAGA, 1, 1 +10, 46, 0, 10, CT, 1, 1 +23, 44, 0, 4, GGGGC, 1, 1 +19, 48, 0, 5, CCGTCG, 1, 1 +35, 9, 6, 3, CCACACCTA, 1, 1 +10, 39, 0, 1, GGGG, 1, 1 +3, 10, 1, 5, GGATACAG, 1, 1 +30, 34, 8, 0, CGGTT, 1, 1 +13, 24, 0, 3, GGGGGTAGG, 1, 1 +23, 41, 1, 5, TCGGGACGC, 1, 1 diff --git a/examples/expected_output/dcr_TINY_1_alpha.n12 b/examples/expected_output/dcr_TINY_1_alpha.n12 new file mode 100644 index 0000000..13f58e8 --- /dev/null +++ b/examples/expected_output/dcr_TINY_1_alpha.n12 @@ -0,0 +1,35 @@ +3, 7, 4, 5, GCT, A00261:687:HMMCCDSX5:3:1101:5050:1094, AAGCTCAGTGATTCAGCCACCTACCTCTGTGTGGTGGCTAAAGCTGCAGGCAACAAGCTAACTTTTGGAGGAGG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTTCCGGTCGTGATAAAGTG, FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF,FFFFF:FFF +20, 11, 5, 4, GGGCGCCTT, A00261:687:HMMCCDSX5:3:1101:18331:1110, ACGCTGAGAGACACTGCTGTGTACTATTGCATCGTCAGGGGCGCCTTTTCTGGTTCTGCAAGGCAACTGACCTTTGGATCTGG, FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAACACCCTGGTTACCCGTCGTGATGCCTTT, FF,FFF,FFFF:,F,FFFFFF:FFFF::FFFFF:F,F,F::F +3, 7, 4, 5, GCT, A00261:687:HMMCCDSX5:3:1101:4417:1219, AAGCTCAGTGATTCAGCCACCTACCTCTGTGTGGTGGCTAAAGCTGCAGGCAACAAGCTAACTTTTGGAGGAGG, FFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTTCCGGTCGTGATAAAGTG, FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +28, 26, 1, 1, ACCTTA, A00261:687:HMMCCDSX5:3:1101:3794:1235, CAGACCGGAGACTCGGCCATCTACCTCTGTGCTGTGGAGACCTTAAATGCTGGCAACAACCGTAAGCTGATTTGGGGATTGGGA, F::,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFF,FFFFF:FFFFFFF:FFFFF, GTCGTGACTGGGAAAACCCTGGTTAGCCGTCGTGATGTGCTA, FFFFFFFFFFF:FF,FFFFF:FFFFFF:FFFF,FFF,FFFFF +15, 37, 4, 12, TT, A00261:687:HMMCCDSX5:3:1101:19027:1313, TCAGCCTGGTGACTCAGCCACCTACCTCTGTGCTGTTTGAAATGAGAAATTAACCTTTGGGACTGG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCACTTTGTCGTGATACGATC, FFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFF +1, 22, 4, 0, , A00261:687:HMMCCDSX5:3:1101:29170:1344, GAAGGAGCTCCAGATGAAAGACTCTGCCTCTTACCTCTGTGCTGTGATGGATAGCAACTATCAGTTAATCTGGGGCGCTGGGACCAA, FFFFFFFFFFFFF,FFF,FFFFFF,FFFFFFF,,FFFFFFFFF:FFFFFFFF:FFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTGCCTGTCGTGATTTCTAT, FFFFFFFFFFFFFFFFFFFFF:FFFFF:FFFF::FFFFF,FF +1, 17, 2, 7, CGACTCGG, A00261:687:HMMCCDSX5:3:1101:21260:1360, GAAGGAGCTCCAGATGAAAGACTCTGCCTCTTACCTCTGTGCTGTGAGACGACTCGGAGGAAACACACCTCTTGTCTTTGGAAAGGGCAC, FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCGCGCGGTCGTGATAATCAT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF +8, 2, 0, 7, CGGGGG, A00261:687:HMMCCDSX5:3:1101:26585:1407, TGGGGGACTCCGCAATGTACTTCTGTGCAATGAGAGAGGGCGGGGGTAGCAGCTATAAATTGATCTTCGGGAGTGGGAC, FFFFFFFFFF,FFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCTCAAAGTCGTGATCTCCTG, FFFFFFFFFFFFF:FFFFF:FF:FFFFFF:FFFF,:FFFFFF +17, 12, 0, 8, TG, A00261:687:HMMCCDSX5:3:1101:30653:1407, CAGCCTGGAGACTCAGCCACCTACTTCTGTGCAGCAAGCATGAACCAGGGAGGAAAGCTTATCTTCGGACAGGG, FFFFFF:FFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTGCCAGTCGTGATGTCTAT, :FFF,FFFFFFFFFFFFFFFFFFF,FFFFFF:F,FFFFFFFF +39, 14, 3, 5, CCGG, A00261:687:HMMCCDSX5:3:1101:30309:1438, ATATGAGCGACGCGGCTGAGTACTTCTGTGTTGTGAGCCGGAACTATGGTCAGAATTTTGTCTTTGGTCCCGGAAC, FFFF::FFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGGTATTCGTCGTGATGGTCCG, ,FF,FF:FFFFFFFFF,FFFFF:FFFFFFFF:FFFF:FFFF, +13, 7, 0, 1, TCGGGCGAGGGG, A00261:687:HMMCCDSX5:3:1101:24532:1454, GGGAGGCAGATGCTGCTGTTTACTACTGTGCTGTGGAGGATCGGGCGAGGGGGATCAAAGCTGCAGGCAACAAGCTAACTTTTGGAGGAGG, F:FF:FFFFFF:FFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTCCTCGTCGTGATGGCTAC, ::FFFFFFFFFFFFFFFFFFF:FFF:FFFFFFFFFF:F:FFF +43, 32, 0, 7, GCT, A00261:687:HMMCCDSX5:3:1101:28890:1454, AAGTGTCAGACTCAGCGGTGTACTTCTGTGCTCTGAGTGAGCTCAATGACATGCGCTTTGGAGCAGGGACCAGACT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCCTCAGTCGTGATATCTAG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFF +13, 7, 0, 1, TCGGGCGAGGGG, A00261:687:HMMCCDSX5:3:1101:24542:1470, GGGAGGCAGATGCTGCTGTTTACTACTGTGCTGTGGAGGATCGGGCGAGGGGGATCAAAGCTGCAGGCAACAAGCTAACTTTTGGAGGAGG, FFFFFFFFFFF:FFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:F:FFFF, GTCGTGACTGGGAAAACCCTGGTTCCTCGTCGTGATGGCTAC, ,FF:FFFFF,FFFFFF:F:FFF:FFF:,FFFFFFFFFF:FF: +18, 16, 4, 6, CCTCCGGCCC, A00261:687:HMMCCDSX5:3:1101:10854:1532, TCCCAGCCTGAAGACTCAGCCACATACCTCTGTGCCCCTCCGGCCCCTGGGGCTGGGAGTTACCAACTCACTTTCGGGAA, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGAGTGACGTCGTGATTTGTAT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +32, 47, 2, 11, TCCTT, A00261:687:HMMCCDSX5:3:1101:13964:1564, TCCCTGAGCGACACTGCTGTGTACTACTGCCTCGTGGGTGATCCTTACAACAGACTCGCTTTTGGGAAGGGGAAC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTGTGCGTCGTGATTACTTT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +1, 23, 3, 3, G, A00261:687:HMMCCDSX5:3:1101:27597:1658, GAAGGAGCTCCAGATGAAAGACTCTGCCTCTTACCTCTGTGCTGTGAGGTATAACACCGACAAGCTCATCTTTGGGACTGGGACCA, FFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCTTTTAGTCGTGATATTCAA, FFFF::FFFFFFFFF::FFFFFFFFFFFFFFFFFFFFF:FFF +39, 14, 3, 5, CCGG, A00261:687:HMMCCDSX5:3:1101:31901:1658, ATATGAGCGACGCGGCTGAGTACTTCTGTGTTGTGAGCCGGAACTATGGTCAGAATTTTGTCTTTGGTCCCGGAAC, FFF,F,F:FFF:FF,FF:FFFF,,F,FFFF:FFFFFF,FFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:F, GTCGTGACTGGGAAAACCCTGGGTATTCGTCGTGATGGTCCG, FFFFFFFFFF:FFFF,FF:FFFFFFFFF:FFFFFFFFFFF,F +15, 36, 5, 0, CGG, A00261:687:HMMCCDSX5:3:1101:3332:1720, TCAGCCTGGTGACTCAGCCACCTACCTCTGTGCTGCGGTGGAATATGGAAACAAGCTGGTCTTTGGCGCAGGAACCAT, FFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCTACCGTCGTGATTCTACG, FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFF:F,FFFF +7, 15, 3, 11, TTC, A00261:687:HMMCCDSX5:3:1101:18041:1736, CAACCTGGAGACTCAGCTGTCTACTTTTGTGCAGAGATTCCAGGCAAATCAACCTTTGGGGATGGGACT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF, GTCGTGACTGGGAAAACCCTGGAGCGTAGTCGTGATCTGTAG, FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFF:,F,:F +4, 7, 3, 8, CAGGG, A00261:687:HMMCCDSX5:3:1101:26350:1752, CTCATCAGAGACTCCCAGCCCAGTGATTCAGCCACCTACCTCTGTGCCGTGACAGGGGCTGCAGGCAACAAGCTAACTTTTGGAGGAGG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCTCTTGTCGTGATATGCAA, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +6, 49, 0, 9, TA, A00261:687:HMMCCDSX5:3:1101:3287:1767, CAACCTGAAGACTCGGCTGTCTACTTCTGTGCAGCAAGTATAGGAGGCTTCAAAACTATCTTTGGAGCAGGAA, FFFFFFFFFFFFFFFF:FFFFF:FFFFFFF:FFFFFFFFFF:FFFFFFFF::FFFFFFFF:FFFFF:FFFFFF, GTCGTGACTGGGAAAACCCTGGTTTCATGTCGTGATAACAAG, FFFFFFFFF:FFFFFFFFFFF,FFFFFFFFFF,FFFF,:F,F +43, 33, 4, 5, CCAGA, A00261:687:HMMCCDSX5:3:1101:2302:1814, AAGTGTCAGACTCAGCGGTGTACTTCTGTGCTCTGACCAGAACCGGCACTGCCAGTAAACTCACCTTTGGGACTGG, FFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTATTGCGTCGTGATCCATAT, FFFFFFFFFFFFF:FFFFF,FFF:FFFFFFFFFFFFFFFFF: +10, 46, 0, 10, CT, A00261:687:HMMCCDSX5:3:1101:6189:1814, CGGGCAGCAGACACTGCTTCTTACTTCTGTGCTACGGACGCTGGAAGCTACATACCTACATTTGGAAGAGGA, FFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGGAACATGTCGGATCCTCATG, FFFFFFFFFFFFFFFFFFFF:FFFFFF:FFFFFFFFFFFFFF +1, 22, 4, 0, , A00261:687:HMMCCDSX5:3:1101:29243:1814, GAAGGAGCTCCAGATGAAAGACTCTGCCTCTTACCTCTGTGCTGTGATGGATAGCAACTATCAGTTAATCTGGGGCGCTGGGACCAA, F:FFFFFFF,FFFFFFFFFFFFFFFF:FFFFFF::FFFFFFFFFF:FFFFFFF:FFFFFFF:FFFFF:F:FFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTGCCTGTCGTGATTTCTAT, FFFFFFFFFFFFFFFF,:FFF,:::F:FFFFF:FFF,FFFFF +32, 9, 0, 7, TAGG, A00261:687:HMMCCDSX5:3:1101:3296:1846, TCCCTGAGCGACACTGCTGTGTACTACTGCCTCGTGGGTGACATAGGCGACTACAAGCTCAGCTTTGGAGCCGGAACCAC, ,FF,,FFFFFFF,FFF:FFFFFFFF:FFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFF,FFFFFFF:F,FFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGCACAGTCGTGATTCCGCG, ,,,F,:FF,:,:FFFFF,FF:,F:FFFF:FF:F,F::F,F:, +22, 19, 0, 11, CTGAA, A00261:687:HMMCCDSX5:3:1101:3314:1846, GCCCAGCCTGGTGATACAGGCCTCTACCTCTGTGCAGGAGCTGAAGACAAGATCATCTTTGGAAAAGGGACACG, FFFFFFFFFFFFFFFF:FFF:FFFFFFFF:FFF:F:F::F,F,F,FFFFFFFF:,F:FFFF,FFFFFF:FFFFF, GTCGTGACTGGGAAAACCCTGGCCACTAGTCGTGATTTAATG, FFFFFF:FFFFFFFFFFFFFFFFFF,FF:FF:FFFF,FFFFF +23, 44, 0, 4, GGGGC, A00261:687:HMMCCDSX5:3:1101:9833:1892, CAGCCTGGAGACTCTGCAGTGTACTTCTGTGCAGCAAGCGGGGGCACTGGAGCCAATAGTAAGCTGACATTTGGAAAAGGA, FFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFF:FFFFFF, GTCGTGACTGGGAAAACCCTGGCGTGACGTCGTGATTTTAAG, FFFFFFFFF:F:FFFFFFF:F,:FFFF:F,FFFFFFFFFFFF +19, 48, 0, 5, CCGTCG, A00261:687:HMMCCDSX5:3:1101:16848:1924, CACCCAGACTACAGATGTAGGAACCTACTTCTGTGCAGGGCCGTCGACAGGCTTTCAGAAACTTGTATTTGGAACTGGCAC, FFFFFFFFFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCACATGGTCGTGATGAGACC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +35, 9, 6, 3, CCACACCTA, A00261:687:HMMCCDSX5:3:1101:19090:1924, CAGACTGGGGACTCAGCTATCTACTTCTGTGCAGCCACACCTACTAACGACTACAAGCTCAGCTTTGGAGCCGGAACCAC, ::,FFFFFFFFFF::FFFFFFFFFFF:FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGACTCAAGTCGTGATTTATTC, FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFF:FFFF +10, 39, 0, 1, GGGG, A00261:687:HMMCCDSX5:3:1101:20871:1939, CGGGCAGCAGACACTGCTTCTTACTTCTGTGCTACGGACGGGGGGGACACGGGCAGGAGAGCACTTACTTTTGGGAGTGGAAC, FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF, GTCGTGACTGGGAAAACCCTGGTACACCGTCGTGATAAATAT, :FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFF:FFFFFFF +3, 10, 1, 5, GGATACAG, A00261:687:HMMCCDSX5:3:1101:10881:1986, AAGCTCAGTGATTCAGCCACCTACCTCTGTGTGGTGAACGGATACAGCTTCAACAAATTTTACTTTGGATCTGGGACCAAAC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTCCTAGTCGTGATTTTAGT, FFFFFFFFFFFFFFFFFFFFFF,FFFFFF:FFFFFFFFFFFF +30, 34, 8, 0, CGGTT, A00261:687:HMMCCDSX5:3:1101:12165:1986, CTGGGGGATGCCGCGATGTATTTCTGTGCTTACGGTTTGTATTCAGGAGGAGGTGCTGACGGACTCACCTTTGGCAA, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF, GTCGTGACTGGGAAAACCCTGGAGATCTGTCGTGATTCGGGC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +13, 24, 0, 3, GGGGGTAGG, A00261:687:HMMCCDSX5:3:1101:24008:1986, GGGAGGCAGATGCTGCTGTTTACTACTGTGCTGTGGAGGAGGGGGTAGGAACTGGGGCAAACAACCTCTTCTTTGGGACTGGAACG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTCAGCGTCGTGATTTGTGA, FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF:F,FFFF:FF +39, 28, 10, 8, CTGTGAGTGAA, A00261:687:HMMCCDSX5:3:1101:6144:2049, ATATGAGCGACGCGGCTGAGTACTTCTGTGCTGTGAGTGAAGGTGGCTACAATAAGCTGATTTTTGGAGCAGG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTATTTGTCGTGATTTCTTT, F:FFFF:FFFFFFFFFF:FFF:FFF,FFFFFFFF:FFFFF,F +23, 41, 1, 5, TCGGGACGC, A00261:687:HMMCCDSX5:3:1101:15727:2080, CAGCCTGGAGACTCTGCAGTGTACTTCTGTGCAGCAAGCTCGGGACGCGCTGGTGGTACTAGCTATGGAAAGCTGACATTTGG, FF:FFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCTATCGTCGTGATACCCCT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF diff --git a/examples/expected_output/dcr_TINY_1_alpha.tsv b/examples/expected_output/dcr_TINY_1_alpha.tsv new file mode 100644 index 0000000..208252b --- /dev/null +++ b/examples/expected_output/dcr_TINY_1_alpha.tsv @@ -0,0 +1,27 @@ +sequence_id v_call d_call j_call junction_aa duplicate_count sequence junction decombinator_id rev_comp productive sequence_aa cdr1_aa cdr2_aa vj_in_frame stop_codon conserved_c conserved_f sequence_alignment germline_alignment v_cigar d_cigar j_cigar av_UMI_cluster_size +1 TRAV12-1 TRAJ17 CVVAKAAGNKLTF 1 CGGAAGGAGGTGGAGCAGGATCCTGGACCCTTCAATGTTCCAGAGGGAGCCACTGTCGCTTTCAACTGTACTTACAGCAACAGTGCTTCTCAGTCTTTCTTCTGGTACAGACAGGATTGCAGGAAAGAACCTAAGTTGCTGATGTCCGTATACTCCAGTGGTAATGAAGATGGAAGGTTTACAGCACAGCTCAATAGAGCCAGCCAGTATATTTCCCTGCTCATCAGAGACTCCAAGCTCAGTGATTCAGCCACCTACCTCTGTGTGGTGGCTAAAGCTGCAGGCAACAAGCTAACTTTTGGAGGAGGAACCAGGGTGCTAGTTAAACCAA TGTGTGGTGGCTAAAGCTGCAGGCAACAAGCTAACTTTT 3, 7, 4, 5, GCT F T RKEVEQDPGPFNVPEGATVAFNCTYSNSASQSFFWYRQDCRKEPKLLMSVYSSGNEDGRFTAQLNRASQYISLLIRDSKLSDSATYLCVVAKAAGNKLTFGGGTRVLVKP NSASQS VYSSGN T F T T 2 +2 TRAV26-1 TRAJ22 CIVRGAFSGSARQLTF 1 GATGCTAAGACCACCCAGCCCCCCTCCATGGATTGCGCTGAAGGAAGAGCTGCAAACCTGCCTTGTAATCACTCTACCATCAGTGGAAATGAGTATGTGTATTGGTATCGACAGATTCACTCCCAGGGGCCACAGTATATCATTCATGGTCTAAAAAACAATGAAACCAATGAAATGGCCTCTCTGATCATCACAGAAGACAGAAAGTCCAGCACCTTGATCCTGCCCCACGCTACGCTGAGAGACACTGCTGTGTACTATTGCATCGTCAGGGGCGCCTTTTCTGGTTCTGCAAGGCAACTGACCTTTGGATCTGGGACACAATTGACTGTTTTACCTG TGCATCGTCAGGGGCGCCTTTTCTGGTTCTGCAAGGCAACTGACCTTT 20, 11, 5, 4, GGGCGCCTT F T DAKTTQPPSMDCAEGRAANLPCNHSTISGNEYVYWYRQIHSQGPQYIIHGLKNNETNEMASLIITEDRKSSTLILPHATLRDTAVYYCIVRGAFSGSARQLTFGSGTQLTVLP TISGNEY GLKNN T F T T 1 +3 TRAV36/DV7 TRAJ38 CAVETLNAGNNRKLIW 1 GAAGACAAGGTGGTACAAAGCCCTCTATCTCTGGTTGTCCACGAGGGAGACACCGTAACTCTCAATTGCAGTTATGAAGTGACTAACTTTCGAAGCCTACTATGGTACAAGCAGGAAAAGAAAGCTCCCACATTTCTATTTATGCTAACTTCAAGTGGAATTGAAAAGAAGTCAGGAAGACTAAGTAGCATATTAGATAAGAAAGAACTTTCCAGCATCCTGAACATCACAGCCACCCAGACCGGAGACTCGGCCATCTACCTCTGTGCTGTGGAGACCTTAAATGCTGGCAACAACCGTAAGCTGATTTGGGGATTGGGAACAAGCCTGGCAGTAAATCCGA TGTGCTGTGGAGACCTTAAATGCTGGCAACAACCGTAAGCTGATTTGG 28, 26, 1, 1, ACCTTA F T EDKVVQSPLSLVVHEGDTVTLNCSYEVTNFRSLLWYKQEKKAPTFLFMLTSSGIEKKSGRLSSILDKKELSSILNITATQTGDSAIYLCAVETLNAGNNRKLIWGLGTSLAVNP VTNFRS LTSSGIE T F T T 1 +4 TRAV21 TRAJ48 1 AAACAGGAGGTGACGCAGATTCCTGCAGCTCTGAGTGTCCCAGAAGGAGAAAACTTGGTTCTCAACTGCAGTTTCACTGATAGCGCTATTTACAACCTCCAGTGGTTTAGGCAGGACCCTGGGAAAGGTCTCACATCTCTGTTGCTTATTCAGTCAAGTCAGAGAGAGCAAACAAGTGGAAGACTTAATGCCTCGCTGGATAAATCATCAGGACGTAGTACTTTATACATTGCAGCTTCTCAGCCTGGTGACTCAGCCACCTACCTCTGTGCTGTTTGAAATGAGAAATTAACCTTTGGGACTGGAACAAGACTCACCATCATACCCA 15, 37, 4, 12, TT F F KQEVTQIPAALSVPEGENLVLNCSFTDSAIYNLQWFRQDPGKGLTSLLLIQSSQREQTSGRLNASLDKSSGRSTLYIAASQPGDSATYLCAV*NEKLTFGTGTRLTIIP T T T T 1 +5 TRAV1-2 TRAJ33 CAVMDSNYQLIW 1 GGACAAAACATTGACCAGCCCACTGAGATGACAGCTACGGAAGGTGCCATTGTCCAGATCAACTGCACGTACCAGACATCTGGGTTCAACGGGCTGTTCTGGTACCAGCAACATGCTGGCGAAGCACCCACATTTCTGTCTTACAATGTTCTGGATGGTTTGGAGGAGAAAGGTCGTTTTTCTTCATTCCTTAGTCGGTCTAAAGGGTACAGTTACCTCCTTTTGAAGGAGCTCCAGATGAAAGACTCTGCCTCTTACCTCTGTGCTGTGATGGATAGCAACTATCAGTTAATCTGGGGCGCTGGGACCAAGCTAATTATAAAGCCAG TGTGCTGTGATGGATAGCAACTATCAGTTAATCTGG 1, 22, 4, 0, F T GQNIDQPTEMTATEGAIVQINCTYQTSGFNGLFWYQQHAGEAPTFLSYNVLDGLEEKGRFSSFLSRSKGYSYLLLKELQMKDSASYLCAVMDSNYQLIWGAGTKLIIKP TSGFNG NVLDGL T F T T 2 +6 TRAV1-2 TRAJ29 CAVRRLGGNTPLVF 1 GGACAAAACATTGACCAGCCCACTGAGATGACAGCTACGGAAGGTGCCATTGTCCAGATCAACTGCACGTACCAGACATCTGGGTTCAACGGGCTGTTCTGGTACCAGCAACATGCTGGCGAAGCACCCACATTTCTGTCTTACAATGTTCTGGATGGTTTGGAGGAGAAAGGTCGTTTTTCTTCATTCCTTAGTCGGTCTAAAGGGTACAGTTACCTCCTTTTGAAGGAGCTCCAGATGAAAGACTCTGCCTCTTACCTCTGTGCTGTGAGACGACTCGGAGGAAACACACCTCTTGTCTTTGGAAAGGGCACAAGACTTTCTGTGATTGCAA TGTGCTGTGAGACGACTCGGAGGAAACACACCTCTTGTCTTT 1, 17, 2, 7, CGACTCGG F T GQNIDQPTEMTATEGAIVQINCTYQTSGFNGLFWYQQHAGEAPTFLSYNVLDGLEEKGRFSSFLSRSKGYSYLLLKELQMKDSASYLCAVRRLGGNTPLVFGKGTRLSVIA TSGFNG NVLDGL T F T T 1 +7 TRAV14/DV4 TRAJ12 CAMREGGGSSYKLIF 1 GCCCAGAAGATAACTCAAACCCAACCAGGAATGTTCGTGCAGGAAAAGGAGGCTGTGACTCTGGACTGCACATATGACACCAGTGATCCAAGTTATGGTCTATTCTGGTACAAGCAGCCCAGCAGTGGGGAAATGATTTTTCTTATTTATCAGGGGTCTTATGACCAGCAAAATGCAACAGAAGGTCGCTACTCATTGAATTTCCAGAAGGCAAGAAAATCCGCCAACCTTGTCATCTCCGCTTCACAACTGGGGGACTCAGCAATGTACTTCTGTGCAATGAGAGAGGGCGGGGGTAGCAGCTATAAATTGATCTTCGGGAGTGGGACCAGACTGCTGGTCAGGCCTG TGTGCAATGAGAGAGGGCGGGGGTAGCAGCTATAAATTGATCTTC 8, 2, 0, 7, CGGGGG F T AQKITQTQPGMFVQEKEAVTLDCTYDTSDPSYGLFWYKQPSSGEMIFLIYQGSYDQQNATEGRYSLNFQKARKSANLVISASQLGDSAMYFCAMREGGGSSYKLIFGSGTRLLVRP TSDPSYG QGSYDQQN T F T T 1 +8 TRAV23/DV6 TRAJ23 CAASMNQGGKLIF 1 GGCCAACAGAAGGAGAAAAGTGACCAGCAGCAGGTGAAACAAAGTCCTCAATCTTTGATAGTCCAGAAAGGAGGGATTTCAATTATAAACTGTGCTTATGAGAACACTGCGTTTGACTACTTTCCATGGTACCAACAATTCCCTGGGAAAGGCCCTGCATTATTGATAGCCATACGTCCAGATGTGAGTGAAAAGAAAGAAGGAAGATTCACAATCTCCTTCAATAAAAGTGCCAAGCAGTTCTCATTGCATATCATGGATTCCCAGCCTGGAGACTCAGCCACCTACTTCTGTGCAGCAAGCATGAACCAGGGAGGAAAGCTTATCTTCGGACAGGGAACGGAGTTATCTGTGAAACCCA TGTGCAGCAAGCATGAACCAGGGAGGAAAGCTTATCTTC 17, 12, 0, 8, TG F T GQQKEKSDQQQVKQSPQSLIVQKGGISIINCAYENTAFDYFPWYQQFPGKGPALLIAIRPDVSEKKEGRFTISFNKSAKQFSLHIMDSQPGDSATYFCAASMNQGGKLIFGQGTELSVKP NTAFDY IRPDVSE T F T T 1 +9 TRAV8-2 TRAJ26 CVVSRNYGQNFVF 1 GCCCAGTCGGTGACCCAGCTTGACAGCCACGTCTCTGTCTCTGAAGGAACCCCGGTGCTGCTGAGGTGCAACTACTCATCTTCTTATTCACCATCTCTCTTCTGGTATGTGCAACACCCCAACAAAGGACTCCAGCTTCTCCTGAAGTACACATCAGCGGCCACCCTGGTTAAAGGCATCAACGGTTTTGAGGCTGAATTTAAGAAGAGTGAAACCTCCTTCCACCTGACGAAACCCTCAGCCCATATGAGCGACGCGGCTGAGTACTTCTGTGTTGTGAGCCGGAACTATGGTCAGAATTTTGTCTTTGGTCCCGGAACCAGATTGTCCGTGCTGCCCT TGTGTTGTGAGCCGGAACTATGGTCAGAATTTTGTCTTT 39, 14, 3, 5, CCGG F T AQSVTQLDSHVSVSEGTPVLLRCNYSSSYSPSLFWYVQHPNKGLQLLLKYTSAATLVKGINGFEAEFKKSETSFHLTKPSAHMSDAAEYFCVVSRNYGQNFVFGPGTRLSVLP SSYSPS YTSAATLV T F T T 2 +10 TRAV2 TRAJ17 CAVEDRARGIKAAGNKLTF 1 AAGGACCAAGTGTTTCAGCCTTCCACAGTGGCATCTTCAGAGGGAGCTGTGGTGGAAATCTTCTGTAATCACTCTGTGTCCAATGCTTACAACTTCTTCTGGTACCTTCACTTCCCGGGATGTGCACCAAGACTCCTTGTTAAAGGCTCAAAGCCTTCTCAGCAGGGACGATACAACATGACCTATGAACGGTTCTCTTCATCGCTGCTCATCCTCCAGGTGCGGGAGGCAGATGCTGCTGTTTACTACTGTGCTGTGGAGGATCGGGCGAGGGGGATCAAAGCTGCAGGCAACAAGCTAACTTTTGGAGGAGGAACCAGGGTGCTAGTTAAACCAA TGTGCTGTGGAGGATCGGGCGAGGGGGATCAAAGCTGCAGGCAACAAGCTAACTTTT 13, 7, 0, 1, TCGGGCGAGGGG F T KDQVFQPSTVASSEGAVVEIFCNHSVSNAYNFFWYLHFPGCAPRLLVKGSKPSQQGRYNMTYERFSSSLLILQVREADAAVYYCAVEDRARGIKAAGNKLTFGGGTRVLVKP VSNAYN GSKP T F T T 2 +11 TRAV9-2 TRAJ43 CALSELNDMRF 1 GGAAATTCAGTGACCCAGATGGAAGGGCCAGTGACTCTCTCAGAAGAGGCCTTCCTGACTATAAACTGCACGTACACAGCCACAGGATACCCTTCCCTTTTCTGGTATGTCCAATATCCTGGAGAAGGTCTACAGCTCCTCCTGAAAGCCACGAAGGCTGATGACAAGGGAAGCAACAAAGGTTTTGAAGCCACATACCGTAAAGAAACCACTTCTTTCCACTTGGAGAAAGGCTCAGTTCAAGTGTCAGACTCAGCGGTGTACTTCTGTGCTCTGAGTGAGCTCAATGACATGCGCTTTGGAGCAGGGACCAGACTGACAGTAAAACCAA TGTGCTCTGAGTGAGCTCAATGACATGCGCTTT 43, 32, 0, 7, GCT F T GNSVTQMEGPVTLSEEAFLTINCTYTATGYPSLFWYVQYPGEGLQLLLKATKADDKGSNKGFEATYRKETTSFHLEKGSVQVSDSAVYFCALSELNDMRFGAGTRLTVKP ATGYPS ATKADDK T F T T 1 +12 TRAV24 TRAJ28 CAPPAPGAGSYQLTF 1 ATACTGAACGTGGAACAAAGTCCTCAGTCACTGCATGTTCAGGAGGGAGACAGCACCAATTTCACCTGCAGCTTCCCTTCCAGCAATTTTTATGCCTTACACTGGTACAGATGGGAAACTGCAAAAAGCCCCGAGGCCTTGTTTGTAATGACTTTAAATGGGGATGAAAAGAAGAAAGGACGAATAAGTGCCACTCTTAATACCAAGGAGGGTTACAGCTATTTGTACATCAAAGGATCCCAGCCTGAAGACTCAGCCACATACCTCTGTGCCCCTCCGGCCCCTGGGGCTGGGAGTTACCAACTCACTTTCGGGAAGGGGACCAAACTCTCGGTCATACCAA TGTGCCCCTCCGGCCCCTGGGGCTGGGAGTTACCAACTCACTTTC 18, 16, 4, 6, CCTCCGGCCC F T ILNVEQSPQSLHVQEGDSTNFTCSFPSSNFYALHWYRWETAKSPEALFVMTLNGDEKKKGRISATLNTKEGYSYLYIKGSQPEDSATYLCAPPAPGAGSYQLTFGKGTKLSVIP SSNFYA MTLNGDE T F T T 1 +13 TRAV4 TRAJ7 CLVGDPYNRLAF 1 CTTGCTAAGACCACCCAGCCCATCTCCATGGACTCATATGAAGGACAAGAAGTGAACATAACCTGTAGCCACAACAACATTGCTACAAATGATTATATCACGTGGTACCAACAGTTTCCCAGCCAAGGACCACGATTTATTATTCAAGGATACAAGACAAAAGTTACAAACGAAGTGGCCTCCCTGTTTATCCCTGCCGACAGAAAGTCCAGCACTCTGAGCCTGCCCCGGGTTTCCCTGAGCGACACTGCTGTGTACTACTGCCTCGTGGGTGATCCTTACAACAGACTCGCTTTTGGGAAGGGGAACCAAGTGGTGGTCATACCAA TGCCTCGTGGGTGATCCTTACAACAGACTCGCTTTT 32, 47, 2, 11, TCCTT F T LAKTTQPISMDSYEGQEVNITCSHNNIATNDYITWYQQFPSQGPRFIIQGYKTKVTNEVASLFIPADRKSSTLSLPRVSLSDTAVYYCLVGDPYNRLAFGKGNQVVVIP NIATNDY GYKTK T F T T 1 +14 TRAV1-2 TRAJ34 CAVRYNTDKLIF 1 GGACAAAACATTGACCAGCCCACTGAGATGACAGCTACGGAAGGTGCCATTGTCCAGATCAACTGCACGTACCAGACATCTGGGTTCAACGGGCTGTTCTGGTACCAGCAACATGCTGGCGAAGCACCCACATTTCTGTCTTACAATGTTCTGGATGGTTTGGAGGAGAAAGGTCGTTTTTCTTCATTCCTTAGTCGGTCTAAAGGGTACAGTTACCTCCTTTTGAAGGAGCTCCAGATGAAAGACTCTGCCTCTTACCTCTGTGCTGTGAGGTATAACACCGACAAGCTCATCTTTGGGACTGGGACCAGATTACAAGTCTTTCCAA TGTGCTGTGAGGTATAACACCGACAAGCTCATCTTT 1, 23, 3, 3, G F T GQNIDQPTEMTATEGAIVQINCTYQTSGFNGLFWYQQHAGEAPTFLSYNVLDGLEEKGRFSSFLSRSKGYSYLLLKELQMKDSASYLCAVRYNTDKLIFGTGTRLQVFP TSGFNG NVLDGL T F T T 1 +15 TRAV21 TRAJ47 CAAVEYGNKLVF 1 AAACAGGAGGTGACGCAGATTCCTGCAGCTCTGAGTGTCCCAGAAGGAGAAAACTTGGTTCTCAACTGCAGTTTCACTGATAGCGCTATTTACAACCTCCAGTGGTTTAGGCAGGACCCTGGGAAAGGTCTCACATCTCTGTTGCTTATTCAGTCAAGTCAGAGAGAGCAAACAAGTGGAAGACTTAATGCCTCGCTGGATAAATCATCAGGACGTAGTACTTTATACATTGCAGCTTCTCAGCCTGGTGACTCAGCCACCTACCTCTGTGCTGCGGTGGAATATGGAAACAAACTGGTCTTTGGCGCAGGAACCATTCTGAGAGTCAAGTCCT TGTGCTGCGGTGGAATATGGAAACAAACTGGTCTTT 15, 36, 5, 0, CGG F T KQEVTQIPAALSVPEGENLVLNCSFTDSAIYNLQWFRQDPGKGLTSLLLIQSSQREQTSGRLNASLDKSSGRSTLYIAASQPGDSATYLCAAVEYGNKLVFGAGTILRVKS DSAIYN IQSSQRE T F T T 1 +16 TRAV12-2 TRAJ17 CAVTGAAGNKLTF 1 CAGAAGGAGGTGGAGCAGAATTCTGGACCCCTCAGTGTTCCAGAGGGAGCCATTGCCTCTCTCAACTGCACTTACAGTGACCGAGGTTCCCAGTCCTTCTTCTGGTACAGACAATATTCTGGGAAAAGCCCTGAGTTGATAATGTTCATATACTCCAATGGTGACAAAGAAGATGGAAGGTTTACAGCACAGCTCAATAAAGCCAGCCAGTATGTTTCTCTGCTCATCAGAGACTCCCAGCCCAGTGATTCAGCCACCTACCTCTGTGCCGTGACAGGGGCTGCAGGCAACAAGCTAACTTTTGGAGGAGGAACCAGGGTGCTAGTTAAACCAA TGTGCCGTGACAGGGGCTGCAGGCAACAAGCTAACTTTT 4, 7, 3, 8, CAGGG F T QKEVEQNSGPLSVPEGAIASLNCTYSDRGSQSFFWYRQYSGKSPELIMFIYSNGDKEDGRFTAQLNKASQYVSLLIRDSQPSDSATYLCAVTGAAGNKLTFGGGTRVLVKP DRGSQS IYSNGD T F T T 1 +17 TRAV9-2 TRAJ44 CALTRTGTASKLTF 1 GGAAATTCAGTGACCCAGATGGAAGGGCCAGTGACTCTCTCAGAAGAGGCCTTCCTGACTATAAACTGCACGTACACAGCCACAGGATACCCTTCCCTTTTCTGGTATGTCCAATATCCTGGAGAAGGTCTACAGCTCCTCCTGAAAGCCACGAAGGCTGATGACAAGGGAAGCAACAAAGGTTTTGAAGCCACATACCGTAAAGAAACCACTTCTTTCCACTTGGAGAAAGGCTCAGTTCAAGTGTCAGACTCAGCGGTGTACTTCTGTGCTCTGACCAGAACCGGCACTGCCAGTAAACTCACCTTTGGGACTGGAACAAGACTTCAGGTCACGCTCG TGTGCTCTGACCAGAACCGGCACTGCCAGTAAACTCACCTTT 43, 33, 4, 5, CCAGA F T GNSVTQMEGPVTLSEEAFLTINCTYTATGYPSLFWYVQYPGEGLQLLLKATKADDKGSNKGFEATYRKETTSFHLEKGSVQVSDSAVYFCALTRTGTASKLTFGTGTRLQVTL ATGYPS ATKADDK T F T T 1 +18 TRAV17 TRAJ6 CATDAGSYIPTF 1 AGTCAACAGGGAGAAGAGGATCCTCAGGCCTTGAGCATCCAGGAGGGTGAAAATGCCACCATGAACTGCAGTTACAAAACTAGTATAAACAATTTACAGTGGTATAGACAAAATTCAGGTAGAGGCCTTGTCCACCTAATTTTAATACGTTCAAATGAAAGAGAGAAACACAGTGGAAGATTAAGAGTCACGCTTGACACTTCCAAGAAAAGCAGTTCCTTGTTGATCACGGCTTCCCGGGCAGCAGACACTGCTTCTTACTTCTGTGCTACGGACGCTGGAAGCTACATACCTACATTTGGAAGAGGAACCAGCCTTATTGTTCATCCGT TGTGCTACGGACGCTGGAAGCTACATACCTACATTT 10, 46, 0, 10, CT F T SQQGEEDPQALSIQEGENATMNCSYKTSINNLQWYRQNSGRGLVHLILIRSNEREKHSGRLRVTLDTSKKSSSLLITASRAADTASYFCATDAGSYIPTFGRGTSLIVHP TSINN IRSNERE T F T T 1 +19 TRAV29/DV5 TRAJ56 CAASGGTGANSKLTF 1 AGTCAACAGAAGAATGATGACCAGCAAGTTAAGCAAAATTCACCATCCCTGAGCGTCCAGGAAGGAAGAATTTCTATTCTGAACTGTGACTATACTAACAGCATGTTTGATTATTTCCTATGGTACAAAAAATACCCTGCTGAAGGTCCTACATTCCTGATATCTATAAGTTCCATTAAGGATAAAAATGAAGATGGAAGATTCACTGTCTTCTTAAACAAAAGTGCCAAGCACCTCTCTCTGCACATTGTGCCCTCCCAGCCTGGAGACTCTGCAGTGTACTTCTGTGCAGCAAGCGGGGGCACTGGAGCCAATAGTAAGCTGACATTTGGAAAAGGAATAACTCTGAGTGTTAGACCAG TGTGCAGCAAGCGGGGGCACTGGAGCCAATAGTAAGCTGACATTT 23, 44, 0, 4, GGGGC F T SQQKNDDQQVKQNSPSLSVQEGRISILNCDYTNSMFDYFLWYKKYPAEGPTFLISISSIKDKNEDGRFTVFLNKSAKHLSLHIVPSQPGDSAVYFCAASGGTGANSKLTFGKGITLSVRP NSMFDY ISSIKDK T F T T 1 +20 TRAV25 TRAJ8 CAGPSTGFQKLVF 1 GGACAACAGGTAATGCAAATTCCTCAGTACCAGCATGTACAAGAAGGAGAGGACTTCACCACGTACTGCAATTCCTCAACTACTTTAAGCAATATACAGTGGTATAAGCAAAGGCCTGGTGGACATCCCGTTTTTTTGATACAGTTAGTGAAGAGTGGAGAAGTGAAGAAGCAGAAAAGACTGACATTTCAGTTTGGAGAAGCAAAAAAGAACAGCTCCCTGCACATCACAGCCACCCAGACTACAGATGTAGGAACCTACTTCTGTGCAGGGCCGTCGACAGGCTTTCAGAAACTTGTATTTGGAACTGGCACCCGACTTCTGGTCAGTCCAA TGTGCAGGGCCGTCGACAGGCTTTCAGAAACTTGTATTT 19, 48, 0, 5, CCGTCG F T GQQVMQIPQYQHVQEGEDFTTYCNSSTTLSNIQWYKQRPGGHPVFLIQLVKSGEVKKQKRLTFQFGEAKKNSSLHITATQTTDVGTYFCAGPSTGFQKLVFGTGTRLLVSP TTLSN LVKSGEV T F T T 1 +21 TRAV5 TRAJ20 CAATPTNDYKLSF 1 GGAGAGGATGTGGAGCAGAGTCTTTTCCTGAGTGTCCGAGAGGGAGACAGCTCCGTTATAAACTGCACTTACACAGACAGCTCCTCCACCTACTTATACTGGTATAAGCAAGAACCTGGAGCAGGTCTCCAGTTGCTGACGTATATTTTTTCAAATATGGACATGAAACAAGACCAAAGACTCACTGTTCTATTGAATAAAAAGGATAAACATCTGTCTCTGCGCATTGCAGACACCCAGACTGGGGACTCAGCTATCTACTTCTGTGCAGCCACACCTACTAACGACTACAAGCTCAGCTTTGGAGCCGGAACCACAGTAACTGTAAGAGCAA TGTGCAGCCACACCTACTAACGACTACAAGCTCAGCTTT 35, 9, 6, 3, CCACACCTA F T GEDVEQSLFLSVREGDSSVINCTYTDSSSTYLYWYKQEPGAGLQLLTYIFSNMDMKQDQRLTVLLNKKDKHLSLRIADTQTGDSAIYFCAATPTNDYKLSFGAGTTVTVRA DSSSTY IFSNMDM T F T T 1 +22 TRAV17 TRAJ5 CATDGGDTGRRALTF 1 AGTCAACAGGGAGAAGAGGATCCTCAGGCCTTGAGCATCCAGGAGGGTGAAAATGCCACCATGAACTGCAGTTACAAAACTAGTATAAACAATTTACAGTGGTATAGACAAAATTCAGGTAGAGGCCTTGTCCACCTAATTTTAATACGTTCAAATGAAAGAGAGAAACACAGTGGAAGATTAAGAGTCACGCTTGACACTTCCAAGAAAAGCAGTTCCTTGTTGATCACGGCTTCCCGGGCAGCAGACACTGCTTCTTACTTCTGTGCTACGGACGGGGGGGACACGGGCAGGAGAGCACTTACTTTTGGGAGTGGAACAAGACTCCAAGTGCAACCAA TGTGCTACGGACGGGGGGGACACGGGCAGGAGAGCACTTACTTTT 10, 39, 0, 1, GGGG F T SQQGEEDPQALSIQEGENATMNCSYKTSINNLQWYRQNSGRGLVHLILIRSNEREKHSGRLRVTLDTSKKSSSLLITASRAADTASYFCATDGGDTGRRALTFGSGTRLQVQP TSINN IRSNERE T F T T 1 +23 TRAV12-1 TRAJ21 CVVNGYSFNKFYF 1 CGGAAGGAGGTGGAGCAGGATCCTGGACCCTTCAATGTTCCAGAGGGAGCCACTGTCGCTTTCAACTGTACTTACAGCAACAGTGCTTCTCAGTCTTTCTTCTGGTACAGACAGGATTGCAGGAAAGAACCTAAGTTGCTGATGTCCGTATACTCCAGTGGTAATGAAGATGGAAGGTTTACAGCACAGCTCAATAGAGCCAGCCAGTATATTTCCCTGCTCATCAGAGACTCCAAGCTCAGTGATTCAGCCACCTACCTCTGTGTGGTGAACGGATACAGCTTCAACAAATTTTACTTTGGATCTGGGACCAAACTCAATGTAAAACCAA TGTGTGGTGAACGGATACAGCTTCAACAAATTTTACTTT 3, 10, 1, 5, GGATACAG F T RKEVEQDPGPFNVPEGATVAFNCTYSNSASQSFFWYRQDCRKEPKLLMSVYSSGNEDGRFTAQLNRASQYISLLIRDSKLSDSATYLCVVNGYSFNKFYFGSGTKLNVKP NSASQS VYSSGN T F T T 1 +24 TRAV38-2/DV8 TRAJ45 CAYGLYSGGGADGLTF 1 GCTCAGACAGTCACTCAGTCTCAACCAGAGATGTCTGTGCAGGAGGCAGAGACCGTGACCCTGAGCTGCACATATGACACCAGTGAGAGTGATTATTATTTATTCTGGTACAAGCAGCCTCCCAGCAGGCAGATGATTCTCGTTATTCGCCAAGAAGCTTATAAGCAACAGAATGCAACAGAGAATCGTTTCTCTGTGAACTTCCAGAAAGCAGCCAAATCCTTCAGTCTCAAGATCTCAGACTCACAGCTGGGGGATGCCGCGATGTATTTCTGTGCTTACGGTTTGTATTCAGGAGGAGGTGCTGACGGACTCACCTTTGGCAAAGGGACTCATCTAATCATCCAGCCCT TGTGCTTACGGTTTGTATTCAGGAGGAGGTGCTGACGGACTCACCTTT 30, 34, 8, 0, CGGTT F T AQTVTQSQPEMSVQEAETVTLSCTYDTSESDYYLFWYKQPPSRQMILVIRQEAYKQQNATENRFSVNFQKAAKSFSLKISDSQLGDAAMYFCAYGLYSGGGADGLTFGKGTHLIIQP TSESDYY QEAYKQQN T F T T 1 +25 TRAV2 TRAJ36 CAVEEGVGTGANNLFF 1 AAGGACCAAGTGTTTCAGCCTTCCACAGTGGCATCTTCAGAGGGAGCTGTGGTGGAAATCTTCTGTAATCACTCTGTGTCCAATGCTTACAACTTCTTCTGGTACCTTCACTTCCCGGGATGTGCACCAAGACTCCTTGTTAAAGGCTCAAAGCCTTCTCAGCAGGGACGATACAACATGACCTATGAACGGTTCTCTTCATCGCTGCTCATCCTCCAGGTGCGGGAGGCAGATGCTGCTGTTTACTACTGTGCTGTGGAGGAGGGGGTAGGAACTGGGGCAAACAACCTCTTCTTTGGGACTGGAACGAGACTCACCGTTATTCCCT TGTGCTGTGGAGGAGGGGGTAGGAACTGGGGCAAACAACCTCTTCTTT 13, 24, 0, 3, GGGGGTAGG F T KDQVFQPSTVASSEGAVVEIFCNHSVSNAYNFFWYLHFPGCAPRLLVKGSKPSQQGRYNMTYERFSSSLLILQVREADAAVYYCAVEEGVGTGANNLFFGTGTRLTVIP VSNAYN GSKP T F T T 1 +26 TRAV29/DV5 TRAJ52 CAASSGRAGGTSYGKLTF 1 AGTCAACAGAAGAATGATGACCAGCAAGTTAAGCAAAATTCACCATCCCTGAGCGTCCAGGAAGGAAGAATTTCTATTCTGAACTGTGACTATACTAACAGCATGTTTGATTATTTCCTATGGTACAAAAAATACCCTGCTGAAGGTCCTACATTCCTGATATCTATAAGTTCCATTAAGGATAAAAATGAAGATGGAAGATTCACTGTCTTCTTAAACAAAAGTGCCAAGCACCTCTCTCTGCACATTGTGCCCTCCCAGCCTGGAGACTCTGCAGTGTACTTCTGTGCAGCAAGCTCGGGACGCGCTGGTGGTACTAGCTATGGAAAGCTGACATTTGGACAAGGGACCATCTTGACTGTCCATCCAA TGTGCAGCAAGCTCGGGACGCGCTGGTGGTACTAGCTATGGAAAGCTGACATTT 23, 41, 1, 5, TCGGGACGC F T SQQKNDDQQVKQNSPSLSVQEGRISILNCDYTNSMFDYFLWYKKYPAEGPTFLISISSIKDKNEDGRFTVFLNKSAKHLSLHIVPSQPGDSAVYFCAASSGRAGGTSYGKLTFGQGTILTVHP NSMFDY ISSIKDK T F T T 1 diff --git a/examples/expected_output/dcr_TINY_1_beta.freq b/examples/expected_output/dcr_TINY_1_beta.freq new file mode 100644 index 0000000..fa04d9a --- /dev/null +++ b/examples/expected_output/dcr_TINY_1_beta.freq @@ -0,0 +1,36 @@ +15, 10, 4, 1, CTACCCCCGCGGGGAC, 2, 3 +15, 10, 4, 1, CTACCCCCGCAAAGAC, 1, 1 +43, 0, 5, 6, GGAGGGACAG, 1, 2 +15, 6, 3, 9, CCTAGCGGAATACTCCTACAC, 1, 1 +9, 6, 4, 0, CTCACGGGGGGTT, 1, 1 +0, 7, 6, 8, CCGGGACTAGCGGGTGA, 1, 1 +26, 12, 0, 0, AGCCAGGGGT, 1, 1 +24, 8, 6, 2, TCCTAGCGGCACCTTTTC, 1, 1 +20, 12, 2, 0, AGGGGCAG, 1, 1 +27, 11, 4, 2, ACCGATG, 1, 1 +31, 12, 7, 3, CCGAGACGACTAGCGGGAG, 1, 1 +14, 10, 7, 5, AACGGACGGGAGGGAAGCCGAGC, 1, 1 +22, 0, 5, 4, GAACAGGGGGC, 1, 1 +4, 6, 0, 4, TGGGCAGCGGGTGTTTAAT, 1, 1 +37, 3, 2, 12, ATCTGGGGGACAGCAAAGG, 1, 1 +20, 6, 1, 10, CCGCAGGGGGCCGC, 1, 1 +43, 8, 3, 0, GGGGAGGACTAGC, 1, 2 +20, 11, 4, 2, CTCTCCAGG, 1, 2 +15, 12, 5, 3, CGACTAGCGGGAGA, 1, 1 +33, 6, 5, 2, CCCGGGACAGGAG, 1, 1 +44, 3, 4, 0, CCCCGCGCGG, 1, 1 +29, 6, 0, 5, GTCTAGCGGGAC, 1, 1 +14, 7, 3, 9, GGTGACTAGCC, 1, 1 +26, 6, 4, 7, CCAGGACTAGCACCGAGG, 1, 1 +15, 4, 4, 0, ATCTAGGGGCCGA, 1, 1 +9, 10, 5, 2, TCGACCGGGGCG, 1, 1 +21, 12, 9, 0, TACCCTTAGCGGGAG, 1, 1 +43, 10, 1, 2, GCGCCATGCAA, 1, 1 +44, 0, 2, 1, CAGGGAATC, 1, 1 +2, 10, 0, 1, GAGGTTTGGACACAC, 1, 1 +31, 12, 3, 3, GAGCGGGAGGGT, 1, 1 +19, 6, 9, 3, TGGACTAGCGGA, 1, 1 +6, 6, 1, 7, GTCTGCTT, 1, 1 +19, 12, 5, 6, GAACTTACGGGGCGCGG, 1, 1 +40, 6, 5, 8, CAAACCATGAGTGGCTCACAGGCCC, 1, 1 +17, 10, 1, 5, GGGGAGGGGGCCAG, 1, 1 diff --git a/examples/expected_output/dcr_TINY_1_beta.n12 b/examples/expected_output/dcr_TINY_1_beta.n12 new file mode 100644 index 0000000..b5ba10d --- /dev/null +++ b/examples/expected_output/dcr_TINY_1_beta.n12 @@ -0,0 +1,48 @@ +15, 10, 4, 1, CTACCCCCGCGGGGAC, A00261:687:HMMCCDSX5:3:1101:8684:1125, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCGGGGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFF:FFFF:FFFFFF:,FFF:FFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGCTAAGTCGTGATTTTTAT, FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +15, 10, 4, 1, CTACCCCCGCGGGGAC, A00261:687:HMMCCDSX5:4:1101:8684:1125, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCGGGGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFF:FFFF:FFFFFF:,FFF:FFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGCTAAGTCGTGATTTTTAT, FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +15, 10, 4, 1, CTACCCCCGCGGAGAC, A00261:687:HMMCCDSX5:5:1101:8684:1125, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCGGAGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFF:FFFF:FFFFFF:,FFF:FFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGCCAAGTCGTGATTTTTAT, FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +15, 10, 4, 1, CTACCCCCGCGAAGAC, A00261:687:HMMCCDSX5:5:1101:8684:1125, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCGAAGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFF:FFFF:FFFFFF:,FFF:FFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGCACAGTCGTGATTTTTAT, FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +15, 10, 4, 1, CTACCCCCGCAAAGAC, A00261:687:HMMCCDSX5:5:1101:8684:1125, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCAAAGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFF:FFFF:FFFFFF:,FFF:FFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGCGGGGTCGTGATTTTTAT, FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +15, 10, 4, 1, CTACCCCCGCGGGGAC, A00261:687:HMMCCDSX5:6:1101:8684:1125, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCGGGGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFF:FFFF:FFFFFF:,FFF:FFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCGATAGTCGTGATTTTTAT, FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +15, 10, 3, 1, CGGGGACCTACCCCC, A00261:687:HMMCCDSX5:7:1101:8684:1125, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTAGCGGGGACCTACCCCCCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFF:FFFF:FFFFFF:,FFF:FFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGCTAAGTCGTGATTTTTAT, FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +43, 0, 5, 6, GGAGGGACAG, A00261:687:HMMCCDSX5:3:1101:29604:1125, GAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCGGAGGGACAGCTGAAGCTTTCTTTGGACAAGGCACCAGACTCAC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCATTCTGTCGTGATTAAAAT, FFFFFFFFFFFFFFFFF:FFFF:FFFFFFF:FF:F:FFFFFF +15, 6, 3, 9, CCTAGCGGAATACTCCTACAC, A00261:687:HMMCCDSX5:3:1101:17300:1141, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTAGCCTAGCGGAATACTCCTACACTGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFF:FF:FFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FF::FFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGGTTGCCGTCGTGATATTAAT, FFFFFFFF:FFFFFFFFFFFFFF:FF::FFFF,F,FFFFFFF +9, 6, 4, 0, CTCACGGGGGGTT, A00261:687:HMMCCDSX5:3:1101:26793:1141, TGGAGGATTCTGGAGTTTATTTCTGTGCCAGCAGCCCTCACGGGGGGTTCTCCTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGTAATGTCGTGATTTCAAC, FFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFF: +0, 7, 6, 8, CCGGGACTAGCGGGTGA, A00261:687:HMMCCDSX5:3:1101:28854:1141, CCTCCTCCCCGACATCTGTATATTTCTGCGCCAGCAGCCGGGACTAGCGGGTGAGGGGAGCTGTTTTTTGGAGAAGGCTCTAGGCT, FFFFF,FF:,,FF:,,,,FF:FF,,:,FFFFF:FFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCTTTCGTCGTGATATCCAT, FFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFF +26, 12, 0, 0, AGCCAGGGGT, A00261:687:HMMCCDSX5:3:1101:22480:1157, CTGGGGGACTCGGCCCTTTATCTTTGCGCCAGCAGCTTGGAGCCAGGGGTCTCCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCCGTTGTCGTGATCCTTAT, :FFFFFFFFFFFFFFFFFFFFFFFFF::FFFFFF:FFFFFFF +24, 8, 6, 2, TCCTAGCGGCACCTTTTC, A00261:687:HMMCCDSX5:3:1101:28926:1172, AGAAGACTCGGCCCTGTATCTCTGTGCCAGCAGTCCTAGCGGCACCTTTTCCACAGATACGCAGTATTTTGGCCCAGGCACCCGGCTGA, FFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTGGGTGTCGTGATCCTTTG, FFFFFFFFFFF,FFFFFFFF,FFF:FFFFFFFFFFF:F:FFF +15, 3, 3, 7, GGGTTTTGGGGCAGTTAG, A00261:687:HMMCCDSX5:3:1101:23737:1266, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTAGGGGTTTTGGGGCAGTTAGTGAAAAACTGTTTTTTGGCAGTGGAACCCAGCT, ::FFFFF,F:F:FF:FF:FF,FF:FFFFFFF:FF::FFFFFF:FFFFFF:FF:FFFF::::FF::F::F,FFFF,F,:F:FFFFFFFF, GTCGTGACTGGGAAAACCCTGGTAGTCAGTCGTGATTGTCAT, F,FF,FFF:FFFF,FFFFF,FF:,FF,:FFFFFFF,,F:FF: +20, 12, 2, 0, AGGGGCAG, A00261:687:HMMCCDSX5:3:1101:29966:1282, GCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGAAAGGGGCAGCTCCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCCTATGTCGTGATCTTTAC, FFFF:FFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF +27, 11, 4, 2, ACCGATG, A00261:687:HMMCCDSX5:3:1101:20627:1297, CTGGACGACTCGGCCCTGTATCTCTGTGCCAGCAGCACCGATGCTGGGGCCAACGTCCTGACTTTCGGGGCCGGCAGCAGG, FFFFFFFFFFFFFFF:FFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCTCGCCGTCGTGATGTTTAT, :FFFFFFFFFFFFFF:FFFFFFFFFFFFF,FF:FFFF:F:FF +43, 0, 5, 4, AGGGGGAC, A00261:687:HMMCCDSX5:3:1101:24641:1391, GAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCAGGGGGACCACTGAAGCTTTCTATGGACAAGGCACCAGACTCAC, F:FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFF:FFF,FF::FFFF,FF:FF:F:F:,FFF,FFFFFFFFFF,FFF:F:, GTCGTGACTGGGAAAACCCTGGTCCCCAGTCGTGATTTTGTT, FFFFFFFFFFF:FFFFF:F:FF,F:FF,FFFFFF,:FFFFFF +31, 12, 7, 3, CCGAGACGACTAGCGGGAG, A00261:687:HMMCCDSX5:3:1101:5638:1423, GAGTCGGCTGCTCCCTCCCAAACATCTGTGTACTTCTGTGCCAGCACCGAGACGACTAGCGGGAGCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTTCGGTCGTGATATGAAAG, FFF:FFF:FFF:FFFFF,FFFFFFFFFFFFFFFF,FFFFFFF +14, 10, 7, 5, AACGGACGGGAGGGAAGCCGAGC, A00261:687:HMMCCDSX5:3:1101:24713:1423, TGGAGGACTCAGCCATGTACTTCTGTGCCAGCAAACGGACGGGAGGGAAGCCGAGCGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFFFFFFFF:F:FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCTGAGGTCGTGATCGATTG, FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFFFFF:FFF +22, 0, 5, 4, GAACAGGGGGC, A00261:687:HMMCCDSX5:3:1101:24379:1438, TTCTCAGTGACTCTGGCTTCTATCTCTGTGCCTGGGAACAGGGGGCCACTGAAGCTTTCTTTGGACAAGGCACCAGACTCAC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTGTACTGTCGTGATTCAGAT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +4, 6, 0, 4, TGGGCAGCGGGTGTTTAAT, A00261:687:HMMCCDSX5:3:1101:8133:1454, AAGCTTGAGGACTCGGCCGTGTATCTCTGTGCCAGCAGCTTAGATGGGCAGCGGGTGTTTAATTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGAAAACCCTGGCCTTAAGTCGTGATTCATTGG, FFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFF,FFFFFFFF +37, 3, 2, 12, ATCTGGGGGACAGCAAAGG, A00261:687:HMMCCDSX5:3:1101:4869:1470, AGGAGGACTCGGCCGTGTATCTCTGTGCCAGCAGCTTAATCTGGGGGACAGCAAAGGAACTGTTTTTTGGCAGTGGAACCCAGCT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF:FFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCTTTGAGTCGTGATCATAAG, FF:FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFF:FFFFFF +20, 6, 1, 10, CCGCAGGGGGCCGC, A00261:687:HMMCCDSX5:3:1101:6777:1485, GCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGAAGCCGCAGGGGGCCGCGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGGTTATAGTCGTGATCAAATC, FFFFFF:FFF,FFF:FFFFFFFF:F:F::FFF:F,FFFF:FF +43, 8, 3, 0, GGGGAGGACTAGC, A00261:687:HMMCCDSX5:3:1101:31331:1485, GAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCTTGGGGAGGACTAGCAGCACAGATACGCAGTATTTTGGCCCAGGCACCCGGCTGA, FFFFFFFFFFFFFFFFFFFFF:,FFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCAAGGGTCGTGATCTCTAT, FFF:FFFFFFFFFFFFFFFFFFFF:FFFFF:FF:FFF:FF:F +20, 11, 4, 2, CTCTCCAGG, A00261:687:HMMCCDSX5:3:1101:18629:1501, GCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGCTCTCCAGGCTGGGGCCAACGTCCTGACTTTCGGGGCCGGCAGCAGG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCGTCGAGTCGTGATGTCATG, FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF::,:F +43, 8, 3, 0, GGGGAGGACTAGC, A00261:687:HMMCCDSX5:3:1101:31448:1501, GAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCTTGGGGAGGACTAGCAGCACAGATACGCAGTATTTTGGCCCAGGCACCCGGCTGA, FFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCAAGGGTCGTGATCTCTAT, FFFFFFFFFFFFFFFFFF:FFFFFFFF:FF:F,FFFFFF:,F +15, 12, 5, 3, CGACTAGCGGGAGA, A00261:687:HMMCCDSX5:3:1101:1768:1517, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTCGACTAGCGGGAGACTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCGTTGAGTCGTGATGGAGTG, FFFFFFFFFFFFFFFFFF:FFFFFFFF,FFFFFFFFFFFFF: +33, 6, 5, 2, CCCGGGACAGGAG, A00261:687:HMMCCDSX5:3:1101:9607:1532, CTGTCGGCTGCTCCCTCCCAGACATCTGTGTACTTCTGTGCCAGCAGTCCCGGGACAGGAGCCTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCATCTTGTCGTGATACCATA, FFFFFFFFFFFFFFFFFF:FFFFFFF,FFFFFFFFFFFFFFF +44, 3, 4, 0, CCCCGCGCGG, A00261:687:HMMCCDSX5:3:1101:30120:1548, CTGGGGGACTCAGCTTTGTATTTCTGTGCCAGCAGCCCCCGCGCGGCAACTAATGAAAAACTGTTTTTTGGCAGTGGAACCCAGCT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCCAAAGTCGTGATCTCTAG, FFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFF +29, 6, 0, 5, GTCTAGCGGGAC, A00261:687:HMMCCDSX5:3:1101:7238:1564, GGGACTCGGCCCTCTATCTCTGTGCCAGCAGCTTGGGTCTAGCGGGACACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFFFF:FFFFFFFFF,FFFFFFFFFFFFFF:FFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCCTTTGTCGTATACCTACG, FFFF:FFFFFFFFF:FFFFFFFFFFFF:FFFFFFFFFF:FFF +14, 7, 3, 9, GGTGACTAGCC, A00261:687:HMMCCDSX5:3:1101:24777:1564, TGGAGGACTCAGCCATGTACTTCTGTGCCAGCAGTGAGGTGACTAGCCGGGAGCTGTTTTTTGGAGAAGGCTCTAGGCT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTCACTGTCGTGGATTCCTG, ,FFFFFFFF:FFFFF:FFFFFFFFFF:FFFFFF,FFFF:FFF +20, 11, 4, 2, CTCTCCAGG, A00261:687:HMMCCDSX5:3:1101:18484:1626, GCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGCTCTCCAGGCTGGGGCCAACGTCCTGACTTTCGGGGCCGGCAGCAGG, FFFFFFFFFFFFFFFFFFFF,:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCGTCGAGTCGTGATGTCATG, FF,FFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFF: +26, 6, 4, 7, CCAGGACTAGCACCGAGG, A00261:687:HMMCCDSX5:3:1101:11966:1642, CTGGGGGACTCGGCCCTTTATCTTTGCGCCAGCAGCCCAGGACTAGCACCGAGGAATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTCCCCTGTCGTGATGAAGAT, FFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:F:FF,FFFFFF +15, 4, 4, 0, ATCTAGGGGCCGA, A00261:687:HMMCCDSX5:3:1101:7283:1705, ATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTAATCTAGGGGCCGATAGCAATCAGCCCCAGCATTTTGGTGATGGGACTCGACTC, FFFFFFFF:FFFFFFFFFFFFFF:FFFFFFF:FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGAATCTCGTCGTGATGGAGGC, :FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF +9, 10, 5, 2, TCGACCGGGGCG, A00261:687:HMMCCDSX5:3:1101:24975:1720, TGGAGGATTCTGGAGTTTATTTCTGTGCCAGCAGCTCGACCGGGGCGCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, F,FFF:FFF:FFFFFFFFFFFF:FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF::FFFFFFFFFFFFFF:FFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTATTACGTCGTGATTTGCTT, FFFFFFFFFFFFFF::FFFFFFF:FFF::FFFF:FFFFF,FF +21, 12, 9, 0, TACCCTTAGCGGGAG, A00261:687:HMMCCDSX5:3:1101:26115:1783, TTGGTGACTCTGCTGTGTATTTCTGTGCCAGTACCCTTAGCGGGAGCTCCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF, GTCGTGACTGGGAAAACCCTGGCGTCATGTCGTGATTAAATG, FFFFFFFF,:FFFFFFFFFFF:FFFFF:FFFFFFFFFF:FFF +43, 10, 1, 2, GCGCCATGCAA, A00261:687:HMMCCDSX5:3:1101:26874:1783, GAGCAGGGGGACTCGGTCATGTATCTCTGTGCCAGCAGCTTAGGCGCCATGCAACAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF, GTCATGACTGGGAAAACCCTGGCGTGGTCGTGATGTACAACC, FFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, +44, 0, 2, 1, CAGGGAATC, A00261:687:HMMCCDSX5:3:1101:7862:1799, CTGGGGGACTCAGCTTTGTATTTCTGTGCCAGCAGCGTCAGGGAATCGAACACTGAAGCTTTCTTTGGACAAGGCACCAGACTCAC, FFFFFFFF:FFFFF:FFFFFFFFFFFFFFFFF:FFFFFF:FFFFFF:FF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F, GTCGTGACTGGGAAAACCCTGGCTCGTTGTCGTGATCGGAAC, FFF::FFFFFFFFFFFF,FFFFF:F::FFFFFFF:FFFFFFF +2, 10, 0, 1, GAGGTTTGGACACAC, A00261:687:HMMCCDSX5:3:1101:9579:1861, CCAGCTCCCAGACATCTGTGTACTTCTGTGCCATCAGTGAGTCGAGGTTTGGACACACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFFFFFFFFFFFFF,FFFF,FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTAACCGGTCGTGATGGAGAG, F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF +31, 12, 3, 3, GAGCGGGAGGGT, A00261:687:HMMCCDSX5:3:1101:16396:1861, GAGTCGGCTGCTCCCTCCCAAACATCTGTGTACTTCTGTGCCAGCAGTGAGAGCGGGAGGGTCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTTGACGTCGTGATTGAATG, :FFFF:FFFFF:FFFFFFFF::FFFFFFFFFFF:FFFFFFFF +19, 6, 9, 3, TGGACTAGCGGA, A00261:687:HMMCCDSX5:3:1101:19027:1877, CCAACCAGACATCTATGTACCTCTGTGCCAGTGGACTAGCGGACTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFF:FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCACTCGTCGTGATTCTCTC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F +33, 12, 0, 0, GAGAGGCGTCGTAAG, A00261:687:HMMCCDSX5:3:1101:23095:1877, CTGTCGGCTGCTCCCTCCCAGACATCTGTGTACTTCTGTGCCAGCAGTTACTCGAGAGGCGTCGTAAGCTCCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FF:FFF:FF:FF:FFFFF:FFF,F::FFFFF:FFFFFFF:FFFFFFFFFFFF:FFF:F:F:FFFFFF:FF:FFFFFFFF:FFFFFF,FFF:FFFFFFFF:::FFF::F, GTCGTGACTGGGAAAACCCTGGGTTGTGGTGGCTCCCCTTTC, FFFF:,FF:FFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFF +6, 6, 1, 7, GTCTGCTT, A00261:687:HMMCCDSX5:3:1101:2546:1892, CCAGGGACTCAGCTGTGTACTTCTGTGCCAGCAGTTTAGGTCTGCTTAATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGGTGAAAGTCGTGATGTTAAT, :FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF +19, 12, 5, 6, GAACTTACGGGGCGCGG, A00261:687:HMMCCDSX5:3:1101:20437:1908, CCAACCAGACATCTATGTACCTCTGTGCCAGCAGTGAACTTACGGGGCGCGGCGAGCAGTACTTCGGGCCGGGCACCAGGCTCACG, FFFFFFFFFFFF:FFFFFFFFFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCTTCGCGTCGTGATTTATGG, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +40, 6, 5, 8, CAAACCATGAGTGGCTCACAGGCCC, A00261:687:HMMCCDSX5:3:1101:5394:1971, GAGCAGCGGGACTCGGCCATGTATCGCTGTGCCAGCAGCCAAACCATGAGTGGCTCACAGGCCCATGAGCAGTTCTTCGGGCCAGGGACACGGCTC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F, GTCGTGACTGGGAAAACCCTGGCCCCCTGTCGTGATCAAATT, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF +17, 10, 1, 5, GGGGAGGGGGCCAG, A00261:687:HMMCCDSX5:3:1101:24542:1971, CCTCACATACCTCTCAGTACCTCTGTGCCAGCAGTGAATGGGGAGGGGGCCAGGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFFFF:F,FFF,F,:FFFFFFFFFFFFFFFFFFFF:FFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCCACATGTCGTGATTTCTTT, FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF: +43, 0, 5, 6, GGAGGGACAG, A00261:687:HMMCCDSX5:3:1101:30942:2002, GAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCGGAGGGACAGCTGAAGCTTTCTTTGGACAAGGCACCAGACTCAC, FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGCATTCTGTCGTGATTAAAAT, F,:FF:FFFFFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFF +20, 10, 3, 9, GTCCGACGTGGGACGGG, A00261:687:HMMCCDSX5:3:1101:26151:2033, GCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGAGTCCGACGTGGGACGGGCCCAGTACTTCGGGCCAGGCACGCGGCTCCT, FFFFFFFFFFFFFFFFF,FFFFFFF:FFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFF, GTCGTGACTGGGAAAACCCTGGTTTTATGTCGTGATTTTAAC, FFFFFFFFFFFFF:FFFFFFFFF,F:,FFFFFFFFFFFFFFF diff --git a/examples/expected_output/dcr_TINY_1_beta.tsv b/examples/expected_output/dcr_TINY_1_beta.tsv new file mode 100644 index 0000000..88e78a7 --- /dev/null +++ b/examples/expected_output/dcr_TINY_1_beta.tsv @@ -0,0 +1,37 @@ +sequence_id v_call d_call j_call junction_aa duplicate_count sequence junction decombinator_id rev_comp productive sequence_aa cdr1_aa cdr2_aa vj_in_frame stop_codon conserved_c conserved_f sequence_alignment germline_alignment v_cigar d_cigar j_cigar av_UMI_cluster_size +1 TRBV20-1 TRBJ2-5 CSATTPAGTQETQYF 2 GGTGCTGTCGTCTCTCAACATCCGAGCTGGGTTATCTGTAAGAGTGGAACCTCTGTGAAGATCGAGTGCCGTTCCCTGGACTTTCAGGCCACAACTATGTTTTGGTATCGTCAGTTCCCGAAACAGAGTCTCATGCTGATGGCAACTTCCAATGAGGGCTCCAAGGCCACATACGAGCAAGGCGTCGAGAAGGACAAGTTTCTCATCAACCATGCAAGCCTGACCTTGTCCACTCTGACAGTGACCAGTGCCCATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCGGGGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCTGGTGCTCG TGCAGTGCTACTACCCCCGCGGGGACCCAAGAGACCCAGTACTTC 15, 10, 4, 1, CTACCCCCGCGGGGAC F T GAVVSQHPSWVICKSGTSVKIECRSLDFQATTMFWYRQFPKQSLMLMATSNEGSKATYEQGVEKDKFLINHASLTLSTLTVTSAHPEDSSFYICSATTPAGTQETQYFGPGTRLLVL DFQATT SNEGSKA T F T T 3 +2 TRBV20-1 TRBJ2-5 CSATTPAKTQETQYF 1 GGTGCTGTCGTCTCTCAACATCCGAGCTGGGTTATCTGTAAGAGTGGAACCTCTGTGAAGATCGAGTGCCGTTCCCTGGACTTTCAGGCCACAACTATGTTTTGGTATCGTCAGTTCCCGAAACAGAGTCTCATGCTGATGGCAACTTCCAATGAGGGCTCCAAGGCCACATACGAGCAAGGCGTCGAGAAGGACAAGTTTCTCATCAACCATGCAAGCCTGACCTTGTCCACTCTGACAGTGACCAGTGCCCATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTACTACCCCCGCAAAGACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCTGGTGCTCG TGCAGTGCTACTACCCCCGCAAAGACCCAAGAGACCCAGTACTTC 15, 10, 4, 1, CTACCCCCGCAAAGAC F T GAVVSQHPSWVICKSGTSVKIECRSLDFQATTMFWYRQFPKQSLMLMATSNEGSKATYEQGVEKDKFLINHASLTLSTLTVTSAHPEDSSFYICSATTPAKTQETQYFGPGTRLLVL DFQATT SNEGSKA T F T T 1 +3 TRBV7-9 TRBJ1-1 CASSGGTAEAFF 1 GATACTGGAGTCTCCCAGAACCCCAGACACAAGATCACAAAGAGGGGACAGAATGTAACTTTCAGGTGTGATCCAATTTCTGAACACAACCGCCTTTATTGGTACCGACAGACCCTGGGGCAGGGCCCAGAGTTTCTGACTTACTTCCAGAATGAAGCTCAACTAGAAAAATCAAGGCTGCTCAGTGATCGGTTCTCTGCAGAGAGGCCTAAGGGATCTTTCTCCACCTTGGAGATCCAGCGCACAGAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCGGAGGGACAGCTGAAGCTTTCTTTGGACAAGGCACCAGACTCACAGTTGTAG TGTGCCAGCAGCGGAGGGACAGCTGAAGCTTTCTTT 43, 0, 5, 6, GGAGGGACAG F T DTGVSQNPRHKITKRGQNVTFRCDPISEHNRLYWYRQTLGQGPEFLTYFQNEAQLEKSRLLSDRFSAERPKGSFSTLEIQRTEQGDSAMYLCASSGGTAEAFFGQGTRLTVV SEHNR FQNEAQ T F T T 2 +4 TRBV20-1 TRBJ2-1 CSASLAEYSYTEQFF 1 GGTGCTGTCGTCTCTCAACATCCGAGCTGGGTTATCTGTAAGAGTGGAACCTCTGTGAAGATCGAGTGCCGTTCCCTGGACTTTCAGGCCACAACTATGTTTTGGTATCGTCAGTTCCCGAAACAGAGTCTCATGCTGATGGCAACTTCCAATGAGGGCTCCAAGGCCACATACGAGCAAGGCGTCGAGAAGGACAAGTTTCTCATCAACCATGCAAGCCTGACCTTGTCCACTCTGACAGTGACCAGTGCCCATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTAGCCTAGCGGAATACTCCTACACTGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGCAGTGCTAGCCTAGCGGAATACTCCTACACTGAGCAGTTCTTC 15, 6, 3, 9, CCTAGCGGAATACTCCTACAC F T GAVVSQHPSWVICKSGTSVKIECRSLDFQATTMFWYRQFPKQSLMLMATSNEGSKATYEQGVEKDKFLINHASLTLSTLTVTSAHPEDSSFYICSASLAEYSYTEQFFGPGTRLTVL DFQATT SNEGSKA T F T T 1 +5 TRBV14 TRBJ2-1 CASSPHGGFSYNEQFF 1 GAAGCTGGAGTTACTCAGTTCCCCAGCCACAGCGTAATAGAGAAGGGCCAGACTGTGACTCTGAGATGTGACCCAATTTCTGGACATGATAATCTTTATTGGTATCGACGTGTTATGGGAAAAGAAATAAAATTTCTGTTACATTTTGTGAAAGAGTCTAAACAGGATGAGTCCGGTATGCCCAACAATCGATTCTTAGCTGAAAGGACTGGAGGGACGTATTCTACTCTGAAGGTGCAGCCTGCAGAACTGGAGGATTCTGGAGTTTATTTCTGTGCCAGCAGCCCTCACGGGGGGTTCTCCTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGTGCCAGCAGCCCTCACGGGGGGTTCTCCTACAATGAGCAGTTCTTC 9, 6, 4, 0, CTCACGGGGGGTT F T EAGVTQFPSHSVIEKGQTVTLRCDPISGHDNLYWYRRVMGKEIKFLLHFVKESKQDESGMPNNRFLAERTGGTYSTLKVQPAELEDSGVYFCASSPHGGFSYNEQFFGPGTRLTVL SGHDN FVKESK T F T T 1 +6 TRBV10-1 TRBJ2-2 1 GATGCTGAAATCACCCAGAGCCCAAGACACAAGATCACAGAGACAGGAAGGCAGGTGACCTTGGCGTGTCACCAGACTTGGAACCACAACAATATGTTCTGGTATCGACAAGACCTGGGACATGGGCTGAGGCTGATCCATTACTCATATGGTGTTCAAGACACTAACAAAGGAGAAGTCTCAGATGGCTACAGTGTCTCTAGATCAAACACAGAGGACCTCCCCCTCACTCTGGAGTCTGCTGCCTCCTCCCAGACATCTGTATATTTCTGCGCCAGCAGCCGGGACTAGCGGGTGAGGGGAGCTGTTTTTTGGAGAAGGCTCTAGGCTGACCGTACTGG 0, 7, 6, 8, CCGGGACTAGCGGGTGA F F DAEITQSPRHKITETGRQVTLACHQTWNHNNMFWYRQDLGHGLRLIHYSYGVQDTNKGEVSDGYSVSRSNTEDLPLTLESAASSQTSVYFCASSRD*RVRGAVFWRRL*ADRT F T T F 1 +7 TRBV5-1 TRBJ2-7 CASSLEPGVSYEQYF 1 AAGGCTGGAGTCACTCAAACTCCAAGATATCTGATCAAAACGAGAGGACAGCAAGTGACACTGAGCTGCTCCCCTATCTCTGGGCATAGGAGTGTATCCTGGTACCAACAGACCCCAGGACAGGGCCTTCAGTTCCTCTTTGAATACTTCAGTGAGACACAGAGAAACAAAGGAAACTTCCCTGGTCGATTCTCAGGGCGCCAGTTCTCTAACTCTCGCTCTGAGATGAATGTGAGCACCTTGGAGCTGGGGGACTCGGCCCTTTATCTTTGCGCCAGCAGCTTGGAGCCAGGGGTCTCCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACGGTCACAG TGCGCCAGCAGCTTGGAGCCAGGGGTCTCCTACGAGCAGTACTTC 26, 12, 0, 0, AGCCAGGGGT F T KAGVTQTPRYLIKTRGQQVTLSCSPISGHRSVSWYQQTPGQGLQFLFEYFSETQRNKGNFPGRFSGRQFSNSRSEMNVSTLELGDSALYLCASSLEPGVSYEQYFGPGTRLTVT SGHRS YFSETQ T F T T 1 +8 TRBV4-2 TRBJ2-3 CASSPSGTFSTDTQYF 1 ACGGGAGTTACGCAGACACCAAGACACCTGGTCATGGGAATGACAAATAAGAAGTCTTTGAAATGTGAACAACATCTGGGGCATAACGCTATGTATTGGTACAAGCAAAGTGCTAAGAAGCCACTGGAGCTCATGTTTGTCTACAACTTTAAAGAACAGACTGAAAACAACAGTGTGCCAAGTCGCTTCTCACCTGAATGCCCCAACAGCTCTCACTTATTCCTTCACCTACACACCCTGCAGCCAGAAGACTCGGCCCTGTATCTCTGTGCCAGCAGTCCTAGCGGCACCTTTTCCACAGATACGCAGTATTTTGGCCCAGGCACCCGGCTGACAGTGCTCG TGTGCCAGCAGTCCTAGCGGCACCTTTTCCACAGATACGCAGTATTTT 24, 8, 6, 2, TCCTAGCGGCACCTTTTC F T TGVTQTPRHLVMGMTNKKSLKCEQHLGHNAMYWYKQSAKKPLELMFVYNFKEQTENNSVPSRFSPECPNSSHLFLHLHTLQPEDSALYLCASSPSGTFSTDTQYFGPGTRLTVL LGHNA YNFKEQ T F T T 1 +9 TRBV29-1 TRBJ2-7 CSVERGSSYEQYF 1 AGTGCTGTCATCTCTCAAAAGCCAAGCAGGGATATCTGTCAACGTGGAACCTCCCTGACGATCCAGTGTCAAGTCGATAGCCAAGTCACCATGATGTTCTGGTACCGTCAGCAACCTGGACAGAGCCTGACACTGATCGCAACTGCAAATCAGGGCTCTGAGGCCACATATGAGAGTGGATTTGTCATTGACAAGTTTCCCATCAGCCGCCCAAACCTAACATTCTCAACTCTGACTGTGAGCAACATGAGCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGAAAGGGGCAGCTCCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACGGTCACAG TGCAGCGTTGAAAGGGGCAGCTCCTACGAGCAGTACTTC 20, 12, 2, 0, AGGGGCAG F T SAVISQKPSRDICQRGTSLTIQCQVDSQVTMMFWYRQQPGQSLTLIATANQGSEATYESGFVIDKFPISRPNLTFSTLTVSNMSPEDSSIYLCSVERGSSYEQYFGPGTRLTVT SQVTM ANQGSEA T F T T 1 +10 TRBV5-4 TRBJ2-6 CASSTDAGANVLTF 1 GAGACTGGAGTCACCCAAAGTCCCACACACCTGATCAAAACGAGAGGACAGCAAGTGACTCTGAGATGCTCTTCTCAGTCTGGGCACAACACTGTGTCCTGGTACCAACAGGCCCTGGGTCAGGGGCCCCAGTTTATCTTTCAGTATTATAGGGAGGAAGAGAATGGCAGAGGAAACTTCCCTCCTAGATTCTCAGGTCTCCAGTTCCCTAATTATAGCTCTGAGCTGAATGTGAACGCCTTGGAGCTGGACGACTCGGCCCTGTATCTCTGTGCCAGCAGCACCGATGCTGGGGCCAACGTCCTGACTTTCGGGGCCGGCAGCAGGCTGACCGTGCTGG TGTGCCAGCAGCACCGATGCTGGGGCCAACGTCCTGACTTTC 27, 11, 4, 2, ACCGATG F T ETGVTQSPTHLIKTRGQQVTLRCSSQSGHNTVSWYQQALGQGPQFIFQYYREEENGRGNFPPRFSGLQFPNYSSELNVNALELDDSALYLCASSTDAGANVLTFGAGSRLTVL SGHNT YYREEE T F T T 1 +11 TRBV6-1 TRBJ2-7 CASTETTSGSYEQYF 1 AATGCTGGTGTCACTCAGACCCCAAAATTCCAGGTCCTGAAGACAGGACAGAGCATGACACTGCAGTGTGCCCAGGATATGAACCATAACTCCATGTACTGGTATCGACAAGACCCAGGCATGGGACTGAGGCTGATTTATTACTCAGCTTCTGAGGGTACCACTGACAAAGGAGAAGTCCCCAATGGCTACAATGTCTCCAGATTAAACAAACGGGAGTTCTCGCTCAGGCTGGAGTCGGCTGCTCCCTCCCAGACATCTGTGTACTTCTGTGCCAGCACCGAGACGACTAGCGGGAGCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACGGTCACAG TGTGCCAGCACCGAGACGACTAGCGGGAGCTACGAGCAGTACTTC 31, 12, 7, 3, CCGAGACGACTAGCGGGAG F T NAGVTQTPKFQVLKTGQSMTLQCAQDMNHNSMYWYRQDPGMGLRLIYYSASEGTTDKGEVPNGYNVSRLNKREFSLRLESAAPSQTSVYFCASTETTSGSYEQYFGPGTRLTVT MNHNS SASEGT T F T T 1 +12 TRBV2 TRBJ2-5 CASKRTGGKPSETQYF 1 GAACCTGAAGTCACCCAGACTCCCAGCCATCAGGTCACACAGATGGGACAGGAAGTGATCTTGCGCTGTGTCCCCATCTCTAATCACTTATACTTCTATTGGTACAGACAAATCTTGGGGCAGAAAGTCGAGTTTCTGGTTTCCTTTTATAATAATGAAATCTCAGAGAAGTCTGAAATATTCGATGATCAATTCTCAGTTGAAAGGCCTGATGGATCAAATTTCACTCTGAAGATCCGGTCCACAAAGCTGGAGGACTCAGCCATGTACTTCTGTGCCAGCAAACGGACGGGAGGGAAGCCGAGCGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCTGGTGCTCG TGTGCCAGCAAACGGACGGGAGGGAAGCCGAGCGAGACCCAGTACTTC 14, 10, 7, 5, AACGGACGGGAGGGAAGCCGAGC F T EPEVTQTPSHQVTQMGQEVILRCVPISNHLYFYWYRQILGQKVEFLVSFYNNEISEKSEIFDDQFSVERPDGSNFTLKIRSTKLEDSAMYFCASKRTGGKPSETQYFGPGTRLLVL SNHLY FYNNEI T F T T 1 +13 TRBV30 TRBJ1-1 CAWEQGATEAFF 1 TCTCAGACTATTCATCAATGGCCAGCGACCCTGGTGCAGCCTGTGGGCAGCCCGCTCTCTCTGGAGTGCACTGTGGAGGGAACATCAAACCCCAACCTATACTGGTACCGACAGGCTGCAGGCAGGGGCCTCCAGCTGCTCTTCTACTCCGTTGGTATTGGCCAGATCAGCTCTGAGGTGCCCCAGAATCTCTCAGCCTCCAGACCCCAGGACCGGCAGTTCATCCTGAGTTCTAAGAAGCTCCTTCTCAGTGACTCTGGCTTCTATCTCTGTGCCTGGGAACAGGGGGCCACTGAAGCTTTCTTTGGACAAGGCACCAGACTCACAGTTGTAG TGTGCCTGGGAACAGGGGGCCACTGAAGCTTTCTTT 22, 0, 5, 4, GAACAGGGGGC F T SQTIHQWPATLVQPVGSPLSLECTVEGTSNPNLYWYRQAAGRGLQLLFYSVGIGQISSEVPQNLSASRPQDRQFILSSKKLLLSDSGFYLCAWEQGATEAFFGQGTRLTVV GTSNPN SVGIG T F T T 1 +14 TRBV11-2 TRBJ2-1 CASSLDGQRVFNYNEQFF 1 GAAGCTGGAGTTGCCCAGTCTCCCAGATATAAGATTATAGAGAAAAGGCAGAGTGTGGCTTTTTGGTGCAATCCTATATCTGGCCATGCTACCCTTTACTGGTACCAGCAGATCCTGGGACAGGGCCCAAAGCTTCTGATTCAGTTTCAGAATAACGGTGTAGTGGATGATTCACAGTTGCCTAAGGATCGATTTTCTGCAGAGAGGCTCAAAGGAGTAGACTCCACTCTCAAGATCCAGCCTGCAAAGCTTGAGGACTCGGCCGTGTATCTCTGTGCCAGCAGCTTAGATGGGCAGCGGGTGTTTAATTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGTGCCAGCAGCTTAGATGGGCAGCGGGTGTTTAATTACAATGAGCAGTTCTTC 4, 6, 0, 4, TGGGCAGCGGGTGTTTAAT F T EAGVAQSPRYKIIEKRQSVAFWCNPISGHATLYWYQQILGQGPKLLIQFQNNGVVDDSQLPKDRFSAERLKGVDSTLKIQPAKLEDSAVYLCASSLDGQRVFNYNEQFFGPGTRLTVL SGHAT FQNNGV T F T T 1 +15 TRBV7-2 TRBJ1-4 CASSLIWGTAKELFF 1 GGAGCTGGAGTCTCCCAGTCCCCCAGTAACAAGGTCACAGAGAAGGGAAAGGATGTAGAGCTCAGGTGTGATCCAATTTCAGGTCATACTGCCCTTTACTGGTACCGACAGAGCCTGGGGCAGGGCCTGGAGTTTTTAATTTACTTCCAAGGCAACAGTGCACCAGACAAATCAGGGCTGCCCAGTGATCGCTTCTCTGCAGAGAGGACTGGGGGATCCGTCTCCACTCTGACGATCCAGCGCACACAGCAGGAGGACTCGGCCGTGTATCTCTGTGCCAGCAGCTTAATCTGGGGGACAGCAAAGGAACTGTTTTTTGGCAGTGGAACCCAGCTCTCTGTCTTGG TGTGCCAGCAGCTTAATCTGGGGGACAGCAAAGGAACTGTTTTTT 37, 3, 2, 12, ATCTGGGGGACAGCAAAGG F T GAGVSQSPSNKVTEKGKDVELRCDPISGHTALYWYRQSLGQGLEFLIYFQGNSAPDKSGLPSDRFSAERTGGSVSTLTIQRTQQEDSAVYLCASSLIWGTAKELFFGSGTQLSVL SGHTA FQGNSA T F T T 1 +16 TRBV29-1 TRBJ2-1 CSVEAAGGREQFF 1 AGTGCTGTCATCTCTCAAAAGCCAAGCAGGGATATCTGTCAACGTGGAACCTCCCTGACGATCCAGTGTCAAGTCGATAGCCAAGTCACCATGATGTTCTGGTACCGTCAGCAACCTGGACAGAGCCTGACACTGATCGCAACTGCAAATCAGGGCTCTGAGGCCACATATGAGAGTGGATTTGTCATTGACAAGTTTCCCATCAGCCGCCCAAACCTAACATTCTCAACTCTGACTGTGAGCAACATGAGCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGAAGCCGCAGGGGGCCGCGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGCAGCGTTGAAGCCGCAGGGGGCCGCGAGCAGTTCTTC 20, 6, 1, 10, CCGCAGGGGGCCGC F T SAVISQKPSRDICQRGTSLTIQCQVDSQVTMMFWYRQQPGQSLTLIATANQGSEATYESGFVIDKFPISRPNLTFSTLTVSNMSPEDSSIYLCSVEAAGGREQFFGPGTRLTVL SQVTM ANQGSEA T F T T 1 +17 TRBV7-9 TRBJ2-3 CASSLGRTSSTDTQYF 1 GATACTGGAGTCTCCCAGAACCCCAGACACAAGATCACAAAGAGGGGACAGAATGTAACTTTCAGGTGTGATCCAATTTCTGAACACAACCGCCTTTATTGGTACCGACAGACCCTGGGGCAGGGCCCAGAGTTTCTGACTTACTTCCAGAATGAAGCTCAACTAGAAAAATCAAGGCTGCTCAGTGATCGGTTCTCTGCAGAGAGGCCTAAGGGATCTTTCTCCACCTTGGAGATCCAGCGCACAGAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCTTGGGGAGGACTAGCAGCACAGATACGCAGTATTTTGGCCCAGGCACCCGGCTGACAGTGCTCG TGTGCCAGCAGCTTGGGGAGGACTAGCAGCACAGATACGCAGTATTTT 43, 8, 3, 0, GGGGAGGACTAGC F T DTGVSQNPRHKITKRGQNVTFRCDPISEHNRLYWYRQTLGQGPEFLTYFQNEAQLEKSRLLSDRFSAERPKGSFSTLEIQRTEQGDSAMYLCASSLGRTSSTDTQYFGPGTRLTVL SEHNR FQNEAQ T F T T 2 +18 TRBV29-1 TRBJ2-6 CSVALQAGANVLTF 1 AGTGCTGTCATCTCTCAAAAGCCAAGCAGGGATATCTGTCAACGTGGAACCTCCCTGACGATCCAGTGTCAAGTCGATAGCCAAGTCACCATGATGTTCTGGTACCGTCAGCAACCTGGACAGAGCCTGACACTGATCGCAACTGCAAATCAGGGCTCTGAGGCCACATATGAGAGTGGATTTGTCATTGACAAGTTTCCCATCAGCCGCCCAAACCTAACATTCTCAACTCTGACTGTGAGCAACATGAGCCCTGAAGACAGCAGCATATATCTCTGCAGCGTTGCTCTCCAGGCTGGGGCCAACGTCCTGACTTTCGGGGCCGGCAGCAGGCTGACCGTGCTGG TGCAGCGTTGCTCTCCAGGCTGGGGCCAACGTCCTGACTTTC 20, 11, 4, 2, CTCTCCAGG F T SAVISQKPSRDICQRGTSLTIQCQVDSQVTMMFWYRQQPGQSLTLIATANQGSEATYESGFVIDKFPISRPNLTFSTLTVSNMSPEDSSIYLCSVALQAGANVLTFGAGSRLTVL SQVTM ANQGSEA T F T T 2 +19 TRBV20-1 TRBJ2-7 CSARLAGDYEQYF 1 GGTGCTGTCGTCTCTCAACATCCGAGCTGGGTTATCTGTAAGAGTGGAACCTCTGTGAAGATCGAGTGCCGTTCCCTGGACTTTCAGGCCACAACTATGTTTTGGTATCGTCAGTTCCCGAAACAGAGTCTCATGCTGATGGCAACTTCCAATGAGGGCTCCAAGGCCACATACGAGCAAGGCGTCGAGAAGGACAAGTTTCTCATCAACCATGCAAGCCTGACCTTGTCCACTCTGACAGTGACCAGTGCCCATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTCGACTAGCGGGAGACTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACGGTCACAG TGCAGTGCTCGACTAGCGGGAGACTACGAGCAGTACTTC 15, 12, 5, 3, CGACTAGCGGGAGA F T GAVVSQHPSWVICKSGTSVKIECRSLDFQATTMFWYRQFPKQSLMLMATSNEGSKATYEQGVEKDKFLINHASLTLSTLTVTSAHPEDSSFYICSARLAGDYEQYFGPGTRLTVT DFQATT SNEGSKA T F T T 1 +20 TRBV6-5 TRBJ2-1 CASSPGTGAYNEQFF 1 AATGCTGGTGTCACTCAGACCCCAAAATTCCAGGTCCTGAAGACAGGACAGAGCATGACACTGCAGTGTGCCCAGGATATGAACCATGAATACATGTCCTGGTATCGACAAGACCCAGGCATGGGGCTGAGGCTGATTCATTACTCAGTTGGTGCTGGTATCACTGACCAAGGAGAAGTCCCCAATGGCTACAATGTCTCCAGATCAACCACAGAGGATTTCCCGCTCAGGCTGCTGTCGGCTGCTCCCTCCCAGACATCTGTGTACTTCTGTGCCAGCAGTCCCGGGACAGGAGCCTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGTGCCAGCAGTCCCGGGACAGGAGCCTACAATGAGCAGTTCTTC 33, 6, 5, 2, CCCGGGACAGGAG F T NAGVTQTPKFQVLKTGQSMTLQCAQDMNHEYMSWYRQDPGMGLRLIHYSVGAGITDQGEVPNGYNVSRSTTEDFPLRLLSAAPSQTSVYFCASSPGTGAYNEQFFGPGTRLTVL MNHEY SVGAGI T F T T 1 +21 TRBV9 TRBJ1-4 CASSPRAATNEKLFF 1 GATTCTGGAGTCACACAAACCCCAAAGCACCTGATCACAGCAACTGGACAGCGAGTGACGCTGAGATGCTCCCCTAGGTCTGGAGACCTCTCTGTGTACTGGTACCAACAGAGCCTGGACCAGGGCCTCCAGTTCCTCATTCAGTATTATAATGGAGAAGAGAGAGCAAAAGGAAACATTCTTGAACGATTCTCCGCACAACAGTTCCCTGACTTGCACTCTGAACTAAACCTGAGCTCTCTGGAGCTGGGGGACTCAGCTTTGTATTTCTGTGCCAGCAGCCCCCGCGCGGCAACTAATGAAAAACTGTTTTTTGGCAGTGGAACCCAGCTCTCTGTCTTGG TGTGCCAGCAGCCCCCGCGCGGCAACTAATGAAAAACTGTTTTTT 44, 3, 4, 0, CCCCGCGCGG F T DSGVTQTPKHLITATGQRVTLRCSPRSGDLSVYWYQQSLDQGLQFLIQYYNGEERAKGNILERFSAQQFPDLHSELNLSSLELGDSALYFCASSPRAATNEKLFFGSGTQLSVL SGDLS YYNGEE T F T T 1 +22 TRBV5-6 TRBJ2-1 CASSLGLAGHNEQFF 1 GACGCTGGAGTCACCCAAAGTCCCACACACCTGATCAAAACGAGAGGACAGCAAGTGACTCTGAGATGCTCTCCTAAGTCTGGGCATGACACTGTGTCCTGGTACCAACAGGCCCTGGGTCAGGGGCCCCAGTTTATCTTTCAGTATTATGAGGAGGAAGAGAGACAGAGAGGCAACTTCCCTGATCGATTCTCAGGTCACCAGTTCCCTAACTATAGCTCTGAGCTGAATGTGAACGCCTTGTTGCTGGGGGACTCGGCCCTCTATCTCTGTGCCAGCAGCTTGGGTCTAGCGGGACACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGTGCCAGCAGCTTGGGTCTAGCGGGACACAATGAGCAGTTCTTC 29, 6, 0, 5, GTCTAGCGGGAC F T DAGVTQSPTHLIKTRGQQVTLRCSPKSGHDTVSWYQQALGQGPQFIFQYYEEEERQRGNFPDRFSGHQFPNYSSELNVNALLLGDSALYLCASSLGLAGHNEQFFGPGTRLTVL SGHDT YYEEEE T F T T 1 +23 TRBV2 TRBJ2-2 CASSEVTSRELFF 1 GAACCTGAAGTCACCCAGACTCCCAGCCATCAGGTCACACAGATGGGACAGGAAGTGATCTTGCGCTGTGTCCCCATCTCTAATCACTTATACTTCTATTGGTACAGACAAATCTTGGGGCAGAAAGTCGAGTTTCTGGTTTCCTTTTATAATAATGAAATCTCAGAGAAGTCTGAAATATTCGATGATCAATTCTCAGTTGAAAGGCCTGATGGATCAAATTTCACTCTGAAGATCCGGTCCACAAAGCTGGAGGACTCAGCCATGTACTTCTGTGCCAGCAGTGAGGTGACTAGCCGGGAGCTGTTTTTTGGAGAAGGCTCTAGGCTGACCGTACTGG TGTGCCAGCAGTGAGGTGACTAGCCGGGAGCTGTTTTTT 14, 7, 3, 9, GGTGACTAGCC F T EPEVTQTPSHQVTQMGQEVILRCVPISNHLYFYWYRQILGQKVEFLVSFYNNEISEKSEIFDDQFSVERPDGSNFTLKIRSTKLEDSAMYFCASSEVTSRELFFGEGSRLTVL SNHLY FYNNEI T F T T 1 +24 TRBV5-1 TRBJ2-1 CASSPGLAPRNEQFF 1 AAGGCTGGAGTCACTCAAACTCCAAGATATCTGATCAAAACGAGAGGACAGCAAGTGACACTGAGCTGCTCCCCTATCTCTGGGCATAGGAGTGTATCCTGGTACCAACAGACCCCAGGACAGGGCCTTCAGTTCCTCTTTGAATACTTCAGTGAGACACAGAGAAACAAAGGAAACTTCCCTGGTCGATTCTCAGGGCGCCAGTTCTCTAACTCTCGCTCTGAGATGAATGTGAGCACCTTGGAGCTGGGGGACTCGGCCCTTTATCTTTGCGCCAGCAGCCCAGGACTAGCACCGAGGAATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGCGCCAGCAGCCCAGGACTAGCACCGAGGAATGAGCAGTTCTTC 26, 6, 4, 7, CCAGGACTAGCACCGAGG F T KAGVTQTPRYLIKTRGQQVTLSCSPISGHRSVSWYQQTPGQGLQFLFEYFSETQRNKGNFPGRFSGRQFSNSRSEMNVSTLELGDSALYLCASSPGLAPRNEQFFGPGTRLTVL SGHRS YFSETQ T F T T 1 +25 TRBV20-1 TRBJ1-5 CSANLGADSNQPQHF 1 GGTGCTGTCGTCTCTCAACATCCGAGCTGGGTTATCTGTAAGAGTGGAACCTCTGTGAAGATCGAGTGCCGTTCCCTGGACTTTCAGGCCACAACTATGTTTTGGTATCGTCAGTTCCCGAAACAGAGTCTCATGCTGATGGCAACTTCCAATGAGGGCTCCAAGGCCACATACGAGCAAGGCGTCGAGAAGGACAAGTTTCTCATCAACCATGCAAGCCTGACCTTGTCCACTCTGACAGTGACCAGTGCCCATCCTGAAGACAGCAGCTTCTACATCTGCAGTGCTAATCTAGGGGCCGATAGCAATCAGCCCCAGCATTTTGGTGATGGGACTCGACTCTCCATCCTAG TGCAGTGCTAATCTAGGGGCCGATAGCAATCAGCCCCAGCATTTT 15, 4, 4, 0, ATCTAGGGGCCGA F T GAVVSQHPSWVICKSGTSVKIECRSLDFQATTMFWYRQFPKQSLMLMATSNEGSKATYEQGVEKDKFLINHASLTLSTLTVTSAHPEDSSFYICSANLGADSNQPQHFGDGTRLSIL DFQATT SNEGSKA T F T T 1 +26 TRBV14 TRBJ2-5 CASSSTGAQETQYF 1 GAAGCTGGAGTTACTCAGTTCCCCAGCCACAGCGTAATAGAGAAGGGCCAGACTGTGACTCTGAGATGTGACCCAATTTCTGGACATGATAATCTTTATTGGTATCGACGTGTTATGGGAAAAGAAATAAAATTTCTGTTACATTTTGTGAAAGAGTCTAAACAGGATGAGTCCGGTATGCCCAACAATCGATTCTTAGCTGAAAGGACTGGAGGGACGTATTCTACTCTGAAGGTGCAGCCTGCAGAACTGGAGGATTCTGGAGTTTATTTCTGTGCCAGCAGCTCGACCGGGGCGCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCTGGTGCTCG TGTGCCAGCAGCTCGACCGGGGCGCAAGAGACCCAGTACTTC 9, 10, 5, 2, TCGACCGGGGCG F T EAGVTQFPSHSVIEKGQTVTLRCDPISGHDNLYWYRRVMGKEIKFLLHFVKESKQDESGMPNNRFLAERTGGTYSTLKVQPAELEDSGVYFCASSSTGAQETQYFGPGTRLLVL SGHDN FVKESK T F T T 1 +27 TRBV3-1 TRBJ2-7 CASTLSGSSYEQYF 1 ACAGCTGTTTCCCAGACTCCAAAATACCTGGTCACACAGATGGGAAACGACAAGTCCATTAAATGTGAACAAAATCTGGGCCATGATACTATGTATTGGTATAAACAGGACTCTAAGAAATTTCTGAAGATAATGTTTAGCTACAATAATAAGGAGCTCATTATAAATGAAACAGTTCCAAATCGCTTCTCACCTAAATCTCCAGACAAAGCTCACTTAAATCTTCACATCAATTCCCTGGAGCTTGGTGACTCTGCTGTGTATTTCTGTGCCAGTACCCTTAGCGGGAGCTCCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACGGTCACAG TGTGCCAGTACCCTTAGCGGGAGCTCCTACGAGCAGTACTTC 21, 12, 9, 0, TACCCTTAGCGGGAG F T TAVSQTPKYLVTQMGNDKSIKCEQNLGHDTMYWYKQDSKKFLKIMFSYNNKELIINETVPNRFSPKSPDKAHLNLHINSLELGDSAVYFCASTLSGSSYEQYFGPGTRLTVT LGHDT YNNKEL T F T T 1 +28 TRBV7-9 TRBJ2-5 CASSLGAMQQETQYF 1 GATACTGGAGTCTCCCAGAACCCCAGACACAAGATCACAAAGAGGGGACAGAATGTAACTTTCAGGTGTGATCCAATTTCTGAACACAACCGCCTTTATTGGTACCGACAGACCCTGGGGCAGGGCCCAGAGTTTCTGACTTACTTCCAGAATGAAGCTCAACTAGAAAAATCAAGGCTGCTCAGTGATCGGTTCTCTGCAGAGAGGCCTAAGGGATCTTTCTCCACCTTGGAGATCCAGCGCACAGAGCAGGGGGACTCGGCCATGTATCTCTGTGCCAGCAGCTTAGGCGCCATGCAACAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCTGGTGCTCG TGTGCCAGCAGCTTAGGCGCCATGCAACAAGAGACCCAGTACTTC 43, 10, 1, 2, GCGCCATGCAA F T DTGVSQNPRHKITKRGQNVTFRCDPISEHNRLYWYRQTLGQGPEFLTYFQNEAQLEKSRLLSDRFSAERPKGSFSTLEIQRTEQGDSAMYLCASSLGAMQQETQYFGPGTRLLVL SEHNR FQNEAQ T F T T 1 +29 TRBV9 TRBJ1-1 CASSVRESNTEAFF 1 GATTCTGGAGTCACACAAACCCCAAAGCACCTGATCACAGCAACTGGACAGCGAGTGACGCTGAGATGCTCCCCTAGGTCTGGAGACCTCTCTGTGTACTGGTACCAACAGAGCCTGGACCAGGGCCTCCAGTTCCTCATTCAGTATTATAATGGAGAAGAGAGAGCAAAAGGAAACATTCTTGAACGATTCTCCGCACAACAGTTCCCTGACTTGCACTCTGAACTAAACCTGAGCTCTCTGGAGCTGGGGGACTCAGCTTTGTATTTCTGTGCCAGCAGCGTCAGGGAATCGAACACTGAAGCTTTCTTTGGACAAGGCACCAGACTCACAGTTGTAG TGTGCCAGCAGCGTCAGGGAATCGAACACTGAAGCTTTCTTT 44, 0, 2, 1, CAGGGAATC F T DSGVTQTPKHLITATGQRVTLRCSPRSGDLSVYWYQQSLDQGLQFLIQYYNGEERAKGNILERFSAQQFPDLHSELNLSSLELGDSALYFCASSVRESNTEAFFGQGTRLTVV SGDLS YYNGEE T F T T 1 +30 TRBV10-3 TRBJ2-5 CAISESRFGHTQETQYF 1 GATGCTGGAATCACCCAGAGCCCAAGACACAAGGTCACAGAGACAGGAACACCAGTGACTCTGAGATGTCACCAGACTGAGAACCACCGCTATATGTACTGGTATCGACAAGACCCGGGGCATGGGCTGAGGCTGATCCATTACTCATATGGTGTTAAAGATACTGACAAAGGAGAAGTCTCAGATGGCTATAGTGTCTCTAGATCAAAGACAGAGGATTTCCTCCTCACTCTGGAGTCCGCTACCAGCTCCCAGACATCTGTGTACTTCTGTGCCATCAGTGAGTCGAGGTTTGGACACACCCAAGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCTGGTGCTCG TGTGCCATCAGTGAGTCGAGGTTTGGACACACCCAAGAGACCCAGTACTTC 2, 10, 0, 1, GAGGTTTGGACACAC F T DAGITQSPRHKVTETGTPVTLRCHQTENHRYMYWYRQDPGHGLRLIHYSYGVKDTDKGEVSDGYSVSRSKTEDFLLTLESATSSQTSVYFCAISESRFGHTQETQYFGPGTRLLVL ENHRY SYGVKD T F T T 1 +31 TRBV6-1 TRBJ2-7 CASSESGRVYEQYF 1 AATGCTGGTGTCACTCAGACCCCAAAATTCCAGGTCCTGAAGACAGGACAGAGCATGACACTGCAGTGTGCCCAGGATATGAACCATAACTCCATGTACTGGTATCGACAAGACCCAGGCATGGGACTGAGGCTGATTTATTACTCAGCTTCTGAGGGTACCACTGACAAAGGAGAAGTCCCCAATGGCTACAATGTCTCCAGATTAAACAAACGGGAGTTCTCGCTCAGGCTGGAGTCGGCTGCTCCCTCCCAGACATCTGTGTACTTCTGTGCCAGCAGTGAGAGCGGGAGGGTCTACGAGCAGTACTTCGGGCCGGGCACCAGGCTCACGGTCACAG TGTGCCAGCAGTGAGAGCGGGAGGGTCTACGAGCAGTACTTC 31, 12, 3, 3, GAGCGGGAGGGT F T NAGVTQTPKFQVLKTGQSMTLQCAQDMNHNSMYWYRQDPGMGLRLIYYSASEGTTDKGEVPNGYNVSRLNKREFSLRLESAAPSQTSVYFCASSESGRVYEQYFGPGTRLTVT MNHNS SASEGT T F T T 1 +32 TRBV28 TRBJ2-1 CASGLADYNEQFF 1 GATGTGAAAGTAACCCAGAGCTCGAGATATCTAGTCAAAAGGACGGGAGAGAAAGTTTTTCTGGAATGTGTCCAGGATATGGACCATGAAAATATGTTCTGGTATCGACAAGACCCAGGTCTGGGGCTACGGCTGATCTATTTCTCATATGATGTTAAAATGAAAGAAAAAGGAGATATTCCTGAGGGGTACAGTGTCTCTAGAGAGAAGAAGGAGCGCTTCTCCCTGATTCTGGAGTCCGCCAGCACCAACCAGACATCTATGTACCTCTGTGCCAGTGGACTAGCGGACTACAATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGTGCCAGTGGACTAGCGGACTACAATGAGCAGTTCTTC 19, 6, 9, 3, TGGACTAGCGGA F T DVKVTQSSRYLVKRTGEKVFLECVQDMDHENMFWYRQDPGLGLRLIYFSYDVKMKEKGDIPEGYSVSREKKERFSLILESASTNQTSMYLCASGLADYNEQFFGPGTRLTVL MDHEN SYDVKM T F T T 1 +33 TRBV12-4 TRBJ2-1 CASSLGLLNEQFF 1 GATGCTGGAGTTATCCAGTCACCCCGGCACGAGGTGACAGAGATGGGACAAGAAGTGACTCTGAGATGTAAACCAATTTCAGGACACGACTACCTTTTCTGGTACAGACAGACCATGATGCGGGGACTGGAGTTGCTCATTTACTTTAACAACAACGTTCCGATAGATGATTCAGGGATGCCCGAGGATCGATTCTCAGCTAAGATGCCTAATGCATCATTCTCCACTCTGAAGATCCAGCCCTCAGAACCCAGGGACTCAGCTGTGTACTTCTGTGCCAGCAGTTTAGGTCTGCTTAATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGTGCCAGCAGTTTAGGTCTGCTTAATGAGCAGTTCTTC 6, 6, 1, 7, GTCTGCTT F T DAGVIQSPRHEVTEMGQEVTLRCKPISGHDYLFWYRQTMMRGLELLIYFNNNVPIDDSGMPEDRFSAKMPNASFSTLKIQPSEPRDSAVYFCASSLGLLNEQFFGPGTRLTVL SGHDY FNNNVP T F T T 1 +34 TRBV28 TRBJ2-7 CASSELTGRGEQYF 1 GATGTGAAAGTAACCCAGAGCTCGAGATATCTAGTCAAAAGGACGGGAGAGAAAGTTTTTCTGGAATGTGTCCAGGATATGGACCATGAAAATATGTTCTGGTATCGACAAGACCCAGGTCTGGGGCTACGGCTGATCTATTTCTCATATGATGTTAAAATGAAAGAAAAAGGAGATATTCCTGAGGGGTACAGTGTCTCTAGAGAGAAGAAGGAGCGCTTCTCCCTGATTCTGGAGTCCGCCAGCACCAACCAGACATCTATGTACCTCTGTGCCAGCAGTGAACTTACGGGGCGCGGCGAGCAGTACTTCGGGCCGGGCACCAGGCTCACGGTCACAG TGTGCCAGCAGTGAACTTACGGGGCGCGGCGAGCAGTACTTC 19, 12, 5, 6, GAACTTACGGGGCGCGG F T DVKVTQSSRYLVKRTGEKVFLECVQDMDHENMFWYRQDPGLGLRLIYFSYDVKMKEKGDIPEGYSVSREKKERFSLILESASTNQTSMYLCASSELTGRGEQYFGPGTRLTVT MDHEN SYDVKM T F T T 1 +35 TRBV7-6 TRBJ2-1 CASSQTMSGSQAHEQFF 1 GGTGCTGGAGTCTCCCAGTCTCCCAGGTACAAAGTCACAAAGAGGGGACAGGATGTAGCTCTCAGGTGTGATCCAATTTCGGGTCATGTATCCCTTTATTGGTACCGACAGGCCCTGGGGCAGGGCCCAGAGTTTCTGACTTACTTCAATTATGAAGCCCAACAAGACAAATCAGGGCTGCCCAATGATCGGTTCTCTGCAGAGAGGCCTGAGGGATCCATCTCCACTCTGACGATCCAGCGCACAGAGCAGCGGGACTCGGCCATGTATCGCTGTGCCAGCAGCCAAACCATGAGTGGCTCACAGGCCCATGAGCAGTTCTTCGGGCCAGGGACACGGCTCACCGTGCTAG TGTGCCAGCAGCCAAACCATGAGTGGCTCACAGGCCCATGAGCAGTTCTTC 40, 6, 5, 8, CAAACCATGAGTGGCTCACAGGCCC F T GAGVSQSPRYKVTKRGQDVALRCDPISGHVSLYWYRQALGQGPEFLTYFNYEAQQDKSGLPNDRFSAERPEGSISTLTIQRTEQRDSAMYRCASSQTMSGSQAHEQFFGPGTRLTVL SGHVS FNYEAQ T F T T 1 +36 TRBV25-1 TRBJ2-5 CASSEWGGGQETQYF 1 GAAGCTGACATCTACCAGACCCCAAGATACCTTGTTATAGGGACAGGAAAGAAGATCACTCTGGAATGTTCTCAAACCATGGGCCATGACAAAATGTACTGGTATCAACAAGATCCAGGAATGGAACTACACCTCATCCACTATTCCTATGGAGTTAATTCCACAGAGAAGGGAGATCTTTCCTCTGAGTCAACAGTCTCCAGAATAAGGACGGAGCATTTTCCCCTGACCCTGGAGTCTGCCAGGCCCTCACATACCTCTCAGTACCTCTGTGCCAGCAGTGAATGGGGAGGGGGCCAGGAGACCCAGTACTTCGGGCCAGGCACGCGGCTCCTGGTGCTCG TGTGCCAGCAGTGAATGGGGAGGGGGCCAGGAGACCCAGTACTTC 17, 10, 1, 5, GGGGAGGGGGCCAG F T EADIYQTPRYLVIGTGKKITLECSQTMGHDKMYWYQQDPGMELHLIHYSYGVNSTEKGDLSSESTVSRIRTEHFPLTLESARPSHTSQYLCASSEWGGGQETQYFGPGTRLLVL MGHDK SYGVNS T F T T 1 diff --git a/examples/input/TINY_1.fq b/examples/input/TINY_1.fq new file mode 100644 index 0000000..fe19069 --- /dev/null +++ b/examples/input/TINY_1.fq @@ -0,0 +1,424 @@ +@A00261:687:HMMCCDSX5:3:1101:24948:1047 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTCTCTGGTGGGAACACGTTTTTCAGGTCCTCTGGAAAGGGAAGAGGTAATGGGGCTAGGGTTGCTCTAAGAGCTGTCTGGTCCTGGTAGGGGCTCTGTGTATGTGTGAGAGAGAAGGCCGGGAAAGGACCATAATGAAGCACAATTATG ++ +FFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:27272:1063 1:N:0:AACCGCGGAT+CTAGCGCTGT +TCTTTCGGTGGGAACACCTTGTTCAGGTCCTCTACAACGGTTAACCTGGTCCCCGAACCGAAGGTGTAGCCATAGTCCGTTGTCCCCACAATGTTACAGCTTTGTACAAAAACAGCCCCTCCCATGGGTCCGCCCCCAGAGCCTGGGAGA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:5050:1094 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAAACACGGCAGGGTCAGGGTTCTGGATATTTGGTTTAACTAGCACCCTGGTTCCTCCTCCAAAAGTTAGCTTGTTGCCTGCAGCTTTAGCCACCACACAGAGGTAGGTGGCTGAATCACTGAGCTTGGAGTCTCTGATGAGCAGGGAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF::FFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18331:1110 1:N:0:ACCCGCGGAT+CTAGCGCTGT +TTATTTACGGCAGGGTCAGGGTTCTGGATATCAGGTAAAACAGTCAATTGTGTCCCAGATCCAAAGGTCAGTTGCCTTGCAGAACCAGAAAAGGCGCCCCTGACGATGCAATAGTACACAGCAGTGTCTCTCAGCGTAGCGTGGGGCAGG ++ +FFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:8684:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGCGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGTCCCCGCGGGGGTAGTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFF:FFF,:FFFFFF:FFFF:FFFFF:FFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:4:1101:8684:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGCGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGTCCCCGCGGGGGTAGTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFF:FFF,:FFFFFF:FFFF:FFFFF:FFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:5:1101:8684:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGCGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGTCTCCGCGGGGGTAGTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFF:FFF,:FFFFFF:FFFF:FFFFF:FFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:5:1101:8684:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGCGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGTCTTCGCGGGGGTAGTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFF:FFF,:FFFFFF:FFFF:FFFFF:FFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:5:1101:8684:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGCGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGTCTTTGCGGGGGTAGTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFF:FFF,:FFFFFF:FFFF:FFFFF:FFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:6:1101:8684:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGCGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGTCCCCGCGGGGGTAGTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFF:FFF,:FFFFFF:FFFF:FFFFF:FFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:7:1101:8684:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGCGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGGGGGTAGGTCCCCGCTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF,FFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFF:FFF,:FFFFFF:FFFF:FFFFF:FFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:3:1101:29604:1125 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAACCGGTGGGAACACCTTGTTCAGGTCCTCTACAACTGTGAGTCTGGTGCCTTGTCCAAAGAAAGCTTCAGCTGTCCCTCCGCTGCTGGCACAGAGATACATGGCCGAGTCCCCCTGCTCTGTGCGCTGGATCTCCAAGGTGGAGAAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:17300:1141 1:N:0:AACCGCGGAT+CTAGCGCTGT +CCCCAAGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCAGTGTAGGAGTATTCCGCTAGGCTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF::FF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26793:1141 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAGGCGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTGTAGGAGAACCCCCCGTGAGGGCTGCTGGCACAGAAATAAACTCCAGAATCCTCCAGTTCTGCAGGCTGCACCTTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:28854:1141 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAGAAGGTGGGAACACGTTTTTCAGGTCCTCCAGTACGGTCAGCCTAGAGCCTTCTCCAAAAAACAGCTCCCCTCACCCGCTAGTCCCGGCTGCTGGCGCAGAAATATACAGATGTCGGGGAGGAGGCAGCAGACTCCAGAGTGAGGGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF:FFFFF,:,,FF:FF,,,,:FF,,:FF,FFFFFF,:F:,FFFFF:,,,FF,FFFF +@A00261:687:HMMCCDSX5:3:1101:22480:1157 1:N:0:CACCGCGGAT+CTAGCGCTGT +AGAGGAGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGTAGGAGACCCCTGGCTCCAAGCTGCTGGCGCAAAGATAAAGGGCCGAGTCCCCCAGCTCCAAGGTGCTCACATTCATC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF::FFFFF +@A00261:687:HMMCCDSX5:3:1101:28926:1172 1:N:0:CACCGCGGAT+CTAGCGCTGT +AGTCGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACTGTCAGCCGGGTGCCTGGGCCAAAATACTGCGTATCTGTGGAAAAGGTGCCGCTAGGACTGCTGGCACAGAGATACAGGGCCGAGTCTTCTGGCTGCAGGGTGTGTAGGTGA ++ +,FFFF:FFFFFFFFFFFFFFFFFF:FF:FFF:FFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF:FFFFFF +@A00261:687:HMMCCDSX5:3:1101:28592:1188 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGAAAAACGGCAGGGTCAGGGTTCTGGATATCAGAAGAGATTAACAGTTCCACCATAGGAATAGTTGCATCCCAACTCCAGTGAGGCTGCTTCAGAGAGAATTACGTGGTGGTTATGCTGGCTCACAGACTGGGCTCTGGCATCTCTCAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:4417:1219 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAAACACGGCAGGGTCAGGGTTCTGGATATTTGGTTTAACTAGCACCCTGGTTCCTCCTCCAAAAGTTAGCTTGTTGCCTGCAGCTTTAGCCACCACACAGAGGTAGGTGGCTGAATCACTGAGCTTGGAGTCTCTGATGAGCAGGGAA ++ +FFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFF, +@A00261:687:HMMCCDSX5:3:1101:3794:1235 1:N:0:ACCCGCGGAT+CTAGCGCTGT +CTCGGTACGGCAGGGTCAGGGTTCTGGATATTCGGATTTACTGCCAGGCTTGTTCCCAATCCCCAAATCAGCTTACGGTTGTTGCCAGCATTTAAGGTCTCCACAGCACAGAGGTAGATGGCCGAGTCTCCGGTCTGGGTGGCTGTGATG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF:FFFFF,FFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,::FFF:FFF:FFFF:F +@A00261:687:HMMCCDSX5:3:1101:23737:1266 1:N:0:AACCGCGGAT+CGAGCGCTGT +CCTCACGGTGGGAACACCTTGTTCAGGTCCTCCAAGACAGAGAGCTGGGTTCCACTGCCAAAAAACAGTTTTTCACTAACTGCCCCAAAACCCCTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFF,FFFF:FFFF,FF,:F,FF,F:FFFFFFFF,FFFFFFFFFFF:F:,F,FFFF,F::F::FF::::FFFF:FF:FFFFFF:FFFFFF::FF:FFFFFFF:FF,FF:FF:FF:F:F,FFFFF::FF:F,FFFF:F:FFF:F:FF +@A00261:687:HMMCCDSX5:3:1101:29966:1282 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTGCCGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGTAGGAGCTGCCCCTTTCAACGCTGCAGAGATATATGCTGCTGTCTTCAGGGCTCATGTTGCTCACAGTCAGAGTTGAG ++ +FFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFF +@A00261:687:HMMCCDSX5:3:1101:20627:1297 1:N:0:ACCCGCGGAT+CTAGCGCTGT +CGCGTGGTGGGAACACGTTTTTCAGGTCCTCCAGCACGGTCAGCCTGCTGCCGGCCCCGAAAGTCAGGACGTTGGCCCCAGCATCGGTGCTGCTGGCACAGAGATACAGGGCCGAGTCGTCCAGCTCCAAGGCGTTCACATTCAGCTCAG ++ +FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:19027:1313 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGTGACACGGCAGGGTCAGGGTTCTGGATATTGGGTATGATGGTGAGTCTTGTTCCAGTCCCAAAGGTTAATTTCTCATTTCAAACAGCACAGAGGTAGGTGGCTGAGTCACCAGGCTGAGAAGCTGCAATGTATAAAGTACTACGTCCT ++ +FFFFFFFFFFFFFFFFFFFFFFFF:FF:FF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:6876:1344 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTATCTACGGCAGGGTCAGGGTTCTGGATATTCGGATTTACTGCCAGGCTTGTTCCCAATCCCCAAATCAGCTTACGGTTGTTGCCAGCATTACACAGTCATAGAAAGCTTTACCAAAACCGACCAGGCCAATGTGTATCAGGAAATCCC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:29170:1344 1:N:0:ACCCGCGGAT+CTAGCGCTGT +TTTTTTACGGCAGGGTCAGGGTTCTGGATATCTGGCTTTATAATTAGCTTGGTCCCAGCGCCCCAGATTAACTGATAGTTGCTATCCATCACAGCACAGAGGTAAGAGGCAGAGTCTTTCATCTGGAGCTCCTTCAAAAGGAGGTAACTG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFF:FFFFFFFF:FFFFFFFFF,,FFFFFFF,FFFFFF,FFF,FFFFFFFFFFFFFF:F,FF:FFFFFF:F +@A00261:687:HMMCCDSX5:3:1101:21260:1360 1:N:0:AACCGCGGAT+CTAGCGCTGT +ATTTTTACGGCAGGGTCAGGGTTCTGGATATTTGCAATCACAGAAAGTCTTGTGCCCTTTCCAAAGACAAGAGGTGTGTTTCCTCCGAGTCGTCTCACAGCACAGAGGTAAGAGGCAGAGTCTTTCATCTGGAGCTCCTTCAAAAGGAGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18611:1376 1:N:0:ACCCGCGGAT+CTAGCGCTGT +CTTTTTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACTGTCAGCCGGGTGCCTGGGCCAAAATACTGCGTATCTGTGCTCACAGCCTGGAGGCCCAGGACAAAAACCGGTGCTGCCCGCTCTGCTGTCCCAGACTCAGCTCGGGTCCTT ++ +:FFFFFFFF,,FFF:FFF:F:FFFF,,,F:FF::FF:,:F:FFF,,,F,F,,F,::FFFFFFF:FF:FF,FFF:FFF,FFFFFF:FFFF,FF,FFFF:,FFFFFFFFFF:FFFFF:FFFFFFFFFF,:F:,F:F:FFFF:FF,FFFFF:F +@A00261:687:HMMCCDSX5:3:1101:7446:1391 1:N:0:AACCGCGGAT+CTAGCGCTGT +CCCCCTGGTGGGAAGGAGGGCCCGTTGGGAGGCCCAGCGGGCAGGAGGAACGGCTACCGAGGCTCCAGCTTAACGGTATTTGGAGGTCAGCACGGTGCTCACAGAAGCCAGGAACTTGTCCAGGGAGGCGTGCACCGCAGGGGTGAACTC ++ +,FF,F:FFFFF:F,FFFFFFFFFF:FFFFFFFFF::FFFF:FF:FF:FFFF:FFF:FFFF,FFFFFFF:FF:F,FFF:FF:FFF::FFFFFFFFF:FFF::FFFFFFFFFFF::F:FF:FFF:FF,FF:FFFFFFFFFFF:F:F:FFFFF +@A00261:687:HMMCCDSX5:3:1101:8024:1391 1:N:0:AACCGCGGAT+CTAGCGCTGT +ATAAGTGGTGGGAACACCTTGTTCAGGTCCTCGAGCACTGTCAGCCGGGTGCCTGGGCCAAAATACTGCGTATCTGCCCCCCCGAGTTGGCGGCTGGGGCCGAGATAAAGGGCTGTGTGTTCTTGCTGCAGGGCGGGTAGGTGGTGGTTT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF,,,,,:F,:,,:,,:,,,:,,F,,,,,:,,F,:,F,F,,,,FF,F,F:,,F,:,,:::,F +@A00261:687:HMMCCDSX5:3:1101:24641:1391 1:N:0:ACCCGCGGAT+CTAGCGCGGT +TGTGCAGGTGGGAACACCTTGTTCAGGTCCTCTACAACTGTGAGTCTGGTGCCTTGTCCATAGAAAGCTTCAGTGGTCCCCCTGCTGCTGGCACAGAGATACATGGCCGAGTCCCCCTGCTCTGTGCGCTGGATCTCCAAGGTGGAGAAA ++ +F:FFFFFF::FFFFFFF:FFF:FFF,FFFFF:FFFFFFF:F:FFF,FFFFFFFFFF,FFF,:F:F:FF:FF,FFFF::FF,FFF:FFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF:F:FFFFFFFFFF:FFF,F,FFFFFFFFF: +@A00261:687:HMMCCDSX5:3:1101:26585:1407 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTGTTACGGCAGGGTCAGGGTTCTGGATATCAGGCCTGACCAGCAGTCTGGTCCCACTCCCGAAGATCAATTTATAGCTGCTACCCCCGCCCTCTCTCATTGCACAGAAGTACATTGCGGAGTCCCCCAGTTGTGAAGCGGAGATGACA ++ +:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFF,FFFFFFFFFFF:FFFFFFFFFFFFFFFFF, +@A00261:687:HMMCCDSX5:3:1101:30653:1407 1:N:0:AACCGCGGAT+CTACCGCTGT +CAATCAACGGCAGGGTCAGGGTTCTGGATATTGGGTTTCACAGATAACTCCGTTCCCTGTCCGAAGATAAGCTTTCCTCCCTGGTTCATGCTTGCTGCACAGAAGTAGGTGGCTGAGTCTCCAGGCTGGGAATCCATGATATGCAATGAG ++ +FFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFF:FFFFF +@A00261:687:HMMCCDSX5:3:1101:5638:1423 1:N:0:CACCGCGGAT+CTAGCGCTGT +AATCTTGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGTAGCTCCCGCTAGTCGTCTCGGTGCTGGCACAGAAGTACACAGATGTTTGGGAGGGAGCAGCCGACTCCAACCCCAGC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24713:1423 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTTGTGGTGGGAACACCTTGTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCGCTCGGCTTCCCTCCCGTCCGTTTGCTGGCACAGAAGTACATGGCTGAGTCCTCCAGCTTTGTGGACCGGATCTTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24379:1438 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGACGTGGTGGGAACACCTTGTTCAGGTCCTCTACAACTGTGAGTCTGGTGCCTTGTCCAAAGAAAGCTTCAGTGGCCCCCTGTTCCCAGGCACAGAGATAGAAGCCAGAGTCACTGAGAAGGAGCTTCTTAGAACTCAGGATGAACTGC ++ +FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:30309:1438 1:N:0:ACCCGCGGAT+CTAGCGCTGT +TTTTATACGGCAGGGTCAGGGTTCTGGATATAGGGCAGCACGGACAATCTGGTTCCGGGACCAAAGACAAAATTCTGACCATAGTTCCGGCTCACAACACAGAAGTACTCAGCCGCGTCGCTCATATGGGCTGAGGGTTTCGTCAGGTGG ++ +FFFFFFFFFFFFFFFFFF:FFFFFFFFFF:FFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF::FFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:8133:1454 1:N:0:AACCGCGGAT+CTAGCGCTGT +CACACCGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTGTAATTAAACACCCGCTGCCCATCTAAGCTGCTGGCACAGAGATACACGGCCGAGTCCTCAAGCTTTGCAGGTTGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24532:1454 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGTCTTACGGCAGGGTCAGGGTTCTGGATATTTGGTTTAACTAGCACCCTGGTTCCTCCTCCAAAAGTTAGCTTGTTGCCTGCAGCTTTGATCCCCCTCGCCCGATCCTCCACAGCACAGTAGTAAACAGCAGCATCTGCCTCCCGCACC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFF:F,FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF:FFFFFF:FF:FFFFFF +@A00261:687:HMMCCDSX5:3:1101:28890:1454 1:N:0:AACCGCGGAT+CTAGCGCTGT +TGAGGTACGGCAGGGTCAGGGTTCTGGATATTTGGTTTTACTGTCAGTCTGGTCCCTGCTCCAAAGCGCATGTCATTGAGCTCACTCAGAGCACAGAAGTACACCGCTGAGTCTGACACTTGAACTGAGCCTTTCTCCAAGTGGAAAGAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:4869:1470 1:N:0:AACCGCGGAT+CTAGCGCTGT +CATGGTGGTGGGAACACCTTGTTCAGGTCCTCCAAGACAGAGAGCTGGGTTCCACTGCCAAAAAACAGTTCCTTTGCTGTCCCCCAGATTAAGCTGCTGGCACAGAGATACACGGCCGAGTCCTCCTGCTGTGTGCGCTGGATCGTCAGA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24542:1470 1:N:0:AACCGCGGAT+GTAGCGCTGT +AGTCTTACGGCAGGGTCAGGGTTCTGGATATTTGGTTTAACTAGCACCCTGGTTCCTCCTCCAAAAGTTAGCTTGTTGCCTGCAGCTTTGATCCCCCTCGCCCGATCCTCCACAGCACAGTAGTAAACAGCAGCATCTGCCTCCCGCACC ++ +F,,F::FF::FF::FFFFFFFFFFFF:F,::FFF:FFFFFFFFFFFFFFFFFFFFFFF:F:FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFF:FFFFFFFFFFFFFFF: +@A00261:687:HMMCCDSX5:3:1101:6777:1485 1:N:0:AACCGCGGAT+CTAGCGCTGT +AACGACGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCGCGGCCCCCTGCGGCTTCAACGCTGCAGAGATATATGCTGCTGTCTTCAGGGCTCATGTTGCTCACAGTCAGAGTTGAG ++ +FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFF,FFFF:F:FFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:31331:1485 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACTGTCAGCCGGGTGCCTGGGCCAAAATACTGCGTATCTGTGCTGCTAGTCCTCCCCAAGCTGCTGGCACAGAGATACATGGCCGAGTCCCCCTGCTCTGTGCGCTGGATCTCC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF,:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18629:1501 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGAGAAGGTGGGAACACGTTTTTCAGGTCCTCCAGCACGGTCAGCCTGCTGCCGGCCCCGAAAGTCAGGACGTTGGCCCCAGCCTGGAGAGCAACGCTGCAGAGATATATGCTGCTGTCTTCAGGGCTCATGTTGCTCACAGTCAGAGTT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FF:FFFFFFFFFFF:F +@A00261:687:HMMCCDSX5:3:1101:31448:1501 1:N:0:ACCCGCGGAT+CTAGCGCTGT +AAAGGTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACTGTCAGCCGGGTGCCTGGGCCAAAATACTGCGTATCTGTGCTGCTAGTCCTCCCCAAGCTGCTGGCACAGAGATACATGGCCGAGTCCCCCTGCTCTGTGCGCTGGATCTCC ++ +F,FFFFF:FFFFF,FFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:1768:1517 1:N:0:AACTGCGGAT+CTAGCGCTGT +ATTGAAGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGTAGTCTCCCGCTAGTCGAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGAGTGGAC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFF,FFFF +@A00261:687:HMMCCDSX5:3:1101:9607:1532 1:N:0:AACCGCGGAT+CTAGCGCTGT +CTCGCTGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTGTAGGCTCCTGTCCCGGGACTGCTGGCACAGAAGTACACAGATGTCTGGGAGGGAGCAGCCGACAGCAGCCTGAGC ++ +FFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:10854:1532 1:N:0:ACCCGCGGAT+CTAGCGCTGT +CAAACAACGGCAGGGTCAGGGTTCTGGATATTTGGTATGACCGAGAGTTTGGTCCCCTTCCCGAAAGTGAGTTGGTAACTCCCAGCCCCAGGGGCCGGAGGGGCACAGAGGTATGTGGCTGAGTCTTCAGGCTGGGATCCTTTGATGTAC ++ +FFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:16866:1548 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGTGTTACGGCAGGGTCAGGGTTCTGGATATGGTATGTGGCTTCAAAACCTTTGTTGCTTCCCTTGTCATCAGCCTTCGTGGCTTTCAGGAGGAGCTGTAGACCTTCTCCAGGATATTGGACATACCAGAAAAGGGAAGGGTATCCTGTG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FF:FFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26214:1548 1:N:0:AACCGCGGAT+CTAGCGCTGT +ACTTGAGGTGGGAACACCTTGTTCAGGTCCTCTGGAAAGGGAAGAGGGGTTGGAGCCAGGGTTGCTCTGAGAGCTGTCTGGTTCTGGTAGGGGCTCTGTGTGTGTGAGAGAGAGCCCTGGGAATGGGCCGTAATGAAGCACAACTATGTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:30120:1548 1:N:0:AACCGCAGAT+CTAGCGCTGT +CCAGTTGGTGGGAACACCTTGTTCAGGTCCTCCAAGACAGAGAGCTGGGTTCCACTGCCAAAAAACAGTTTTTCATTAGTTGCCGCGCGGGGGCTGCTGGCACAGAAATACAAAGCTGAGTCCCCCAGCTCCAGAGAGCTCAGGTTTAGT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:7238:1564 1:N:0:AACCGCGGAT+CTAGCGCTGT +TACTCCGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTGTGTCCCGCTAGACCCAAGCTGCTGGCACAGAGATAGAGGGCCGAGTCCCCCAGCAACAAGGCGTTCACATTCAGC ++ +FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF:FFFFFFFFFFFFFF,FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:13964:1564 1:N:0:CACCGCGGAT+CTAGCGCTGT +AGTGGTACGGCAGGGTCAGGGTTCTGGATATTTGGTATGACCACCACTTGGTTCCCCTTCCCAAAAGCGAGTCTGTTGTAAGGATCACCCACGAGGCAGTAGTACACAGCAGTGTCGCTCAGGGAAACCCGGGGCAGGCTCAGAGTGCTG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24777:1564 1:N:0:AACCGCGGAT+CTAGCGCTGT +CCCCGTGGTGGGAACACGTTTTTCAGGTCCTCCAGTACGGTCAGCCTAGAGCCTTCTCCAAAAAACAGCTCCCGGCTAGTCACCTCACTGCTGGCACAGAAGTACATGGCTGAGTCCTCCAGCTTTGTGGACCGGATCTTCAGAGTGAAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF,FFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18195:1595 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGAGCAACGGCAGGGTCAGGGTTCGGGATATCAGAAGAGATTAACAGTTCCACCATAGGAATAGTTGCATCCCAACTCCAGTGAGGCTGCTTCAGAGAGAATTACGTGGTGGTTATGCTGGCTCACAGACTGGGCTCTGGCATCTCTCAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18484:1626 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGAGAAGGTGGGAACACGTTTTTCAGGTCCTCCAGCACGGTCAGCCTGCTGCCGGCCCCGAAAGTCAGGACGTTGGCCCCAGCCTGGAGAGCAACGCTGCAGAGATATATGCTGCTGTCTTCAGGGCTCATGTTGCTCACAGTCAGAGTT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF:,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF +@A00261:687:HMMCCDSX5:3:1101:11966:1642 1:N:0:AACCGCGGAT+CTAGCGCTGT +TGCTGTGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTCCTCGGTGCTAGTCCTGGGCTGCTGGCGCAAAGATAAAGGGCCGAGTCCCCCAGCTCCAAGGTGCTCACATTCATC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26874:1658 1:N:0:AACCGCGGAT+CTAGCGCTGT +ATATAAACGGCAGGGTCAGGGTTCTGGATATAGGCAGTTCAGGGTGACTGCTTTCCTCACTGGCATGGATACTGATGACTGGGCTTGAGTAACCTTCTGGGCCACACTTGATCCAGAGTAGCTGAAGGCCACAAATACACACAGCAGGCT ++ +FF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF,FFFFFFFFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:3:1101:27597:1658 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTAGTAACGGCAGGGTCAGGGTTCTGGATATTTGGAAAGACTTGTAATCTGGTCCCAGTCCCAAAGATGAGCTTGTCGGTGTTATACCTCACAGCACAGAGGTAAGAGGCAGAGTCTTTCATCTGGAGCTCCTTCAAAAGGAGGTAACTG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF,FFFFFFF,,FFF +@A00261:687:HMMCCDSX5:3:1101:31901:1658 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTTATACGGCAGGGTCAGGGTTCTGGATATAGGGCAGCACGGACAATCTGGTTCCGGGACCAAAGACAAAATTCTGACCATAGTTCCGGCTCACAACACAGAAGTACTCAGCCGCGTCGCTCATATGGGCTGAGGGTTTCGTCAGGTGG ++ +,FFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFF:FF:FFFFFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFF,FFFFFF:FFFF,F,,FFFF:FF,FF:FFF:F,F,FFF,FFFFFF,FFFFF:FFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:27118:1673 1:N:0:AACCGCGGAT+CTAGCGCTGT +ACTTGAGGTGGGAACACCTTGTTCAGGTCCTCTGGAAAGGGAAGAGGGGTTGGAGCCAGGGTTGCTCTGAGAGCTGTCTGGTTCTGGTAGGGGCTCTGTGTGTGTGAGAGAGAGCCCTGGGAATGGGCCGTAATGAAGCACAACTATGTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFF,FFF:FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:7283:1705 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTTTTGGTGGGAACACCTTGTTCAGGTCCTCTAGGATGGAGAGTCGAGTCCCATCACCAAAATGCTGGGGCTGATTGCTATCGGCCCCTAGATTAGCACTGCAGATGTAGAAGCTGCTGTCTTCAGGATGGGCACTGGTCACTGTCAGA ++ +FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF:FFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:3332:1720 1:N:0:AACCGCGGAT+CTAGCGCTGT +TGAAATACGGCAGGGTCAGGGTTCTGGATATAGGACTTGACTCTCAGAATGGTTCCTGCGCCAAAGACCAGCTTGTTTCCATATTCCACCGCAGCACAGAGGTAGGTGGCTGAGTCACCAGGCTGAGAAGCTGCAATGTATAAAGTACTA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFF:FFFFFFFFFFF:FFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24722:1720 1:N:0:AACCGCGGAT+CTAGCGCTGT +TGTGTTACGGCAGGGTCAGGGTTCTGGATATCTGGTTGGACCAAGACTGTTCTTGCTACAGGCTGCTCCAGCGGTGGCAACAACAACACCATGGAGTGTCCCTTCTTCCTCACAGTTTTCCTTCAAACGCTGGCTTTGTGGAGGGCATTG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24975:1720 1:N:0:AACCGCGGAT+CTAGCGCTGT +AACGGGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGCGCCCCGGTCGAGCTGCTGGCACAGAAATAAACTCCAGAATCCTCCAGTTCTGCAGGCTGCACCTTCAGAGTAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFF::FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF:FFFFFFFFFFFF:FFF:FFF,FFFF,FF,FFF,FFFFFFFFFFFFFF:F +@A00261:687:HMMCCDSX5:3:1101:14172:1736 1:N:0:AACCGCGGAT+CTAGCGCTGT +CCTGGTACGGCAGGGTCAGGGTTCTGGATATTGAGTTCCACTTTTAGCTGAGTGCCTGTCCCAAAGGTGAGTTTGTTTCCTCCTCCCGTGAGTATCACAGTGTTTGATGCCTCACAATAAACTCCCTTTTGTCACTTCACCCCTCTTCTT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18041:1736 1:N:0:AACCGCGGAT+CTAGCGCTGT +ATTCTTACGGCAGGGTCAGGGTTCTGGATATTTGGCTTCACAGTGAGCGTAGTCCCATCCCCAAAGGTTGATTTGCCTGGAATCTCTGCACAAAAGTAGACAGCTGAGTCTCCAGGTTGAGTAGCTGCAATTTGCAGAGAGAGATGTTTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:28800:1736 1:N:0:AACCGCGGAT+CTAGCGCTGT +ACAACAGGTGGGAACACGTTTTTCAGGTCCTCGAGCACTGTCAGCCGGGTGCCTGGGCCAAAATACTGCGTATCTGTGCTCACAGCCTGGAGGCCCAGGACAAAAACCGGTGCTGCCCGCTCTGCTGTCCCAGACTCAGCTCGGGTCCTT ++ +:FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26350:1752 1:N:0:AACCGCGGAT+CTAGCGCTGT +ATAAAAACGGCAGGGTCAGGGTTCTGGATATTTGGTTTAACTAGCACCCTGGTTCCTCCTCCAAAAGTTAGCTTGTTGCCTGCAGCCCCTGTCACGGCACAGAGGTAGGTGGCTGAATCACTGGGCTGGGAGTCTCTGATGAGCAGAGAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:3287:1767 1:N:0:AACCGCGGAT+CTAGCGCTGT +TAAAACACGGCAGGGTCAGGGTTCTGGATATTTGCTTTAACAAATAGTCTTGTTCCTGCTCCAAAGATAGTTTTGAAGCCTCCTATACTTGCTGCACAGAAGTAGACAGCCGAGTCTTCAGGTTGGGTCTCTGTGATGTGCAGGGAGAAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF:F:FFFFFFFFFFF:FFFFF:FFFFFFFF::FFFFFFFF:FFFFFFFFFF:FFFFFFF:FFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26115:1783 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAAAGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGTAGGAGCTCCCGCTAAGGGTACTGGCACAGAAATACACAGCAGAGTCACCAAGCTCCAGGGAATTGATGTGAAGATTTA ++ +:FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:,FFFFFFFF:FF:FFF +@A00261:687:HMMCCDSX5:3:1101:26295:1783 1:N:0:AACCGCGGAT+CTAGCGCTGT +ATAGAAACGGCAGGGTCAGGGTTCTGGATATCAGGCCAGACAGTCAACTGAGTTCCTCTTCCAAAGTATAGCCTCCCCAGGGTTGAGCCTCTGTCGGCACAATGCTGGTCCCCTTTACATGAACCCAGCCCTAATGTTTCCCATCCCCTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26874:1783 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAAGAGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGTTGCATGGCGCCTAAGCTGCTGGCACAGAGATACATGACCGAGTCCCCCTGCTCTGTGCGCTGGATCTCCAAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF:FFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:7862:1799 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTGTTGGTGGGAACACCTTGTTCAGGTCCTCTACAACTGTGAGTCTGGTGCCTTGTCCAAAGAAAGCTTCAGTGTTCGATTCCCTGACGCTGCTGGCACAGAAATACAAAGCTGAGTCCCCCAGCTCCAGAGAGCTCAGGTTTAGTTCA ++ +FFFFFFFFFF:FFFFFFFFFFFFFFF:FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FF:FFFFFF:FFFFFF:FFFFFFFFFFFFFFFFF:FFFFF:FFFFFFFFFFF,FF:FFFFFFFFFFFF:FFFFF +@A00261:687:HMMCCDSX5:3:1101:2302:1814 1:N:0:AACCGCAGAT+CTAGCGCTGT +ACAAAAACGGCAGGGTCAGGGTTCTGGATATCGAGCGTGACCTGAAGTCTTGTTCCAGTCCCAAAGGTGAGTTTACTGGCAGTGCCGGTTCTGGTCAGAGCACAGAAGTACACCGCTGAGTCTGACACTTGAACTGAGCCTTTCTCCAAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:6189:1814 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAGGAACGGCAGGGTCAGGGTTCTGGATATACGGATGAACAATAAGGCTGGTTCCTCTTCCAAATGTAGGTATGTAGCTTCCAGCGTCCGTAGCACAGAAGTAAGAAGCAGTGTCTGCTGCCCGGGAAGCCGTGATCAACAAGGAACTG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:3:1101:13096:1814 1:N:0:AACCGCGGAT+CTAGCGCTGT +CCCCCTGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTGTATTCCGTCCCCACAATGTTACAGCTTTGTACAAAAACAGCCCCTCCCATGGGTCCGCCCCCAGAGCCTGGGAGA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:29243:1814 1:N:0:ACCCGCGGAT+GTAGCGCTGT +TTTTTTACGGCAGGGTCAGGGTTCTGGATATCTGGCTTTATAATTAGCTTGGTCCCAGCGCCCCAGATTAACTGATAGTTGCTATCCATCACAGCACAGAGGTAAGAGGCAGAGTCTTTCATCTGGAGCTCCTTCAAAAGGAGGTAACTG ++ +F,FFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFF:FFFFFFF:FFFFFFF:FFFFFFFFFF::FFFFFF:FFFFFFFFFFFFFFFF,FFFFFFF:FFFFFFFFF,FFFFFF +@A00261:687:HMMCCDSX5:3:1101:29993:1830 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTTTTACGGCAGGGTCAGGGTTCTGGATATCTGCTGCGTCCCCCGAGGTCCGGCTTCCTGCCTCCCATCTCTCTCCATTTCCTCTGGTTTTTGTCTTGAGGCTTACTGATTTCTCATGTTCAGCATCTCTCTTTCTCAAGGTCTCATTT ++ +F:FFFFFFFFFFFFFFFFFFFFFFFFFF::,FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFF,FFFFFFF:FFFFFFF,FFFFFFFFF:FF +@A00261:687:HMMCCDSX5:3:1101:3296:1846 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTACCACGGCAGGGTCAGGGTTCTGGATATTTGCTCTTACAGTTACTGTGGTTCCGGCTCCAAAGCTGAGCTTGTAGTCGCCTATGTCACCCACGAGGCAGTAGTACACAGCAGTGTCGCTCAGGGAAACCCCGGGCAGGCTCAGAGTG ++ +FF:FFFF:FFFFFF:FFFF:F:,FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,F:FFFFFFF,FFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:FFF:FFFFFFFF:FFF,FFFFFFF,,FF,,,FFF,F:FFF,FFF,:F:FF: +@A00261:687:HMMCCDSX5:3:1101:3314:1846 1:N:0:ACCCGCGGAT+CTAGCGCTGT +TAAATAACGGCAGGGTCAGGGTTCTGGATATTGGGGAGAATATGAAGTCGTGTCCCTTTTCCAAAGATGATCTTGTCTTCAGCTCCTGCACAGAGGTAGAGGCCTGTATCACCAGGCTGGGCTGCAGTGATGTGGAGAGAACTGTCCTTT ++ +FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF:FFFFFF,FFFF:F,:FFFFFFFF,F,F,F::F:F:FFF:FFFFFFFF:FFF:FFFFFFFFFFFFFFFFFFFFF,F:FFFFFFFFFFF,F,FFFFFF +@A00261:687:HMMCCDSX5:3:1101:9579:1861 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAGGAGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCTTGGGTGTGTCCAAACCTCGACTCACTGATGGCACAGAAGTACACAGATGTCTGGGAGCTGGTAGCGGACTCCAGA ++ +:FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF,FFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:16396:1861 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGACAAGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGTAGACCCTCCCGCTCTCACTGCTGGCACAGAAGTACACAGATGTTTGGGAGGGAGCAGCCGACTCCAACCCCAGCAGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F +@A00261:687:HMMCCDSX5:3:1101:19027:1877 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAATAGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTGTAGTCCGCTAGTCCACTGGCACAGAGGTACATAGATGTCTGGTTGGTGCTGGCGGACTCCAGAATCAGGGAGAAG ++ +FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF:FFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:23095:1877 1:N:0:AACCGCGGAT+CTAGCGCTGT +CTTTGTGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGTAGGAGCTTACGACGCCTCTCGAGTAACTGCTGGCACAGAAGTACACAGATGTCTGGGAGGGAGCAGCCGACAGCAGC ++ +F:FFF,:FFF:FF:FFF:FFFFF,FFFF,FFFFFF:FFF::FFF:::FFFFFFFF:FFF,FFFFFF:FFFFFFFF:FF:FFFFFF:F:F:FFF:FFFFFFFFFFFF:FFFFFFF:FFFFF::F,FFF:FFFFF:FF:FF:FFF:FFFFFF +@A00261:687:HMMCCDSX5:3:1101:1099:1892 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGAGGAGGTGGGAACACGTTTTTCAGGTCCTCGAGCACTGTCAGCCGGGTGCCTGGGCCAAAATACTGCGTATCTGTGCTCACAGCCTGGAGGCCCAGGACAAAAACCGGTGCTGCCCGCTCTGCTGTCCCAGACTCAGCTCGGGTCCTT ++ +:FFFFFFFFFFF:FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:2546:1892 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTCTCTGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTAAGCAGACCTAAACTGCTGGCACAGAAGTACACAGCTGAGTCCCTGGGTTCTGAGGGCTGGATCTTCAGAGTGGAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:9833:1892 1:N:0:AACCGCGGAT+CTAGCGCTGT +AATGGAACGGCAGGGTCAGGGTTCTGGATATCTGGTCTAACACTCAGAGTTATTCCTTTTCCAAATGTCAGCTTACTATTGGCTCCAGTGCCCCCGCTTGCTGCACAGAAGTACACTGCAGAGTCTCCAGGCTGGGAGGGCACAATGTGC ++ +:,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:10999:1908 1:N:0:AACCGCGGAT+CTAGCGCTGT +ATTTTTACGGCAGGGTCAGGGTTCTGGATATCAGAAGAGATACACTGAAACAGACGATGAGTAGTTGCACCTCAGCTCCACAGGGGCTTCTTCAAAGACAGGGACTTGGCTGTCAAGCTGGGTCACAGACTGGGCTCTGGTTCCTCCCAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:20437:1908 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTCTTAGGTGGGAACACGTTTTTCAGGTCCTCTGTGACCGTGAGCCTGGTGCCCGGCCCGAAGTACTGCTCGCCGCGCCCCGTAAGTTCACTGCTGGCACAGAGGTACATAGATGTCTGGTTGGTGCTGGCGGACTCCAGAATCAGGGAG ++ +FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:8314:1924 1:N:0:AACCGCGGAT+CTAGCGCTGT +TGAATTGGTGGGAACACCTTGTTCAGGTCCTCTACAACTGTGAGTCTGGTGCCTTGTCCAAAGAAAGCTTCAGTGTTCACACAGTGACAGGGGTCAAGGTGAAAATCACATTGAAGGCATAGAGAGGGGAGAGGGCCCTGGCTAGGATGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,:FFFFF +@A00261:687:HMMCCDSX5:3:1101:16848:1924 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGTGTTACGGCAGGGTCAGGGTTCTGGATATTTGGACTGACCAGAAGTCGGGTGCCAGTTCCAAATACAAGTTTCTGAAAGCCTGTCGACGGCCCTGCACAGAAGTAGGTTCCTACATCTGTAGTCTGGGTGGCTGTGATGTGCAGGGAG ++ +FFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:19090:1924 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTTTTACGGCAGGGTCAGGGTTCTGGATATTTGCTCTTACAGTTACTGTGGTTCCGGCTCCAAAGCTGAGCTTGTAGTCGTTAGTAGGTGTGGCTGCACAGAAGTAGATAGCTGAGTCCCCAGTCTGGGTGTCTGCAATGCGCAGAGAC ++ +FFFFFFFFFFF:FFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF:FFFFFFFFFFF::FFFFFFFFFF,::FF,FFFFFF:F:FFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:29903:1924 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTGTTGGTGGGAACACCTTGTTCAGGTCCTCTAGGATGGAGAGTCGAGTCCCATCACCAAAATGCTGGGGCTGATTGGTCCCTGTCCCCACAATGTTACAGCTTTGTACAAAAACAGCCCCTCCCATGGGTCCGCCCCCAGAGCCTGGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFF:FF,FFF +@A00261:687:HMMCCDSX5:3:1101:20871:1939 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTTTTACGGCAGGGTCAGGGTTCTGGATATTTGGTTGCACTTGGAGTCTTGTTCCACTCCCAAAAGTAAGTGCTCTCCTGCCCGTGTCCCCCCCGTCCGTAGCACAGAAGTAAGAAGCAGTGTCTGCTGCCCGGGAAGCCGTGATCAAC ++ +FFF:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:5394:1971 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGAGGTGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATGGGCCTGTGAGCCACTCATGGTTTGGCTGCTGGCACAGCGATACATGGCCGAGTCCCGCTGCTCTGTGCGCTGGATC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24542:1971 1:N:0:AACCGCGGAT+CTAGCGCTGT +AATCTTGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGTCTCCTGGCCCCCTCCCCATTCACTGCTGGCACAGAGGTACTGAGAGGTATGTGAGGGCCTGGCAGACTCCAGGGTCAGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFF:FFFFFFFFFFFFFFFFFFFF:,F,FFF,F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26205:1971 1:N:0:AACCGCGGAT+CTAGCGCTGT +CTCTCAGGTGGGAACACGTTTTTCAGGTCCTCTAGCACGGTGAGCCGTGTCCCTGGCCCGAAGAACTGCTCATTGTAGGAGCACAGTGGGAAGGGGCTGCCCAGAATTCCTTCTCCTTCCCCACCTCCCTTAGCACAGACAGAGGGCACA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:10881:1986 1:N:0:AACCGCGGAT+CTAGCGCTGT +TTTTTTACGGCAGGGTCAGGGTTCTGGATATTTGGTTTTACATTGAGTTTGGTCCCAGATCCAAAGTAAAATTTGTTGAAGCTGTATCCGTTCACCACACAGAGGTAGGTGGCTGAATCACTGAGCTTGGAGTCTCTGATGAGCAGGGAA ++ +FFFFFFFFFFFFFFFFFFFFFFF:FFF:F:FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:12165:1986 1:N:0:AACCGCGGAT+CTAGCGCTGT +CGAATCACGGCAGGGTCAGGGTTCTGGATATAGGGCTGGATGATTAGATGAGTCCCTTTGCCAAAGGTGAGTCCGTCAGCACCTCCTCCTGAATACAAACCGTAAGCACAGAAATACATCGCGGCATCCCCCAGCTGTGAGTCTGAGATC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFF:F +@A00261:687:HMMCCDSX5:3:1101:24008:1986 1:N:0:ACCCGCGGAT+CTAGCGCTGT +TAAAACACGGCAGGGTCAGGGTTCTGGATATAGGGAATAACGGTGAGTCTCGTTCCAGTCCCAAAGAAGAGGTTGTTTGCCCCAGTTCCTACCCCCTCCTCCACAGCACAGTAGTAAACAGCAGCATCTGCCTCCCGCACCTGGAGGATG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:30942:2002 1:N:0:AACCGCGGAT+CTAGCGCTGT +AAAACCGGTGGGAACACCTTGTTCAGGTCCTCTACAACTGTGAGTCTGGTGCCTTGTCCAAAGAAAGCTTCAGCTGTCCCTCCGCTGCTGGCACAGAGATACATGGCCGAGTCCCCCTGCTCTGTGCGCTGGATCTCCAAGGTGGAGAAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26151:2033 1:N:0:AACCGCGGAT+CTAGCGCTGT +AGTGTAGGTGGGAACACGTTTTTCAGGTCCTCGAGCACCAGGAGCCGCGTGCCTGGCCCGAAGTACTGGGCCCGTCCCACGTCGGACTCAACGCTGCAGAGATATATGCTGCTGTCTTCAGGGCTCATGTTGCTCACAGTCAGAGTTGAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFF:FFFFFFF,FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF:FFFF:FF:FFF +@A00261:687:HMMCCDSX5:3:1101:6144:2049 1:N:0:AACCGCGGAT+CTAGCGCTGG +TTTTTTACGGCAGGGTCAGGGTTCTGGATATATGGGTGTACAGCCAGCCTGGTCCCTGCTCCAAAAATCAGCTTATTGTAGCCACCTTCACTCACAGCACAGAAGTACTCAGCCGCGTCGCTCATATGGGCTGAGGGTTTCGTCAGGTGG ++ +FFFF:FFFFFFFFFFFFFFFFFFFFFFF,FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF,FFF +@A00261:687:HMMCCDSX5:3:1101:15727:2080 1:N:0:AACCGCGGAT+CTAGCGCTGT +TAAATCACGGCAGGGTCAGGGTTCTGGATATTTGGATGGACAGTCAAGATGGTCCCTTGTCCAAATGTCAGCTTTCCATAGCTAGTACCACCAGCGCGTCCCGAGCTTGCTGCACAGAAGTACACTGCAGAGTCTCCAGGCTGGGAGGGC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF diff --git a/examples/input/TINY_2.fq b/examples/input/TINY_2.fq new file mode 100644 index 0000000..c6ba360 --- /dev/null +++ b/examples/input/TINY_2.fq @@ -0,0 +1,424 @@ +@A00261:687:HMMCCDSX5:3:1101:24948:1047 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCATATGTCGTGATCTGTTCAAATGGCAATGAGGAGACATCACCTGGAATGTTAGGCAGTGCCTAACTGGGGGATGGACAGACAATGGGCAGTGCCAACCCATAGGGTGGATACAAAAGACAGGCAAG ++ +:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:27272:1063 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGATCCTGAGGACAGTGCCTGGAGGTGAGAAGGAAGCCCCCGGCCTGGTCCATACCCCACCACCAACTTGCATAATGGGGGGTGATGTCACCCACCCTCCACTCCCCTCAAAGGAGCAGCTGCTCTGGTGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFF,FF,FFF:FF,FFFFFF,,FFFFF +@A00261:687:HMMCCDSX5:3:1101:5050:1094 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTTCCGGTCGTGATAAAGTGGAGTATCTCTGCTATGTTCATTTCTTTTTGGATTGAAAATTTTAATCCTCAGTGAACCAGGGCAGAAAAGAATGATGATATCCTTGAGAGTTTTACTGGTGATCCTGT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF,FFFFF:FFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFF:FFFFFFF::F:FFFF,FFFFFFFFFF:FFFFFFF,F:FFFFFFFF:F:FFFFFFFF:FFF +@A00261:687:HMMCCDSX5:3:1101:18331:1110 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAACACCCTGGTTACCCGTCGTGATGCCTTTGAAACCTGGGGAGGAAGACATAGAAAAGCAACTGCTCTATCATTGGCTGACAAGATCCTGACAACAACACAATGAGAGCTGCACTCCTTGCAAGTCTGTACTAGTAGC ++ +FF,FFF,FFFF:,F,FFFFFF:FFFF::FFFFF:F,F,F::FFFF,FFFFF:F,FFFFFF:FFFFF:FF,F:FFFFFFFFFFF,FFFFF:FF:FFFFFFFFFFF:FFFF:F:FFFF,:FFFFFF::::FF,,F:FFFF,,FFFFFFF,:: +@A00261:687:HMMCCDSX5:3:1101:8684:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCTAAGTCGTGATTTTTATGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTC ++ +FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF::FFF,FFFFFFFFFFFFFFF,FFF:F:FFFFF:::FFF:FF,F,FFF:F,FFFFF,FF,FF:FFF:FFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:4:1101:8684:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCTAAGTCGTGATTTTTATGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTC ++ +FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF::FFF,FFFFFFFFFFFFFFF,FFF:F:FFFFF:::FFF:FF,F,FFF:F,FFFFF,FF,FF:FFF:FFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:5:1101:8684:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCCAAGTCGTGATTTTTATGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTC ++ +FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF::FFF,FFFFFFFFFFFFFFF,FFF:F:FFFFF:::FFF:FF,F,FFF:F,FFFFF,FF,FF:FFF:FFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:5:1101:8684:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCACAGTCGTGATTTTTATGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTC ++ +FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF::FFF,FFFFFFFFFFFFFFF,FFF:F:FFFFF:::FFF:FF,F,FFF:F,FFFFF,FF,FF:FFF:FFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:5:1101:8684:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCGGGGTCGTGATTTTTATGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTC ++ +FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF::FFF,FFFFFFFFFFFFFFF,FFF:F:FFFFF:::FFF:FF,F,FFF:F,FFFFF,FF,FF:FFF:FFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:6:1101:8684:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCGATAGTCGTGATTTTTATGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTC ++ +FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF::FFF,FFFFFFFFFFFFFFF,FFF:F:FFFFF:::FFF:FF,F,FFF:F,FFFFF,FF,FF:FFF:FFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:7:1101:8684:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCTAAGTCGTGATTTTTATGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTC ++ +FFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF::FFF,FFFFFFFFFFFFFFF,FFF:F:FFFFF:::FFF:FF,F,FFF:F,FFFFF,FF,FF:FFF:FFF:FFFFFFFFFFFFFFFFF,FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:29604:1125 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCATTCTGTCGTGATTAAAATTGACAGTGACCCTGATCTGGTAAAGCTCCCATCCTGCCCTGACCCTGCCATGGGCACCAGCCTCCTCTGCTGGATGGCCCTGTGTCTCCTGGGGGCAGATCACGCAGA ++ +FFFFFFFFFFFFFFFFF:FFFF:FFFFFFF:FF:F:FFFFFFFFFFFFF:FFFFFFFFF,FFFFFF:FFFFFF:FFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFF +@A00261:687:HMMCCDSX5:3:1101:17300:1141 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGTTGCCGTCGTGATATTAATAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTCTCTCAACATCCGAGCTGGGTTATCTGTAA ++ +FFFFFFFF:FFFFFFFFFFFFFF:FF::FFFF,F,FFFFFFFFFFFFFFFFFF:FFF,FFFFFFFF:FF:FFFFFFFFFFF:F:F:,:FFFFFFFFF:FFFFFFFFF::FFFF:FFFFFFFFFFFF,,FFFFFFFF,FFF,:FF:FF,FF +@A00261:687:HMMCCDSX5:3:1101:26793:1141 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGTAATGTCGTGATTTCAACACTCTGAAGGTGCAGCCTGCAGAACTGGAGGATTCTGGAGTTTATTTCTGTGCCAGCAGCCCTCACGGGGGGTTCTCCTACAATGAGCAGTTCTTCGGGCCAGGGACA ++ +FFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:28854:1141 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCTTTCGTCGTGATATCCATGATGGGCTGAGGCTGATCCATTACTCATATGGTGTTCAAGACACTAACAAAGGAGAAGTCTCAGATGGCTACAGTGTCTCTAGATCAAACACAGAGGACCTCCCCCTC ++ +FFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFF,FFFFF:FFFFFFFFFFF,FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF:,FFFFFF:FFFF +@A00261:687:HMMCCDSX5:3:1101:22480:1157 2:N:0:CACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCCGTTGTCGTGATCCTTATGCCCTTTGTGCCATGGGCTCCAGGCTGCTCTGTTGGGTGCTGCGTTGTCTCCTGGGAGCAGGCCCAGTAAAGGCTGGAGTCACTCAAACTCCAAGATATCTGATCAAA ++ +:FFFFFFFFFFFFFFFFFFFFFFFFF::FFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFF,FFFFFFFFFFFF,:FFFFFFFFFFFFFFFFFFFFFF:,:FFFFFFFFFFFFFFFFFFFF:F:FF,F,,:FFF:F,,F +@A00261:687:HMMCCDSX5:3:1101:28926:1172 2:N:0:CACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTGGGTGTCGTGATCCTTTGGTCAGGAGCAGTGACATCACAGGAAAAACCACCAACCAAGGCCAAGGAGACCAGAGCCCAGCACCTCACCCAGAGGACCCCAGTCAGAGGCCCCATCTCAGACCCGAG ++ +FFFFFFFFFFF,FFFFFFFF,FFF:FFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFF:FFFFFFFF,FFFFFFFFFF,FFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:28592:1188 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCTATGTCGTGATGAATATGAAGAGACTTGCCTAGCCCAACCTTCCTCACGCTCACTATTCTCAAGACCTGGGTTCCAGCCACTTTCCTACTGGCCCCGAGGAGAATTTCCAAAGAGACGCCTGCAG ++ +FFFFFF:F:FFFFFFFFFFFFFFFFFFFF:FF,FFFFFFFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFF:FFFFFF:FF,FFFFFF:FFFFFFF:FFFF:FFFFFF:FF,FFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:4417:1219 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTTCCGGTCGTGATAAAGTGGAGTATCTCTGCTATGTTCATTTCTTTTTGGATTGAAAATTTTAATCCTCAGTGAACCAGGGCAGAAAAGAATGATGATATCCTTGAGAGTTTTACTGGTGATCCTGT ++ +FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF:FF:FF:F,FFFFFFFFFFFFF,FFFF,FFFF,FFFFFFFF::FFFFFF::FFFFF:F:FFFFFFF,F:FFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:3794:1235 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTAGCCGTCGTGATGTGCTATGACGTTTTGCCCAGGAAAACACACTTGATAACTGAAGGATGATGAAGTGTCCACAGGCTTTACTAGCTATCTTTTGGCTTCTACTGAGCTGGGTGAGCAGTGAAGAC ++ +FFFFFFFFFFF:FF,FFFFF:FFFFFF:FFFF,FFF,FFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFF,FFF,FFFF,FFFFFFFFF:FFFFFFFFFFFFFFFF:F:FFFFFF +@A00261:687:HMMCCDSX5:3:1101:23737:1266 2:N:0:AACCGCGGAT+CGAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTAGTCAGTCGTGATTGTCATTGAGTCATCCCTCCTCGCTGGTGAATGGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGACGATGCTGCTGCTTCTGCTGCTTCTGGGG ++ +F,FF,FFF:FFFF,FFFFF,FF:,FF,:FFFFFFF,,F:FF::FFFFFFFFFFFFF,FFF,FFFFF,FFF:FF:,FFFFFF:::,FF:FFFFFFFFFF,FFF:,FF:FFFF::,FFFFFFF,,:FFFF:FFFFF,F,:FFFF:,FFFFF, +@A00261:687:HMMCCDSX5:3:1101:29966:1282 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCCTATGTCGTGATCTTTACGCTTTTCTCAGGGGAGAGGCCATCACTTGAAGATGCTGAGTCTTCTGCTCCTTCTCCTGGGACTAGGCTCTGTGTTCAGTGCTGTCATCTCTCAAAAGCCAAGCAGGG ++ +FFFF:FFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:20627:1297 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCTCGCCGTCGTGATGTTTATGAAGCCCTGAATCAGGTGCAGTGCTGCCTGCCCCACTGTGCCATGGGCCCTGGGCTCCTCTGCTGGGTGCTGCTTTGTCTCCTGGGAGCAGGCTCAGTGGAGACTGGA ++ +:FFFFFFFFFFFFFF:FFFFFFFFFFFFF,FF:FFFF:F:FFFFFFFFFFFFFFFF:FFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFFFFF:,FFFFFF +@A00261:687:HMMCCDSX5:3:1101:19027:1313 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCACTTTGTCGTGATACGATCTTGTGAGCCATTCTCCATATTTCAGATATAAGATTTCAGTTCTCAGTGAGTCTAAGTGACAGAAGGAATGGAGACCCTCTTGGGCCTGCTTATCCTTTGGCTGCAGCT ++ +FFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF,:FFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFF,FF:FFF:FFF:F:FFFFFFF,FFFFFF:FFFFFFFFFFF:FFFFF:FF +@A00261:687:HMMCCDSX5:3:1101:6876:1344 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTATCAGTCGTGATGCATGTGGCTGGCCTGGTCAGAAGGCAGCTTTGTCAGGAGGTGGAAGGATGAACTAGTTAACCTCAGAAGAGCCTTTCCACTCCGACCATGCTGGGAATAGACCGCCATAATAA ++ +FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFF +@A00261:687:HMMCCDSX5:3:1101:29170:1344 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTGCCTGTCGTGATTTCTATGACCTGAAGTGTCTACCTTCTGCAGACTCCAATGGCTCAGGAACTGGGAATGCAGTGCCAGGCTCGTGGTATCCTGCAGCAGATGTGGGGAGTTTTCCTTCTTTATGT ++ +FFFFFFFFFFFFFFFFFFFFF:FFFFF:FFFF::FFFFF,FF,FF,:FF,FFFFFFFFFFFFFFFFFFFFFF,F:FFFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFF:FFFFFFFF,FF,FFFF:FFF:FF +@A00261:687:HMMCCDSX5:3:1101:21260:1360 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGCGCGGTCGTGATAATCATGTTCCTTCCCTCACCCACATGAAGTGTCTACCTTCTGCAGACTCCAATGGCTCAGGAACTGGGAATGCAGTGCCAGGCTCGTGGTATCCTGCAGCAGATGTGGGGAGT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF:F,FFFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFF,FF:FF,FFFF,FFFF,F,FFFF,F,FF:F,F,FFFFFFF,FFFFF,FFFFFFFFFFF:F,FFFFFF +@A00261:687:HMMCCDSX5:3:1101:18611:1376 2:N:0:ACCCGCGGAT+CTAGCGCTGT +TTCGTGACAGGGAAAACCCTGGGGGCCGGTCGTGATCCCTAATCTCTGAGAGGGGCTGTGCTGAGAGGCGCTGCTTGTCGTCTGGGCGGAGGACTCCTGGTTCTGGGTGCTGGGAGAGCGATGGGGCTCTCAGCGGTGGGAAGGACCCGA ++ +,F,,FFF,,:,,FFFF,,,F,,F,,,F,FFF,FFF,,:,F:F,,,F:,F,:,,FFFFF,F,,:F:F,F,,:,,FF,F,,,F,FF,,FF:F,FFFF:,F,:FF,F,FFFFF,,FF,,:F::FF::FFFFF:,FFFFFFF,:F,:,FF:FFF +@A00261:687:HMMCCDSX5:3:1101:7446:1391 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGTTAGTCGTGATTTTTGGACTCTTCTGGTCCCCACAGACTCAGAGAGAACCCACCATGGTGCTGTCTCCTGCCGACAAGACCCACGTCCAGGCCGCCTGGGGTAAGGTCGGCGCGCCCGCGGGCGAG ++ +,FFF:FFF::F:F:F::FFFFFFF::,,,::::F,,,FFF:FF,F:,:FF,FFFFFFFF,FFFFFFFFFFFFF,F,FFFFFFFFF:FFFF:F,,FFF,FFF,::F,FF,,FF:,FFFFF:,,FFFF,F:F,,:FF,FFFFF,:,,FF,,F +@A00261:687:HMMCCDSX5:3:1101:8024:1391 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGACTTACGTCGTGATCGAATGGCTCTGCTGTGCGGTTCTCTGTCTCCTGGGAGCAGTTCCCATAGACACTGAAGTTACCCAGACACCAAAACACCTGGTCATGGGAATGACAAATAAGAAGTCTTTGAA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFFF +@A00261:687:HMMCCDSX5:3:1101:24641:1391 2:N:0:ACCCGCGGAT+CTAGCGCGGT +GTCGTGACTGGGAAAACCCTGGTCCCCAGTCGTGATTTTGTTGAAAGATCCCATACTGCCCTGACCCTGCCATGGGCACCAGCCTCCTCTGCTGGCTTGCCCTGGGTCTCGTGGGGGCAGATCACGCAGATACTGGAGTCTCCCAGGACC ++ +FFFFFFFFFFF:FFFFF:F:FF,F:FF,FFFFFF,:FFFFFF:FFF,,FF:FF:,F,F,:F,,,F,FF,:,,F,F,F,,:F:::F:F,,:FF,,:,F,F,,F,,,,:F:F,:,:,FFFFF,F,F,,F,F,FF,F,FF:FFFF,::,F:FF +@A00261:687:HMMCCDSX5:3:1101:26585:1407 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCTCAAAGTCGTGATCTCCTGGCTTCTCTCACCATGCCAGGTTCACTTCACAGTACAGAGTCCTGAAAATAAAGAAGAAAATTTTTTTTTATCTAGACAAAGAACCAAACATGTCACTGTCTAGCCTGC ++ +FFFFFFFFFFFFF:FFFFF:FF:FFFFFF:FFFF,:FFFFFFFFF,FFF::FFFFFFFFFFF:FF:F,FFFFFFFFFFFFFFFF,F:FFFFF,:F::FFFFFFFFFF:FFF,:F:,FF,FFFFFF:FF:FFFF:FF:F:,FFF,:::FFF +@A00261:687:HMMCCDSX5:3:1101:30653:1407 2:N:0:AACCGCGGAT+CTACCGCTGT +GTCGTGACTGGGAAAACCCTGGTTGCCAGTCGTGATGTCTATGAGCTGGTTGGGAAGACTGGAAGACCACCTGGGCTGTCATTGAGCTCTGGTGCCAGGAGGAATGGACAAGATCTTAGGAGCATCATTTTTAGTTCTGTGGCTTCAACT ++ +:FFF,FFFFFFFFFFFFFFFFFFF,FFFFFF:F,FFFFFFFF:FFFFFFF:,FFFF:FFFFFFFFFF,FFFF,,FFFF:FFFFFFF:FFFFFFFF,FFFF:FFFFFFFFFF:FF:FFFFFFF:FFFF,FFFFFFFFFFF:FF:FF:FFFF +@A00261:687:HMMCCDSX5:3:1101:5638:1423 2:N:0:CACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTTCGGTCGTGATATGAAAGACCCCTCCATCCTGCAGTTCCTGCCATGAGCCTCGGGCTCCTGTGCTGTGGGGCCTTTTCTCTCCTGTGGGCAGGTCCAGTGAATGCTGGTGTCACTCAGACCCCACA ++ +FFF:FFF:FFF:FFFFF,FFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFF::FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFF,FF::FFFF:F::FFFFFFF:F,F +@A00261:687:HMMCCDSX5:3:1101:24713:1423 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCTGAGGTCGTGATCGATTGGCAATTTTTAGTCTCTTGAAAGCAGGACTCACAGAACCTGAAGTCACCCAGACTCCCAGCCATCAGGTCACACAGATGGGACAGGAAGTGATCTTGCGCTGTGTCCCC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFF,FFFFFFFFFFFFF: +@A00261:687:HMMCCDSX5:3:1101:24379:1438 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGTACTGTCGTGATTCAGATGCTCAGCCTCCAGACCCCAGGACCGGCAGTTCATCCTGAGTTCTAAGAAGCTCCTTCTCAGTGACTCTGGCTTCTATCTCTGTGCCTGGGAACAGGGGGCCACTGAAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFF,FFFFFF +@A00261:687:HMMCCDSX5:3:1101:30309:1438 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGTATTCGTCGTGATGGTCCGGTTTTTCAAAGCCAGGGACTTGTCCAGCCCAACCTCCCCATTGCTCCTAGCTCCTGAGGCTCAGGACCCCTGGCTTCTGTCCTCCCTGCTCAGGGTCCTGAAGCGTTG ++ +,FF,FF:FFFFFFFFF,FFFFF:FFFFFFFF:FFFF:FFFF,FF::FFFFFFFFFFFFFFF:,FFFFFF:FFF,F:FFFFFFF,FFFFFFFFFFF:FFFF:FFFFFF,FFFFFFF,:F:FFFFFFFFFFFFFFFFF,FF,FF,:FFFFFF +@A00261:687:HMMCCDSX5:3:1101:8133:1454 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGAAAACCCTGGCCTTAAGTCGTGATTCATTGGATCCTTCCCTGACCCTGCCATGGGCACCAGGCTCCTCTGCTGGGCGGCCCTCTGTCTCCTGGGAGCAGAACTCACAGAAGCTGGAGTTGCCCAGTCTCCCAGATATAA ++ +FFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFF,FFFFFFFFFFF::FFFFFFFFFFFFFFFFFFFFFFFFF:FFFF:FFF:FFFFFFFFFFFFFFFFFFFF:F:FFFF,FFFFFFFFFFFFFFF,FFFFFF,FFFFFFFFFFFFF,FFF +@A00261:687:HMMCCDSX5:3:1101:24532:1454 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTCCTCGTCGTGATGGCTACGGCATTTTGGCCATGGCTTTGCAGAGCACTCTGGGGGCGGTGTGGCTAGGGCTTCTCCTCAACTCTCTCTGGAAGGTTGCAGAAAGCAAGGACCAAGTGTTTCAGCCT ++ +::FFFFFFFFFFFFFFFFFFF:FFF:FFFFFFFFFF:F:FFFF,:FFFFFFFFFF::FFFFFFFF,FFFF:FFFFF:FFFFF:FFFFFFFFFF:FF:FFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFFFFFFFFFF:FF +@A00261:687:HMMCCDSX5:3:1101:28890:1454 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCCTCAGTCGTGATATCTAGGATTATTTGGGTAACACACTAAAGATGAACTATTCTCCAGGCTTAGTATCTCTGATACTCTTACTGCTTGGAAGAACCCGTGGAGATTCAGTGACCCAGATGTAAGGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFF:,:FFFF:FF:FFFFFFFFFFFFFF:FF::FFF,FFFFFFFFFFFFFF,FFFFFFF:FFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFF:F +@A00261:687:HMMCCDSX5:3:1101:4869:1470 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCTTTGAGTCGTGATCATAAGGATCTGGTAAAGCTCCCATCCTGCCCTGACCCTGCCATGGGCACCAGGCTCCTCTTCTGGGTGGCCTTCTGTCTCCTGGGGGCAGATCACACAGGAGCTGGAGTCTCC ++ +FF:FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFF:FFFFFFFF::FFFFFFFFFFFFFFFFFF:FFFFFFFFFF:FFFFF:FFFFFF:FFFFFFF:FFFFFFFFFFF:FFFFFFFF,,FFFFFFFFFFFFFF:F:FFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24542:1470 2:N:0:AACCGCGGAT+GTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTCCTCGTCGTGATGGCTACGGCATTTTGGCCATGGCTTTGCAGAGGACTCTGGGGGCTGTGTGGCTAGGGCTTCTCCTCAACTCTCTCTGGAAGGTTGCAGAAAGCAAGGACCAAGTTTTTCAGCCT ++ +,FF:FFFFF,FFFFFF:F:FFF:FFF:,FFFFFFFFFF:FF:F::FFFFFFFF::F:F:,FFFFFF:F,FF:FFF,F,F:,:FFFFFFFFF,F,F,,F:FFFF:FFFFFFFF:,FFF:FFF::FFFFFFFFFFF:,FF::,FFF::F,FF +@A00261:687:HMMCCDSX5:3:1101:6777:1485 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGTTATAGTCGTGATCAAATCGCTTTTCTCATGGGAGAGGCCCTAACTTGAAGATGCTGAGTCTTCTGCTCCTTCTCCTGGGACTAGGCTCTGTGTTCAGTGCTGTCATCTCTCAACAGCCAAGCAGGG ++ +FFFFFF:FFF,FFF:FFFFFFFF:F:F::FFF:F,FFFF:FFFF:FFFFFFF,FFFF:FFF,F,,,FF:F:,F,FFFFF,F,F:,F:,F,,,F:F:,:,F:,FFF,F:,FF,,FF:FFFFFFFFFFF,F,FF:FF:F,FFFFF:FFF:FF +@A00261:687:HMMCCDSX5:3:1101:31331:1485 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCAAGGGTCGTGATCTCTATGTCACAGTGACCCTGATCTGGTAAAGCTCCCATCCTGCCCTGACCCTGCCATGGGCACCAGCCTCCTCTGCTGGATGGCCCTGTGTCTCCTGGGGGCAGATCACGCAG ++ +FFF:FFFFFFFFFFFFFFFFFFFF:FFFFF:FF:FFF:FF:FFF::FFFFFFF:FFFF:FFFFFF,FFFFFF::FFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF::FFF:FFFFFFFFF,FFFFFFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18629:1501 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGTCGAGTCGTGATGTCATGGCTTTTCTCAGGGGAGAGGCCATCACTTGAAGATGCTGAGTCTTCTGCTCCTTCTCCTGGGACTAGGCTCTGTGTTCAGTGCTGTCATCTCTCAAAAGCCAAGCAGGG ++ +FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF::,:FFFFFFF:FF:FFFFFF,:FFF,:,:FFFFF,FFFFFFFFF:FFFF,FFFF,FFFFFFFFFFFFFFFF::F,F:F:,F:FFF,,F:F:FF:FFFF:F,FFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:31448:1501 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCAAGGGTCGTGATCTCTATGTCACAGTGACCCTGATCTGGTAAAGCTCCCATCCTGCCCTGACCCTGCCATGGGCACCAGCCTCCTCTGCTGGATGGCCCTGTGTCTCCTGGGGGCAGATCACGCAG ++ +FFFFFFFFFFFFFFFFFF:FFFFFFFF:FF:F,FFFFFF:,FFF:FFF::FFFFFFF:FF,F:FF:FFFFF:FFFFFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFFFFFFF::FF:FFFFFFFFFFFF,FFFFFF:F::FFF:FFF +@A00261:687:HMMCCDSX5:3:1101:1768:1517 2:N:0:AACTGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGTTGAGTCGTGATGGAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGATGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGCAGGCTCCGGGCTTGGTGCTGTCGTCTCTCAA ++ +FFFFFFFFFFFFFFFFFF:FFFFFFFF,FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF:FFFFFFFFF::F,:FFFFFFFFFFFFFF:FFFF +@A00261:687:HMMCCDSX5:3:1101:9607:1532 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCATCTTGTCGTGATACCATAGACAGATACAGAAGACCCCTCCGTCATGCAGCATCTGCCATGAGCATCGGCCTCCTGTGCTGTGCAGCCTTGTCTCTCCTGTGGGCAGGTCCAGTGAATGCTGGTGTC ++ +FFFFFFFFFFFFFFFFFF:FFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFF::F,:FF:F:FFFFFFFFFFFF::F::FFFFFFFFFFFFFFFF:FFFFF:FFFFFF:F,F,FFFFFFFFFFF:FFFFFF:FFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:10854:1532 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGAGTGACGTCGTGATTTGTATGGTTTTTCTGCTGTGGGTACGTGAGCAGGAAACATGGAGAAGAATCCTTTGGCAGCCCCATTACTAATCCTCTGGTTTCATCTTGACTGCGTGAGCAGCATACTGAAC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F:FFFFFFFFFFFF:FFF,FFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF::FFFFF +@A00261:687:HMMCCDSX5:3:1101:16866:1548 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTACCTTGTCGTGATCTTTGAGATTATTTGGGTAACACACTAAAGATGAACTATTCTCCAGGCTTAGTATCTCTGATACTCTTACTGCTTGGAAGAACCCGTGGAGATTCAGTGACCCAGATGGAAGGG ++ +FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26214:1548 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCATGCTGTCGTGATGCTGACAGCAAGGCTTTGAAAAATAAAGATACACAGATAAATTATTTGGATAGATGATCAGACAAGCCTCAGTAAAAACAGCCAAGACAATCAGGATATAATGTGACCATAGGA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFF:FFFFF:FFFFFFFFF:FFF:FFF,FFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF::FFFFFF:,FF:FF:FFFFFF,FFFF:FFFFF +@A00261:687:HMMCCDSX5:3:1101:30120:1548 2:N:0:AACCGCAGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCCAAAGTCGTGATCTCTAGGACCAGCCCCAAGCTAGGAGATCCTGCCATGGGCTTCAGGCTCCTCTGCTGTGTGGCCTTTTGTCTCCTGGGAGCAGGCCCAGTGGATTCTGGAGTCACACAAACCCC ++ +FFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:7238:1564 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCCTTTGTCGTATACCTACGGAGAACCAAGCCCTGAATCAGATGCAGTGCTTCCTGTCCCTCTGTGCCATGGGCCCCGGGCTCCTCTGCTGGGCACTGCTTTGTCTCCTGGGAGCAGGCTTAGTGGAC ++ +FFFF:FFFFFFFFF:FFFFFFFFFFFF:FFFFFFFFFF:FFFFFFFFFFFF:FFFF:FFFFFFFFFFF:FFFFFFFFFFF,FFFFFFFFFFFF:FFFFF:FFFFFFFFFFFF:,FFFFFFFF,FFFF:FFFFF:FFFFFFF:FFFFFF:F +@A00261:687:HMMCCDSX5:3:1101:13964:1564 2:N:0:CACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTGTGCGTCGTGATTACTTTTCCCTGCCGACAGAAAGTCCAGCACTCTGAGCCTGCCCCGGGTTTCCCTGAGCGACACTGCTGTGTACTACTGCCTCGTGGGTGATCCTTACAACAGACTCGCTTTTG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F,FFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24777:1564 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTCACTGTCGTGGATTCCTGCTGTGATCCTGCCATGGATACCTGGCTCGTATGCTGGGCAATTTTTAGTCTCTTGAAAGCAGGACTCACAGAACCTGAAGTCACCCAGACTCCCAGCCATCAGGTCAC ++ +,FFFFFFFF:FFFFF:FFFFFFFFFF:FFFFFF,FFFF:FFFF::FFF,:FF,:::FFFF:FF::,,:F::FFFFFFFFFFFFFFFFFFFF,FF,FF:FF,FFFFFFFFFFFFF:FFFFFF::FFFFFFFFFFFFFFFF,FFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:18195:1595 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTACGAGTCGTGATATCTACGACCAAAACAAGAGACTTGCCTAGCCCAACCTTCCTCACGCTCGCTATTCTCAAGACCTGGGTTCCAGCCACTTTCCTACTGGCCCCGAGGAGAATTTCCAAAGAGAC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFF,FFF,,FFFFFFFFFFFF:FFF,FFFFFFFFFFFFFF:FFFF,FFFFFFFFF,, +@A00261:687:HMMCCDSX5:3:1101:18484:1626 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGTCGAGTCGTGATGTCATGGCTTTTCTCAGGGGAGAGGCCATCACTTGAAGATGCTGAGTCTTCTGCTCCTTCTCCTGGGACTAGGCTCTGTGTTCAGTGCTGTCATCTCTCAAAAGCCAAGCAGGG ++ +FF,FFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF,FFFFF,:FFFFFF:F:,:FF::FF:FFF::FFF,:F:FF,FFFFF,:FFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:11966:1642 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCCCCTGTCGTGATGAAGATTGACAGACACAGTGCTGCCTGCCCCTTTGTGCCATGGGCTCCAGGATGCTCTGTTGGGAGCTGCTTTGTCTCCTGGGAGCAGGCCCAGTCAAGGCTGGAGTCACTCAA ++ +FFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:F:FF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFF:FFFFFF,FFFFF,F:FFFFFFFFFFFFFFFFFFF:FF,FFFFFFFFFFFFF:FFF, +@A00261:687:HMMCCDSX5:3:1101:26874:1658 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCTCTCAGTCGTGATGGCACGCTGACTGAACTTTGAGCTTCAGGCAGCACAACTCACATTTGTGCAAAGAGCTACATGCCACATGCTGTTCTCCAGCCTGCTGTGTGTATTTGTGGCCTTCAGCTACTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFF:FFF,FFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:27597:1658 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCTTTTAGTCGTGATATTCAAGGTTTCCTTCCCTCACCCACATGAAGTGTCTACCTTCTTCAGACTCCAATGGCTCAGGAACTGGGAATGCAGAGCCAGGCTCGTGGTATCCTGCAGCAGATGTGGGGA ++ +FFFF::FFFFFFFFF::FFFFFFFFFFFFFFFFFFFFF:FFFFFF,FFFFFFFF,FFF:FFFFFFFF:F,FFFFFFFFFF,FFFFFFFF:F:FFFFFFFFF,FFFFFFFFF,F:,:FFFFFF:FFFFF:FFFF,FFFFFF:FFF::F:FF +@A00261:687:HMMCCDSX5:3:1101:31901:1658 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGTATTCGTCGTGATGGTCCGGTTTTTCAAAGCCAGGGACTTGTCCAGCCCAACCTCCCCATTGCTCCTAGCTCCTGAGGCTCAGGACCCCTGGCTTCTGTCCTCCCTGCTCAGGGTCCTGCAGCGTTG ++ +FFFFFFFFFF:FFFF,FF:FFFFFFFFF:FFFFFFFFFFF,FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFF:FF,FFFFFF,FFFFFFFFFF,FFFFFFF,FFFFFFFF,F:FFFF:FFFFFFFFFFFFFFFFFFFFFFF:F +@A00261:687:HMMCCDSX5:3:1101:27118:1673 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCATGCTGTCGTGATGCTGACAGCAAGGCTTTGAAAAATAAAGATACACAGATAAATTATTTGGATAGATGATCAGACAAGCCTCAGTAAAAACAGCCAAGACAATCAGGATATAATGTGACCATAGGA ++ +FFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:7283:1705 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGAATCTCGTCGTGATGGAGGCAGTGGTCACAACTCTCCCCAGAGAAGGTGGTGTGAGGCCATCACGGAAGATGCTGCTGCTTCTGCTGCTTCTGGGGCCAGGCTCCGGGCTTGGTGCTGTCGTCTCTCA ++ +:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF:FFFFF +@A00261:687:HMMCCDSX5:3:1101:3332:1720 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCTACCGTCGTGATTCTACGGAGAGTGTTTATCTTGTGAGCCATTCTCCATATTTCAGATATAAGATTTCAGTTCTCAGTGAGTCTAAGTGACAGAAGGAATGGAGACCCTCTTGGGCCTGCTTATCC ++ +FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFF:F,FFFFF,F,F,F:,:FF:,FFFF:FFFFFFFF:FF,FF:FFFF,FF::,,,FFFFFF:F::FFF,FF:,:,,::FF:,FF,,FF,FF:F,,:FFF,::FFFFFF:FFF:FF:F +@A00261:687:HMMCCDSX5:3:1101:24722:1720 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTTATCGTCGTGATTTTAAGGAAAAAGACCTTGCTGAGGGCCCCATGGGCAAGAACAAGTGTGAAGAACCCTACTATGGTTTTTTGCCCAATGCCCTCCACAAAGCCAGCGTTTGAAGGAAAACTGTG ++ +:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFF:FF::F:FFFFFFFFFFFFFFFFFFFF,FFF:FF:F:FFFF:FFF:FFFFFFF:FFFFFFF::FFFFFFFFFFFFFFFFF::FFFFFFFFF,FFFF:FF::,:: +@A00261:687:HMMCCDSX5:3:1101:24975:1720 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTATTACGTCGTGATTTGCTTGGGCTTGTAAGCTCCTTCATCTGGAAATGTGATTTACCTGGGTCCTGCCATGGTTTCCAGGCTTCTCAGTTTAGTGTCCCTTTGTCTCCTGGGAGCAAAGCACATAGA ++ +FFFFFFFFFFFFFF::FFFFFFF:FFF::FFFF:FFFFF,FF:FFFF,FFFF:FFFFFFFFFFFF,,FFFFFFFFFFFF:FFFFFFFFF,FFFFFFFF::F:FFFFFFFFF,FFFFFFFFFFF:F:F:FF:FFFF:FFF:FFFFFF:FFF +@A00261:687:HMMCCDSX5:3:1101:14172:1736 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGTTTTGTCGTGATTGCCCCGGTCTACCTTCTGCAGACTACAGTGGCTCAGGAACCGGGGATGCAGTGCCAGGCTCATGGTATCCTGCAGCAGATGTGGGGAGCTTTCCTTCTCTATGTTTCCAGGAA ++ +FFFFFFFFFFFFFF:FFFF,F:FFF::FFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF:FFFFF:FFFFFFFFFFFFFFFFFFFFFF,FF:FFF:FFFFFFFFFFFFFFFFFF:FFF:,F:FFF::FFFFFF:F:FFF::,,FFF +@A00261:687:HMMCCDSX5:3:1101:18041:1736 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGAGCGTAGTCGTGATCTGTAGGAGCTCTTATGGCTGGAGATTGCAGGTTTATGACTGATCCTATTTGGGAAGAACAATGATGGCAGGCATTCGAGCTTTATTTATGTACTTGTGGCTGCAGCTGGACTG ++ +FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFF:,F,:FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF:,FFFFFFFFFFF:FFFFF:F:FFFFFF:F,FFFF:FFFFF:FFFF:FFFFFF:FFFFFFFFFFFFFFFFFFFF::FF +@A00261:687:HMMCCDSX5:3:1101:28800:1736 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTATTAGTCGTGATTTATATGGGGACTCCTGGTTCTGGGTGCTGGGAGAGCGATGGGGCTCTCAGCGGTGGGAAGGACCCGAGCTGAGTCTGGGACAGCAGAGCGGGCAGCACCGGTTTTTGTCCTGG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFF::FFFF:F:FFFF:F:FFFF:FFFF,FF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:FF,FFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26350:1752 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCTCTTGTCGTGATATGCAAGATTCTAGGTGCATTTTTTAAGGGTTTAAAATTTGAATCCTCAGTGAACCAGGGCAGAGAAGAATGATGAAATCCTTGAGAGTTTTACTAGTGATCCTGTGGCTTCAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFF:FFFF,FFFFFFFFFFF,F,FFFFFFFF:FF,,FFF:FFFFFFF,FF,FFFFFF,FFFFFFFFF,F:FF:FFFFFFF,FFF:FFFFF +@A00261:687:HMMCCDSX5:3:1101:3287:1767 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTTCATGTCGTGATAACAAGGATAACTGGAGGTTGCAGGTCAATGACTGATCTTAATTGGGAAGAACAAGGATGACATCCATTCGAGCTGTAGTTATATTCCTGTGGCTGCAGCTGGACTTGGTGAAT ++ +FFFFFFFFF:FFFFFFFFFFF,FFFFFFFFFF,FFFF,:F,FFFFFFF,,,FFFFFF,FFFFFFF:,FF,FF:FF,FF,FF:FFFFFFFFF:FFF:FFFFFF:,:FFFFFFFFF,F:FFFFFF:FFFFFFFFFF:FF:F,F,FF,FFF:F +@A00261:687:HMMCCDSX5:3:1101:26115:1783 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGTCATGTCGTGATTAAATGGACCTCTCAGCGGCAGTGGAAACCACAGCCTAGTCCTCTCACCACTGCAGACCAGAATCCTGCCCTGGGCCTTGCCTGGTCTGCCTCACTCTGCCCTGGGCTGCAGGC ++ +FFFFFFFF,:FFFFFFFFFFF:FFFFF:FFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF:FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF:F,FFF,FFF::: +@A00261:687:HMMCCDSX5:3:1101:26295:1783 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACGGGAAAACCCTGGCTAATAGTCGTGATATCTACGGGTCCCTGGAGATATAGCATCTGAGGACTGACTGTCTAGCTTAGACTGCTGTCTGCACAAGAAGAAGACTTACAGATGCTTAAGGAGGATGGAGGCAAATTTTCAACC ++ +FFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26874:1783 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCATGACTGGGAAAACCCTGGCGTGGTCGTGATGTACAACCCGCCTTTATTGGTACCGACAGACCCTGGGGCAGGGCCCAGAGTTTCTGACTTACTTCCAGAATGAAGCTCAACTAGAAAAATCAAGGCTGCTCAGTGATCGGTTCTCT ++ +FFF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFF,FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FFFFFFFFF:FFFFF:FFFFF:FFF:FF,:FFFFFFFF,FFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:7862:1799 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCTCGTTGTCGTGATCGGAACGACCAGCCCCAAGCTAGGAGATCCTGCCATGGGCTTCAGGCTCCTCTGCTGTGTGGCCTTTTGTCTCCTGGGAGCAGGCCCAGTGGATTCTGGAGTCACACAAACCCC ++ +FFF::FFFFFFFFFFFF,FFFFF:F::FFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF:F,FFFFF:FFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFF:FFFFFFF,:FFFFF,FFFFF +@A00261:687:HMMCCDSX5:3:1101:2302:1814 2:N:0:AACCGCAGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTATTGCGTCGTGATCCATATGATTATTTGGGTAACACACTAAAGATGAACTATTCTCCAGGCTTAGTATCTCTGATACTCTTACTGCTTGGAAGAACCCGTGGAGATTCAGTGACCCAGATGGAAGGG ++ +FFFFFFFFFFFFF:FFFFF,FFF:FFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF:FF:FFFFFFFFFFFF,FFFFF:FFFFFFFFFFFFFFFFFF:F,FFFFFF:FFFFF +@A00261:687:HMMCCDSX5:3:1101:6189:1814 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGAACATGTCGGATCCTCATGAGAAGAGGATCCTCAGGCCTTGAGCATCCAGGAGGGTGAAAATGCCACCATGAACTGCAGTTACAAAACTAGTATAAACAATTTACAGTGGTATAGACAAAATTCAGG ++ +FFFFFFFFFFFFFFFFFFFF:FFFFFF:FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:13096:1814 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGCAATGTCGTGATATATTTGATCCTGAGGACAGTGCCTGGAGGTGAGAAGGAAGCCCCCGGCCTGGTCCATACCCCACCACCAACTTGCATAATGGGGGGTGATGTCACCCACCCTCCACTCCCCTC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF +@A00261:687:HMMCCDSX5:3:1101:29243:1814 2:N:0:ACCCGCGGAT+GTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTGCCTGTCGTGATTTCTATGACATGAAGTGTCTACCTTCTTCAGACTCCAATGGATCAGGAACTGGGAATGCAGTGCCAGGCTCGTGGTATCCTGCAGCAGATGTGGGGAGTTTTCCTTCTTTATGT ++ +FFFFFFFFFFFFFFFF,:FFF,:::F:FFFFF:FFF,FFFFFFFFFFFFF::FFFFFF,F:F,,FF:FF,,FFFFFF,FFFFFFFF,FFFFFFF:FFFFFFFFFFFF:,FF,:FFF,FFFFFF,F,FFF,:F,FF:F:,FFFFFFFFF:F +@A00261:687:HMMCCDSX5:3:1101:29993:1830 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGGTTTGTCGTGATACACGCGAAAAATGGAAATGAGACCTTTAGAAAGAGAGATGCTGAACATGAGAAATCAGTAAGCCTCAAGACAAAAACCAGAGGAAATGGAGAGAGATGGGAGGCAGGAAGCCG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFF::FFFFFFF:FFFFFF,FFFFFFF,FFFF:FFF:FFFFFFFF:FF:FF:FFF:FFFFFFFFFFFFFFFFFFFF:F:FFFFFF:FFFFFFFFFFFFFFFFFFFF:FFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:3296:1846 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTGCACAGTCGTGATTCCGCGGAGACACACCAAAAGAGCCTAGAACCTGGGTCCTAGTTTGCACCTAGAAAATGATGCAAGTGGCGAGAGTGATCGTGTTCCAGACCCAGAGTACTTTGAGCCTTGCTA ++ +,,,F,:FF,:,:FFFFF,FF:,F:FFFF:FF:F,F::F,F:,F:FF,FFF,FF:FFFF,F,FF,:FFFFF,FFFFFFFFFFF:F,F::F:F,F,:F,FFFFF:F:FF:F,:F:F,,F,FFF,F,F::FF::,,FFF,FFFFFFF:FF:FF +@A00261:687:HMMCCDSX5:3:1101:3314:1846 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCACTAGTCGTGATTTAATGTAGCACTCTTCTAGCCCAGAGAAGTCTGTTCCAGGACGGGCTCTTTCAGGAGCAGCTAAAGTCAGGGGCCATGTCCACCATGTGATAGAAAGACAAGATGGTCCTGAA ++ +FFFFFF:FFFFFFFFFFFFFFFFFF,FF:FF:FFFF,FFFFFFFFFFFFFFFFFF:FFFFFFFF,FFFF:F,FFFFFFF:FFFFFFFFF:FFFFF:FFF,FF:FFFFFFFFFFFFFFF::FF,F:F,:FFFFFFFFFFFFFFFFFF::F, +@A00261:687:HMMCCDSX5:3:1101:9579:1861 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTAACCGGTCGTGATGGAGAGATGAGATCCTGGCCTGGACCTGAAATGGGCACAAGGTTGTTCTTCTATGTGGCCCTTTGTCTCCTGTGGACAGGACACATGGATGCTGGAATCACCCAGAGCCCAAGA ++ +F:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFF,FFFFFFFFFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:16396:1861 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTTGACGTCGTGATTGAATGGCTTTTCACCAATGCACAGACCCAGAGGACCCCTCCATCCTGCAGTTCCTGCCATGAGCCTCGGGCTCCTGTGCTGTGGGGCCTTTTCTCTCCTGTGGGCAGGTCCAG ++ +:FFFF:FFFFF:FFFFFFFF::FFFFFFFFFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFF::: +@A00261:687:HMMCCDSX5:3:1101:19027:1877 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCACTCGTCGTGATTCTCTCGATTCTTTCTTCAAAGCAGCCATGGGAATCAGGCTCCTCTGTCGTGTGGCCTTTTGGTTCCTGGCGGTAGGCCTCGTAGATGTGAAAGTAACCCAGAGCTCGAGATAT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:F::FFFFF::,FFFFFFFFFFFFFFFF:F,F,FFFFFFFFFFFFFFFF:F,FF,F:F,FF:FF:F:,F:FFF,:,:F,,FF,FFF,FF,:FF,FFF:FF:FF::FFFFF +@A00261:687:HMMCCDSX5:3:1101:23095:1877 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGTTGTGGTGGCTCCCCTTTCATCAATGCACAGATACAGAAGACCCCTCCGTCATGCAGCATCTGCCATGAGCATCGGCCTCCTGTGCTGTGCAGCCTTGTCTCTCCGGTGGGCAGGTCCAGTGAATGC ++ +FFFF:,FF:FFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFFFF:FFFF:F:FF:FF:FFFFFFFFFFFFFFFFFF:F,FFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF:FFFFFFF,FFFFFFFFFFFFFFFFFFF,F +@A00261:687:HMMCCDSX5:3:1101:1099:1892 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGGGGTTGTCGTGATGGCTAAGGAGGGGCTGTGCTGAGAGGCGCTGCTGGGCGTCTGGGCGGAGGACTCCTGGTTCTGGGTGCTGGGAGAGCGATGGGGCTCTCAGCGGTGGGAAGGACCCGAGCTGAG ++ +FFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:2546:1892 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGTGAAAGTCGTGATGTTAATGATTTCTTTGCTCATGCTCACAGAGGGCCTGGTCTAGAATATTCCACATCTGCTCTCACTCTGCCATGGACTCCTGGACCTTCTGCTGTGTGTCCCTTTGCATCCTGG ++ +:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFF:FF:FFFFFFFFFFFFFF,FFF:FFFFFFFFFFF,FFFFFFFF:F, +@A00261:687:HMMCCDSX5:3:1101:9833:1892 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCGTGACGTCGTGATTTTAAGTGGTTCATATGTAAAATGAAGGGTCTGTGGAAGGACATGAATAAAGCACAGGAGGTTGAAGTCAGATTTGCAGCTTTCTAGGCAGGAGACAAGACAATCTGCATCTTC ++ +FFFFFFFFF:F:FFFFFFF:F,:FFFF:F,FFFFFFFFFFFFFFFFFFFF:FF,FFFF:FF,F:FF:FFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFF:FFFFFFFFFFFFFFF:FFFFFF:FFFFFFFF:FFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:10999:1908 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCCGTGGTCGTGATTACTAGGGTCAGAACACATCCAGGCTCCTTAAGGGAAAGCCTCTTTCTGTTTCTGAAACTTTTCAAAGCCAGGGACTTGTCCAATCCAACCTCCTCACAGTTCCTAGCTCCTGA ++ +FFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFF:FF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFF:F:FFFFFFFFFFFFFFF:F:FF:F,FFFFFFFFFF:F:FFF::FFFF,FFFF:FF +@A00261:687:HMMCCDSX5:3:1101:20437:1908 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCTTCGCGTCGTGATTTATGGGCTATTCTTTCTTCAAAGCAGCCATGGGAATCAGGCTCCTCTGTCGTGTGGCCTTTTGTTTCCTGGCTGTAGGCCTCGTAGATGTGAAAGTAACCCAGAGCTCGAGAT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,FFFFF,:FFFFF:FF:FFFFFFFFFF:FFFFFFF:FF,FFFFFFFFFF:,:FFF:FFF,F,FFFFFFFFFFFFFFFFF,FFFFFFFFF:FFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:8314:1924 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGGCATTGGTCGTGATTCTCTGTGTTAGCTCACTGTCCTCCATCCTAGCCAGGGCCCTCTCCCCTCTCTATGCCTTCAATGTGATTTTCACCTTGACCCCTGTCACTGTGTGAACACTGAAGCTTTCTTT ++ +FFFFFFFFFFFFFFFFFFFFFFFF,FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:16848:1924 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCACATGGTCGTGATGAGACCGCAACTATTTCTATTGATCTGGAAGAAGCTTGTACCAGGCAACCCATTGAGGAGAAGTTGGATGAAGAGGGAGAGGGAGATGCTACTCATCACATCAATGTTGGTCGT ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,F:FFFFFF:FFFFFFFFFFFFFFF,FFFFFFF:FFFFFFFFFFFF:F:FFF:FF,FFF:FFFFFF:FFFFFFFFF::FFFFFFFFFFFF:FFFFFFFFFFFFF,,,F +@A00261:687:HMMCCDSX5:3:1101:19090:1924 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGACTCAAGTCGTGATTTATTCTGCTGCTCTGACAGGAGGCAGAGCTCAGGTATGAGTAGAGGAGAGGATGTGGAGCAGAGTCTTTTCCTGAGTGTCCGAGAGGGAGACAGCTCCGTTATAAACTGCACT ++ +FFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFF:FFFFF:FFFFFFFFFFFFFFFFF::FFF:FFFF,FFFFFF:FFF,F:FFFFFFF:FFFFFFFFFFFFFFFF:,FFFFFFFF,FFFFFFFFFFFF:FFFFF,FFFFF,FFFF: +@A00261:687:HMMCCDSX5:3:1101:29903:1924 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCTTGAGTCGTGATATGTGATAAAAACATCCTGAGGACAGTGCCTGGAGGTGAGAAGGAAGCCCCCGGCCTGGTCCATACCCCACCACCAACTTGCATAATGGGGGGTGATGTCACCCACCCACCACT ++ +FF:F:FFFFFFFFFFFFFFFF:FFFFFFFFFF,FF:FFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFF:::FF:FFFFFFFFFFFF:F:FFFFFFFFFFFFFFFFFFFFFFFF::FFFFFF:F:FFFFFF:FFFFF,F:::FF,FF,,, +@A00261:687:HMMCCDSX5:3:1101:20871:1939 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTACACCGTCGTGATAAATATGAAAGAGAGCTCATCTGTGACTGAGGAGCCTTGCTCCATTTCCGGTCTTCTGTGATTTCAATAAGGAAGAAGAATGGAAACTCTCCTGGGAGTGTCTTTGGTGATTCT ++ +:FFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFF,F::FFFFFFFFFFFF,,F,FFFF,FFFFFFFFF,FF:FFF:,FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFF:F:,FF,FFFFF,FFF +@A00261:687:HMMCCDSX5:3:1101:5394:1971 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCCCCTGTCGTGATCAAATTGATCTGGTAAAGCCCTCATCCTGTCCTGACCCTGCCATGGGCACCAGTCTCCTATGCTGGGTGGTCCTGGGTTTCCTAGGGACAGATCACACAGGTGCTGGAGTCTCC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFF,F,:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF::FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF,FFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24542:1971 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCCACATGTCGTGATTTCTTTGACTACTGGGAGACATCCTCTCTAGCCCCAACTGTGCCATGACTATCAGGCTCCTCTGCTACGTGGGCTTTTATTTTCTGGGGGCAGGCCTCATGGAAGCTGACATCT ++ +FFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFF:FFFFFFFFFFFFFFF,FFFFF:F:FFFFFFFF:FF:FF:FFFFFFFFFFFF:FFFFFFFFF,F:FFFF:FFFFF:FFF,FF,FF:FFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26205:1971 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTCCTCGTCGTGATTTATATGACCAATGCCCCTTCAGATTTCAAATCCACTGTGTTGTCCCCCAGCCAAGTGGATTCTCCTCTGCAAATTGGTGGTGGCCTCATGCAAGATCCAGGTTACCGTGTCCA ++ +FFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFF,FF,FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF,F:FFFFFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFF:FFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFFFFFFFFFF: +@A00261:687:HMMCCDSX5:3:1101:10881:1986 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTCCTAGTCGTGATTTTAGTTCCTGGACCCTTCAATGTTCCAGAGGGAGCCACTGTCGCTTTCAACTGTACTTACAGCAAAAGTTCTTCTCAGTCTTTCTTCTGGTAAAGACAGGATTTCAGGAAAGA ++ +FFFFFFFFFFFFFFFFFFFFFF,FFFFFF:FFFFFFFFFFFFF:FFFFFFFFFF:FFFFFFFFFFF:FFFFFFFFFFFFFF,FF,FFF,:F,,,:F:F,:FF,F,:,:,:,:FF,F,::F,FF,:FFFF,F,F,F,F:F:,,:F:,FF,: +@A00261:687:HMMCCDSX5:3:1101:12165:1986 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGAGATCTGTCGTGATTCGGGCACAAGCAGCCTCCCAGCAGGCAGATGATTCTCGTTATTCGCCAAGAAGCTTATAAGCAACAGAATGCAACAGAGAATCGTTTCTCTGTGAACTTCCAGAAAGCAGCCA ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FF:FF:FFFFFFF:FFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:24008:1986 2:N:0:ACCCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTCAGCGTCGTGATTTGTGAAGGCATTTTGGCCATGGCTTTGCAGAGCACTCTGGGGGCGGTGTGGCTAGGGCTTCTCCTCAACTCTCTCTGGAAGGTTGCAGAAAGCAAGGACCAAGTGTTTCAGCC ++ +FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF:F,FFFF:FFF:FF,:,FFFFFFFFFFFFFFFFFFFFFF,,F:FFFFF,FFFFFFFFFFFFFFFFFFF:FFFFFFFF:F:FFF:F:FF,FFFFFFF:F::FFFFFF:FFFF,FFFFF: +@A00261:687:HMMCCDSX5:3:1101:30942:2002 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGCATTCTGTCGTGATTAAAATTGACAGTGACCCTGATCTGGTAAAGCTCCCATCCTGCCCTGACCCTGCCATGGGCACCAGCCTCCTCTGCTGGATGGCCCTGTGTCTCCTGGGGGCAGATCACGCAGA ++ +F,:FF:FFFFFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFF:FFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF +@A00261:687:HMMCCDSX5:3:1101:26151:2033 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTTTTATGTCGTGATTTTAACGCTTTTCTCAGGGGAGAGGCCATCACTTGAAGATGCTGAGTCTTCTGCTCCTTCTCCTGGGACTAGGCTCTGTGTTCAGTGCTGTCATCTCTCAAAAGCCAAGCAGGG ++ +FFFFFFFFFFFFF:FFFFFFFFF,F:,FFFFFFFFFFFFFFFFFFFFFFFFF:FF:F:FFF:FFFFF,,FFFFFFFFFFFFF:FFFFFFFF:FFFFFFFFFFFFFFF,FFFFFFFFFFFFFFFFFFF,FFFFFFF,FFFFFF:FFFFFFF +@A00261:687:HMMCCDSX5:3:1101:6144:2049 2:N:0:AACCGCGGAT+CTAGCGCTGG +GTCGTGACTGGGAAAACCCTGGTTATTTGTCGTGATTTCTTTAACAATTTTTGAAACCCTTCAAAGGCAGAGACTTGTCCAGCCTAAACTGCCTGCTGCTCCTAGCTCCTGAGGCTCAGGGCCCTTGGCTTCTGTCCCCTCTGCTCAGGG ++ +F:FFFF:FFFFFFFFFF:FFF:FFF,FFFFFFFF:FFFFF,FFF:,FFF,F,FFFFF,FF:,FFFFFFFFFFF::,FFF:FF:FFFF,FFFF:F:FFFFFFFFFFFFFFFFFFFFFFF:,FF,FF,:,,FFFFF,,F,FF:F:F,FFF:F +@A00261:687:HMMCCDSX5:3:1101:15727:2080 2:N:0:AACCGCGGAT+CTAGCGCTGT +GTCGTGACTGGGAAAACCCTGGTCTATCGTCGTGATACCCCTTTGGTACAAAAAATACCCTGCTGAAGGTCCTACATTCCTGATATCTATAAGTTCCATTAAGGATAAAAATGAAGATGGAAGATTCACTGTCTTCTTAAACAAAAGTGC ++ +FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFFFFFFFFFF:FFFFFFFFFFFFFFFFFFFF