# HG changeset patch # User rnateam # Date 1506674176 14400 # Node ID 8cd7fafbbdc4b8be9319e71d295f249d0566ea8b planemo upload for repository https://github.com/bgruening/galaxytools/tree/master/tools/rna_tools/mqc/ commit 27d64fdb96ed232851b1b29cd66250a682df6f18-dirty 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 @@ -0,0 +1,156 @@ + + quality control of ribosome profiling mapping results + + mqc + + + 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 + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +.. 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. + + + + + @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}} + + + 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 @@ -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/. diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/BIOBIX_logo.png Binary file test-data/mQC_test_data_mESC/BIOBIX_logo.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/annotation_coding.png Binary file test-data/mQC_test_data_mESC/annotation_coding.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/annotation_noncoding.png Binary file test-data/mQC_test_data_mESC/annotation_noncoding.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/cumulative.png Binary file test-data/mQC_test_data_mESC/cumulative.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/density.png Binary file test-data/mQC_test_data_mESC/density.png has changed 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 @@ -0,0 +1,350 @@ + + + + Mapping QC Report test_data_mESC + + + + + + + + + + + +
+ +

Analysis information

+

+ + + + + + + + + + + + + + + + + + + + + + + + + + + + +
FeatureValue
Speciesmouse
Ensembl version86
Selected offset sourcestandard
Used only unique alignmentsY
Total mapped genomic sequences10 425
+

+ + +

Offsets overview

+

+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +
RPF lengthOffset
2612
2712
2812
2912
3012
3113
3213
3313
3414
+

+ + + +

Gene distributions

+

+

+ Ranked genes +
+

+

+

+ Cumulative genes +
+

+

+

+ Genes density +
+

+ + +

Metagenic classification

+

+

+ Metagenic classification coding +
+

+

+

+ Noncoding classification +
+

+ + +

Total phase distribution

+

+

+ total phase plot +
+

+ + +

RPF phase distribution

+

+

+ rpf phase plot +
+

+ + +

Phase - relative position distribution

+

+

+ phase relpos distr +
+

+ + +

Triplet identity plots

+

+

+ triplet identity plots +
+

+

+
+ + + + + \ No newline at end of file diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/phase_relpos_distr.png Binary file test-data/mQC_test_data_mESC/phase_relpos_distr.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/rankedgenes.png Binary file test-data/mQC_test_data_mESC/rankedgenes.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/rpf_phase.png Binary file test-data/mQC_test_data_mESC/rpf_phase.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/tot_phase.png Binary file test-data/mQC_test_data_mESC/tot_phase.png has changed diff -r 000000000000 -r 8cd7fafbbdc4 test-data/mQC_test_data_mESC/triplet_id.png Binary file test-data/mQC_test_data_mESC/triplet_id.png has changed 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 @@ -0,0 +1,10429 @@ +@HD VN:1.4 SO:coordinate +@SQ SN:Y LN:91744698 +@PG ID:STAR PN:STAR VN:STAR_2.5.2a CL: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 +@CO user 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 +NS500361:73:HJNVLBGXX:4:23409:11240:5790 16 Y 722631 255 17M301816N15M * 0 0 GTTGTTTTCTCTTCCTGTGGGTACTGAGATGT EEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:1 +NS500361:73:HJNVLBGXX:1:22310:23270:7707 16 Y 722631 255 17M301816N15M * 0 0 GTTGTTTTCTCTTCCTGTGGGTACTGAGATGT EEEEEEEEEEEEEEAEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:1 +NS500361:73:HJNVLBGXX:1:23209:23642:8982 16 Y 722631 255 17M301816N15M * 0 0 GTTGTTTTCTCTTCCTGCGGGTACTGAGATGT EEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:23 nM:i:2 +NS500361:73:HJNVLBGXX:1:12308:24675:14201 16 Y 722631 255 17M301816N15M * 0 0 GTTGTTTTCTCTTCCTATGGGTACTGAGATGT EEEEEEEEEEEEEEEAEEEAEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:23 nM:i:2 +NS500361:73:HJNVLBGXX:1:12206:3848:9150 16 Y 722631 255 17M301816N15M * 0 0 GTTGTTTTCTCTTCCTGTGGGTACTGAGATGT EEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:1 +NS500361:73:HJNVLBGXX:1:22204:25492:6621 16 Y 722631 255 17M301816N14M * 0 0 GTTGTTTTCTCTTCCTGTGGGTACTGAGATG EEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:24 nM:i:1 +NS500361:73:HJNVLBGXX:3:13504:20167:11811 16 Y 722631 255 17M301816N15M * 0 0 GTTGTTTTCTCTTCCTGTGGGTACTGAGATGT EEEEEEEEEA6EE/6EEE/EEEE/EEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:1 +NS500361:73:HJNVLBGXX:1:21103:15762:1449 16 Y 722631 255 17M301816N15M * 0 0 GTTGTTTTCTCTTCCTGTGGGTACTGAGATGT EEEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:1 +NS500361:73:HJNVLBGXX:2:11107:22740:3995 16 Y 722631 255 17M301816N14M * 0 0 GTTGTTTTCTCTTCCTGTGGGTACTGAGATG EEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:24 nM:i:1 +NS500361:73:HJNVLBGXX:1:13202:23464:15291 16 Y 724726 255 26M * 0 0 ATGGACATCTACAAATCTGCAAAACT EEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:0 +NS500361:73:HJNVLBGXX:3:12607:8943:20033 16 Y 725288 255 29M * 0 0 AATAGGCGCAATATCATCTTTATATTCAC AEEE//E/EAAEE6EAAEEEEEE//A6AA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:1:12207:19105:6726 16 Y 725422 255 24M * 0 0 AGGCCAACTTTATGATGTCGCATT EEEEEEEEEEAEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:23 nM:i:0 +NS500361:73:HJNVLBGXX:2:13207:22218:4691 16 Y 725422 255 24M * 0 0 AGGCCAACTTTATGATGTCGCATT EEEEEEEEEEAEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:23 nM:i:0 +NS500361:73:HJNVLBGXX:1:21104:16658:11119 16 Y 725483 255 28M * 0 0 TCTTGCAGAAGTCACAGCTGTGAGGATA EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:21509:8970:13183 16 Y 725483 255 28M * 0 0 TCTTGCAGAAGTCACAGCTGTGAGGATA EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:11305:14257:10940 16 Y 725491 255 28M * 0 0 AAGTCACAGCTGTGAGGATAGTCTTTCG EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:21604:14050:10014 16 Y 725495 255 28M * 0 0 CACAGCTGTGAGGATAGTCTTTCGTATG EAAEEEEEEAA//EEEE6EEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:12608:10057:4412 16 Y 725495 255 27M * 0 0 CACAGCTGTGAGGATAGTCTTTCGTAT EEEEEEEEEEEEEEEEEEEEAEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:4:13603:14845:20304 16 Y 725495 255 28M * 0 0 CACAGCTGTGAGGATAGTCTTTCGTATG EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:12312:9112:13399 16 Y 725499 255 27M * 0 0 GCTGTGAGGATAGTCTTTCGTATGAAT EEEEEEEEEEEEE/EEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:3:12408:22070:13397 16 Y 725512 255 29M * 0 0 TCTTTCGTATGAATGGAGATAACGTGGCG EEEEEEEEEEAEEEEEAEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:11501:15814:18977 16 Y 725513 255 28M * 0 0 CTTTCGTATGAATGGAGATAACGTGGCG EEEEEEEEEEEEEEEEEEAEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:4:21609:13899:4245 16 Y 725533 255 29M * 0 0 ACGTGGCGCTTAAAACCTGAGGCATCTTT EEEEEEEEEEEEEEEEEEEEAEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:4:22506:12772:9978 16 Y 725533 255 29M * 0 0 ACGTGGCGCTTAAAACCTGAGGCATCTTT E/EEEEEEEEEEEEEEEE/EEEE/AAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:1:23110:9540:1760 16 Y 725731 255 28M * 0 0 GGACTATTAAAGTCACAGTGTCTACATT EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:13110:24181:17931 16 Y 725731 255 28M * 0 0 GGACTATTAAAGTCACAGTGTCTACATT EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:11309:18022:10277 16 Y 725959 255 29M * 0 0 GCATGTTGCTGAGCCTCTTTGTTATCTGA E/EEEEEEEEEEA6EEEEEEEAE/AAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:1:22101:10701:14764 16 Y 725959 255 28M * 0 0 GCATGTTGCTGAGCCTCTTTGTTATCTG EEEEAEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:23309:25574:11672 16 Y 725968 255 30M * 0 0 TGAGCCTCTTTGTTATCTGAGAAAGTCAGG EEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:29 nM:i:0 +NS500361:73:HJNVLBGXX:3:12410:24897:2313 16 Y 725968 255 30M * 0 0 TGAGCCTCTTTGTTATCTGAGAAAGTCAGG EEEEE6AEEEEEEEEEEAEEAA/EEAA/AA NH:i:1 HI:i:1 AS:i:29 nM:i:0 +NS500361:73:HJNVLBGXX:2:11212:6567:4627 16 Y 725972 255 27M * 0 0 CCTCTTTGTTATCTGAGAAAGTCAGGA EEAEEEEEEEEEEEEEEEAEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:4:13502:20069:14505 16 Y 725977 255 29M * 0 0 TTGTTATCTGAGAAAGTCAGGAGACAGAT EEEEEEEEEEEEEEEEEEEEE/AEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:13403:5195:13996 16 Y 725977 255 29M * 0 0 TTGTTATCTGAGAAAGTCAGGAGACAGAT EEEEEEEEEEE/AEEAEEAEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:12401:5836:6006 16 Y 725977 255 29M * 0 0 TTGTTATCTGAGAAAGTCAGGAGACAGAT EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:2:13211:25545:8451 16 Y 725979 255 27M * 0 0 GTTATCTGAGAAAGTCAGGAGACAGAT EEEEEEEEEEEAEEEEEEEEE/AAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:2:21105:8227:15970 16 Y 726292 255 30M * 0 0 TTTTCTGTATGCATTGTTTTGTGAGTACAC EEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:29 nM:i:0 +NS500361:73:HJNVLBGXX:2:21206:10055:3128 16 Y 726292 255 30M * 0 0 TTTTCTGTATGCATTGTTTTGTGAGTACAC EEAEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:29 nM:i:0 +NS500361:73:HJNVLBGXX:1:12106:9954:15669 16 Y 726292 255 29M * 0 0 TTTTCTGTATGCATTGTTTTGTGAGTACA EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:1:22101:24586:13939 16 Y 726293 255 27M * 0 0 TTTCTGTATGCATTGTTTTGTGAGTAC EEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:2:23101:16296:1903 16 Y 726311 255 28M * 0 0 TGTGAGTACACAAAGTCCCAGCATGAGA EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:13306:11451:6227 16 Y 726314 255 26M * 0 0 GAGTACACAAAGTCCCAGCATGAGAA EEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:0 +NS500361:73:HJNVLBGXX:2:23104:16803:5540 16 Y 726478 255 29M * 0 0 AGGTATTCAGGATGGTTTTTTGTATGTCT EAEEEEEAEAAEEAAEEEEEEEEEA/AAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:4:22604:22115:20105 16 Y 726478 255 29M * 0 0 AGGTATTCAGGATGGTTTTTTGTATGTCT EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:4:23403:21768:3990 16 Y 726478 255 29M * 0 0 AGGTATTCAGGATGGTTTTTTGTATGTCT E/EEEEEEAEEEEEE6EEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:1:11101:6999:6403 16 Y 726481 255 29M * 0 0 TATTCAGGATGGTTTTTTGTATGTCTTTT AEAEEEEEEEEEAEEEEEE//AEAAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:13604:22363:10245 16 Y 726481 255 28M * 0 0 TATTCAGGATGGTTTTTTGTATGTCTTT EEEEEEEEEEEEEEEEEEEEEEEAAAA/ NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:4:12407:26241:4867 16 Y 726482 255 28M * 0 0 ATTCAGGATGGTTTTTTGTATGTCTTTT EEEEEEEEEEEEEEEEEEEEAEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:11208:12298:14791 16 Y 726482 255 26M * 0 0 ATTCAGGATGGTTTTTTGTATGTCTT EEEEEEEEEEEEEEEEEEEAEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:0 +NS500361:73:HJNVLBGXX:2:23108:15895:1433 16 Y 726482 255 28M * 0 0 ATTCAGGATGGTTTTTTGTATGTCTTTT EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:21503:1309:13837 16 Y 726484 255 29M * 0 0 TCAGGATGGTTTTTTGTATGTCTTTTCAA AEEEEEEEEEEEE/EEEAEEEEEEA//A/ NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:21502:26798:13836 16 Y 726484 255 29M * 0 0 TCAGGATGGTTTTTTGTATGTCTTTTCAA EEEEEEEEEEEEEAEEEEEE/EEAAAAA6 NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:4:21512:1491:12229 16 Y 726493 255 25M * 0 0 TTTTTTGTATGTCTTTTCAAAAACC EEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:24 nM:i:0 +NS500361:73:HJNVLBGXX:1:23302:25970:19144 16 Y 726494 255 28M * 0 0 TTTTTGTATGTCTTTTCAAAAACCTTTT EEEEEEEEEEEEEEEEEEEEEEAAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:21210:1750:11726 16 Y 726494 255 28M * 0 0 TTTTTGTATGTCTTTTCAAAAACCTTTT EEEEEEEEEEEEEEEEAEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:13204:21906:7243 16 Y 726496 255 28M * 0 0 TTTGTATGTCTTTTCAAAAACCTTTTGG EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:13110:3789:15320 16 Y 726496 255 28M * 0 0 TTTGTATGTCTTTTCAAAAACCTTTTGG EEEEEAEEEEEEEEEEEEEAEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:4:21505:11142:14865 16 Y 726496 255 28M * 0 0 TTTGTATGTCTTTTCAAAAACCTTTTGG EAEEEEEEEEE6EEEE/EEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:11211:19517:2361 16 Y 726496 255 27M * 0 0 TTTGTATGTCTTTTCAAAAACCTTTTG EEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:1:11212:19915:9098 16 Y 726496 255 27M * 0 0 TTTGTATGTCTTTTCAAAAACCTTTTG EEEEEEEEEEEEEEEEEAEEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:2:11303:17950:2657 16 Y 727567 255 28M * 0 0 TCTGTCAAGGCCACCAGACTCATCTTGG EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:12305:23810:16594 16 Y 727592 255 29M * 0 0 TGGTTTAAGAGAGCACTGGCAGTGACATT EE/EAEEEEEEE/EEAEEEEEAE/A/AAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:2:12310:11400:6992 16 Y 727592 255 29M * 0 0 TGGTTTAAGAGAGCACTGGCAGTGACATT EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:21508:2091:12222 16 Y 727593 255 28M * 0 0 GGTTTAAGAGAGCACTGGCAGTGACATT EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:21503:7990:1980 16 Y 727593 255 28M * 0 0 GGTTTAAGAGAGCACTGGCAGTGACATT EEEEEEEEEAEEEEEEEEAEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:4:21407:16520:14101 16 Y 727670 255 28M * 0 0 CATGATATAATGTTATGGGATAGGAAGT EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:21203:19260:3759 16 Y 729713 255 24M3191N3M * 0 0 CATCTTCATCTCCTACAATTACTTCAG EEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:2:22105:16050:10141 16 Y 729713 255 24M3191N4M * 0 0 CATCTTCATCTCCTACAATTACTTCAGT EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:11202:21522:13364 16 Y 729719 255 18M3191N10M * 0 0 CATCTCCTACAATTACTTCAGTATCTTC EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:11311:25158:12346 16 Y 729719 255 18M3191N10M * 0 0 CATCTCCTACAATTACTTCAGTATCTTC EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:12405:9825:17669 16 Y 729719 255 18M3191N10M * 0 0 CATCTCCTACAATTACTTCAGTATCTTC EEEEEEEEEEEEEE/EAEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:22303:9789:12799 16 Y 729719 255 18M3191N10M * 0 0 CATCTCCTACAATTACTTCAGTATCTTC EEEEEEEEEEEEEEEEEEEEEEEAAAA6 NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:4:21507:15236:13413 16 Y 729732 255 5M3191N26M * 0 0 TACTTCAGTATCTTCTTCTTCCTTTTGAGAA EEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:30 nM:i:0 +NS500361:73:HJNVLBGXX:4:12406:12329:11728 16 Y 732921 255 27M * 0 0 TACTTACCAGTATCTTCTTCTTCCTTT EEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:2:12106:24462:8540 16 Y 732921 255 27M * 0 0 TACTTACCAGTATCTTCTTCTTCCTTT EEEEEEEEEEEEEEAEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:2:12112:26279:6020 16 Y 732940 255 28M * 0 0 CTTCCTTTTGAGAATCACTGGCTGACAT EEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:21109:5985:3462 16 Y 732943 255 27M * 0 0 CCTTTTGAGAATCACTGGCTGACATAT E6EEEA/E/EEEEEEEEEEEAEAAAA/ NH:i:1 HI:i:1 AS:i:26 nM:i:0 +NS500361:73:HJNVLBGXX:2:21202:24251:6855 16 Y 732995 255 26M * 0 0 GCTAGTGCTCTGATCTGTGACTTCAG EEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:25 nM:i:0 +NS500361:73:HJNVLBGXX:1:11210:22310:7083 16 Y 733001 255 31M * 0 0 GCTCTGATCTGTGACTTCAGCTTCATTGCCA EEEEEEEAEEEEEEEEEEEEEEAEEEAAAAA NH:i:1 HI:i:1 AS:i:30 nM:i:0 +NS500361:73:HJNVLBGXX:2:13312:20748:11434 16 Y 733001 255 29M * 0 0 GCTCTGATCTGTGACTTCAGCTTCATTGC EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:1:22306:15193:19566 16 Y 733001 255 31M * 0 0 GCTCTGATCTGTGACTTCAGCTTCATTGCCA EEEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:30 nM:i:0 +NS500361:73:HJNVLBGXX:3:13408:1920:10276 16 Y 733014 255 29M * 0 0 ACTTCAGCTTCATTGCCATTTTTGGTCTC E6EE/EE6EEAEEEAEEEEEEEE/AA6AA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:23403:21966:2818 16 Y 733017 255 28M * 0 0 TCAGCTTCATTGCCATTTTTGGTCTCAC EEEEEEEEEEEEEEEEEE/EAAEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:13401:20351:7695 16 Y 733017 255 28M * 0 0 TCAGCTTCATTGCCATTTTTGGTCTCAC EEEEEEEEEEEEEEEEEAEEEAAAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:2:11110:11069:4039 16 Y 733020 255 29M * 0 0 GCTTCATTGCCATTTTTGGTCTCACTCTC EEEEEEEEEEEEEEEAE/EEEEEEAAAA6 NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:23409:8631:10130 16 Y 733029 255 29M * 0 0 CCATTTTTGGTCTCACTCTCTACTGCGTG EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:21401:11918:19891 16 Y 733032 255 29M * 0 0 TTTTTGGTCTCACTCTCTACTGCGTGTAT EEEEEEEAEEEE/A/EEEEEAEAE6AAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:4:23405:9031:19035 16 Y 733037 255 28M * 0 0 GGTCTCACTCTCTACTGCGTGTATAGTT EEAEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:3:13608:8602:1843 16 Y 733039 255 28M * 0 0 TCTCACTCTCTACTGCGTGTATAGTTTC EEEEEEEEEEAEEAE/AEEEAEEAAAAA NH:i:1 HI:i:1 AS:i:27 nM:i:0 +NS500361:73:HJNVLBGXX:1:11104:10544:8036 16 Y 735117 255 29M * 0 0 TTCATGGTAACTTCAGAGGAACCTTCAAG EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:2:21204:12349:5997 16 Y 735119 255 30M * 0 0 CATGGTAACTTCAGAGGAACCTTCAAGATC EEEEEEAE6EEEEEEAEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:29 nM:i:0 +NS500361:73:HJNVLBGXX:1:23103:4097:20225 16 Y 735119 255 30M * 0 0 CATGGTAACTTCAGAGGAACCTTCAAGATC EEEEEEEEE/EEEEEEAEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:29 nM:i:0 +NS500361:73:HJNVLBGXX:3:12505:9897:18665 16 Y 735119 255 30M * 0 0 CATGGTAACTTCAGAGGAACCTTCAAGATC EEEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:29 nM:i:0 +NS500361:73:HJNVLBGXX:2:13204:13512:5889 16 Y 735120 255 29M * 0 0 ATGGTAACTTCAGAGGAACCTTCAAGATC EEEEEEEEEEEEEEEEEEEEEEEEAAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:4:13409:6583:10081 16 Y 735120 255 29M * 0 0 ATGGTAACTTCAGAGGAACCTTCAAGATC EAAAEAE6EA//EEEE//E6EA//AAAAA NH:i:1 HI:i:1 AS:i:28 nM:i:0 +NS500361:73:HJNVLBGXX:3:13501:12970:19211 16 Y 735120 255 29M * 0 0 ATGGTAACTTCAGAGGAACCTTCAAGATC EEAEEEA/E/