| Next changeset 1:523dd2a71270 (2017-11-07) |
|
Commit message:
planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/mqc/ commit 27d64fdb96ed232851b1b29cd66250a682df6f18-dirty |
|
added:
mQC.xml readme.rst test-data/mQC_test_data_mESC/BIOBIX_logo.png test-data/mQC_test_data_mESC/annotation_coding.png test-data/mQC_test_data_mESC/annotation_noncoding.png test-data/mQC_test_data_mESC/cumulative.png test-data/mQC_test_data_mESC/density.png test-data/mQC_test_data_mESC/mQC_test_data_mESC.html test-data/mQC_test_data_mESC/phase_relpos_distr.png test-data/mQC_test_data_mESC/rankedgenes.png test-data/mQC_test_data_mESC/rpf_phase.png test-data/mQC_test_data_mESC/tot_phase.png test-data/mQC_test_data_mESC/triplet_id.png test-data/test_data_mESC.sam |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 mQC.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/mQC.xml Fri Sep 29 04:36:16 2017 -0400 |
| b |
| @@ -0,0 +1,156 @@ +<tool id="mqc" name="mQC" version="1.7"> + <description>quality control of ribosome profiling mapping results</description> + <requirements> + <requirement type="package" version="1.7">mqc</requirement> + </requirements> + <command> + mQC.pl + --galaxy Y + --galaxytest $galaxytest + --experiment_name $exp_name + --cores $cores + --species $species + --ens_v $ensv + --unique $uniq + --mapper $mapper + --maxmultimap $maxmultimap + --ens_db "get" + --offset $offsetSelect.offset + --min_length_gd $min_l_gd + --max_length_gd $max_l_gd + --plotrpftool $plotrpftool + --outhtml $MappingQC_output + --outfolder $MappingQC_output.files_path + + #if $inputtypeSelect.inputtype == "sam" + --samfile '$inputtypeSelect.samfile' + --galaxysam Y + #else if $inputtypeSelect.inputtype == "bam" + --samfile '$inputtypeSelect.bamfile' + --galaxysam N + #end if + #if $offsetSelect.offset == "from_file" + --offset_file $offsetSelect.offsetFile + #else if $offsetSelect.offset == "plastid" + #if $inputtypeSelect.inputtype == "sam" + --plastid_bam "convert" + #else if $inputtypeSelect.inputtype == "bam" + --plastid_bam $inputtypeSelect.bamfile + #end if + --min_length_plastid $offsetSelect.min_l_plastid + --max_length_plastid $offsetSelect.max_l_plastid + #end if + + </command> + <inputs> + <param name="galaxytest" type="hidden" value="N"/> + <conditional name="inputtypeSelect"> + <param name="inputtype" type="select" label="Select the format of your input data"> + <option value="sam" selected="true">SAM format</option> + <option value="bam">BAM format</option> + </param> + <when value="sam"> + <param name="samfile" type="data" format="sam" label="Choose the sam file on which you want to do quality control"/> + </when> + <when value="bam"> + <param name="bamfile" type="data" format="bam" label="Choose the bam file on which you want to do quality control"/> + </when> + </conditional> + <param name="exp_name" type="text" label="Enter the experiment name"/> + <param name="cores" type="integer" size="5" value="5" label="Enter the number of cores to run the mapping quality control"/> + <param name="species" type="select" label="Select the species"> + <option value="mouse">Mouse</option> + <option value="human" selected="true">Human</option> + <option value="fruitfly">Fruitfly</option> + </param> + <param name="ensv" type="integer" size="5" value="86" label="Enter the Ensembl database version"/> + <param name="uniq" type="select" label="Perform mQC with unique mapped reads only"> + <option value="Y" selected="true">Yes</option> + <option value="N">No</option> + </param> + <param name="plotrpftool" type="select" label="Choose the style of the RPF phase plot"> + <option value="grouped2D" selected="true">Grouped 2D bar chart</option> + <option value="pyplot3D">3D bar chart (may suffer from Escher effects)</option> + </param> + <param name="mapper" type="select" label="Which mapper is used to generate your SAM file?"> + <option value="STAR" selected="true">STAR</option> + <option value="TopHat2">TopHat2</option> + <option value="HiSat2">HiSat2</option> + </param> + <param name="maxmultimap" value="16" type="integer" label="Only use alignments that mapped fewer times than this value"/> + <param name="min_l_gd" value="26" type="integer" label="Enter the minimum RPF length for gene distribution construction"/> + <param name="max_l_gd" value="34" type="integer" label="Enter the maximum RPF length for gene distribution construction"/> + <conditional name="offsetSelect"> + <param name="offset" type="select" label="The way P-site offsets were determined"> + <option value="standard" selected="true">Standard offsets</option> + <option value="plastid">Calculate offsets with plastid</option> + <option value="from_file">Offsets from custom file</option> + </param> + <when value="standard"/> + <when value="from_file"> + <param name="offsetFile" type="data" format="tabular" label="Choose a text file with the offsets. Each line should be structured as RPF length - offset, seperated by a tab or a space" /> + </when> + <when value="plastid"> + <param name="min_l_plastid" value="22" type="integer" label="Enter the minimum RPF length for plastid offset calculation"/> + <param name="max_l_plastid" value="34" type="integer" label="Enter the maximum RPF length for plastid offset calculation"/> + </when> + </conditional> + </inputs> + <outputs> + <data format="html" name="MappingQC_output" label="HTML output" /> + </outputs> + <tests> + <test> + <param name="galaxytest" value="Y"/> + <param name="inputtype" value="sam"/> + <param name="samfile" value="test_data_mESC.sam" ftype="sam"/> + <param name="exp_name" value="test_data_mESC"/> + <param name="cores" value="10"/> + <param name="species" value="mouse"/> + <param name="ensv" value="86"/> + <param name="uniq" value="Y"/> + <param name="plotrpftool" value="grouped2D"/> + <param name="mapper" value="STAR"/> + <param name="maxmultimap" value="16"/> + <param name="min_l_gd" value="26"/> + <param name="max_l_gd" value="34"/> + <param name="offset" value="standard"/> + <output name="MappingQC_output" file="mQC_test_data_mESC/mQC_test_data_mESC.html" ftype="html"/> + </test> + </tests> + <help> +.. class:: infomark + +**What it does** + +MappingQC is a tool to easily generate some figures which give a nice overview of the quality of the mapping of ribosome profiling data. More specific, it gives an overview of the P site offset calculation, the gene distribution and the metagenic classification. Furthermore, MappingQC does a thorough analysis of the triplet periodicity and the linked triplet phase (typical for ribosome profiling) in the canonical transcript of your data. Especially, the link between the phase distribution and the RPF length, the relative sequence position and the triplet identity are taken into account. + +There is also a version from MappingQC available, which is directly linked to PROTEOFORMER and its SQLite results database structure. This makes its usage in combination with PROTEOFORMER more straight forward and can speed up the total analysis process. + +----- + +.. class:: infomark + +**Input** + +1) SAM file with mapping results + +2) Ensembl DB with annotation for the selected species and Ensembl version + +------ + +.. class:: infomark + +**Output** + +1) HTML file: overview file of the tables and figures. Main output of MappingQC + +2) ZIP file: contains a folder with all necessary figures. In this folder you can also find the overview HTML file. + + </help> + <citations> + <citation type="bibtex"> + @article{verbruggen_menschaert_2017, title={MappingQC: a tool for quality control of ribosome profiling mapping results}, journal={Unpublished}, author={Verbruggen, Steven and Menschaert, Gerben}, year={2017}} + </citation> + </citations> +</tool> |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 readme.rst --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/readme.rst Fri Sep 29 04:36:16 2017 -0400 |
| b |
| @@ -0,0 +1,18 @@ +mQC galaxy wrapper (stand-alone version) +======================================= + +This repository contains all the files needed to embed the stand-alone version of MappingQC in a Galaxy environment. +More information on mQC and the core program files can be found at https://github.com/Biobix/mQC. +The authors of mQC can be contacted at Steven.Verbruggen@UGent.be or Gerben.Menschaert@UGent.be + +========= +Copyright +========= + +Copyright (C) 2017 Verbruggen Steven & Menschaert Gerben + +This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. + +This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. + +You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/. |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/BIOBIX_logo.png |
| b |
| Binary file test-data/mQC_test_data_mESC/BIOBIX_logo.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/annotation_coding.png |
| b |
| Binary file test-data/mQC_test_data_mESC/annotation_coding.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/annotation_noncoding.png |
| b |
| Binary file test-data/mQC_test_data_mESC/annotation_noncoding.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/cumulative.png |
| b |
| Binary file test-data/mQC_test_data_mESC/cumulative.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/density.png |
| b |
| Binary file test-data/mQC_test_data_mESC/density.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/mQC_test_data_mESC.html --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/mQC_test_data_mESC/mQC_test_data_mESC.html Fri Sep 29 04:36:16 2017 -0400 |
| b |
| b'@@ -0,0 +1,350 @@\n+<!DOCTYPE html>\n+<html>\n+<head>\n+ <title>Mapping QC Report test_data_mESC</title>\n+ <meta charset="utf-8"></meta>\n+ <meta name="description" content="Overview HTML of all mappingQC results"></meta>\n+ <link href="https://fonts.googleapis.com/css?family=Indie+Flower" rel="stylesheet">\n+ <style>\n+ *{\n+ box-sizing: border-box;\n+ font-family: \'Helvetica Neue\', Helvetica, Arial, sans-serif;\n+ }\n+\n+\n+ nav{\n+ float:left;\n+ padding: 15px;\n+ width: 17%;\n+ position: fixed;\n+ height: 100%;\n+ overflow: auto;\n+ border-right: ridge;\n+ border-color: lightgrey;\n+ margin-top: 60px;\n+ margin-left: -20px;\n+ padding-left: 20px;\n+ background-color: white;\n+ z-index:1;\n+ }\n+\n+ nav ul {\n+ list-style-type: none;\n+ margin: 0px;\n+ padding: 5px;\n+ padding-top: 15px;\n+\n+ }\n+\n+ nav li{\n+ padding: 8px;\n+ margin-bottom: 8px;\n+ background-color: #33b5e5;\n+ color: #ffffff;\n+ }\n+\n+ nav li:hover {\n+ background-color: #0099cc;\n+ }\n+\n+ #content{\n+ position: absolute;\n+ margin-left:19%;\n+ height: 76%;\n+ }\n+\n+ #rpf_phase{\n+ width:100%;\n+ }\n+\n+ #header{\n+ background-color: grey;\n+ color: white;\n+ position:fixed;\n+ height: 2.7cm;\n+ width:110%;\n+ padding: 15px;\n+ padding-top: 10px;\n+ margin-left: -10px;\n+ margin-top: -30px;\n+ margin-right: -10px;\n+ margin-bottom: 10px;\n+ overflow: visible;\n+ z-index: 2;\n+ }\n+\n+ #mappingqc{\n+ font-family: \'Indie Flower\', cursive;\n+ font-size: 44px;\n+ padding-left: 10px;\n+ position: relative;\n+ z-index: 4;\n+ }\n+ #run_name{\n+ padding-left: 43px;\n+ position: relative;\n+ z-index: 4;\n+ }\n+\n+ #biobix_logo{\n+ height:60%;\n+ position: absolute;\n+ right: 200px;\n+ top: 30px;\n+ }\n+\n+ a {\n+ color: inherit;\n+ text-decoration: none;\n+ }\n+\n+ .anchor{\n+ display: block;\n+ height: 14%; /*same height as header*/\n+ margin-top: -10px; /*same height as header*/\n+ visibility: hidden;\n+ }\n+\n+ #analysis_info_table{\n+ border-style: none;\n+ border-width: 0px;\n+ }\n+\n+ th {\n+ border-style: solid;\n+ border-width: 0px;\n+ border-color: white;\n+ border-collapse: collapse;\n+ padding: 5px;\n+ background-color: #33b5e5;\n+ color: #ffffff;\n+ }\n+\n+ td {\n+ border-style: solid;\n+ border-width: 0px;\n+ border-color: white;\n+ border-collapse: collapse;\n+ background-color: #f2f2f2;\n+ padding: 5px;\n+ }\n+\n+ img {\n+ max-width: 98%;\n+ height: auto;\n+ width: auto\\9; /* ie8 */\n+ }\n+\n+ #ranked_genes, #cumulative, #genes_density, #annotation_coding, #annotation_noncoding {\n+ width: 20cm;\n+ }\n+\n+ #offset_table {\n+ float: left;\n+ display: block;\n+ margin-right: 120px;\n+ }\n+\n+ #plastid_img {\n+ float: left;\n+ display: block;\n+ max-width: 600px;\n+ }\n+\n+ #section3 {\n+ clear: left;\n+ }\n+\n+\n+ #footer{\n+ background-color: grey;\n+ color: white;\n+ position: fixed;\n+ bottom: 0cm;\n+ padding-left: 30px;\n+ '..b'n</td>\n+ <td>86</td>\n+ </tr>\n+ \n+ <tr>\n+ <td>Selected offset source</td>\n+ <td>standard</td>\n+ </tr>\n+ <tr>\n+ <td>Used only unique alignments</td>\n+ <td>Y</td>\n+ </tr>\n+ <tr>\n+ <td>Total mapped genomic sequences</td>\n+ <td>10 425</td>\n+ </tr>\n+ \n+ </table>\n+ </p>\n+\n+ <span class="anchor" id="section2"></span>\n+ <h2 id="plastid">Offsets overview</h2>\n+ <p>\n+ <table id="offset_table">\n+ <tr>\n+ <th id="table_header">RPF length</th>\n+ <th id="table_header">Offset</th>\n+ </tr>\n+ <tr>\n+ <td>26</td>\n+ <td>12</td>\n+ </tr>\n+ <tr>\n+ <td>27</td>\n+ <td>12</td>\n+ </tr>\n+ <tr>\n+ <td>28</td>\n+ <td>12</td>\n+ </tr>\n+ <tr>\n+ <td>29</td>\n+ <td>12</td>\n+ </tr>\n+ <tr>\n+ <td>30</td>\n+ <td>12</td>\n+ </tr>\n+ <tr>\n+ <td>31</td>\n+ <td>13</td>\n+ </tr>\n+ <tr>\n+ <td>32</td>\n+ <td>13</td>\n+ </tr>\n+ <tr>\n+ <td>33</td>\n+ <td>13</td>\n+ </tr>\n+ <tr>\n+ <td>34</td>\n+ <td>14</td>\n+ </tr>\n+ \n+ </table>\n+ </p>\n+ \n+\n+ <span class="anchor" id="section3"></span>\n+ <h2 id="gene_distributions">Gene distributions</h2>\n+ <p>\n+ <div class="img">\n+ <img src="rankedgenes.png" alt="Ranked genes" id="ranked_genes">\n+ </div>\n+ </p>\n+ <p>\n+ <div class="img">\n+ <img src="cumulative.png" alt="Cumulative genes" id="cumulative">\n+ </div>\n+ </p>\n+ <p>\n+ <div class="img">\n+ <img src="density.png" alt="Genes density" id="genes_density">\n+ </div>\n+ </p>\n+\n+ <span class="anchor" id="section4"></span>\n+ <h2 id="metagenic_classification">Metagenic classification</h2>\n+ <p>\n+ <div class="img">\n+ <img src="annotation_coding.png" alt="Metagenic classification coding" id="annotation_coding">\n+ </div>\n+ </p>\n+ <p>\n+ <div class="img">\n+ <img src="annotation_noncoding.png" alt="Noncoding classification" id="annotation_noncoding">\n+ </div>\n+ </p>\n+\n+ <span class="anchor" id="section5"></span>\n+ <h2 id="tot_phase">Total phase distribution</h2>\n+ <p>\n+ <div class="img">\n+ <img src="tot_phase.png" alt="total phase plot" id="tot_phase_img">\n+ </div>\n+ </p>\n+\n+ <span class="anchor" id="section6"></span>\n+ <h2 id="phase_rpf_distr">RPF phase distribution</h2>\n+ <p>\n+ <div class="img">\n+ <img src="rpf_phase.png" alt="rpf phase plot" id="rpf_phase_img">\n+ </div>\n+ </p>\n+\n+ <span class="anchor" id="section7"></span>\n+ <h2 id="phase_relpos_distr">Phase - relative position distribution</h2>\n+ <p>\n+ <div class="img">\n+ <img src="phase_relpos_distr.png" alt="phase relpos distr" id="phase_relpos_distr_img">\n+ </div>\n+ </p>\n+\n+ <span class="anchor" id="section8"></span>\n+ <h2 id="triplet_identity">Triplet identity plots</h2>\n+ <p>\n+ <div class="img">\n+ <img src="triplet_id.png" alt="triplet identity plots" id="triplet_id_img">\n+ </div>\n+ </p>\n+ <br><br>\n+ </div>\n+\n+ <div id="footer">\n+ <p id="footer_content">Generated with mQC - BioBix lab Ghent (Belgium) - Steven Verbruggen</p>\n+ </div>\n+\n+</body>\n+</html>\n\\ No newline at end of file\n' |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/phase_relpos_distr.png |
| b |
| Binary file test-data/mQC_test_data_mESC/phase_relpos_distr.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/rankedgenes.png |
| b |
| Binary file test-data/mQC_test_data_mESC/rankedgenes.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/rpf_phase.png |
| b |
| Binary file test-data/mQC_test_data_mESC/rpf_phase.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/tot_phase.png |
| b |
| Binary file test-data/mQC_test_data_mESC/tot_phase.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/triplet_id.png |
| b |
| Binary file test-data/mQC_test_data_mESC/triplet_id.png has changed |
| b |
| diff -r 000000000000 -r 8cd7fafbbdc4 test-data/test_data_mESC.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/test_data_mESC.sam Fri Sep 29 04:36:16 2017 -0400 |
| b |
| b'@@ -0,0 +1,10429 @@\n+@HD\tVN:1.4\tSO:coordinate\n+@SQ\tSN:Y\tLN:91744698\n+@PG\tID:STAR\tPN:STAR\tVN:STAR_2.5.2a\tCL:STAR --runThreadN 20 --genomeDir /data/igenomes//Mus_musculus/Ensembl/GRCm38/Sequence/STARIndex/mmu.GRCm38.86.genome.35bpOverhang --genomeLoad NoSharedMemory --readFilesIn /data2/steven/mappingQC/proteoformer/mESC/non-uniq_plastid/fastq/fastq1_norrna_nosnrna_notrna.fq --outFileNamePrefix /data2/steven/mappingQC/proteoformer/mESC/non-uniq_plastid/STAR/fastq1/ --outMultimapperOrder Random --outSAMtype BAM SortedByCoordinate --outSAMattributes Standard --outFilterMultimapNmax 16 --outFilterMismatchNmax 2 --seedSearchStartLmaxOverLread 0.5\n+@CO\tuser command line: STAR --outSAMattributes Standard --outSAMtype BAM SortedByCoordinate --genomeLoad NoSharedMemory --seedSearchStartLmaxOverLread .5 --outFilterMultimapNmax 16 --outMultimapperOrder Random --outFilterMismatchNmax 2 --genomeDir /data/igenomes//Mus_musculus/Ensembl/GRCm38/Sequence/STARIndex/mmu.GRCm38.86.genome.35bpOverhang --runThreadN 20 --outFileNamePrefix /data2/steven/mappingQC/proteoformer/mESC/non-uniq_plastid/STAR/fastq1/ --readFilesIn /data2/steven/mappingQC/proteoformer/mESC/non-uniq_plastid/fastq/fastq1_norrna_nosnrna_notrna.fq\n+NS500361:73:HJNVLBGXX:4:23409:11240:5790\t16\tY\t722631\t255\t17M301816N15M\t*\t0\t0\tGTTGTTTTCTCTTCCTGTGGGTACTGAGATGT\tEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:1:22310:23270:7707\t16\tY\t722631\t255\t17M301816N15M\t*\t0\t0\tGTTGTTTTCTCTTCCTGTGGGTACTGAGATGT\tEEEEEEEEEEEEEEAEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:1:23209:23642:8982\t16\tY\t722631\t255\t17M301816N15M\t*\t0\t0\tGTTGTTTTCTCTTCCTGCGGGTACTGAGATGT\tEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:23\tnM:i:2\n+NS500361:73:HJNVLBGXX:1:12308:24675:14201\t16\tY\t722631\t255\t17M301816N15M\t*\t0\t0\tGTTGTTTTCTCTTCCTATGGGTACTGAGATGT\tEEEEEEEEEEEEEEEAEEEAEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:23\tnM:i:2\n+NS500361:73:HJNVLBGXX:1:12206:3848:9150\t16\tY\t722631\t255\t17M301816N15M\t*\t0\t0\tGTTGTTTTCTCTTCCTGTGGGTACTGAGATGT\tEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:1:22204:25492:6621\t16\tY\t722631\t255\t17M301816N14M\t*\t0\t0\tGTTGTTTTCTCTTCCTGTGGGTACTGAGATG\tEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:24\tnM:i:1\n+NS500361:73:HJNVLBGXX:3:13504:20167:11811\t16\tY\t722631\t255\t17M301816N15M\t*\t0\t0\tGTTGTTTTCTCTTCCTGTGGGTACTGAGATGT\tEEEEEEEEEA6EE/6EEE/EEEE/EEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:1:21103:15762:1449\t16\tY\t722631\t255\t17M301816N15M\t*\t0\t0\tGTTGTTTTCTCTTCCTGTGGGTACTGAGATGT\tEEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:2:11107:22740:3995\t16\tY\t722631\t255\t17M301816N14M\t*\t0\t0\tGTTGTTTTCTCTTCCTGTGGGTACTGAGATG\tEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:24\tnM:i:1\n+NS500361:73:HJNVLBGXX:1:13202:23464:15291\t16\tY\t724726\t255\t26M\t*\t0\t0\tATGGACATCTACAAATCTGCAAAACT\tEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:0\n+NS500361:73:HJNVLBGXX:3:12607:8943:20033\t16\tY\t725288\t255\t29M\t*\t0\t0\tAATAGGCGCAATATCATCTTTATATTCAC\tAEEE//E/EAAEE6EAAEEEEEE//A6AA\tNH:i:1\tHI:i:1\tAS:i:28\tnM:i:0\n+NS500361:73:HJNVLBGXX:1:12207:19105:6726\t16\tY\t725422\t255\t24M\t*\t0\t0\tAGGCCAACTTTATGATGTCGCATT\tEEEEEEEEEEAEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:23\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:13207:22218:4691\t16\tY\t725422\t255\t24M\t*\t0\t0\tAGGCCAACTTTATGATGTCGCATT\tEEEEEEEEEEAEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:23\tnM:i:0\n+NS500361:73:HJNVLBGXX:1:21104:16658:11119\t16\tY\t725483\t255\t28M\t*\t0\t0\tTCTTGCAGAAGTCACAGCTGTGAGGATA\tEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:27\tnM:i:0\n+NS500361:73:HJNVLBGXX:3:21509:8970:13183\t16\tY\t725483\t255\t28M\t*\t0\t0\tTCTTGCAGAAGTCACAGCTGTGAGGATA\tEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:27\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:11305:14257:10940\t16\tY\t725491\t255\t28M\t*\t0\t0\tAAGTCACAGCTGTGAGGATAGTCTTTCG\tEEEEEEEEEEEEEEEEEEEEEEEAAAAA\tNH:i:1\tHI:i:1\tAS:i:27\tnM:i:0\n+NS500361:73:HJNVLBGXX:3:21604:14050:10014\t16\tY\t725495\t255\t28M\t*\t0\t0\tCACAGCTGTGAGGATAGTCTTTCGTATG\tEAAEEEEEEAA//EEEE6EEEEEAAAAA\tN'..b'61:73:HJNVLBGXX:2:13207:4683:10961\t0\tY\t90810872\t255\t27M\t*\t0\t0\tAGAGGGTTGCTCTGGTGCTGCACCATT\tAAAAAEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:26\tnM:i:0\n+NS500361:73:HJNVLBGXX:3:22508:25915:12474\t0\tY\t90810934\t255\t26M\t*\t0\t0\tACCAGGACCCAGGGAACTCGCACAGA\t6AAAAEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:21210:24559:18709\t0\tY\t90810938\t255\t23M\t*\t0\t0\tGGACCCAGGGAACTCGCACAGAT\tAA6AAEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:22\tnM:i:0\n+NS500361:73:HJNVLBGXX:1:21201:15802:1362\t0\tY\t90811008\t255\t31M\t*\t0\t0\tTGCGACGCTCACATTTGCACACTGAGGCACT\tAAAAAEEEEEEEEEEEEEEEEEEEEEEEAEE\tNH:i:1\tHI:i:1\tAS:i:30\tnM:i:0\n+NS500361:73:HJNVLBGXX:4:22508:5702:2589\t0\tY\t90811274\t255\t27M\t*\t0\t0\tGACAGATACAGAGTCAGACAGACACAG\tAAAAAEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:26\tnM:i:0\n+NS500361:73:HJNVLBGXX:4:12507:16135:4243\t0\tY\t90811274\t255\t27M\t*\t0\t0\tGACAGATACAGAGTCAGACAGACACAG\tAAAAAEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:26\tnM:i:0\n+NS500361:73:HJNVLBGXX:1:11209:1852:12301\t0\tY\t90811401\t255\t29M\t*\t0\t0\tAGATAGACAGGGAGTCAGACAGACAGACA\tAAAAAEEEEEEEEEEEEEEAEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:28\tnM:i:0\n+NS500361:73:HJNVLBGXX:4:21510:25224:16064\t0\tY\t90811401\t255\t29M\t*\t0\t0\tAGATAGAAAGGGAGTCAGACAGACAGACA\tAAAAAEA/EEAEE6AAEEEEEEE<EE6EE\tNH:i:1\tHI:i:1\tAS:i:26\tnM:i:1\n+NS500361:73:HJNVLBGXX:4:22402:18176:16608\t0\tY\t90811401\t255\t29M\t*\t0\t0\tAGATAGACAGGGAGTCAGACAGACAGACA\tAAAAAEEEEEEEEEEEEEAEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:28\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:12208:22759:14289\t0\tY\t90811401\t255\t29M\t*\t0\t0\tAGATAGACAGGGAGTCAGACAGACAGACA\tAAAAAEEEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:28\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:22307:4582:6010\t0\tY\t90811405\t255\t26M\t*\t0\t0\tAGACAGGGAGTCAGACAGACAGACAC\tAAAAAEAEEAEEEAEEE6EEE/EEEE\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:0\n+NS500361:73:HJNVLBGXX:4:11507:18075:6161\t0\tY\t90811405\t255\t26M\t*\t0\t0\tAGACAGGGAGTCAGACAGACAGACAC\tAAAAAEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:11206:3104:7237\t0\tY\t90812454\t255\t28M\t*\t0\t0\tGACAGAGAGTCAGAAGGACAGATAGCGA\tAAAAAEEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:27\tnM:i:0\n+NS500361:73:HJNVLBGXX:3:12401:10070:17212\t0\tY\t90812454\t255\t28M\t*\t0\t0\tGACAGAGAGTCAGAAGGACAGATAGCGA\tAAAAAEEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:27\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:13308:10193:1797\t0\tY\t90812454\t255\t28M\t*\t0\t0\tGACAGAGAGTCAGAAGGACAGATAGCGA\tAAAAAEEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:27\tnM:i:0\n+NS500361:73:HJNVLBGXX:4:12402:21969:12394\t0\tY\t90812474\t255\t31M\t*\t0\t0\tGATAGCGAGACAGAGACAGACAGAGTCAGAG\tAAAAA/EAEAEEEEE/EEE6EEEE/EEEEEE\tNH:i:1\tHI:i:1\tAS:i:30\tnM:i:0\n+NS500361:73:HJNVLBGXX:2:22308:10582:19050\t0\tY\t90812721\t255\t27M\t*\t0\t0\tACAGAGATAAAGAGACAGACATACAAA\tAAAAAEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:26\tnM:i:0\n+NS500361:73:HJNVLBGXX:3:12508:21210:13067\t0\tY\t90812728\t255\t27M\t*\t0\t0\tTAAAGAGACAGACATACAAACAGGGAG\tAAAAAEEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:26\tnM:i:0\n+NS500361:73:HJNVLBGXX:1:12205:9276:5822\t0\tY\t90813184\t255\t28M\t*\t0\t0\tGTTAGGGTTAGATTAGGGTTAGTGGGGT\tAA/AAAEA/AEEEEEAAA/EEEE6E6EE\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:1:12207:12022:13195\t0\tY\t90813184\t255\t28M\t*\t0\t0\tGTTAGGGTTAGATTAGGGTTAGTGGGGT\tAAA6AEAEEEEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:1:12301:3949:19074\t0\tY\t90813184\t255\t28M\t*\t0\t0\tGTTAGGGTTAGATTAGGGTTAGTGGGGT\tAAAAAEEEEEEEEEEEEEAEEAEEEEEE\tNH:i:1\tHI:i:1\tAS:i:25\tnM:i:1\n+NS500361:73:HJNVLBGXX:3:13508:2581:9033\t0\tY\t90813191\t255\t21M\t*\t0\t0\tTTAGATTAGGGTTAGGGGGGT\tAAAAAEAAEEEAEEEEEEEAE\tNH:i:1\tHI:i:1\tAS:i:20\tnM:i:0\n+NS500361:73:HJNVLBGXX:4:11504:4680:2735\t0\tY\t90828813\t255\t2S20M2S\t*\t0\t0\tGGGCCGCCGACACGGCCGGACCCG\tAAAAAEEEEEEEEEEEEAEEEEEE\tNH:i:1\tHI:i:1\tAS:i:17\tnM:i:1\n+NS500361:73:HJNVLBGXX:2:11203:3765:19149\t0\tY\t90828813\t255\t2S16M2S\t*\t0\t0\tGGGCCGCCGACACGGCCGGA\tAAAAAEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:15\tnM:i:0\n+NS500361:73:HJNVLBGXX:1:11109:3113:2386\t16\tY\t90828813\t255\t20M1S\t*\t0\t0\tGCCGCCGACACGGCCGGACCC\tEEAEEEA/EAEEEEEEAAA/A\tNH:i:1\tHI:i:1\tAS:i:17\tnM:i:1\n+NS500361:73:HJNVLBGXX:4:11610:8616:18356\t0\tY\t90828813\t255\t20M7S\t*\t0\t0\tGCCGCCGACACGGCCGGACCCGCCGCC\tAAAAAEE6EEEEEEEEEEEEEEEEEEE\tNH:i:1\tHI:i:1\tAS:i:17\tnM:i:1\n' |
