# HG changeset patch # User iuc # Date 1515629382 18000 # Node ID a6477bafd52294ff3965a5f4e06974592b56ed52 planemo upload commit eea727c3bdfe36d9d16036d5ab79fb8b27c4e82e diff -r 000000000000 -r a6477bafd522 macros.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/macros.xml Wed Jan 10 19:09:42 2018 -0500 @@ -0,0 +1,23 @@ + + + + + 10.1101/gr.209601.116 + + @misc{githubUMI-tools, + title = {UMI-tools}, + publisher = {GitHub}, + journal = {GitHub repository}, + url = {https://github.com/CGATOxford/UMI-tools}, + } + + + + + + umi_tools + + + + 0.5.3 + diff -r 000000000000 -r a6477bafd522 test-data/dedup_out1.bam Binary file test-data/dedup_out1.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/dedup_out2.bam Binary file test-data/dedup_out2.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/dedup_out3.bam Binary file test-data/dedup_out3.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/dedup_out4.bam Binary file test-data/dedup_out4.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/dedup_out5.bam Binary file test-data/dedup_out5.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/dedup_out6.bam Binary file test-data/dedup_out6.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_in1.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/group_in1.sam Wed Jan 10 19:09:42 2018 -0500 @@ -0,0 +1,666 @@ +@HD VN:1.0 SO:coordinate +@PG ID:Bowtie VN:1.1.2 CL:"bowtie --wrapper basic-0 --threads 4 -v 2 -m 10 -a /ifs/mirror/genomes/bowtie/mm9 /dev/fd/63 --sam" +@SQ SN:chr10 LN:129993255 +@SQ SN:chr11 LN:121843856 +@SQ SN:chr12 LN:121257530 +@SQ SN:chr13 LN:120284312 +@SQ SN:chr14 LN:125194864 +@SQ SN:chr15 LN:103494974 +@SQ SN:chr16 LN:98319150 +@SQ SN:chr17 LN:95272651 +@SQ SN:chr18 LN:90772031 +@SQ SN:chr19 LN:61342430 +@SQ SN:chr1 LN:197195432 +@SQ SN:chr2 LN:181748087 +@SQ SN:chr3 LN:159599783 +@SQ SN:chr4 LN:155630120 +@SQ SN:chr5 LN:152537259 +@SQ SN:chr6 LN:149517037 +@SQ SN:chr7 LN:152524553 +@SQ SN:chr8 LN:131738871 +@SQ SN:chr9 LN:124076172 +@SQ SN:chrM LN:16299 +@SQ SN:chrX LN:166650296 +@SQ SN:chrY LN:15902555 +SRR2057595.10000392_AAGTA 16 chr19 31326011 255 57M * 0 0 * * XA:i:1 MD:Z:23A33 NM:i:1 UG:i:512 BX:Z:AAGTA +SRR2057595.10000468_GCTAT 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:706 BX:Z:GCTAT +SRR2057595.10000704_TGATG 16 chr19 3544146 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:1 BX:Z:TGATG +SRR2057595.10000838_TAGTA 0 chr19 10274608 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:324 BX:Z:TAGTA +SRR2057595.10001118_TCTCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:1 MD:Z:43A6 NM:i:1 UG:i:36 BX:Z:TCTCA +SRR2057595.10001329_CGGTC 16 chr19 28083193 255 60M * 0 0 * * XA:i:0 MD:Z:60 NM:i:0 UG:i:497 BX:Z:CGGTC +SRR2057595.10001382_CTGCC 16 chr19 3957282 255 22M * 0 0 * * XA:i:2 MD:Z:17C1G2 NM:i:2 UG:i:9 BX:Z:CTGCC +SRR2057595.10001382_CTGCC 16 chr19 3970897 255 22M * 0 0 * * XA:i:2 MD:Z:17C1G2 NM:i:2 UG:i:10 BX:Z:CTGCC +SRR2057595.10002047_ATTAT 16 chr19 7546213 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:278 BX:Z:ATTAT +SRR2057595.10002219_ATCAG 16 chr19 5796783 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:234 BX:Z:ATCAG +SRR2057595.10002597_AATTA 16 chr19 50733887 255 64M * 0 0 * * XA:i:2 MD:Z:58C0C4 NM:i:2 UG:i:645 BX:Z:AATTA +SRR2057595.10003084_GTCTA 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:124 BX:Z:GTCTA +SRR2057595.10003230_ACATT 0 chr19 8799370 255 64M * 0 0 * * XA:i:0 MD:Z:64 NM:i:0 UG:i:297 BX:Z:ACATT +SRR2057595.10003317_CGGCC 16 chr19 52641725 255 20M * 0 0 * * XA:i:2 MD:Z:1C13A4 NM:i:2 UG:i:648 BX:Z:CGGCC +SRR2057595.10003567_ATGTC 0 chr19 37052811 255 64M * 0 0 * * XA:i:1 MD:Z:11C52 NM:i:1 UG:i:542 BX:Z:ATGTC +SRR2057595.10003576_AATCG 16 chr19 4078303 255 67M * 0 0 * * XA:i:2 MD:Z:38A21A6 NM:i:2 UG:i:88 BX:Z:AATCG +SRR2057595.10003630_GTGTT 16 chr19 4078298 255 52M * 0 0 * * XA:i:2 MD:Z:42G0A8 NM:i:2 UG:i:45 BX:Z:GTGTT +SRR2057595.10003869_TCTTG 16 chr19 4078423 255 67M * 0 0 * * XA:i:1 MD:Z:19A47 NM:i:1 UG:i:111 BX:Z:TCTTG +SRR2057595.1000404_AACTT 16 chr19 4078420 255 64M * 0 0 * * XA:i:2 MD:Z:2T19A41 NM:i:2 UG:i:55 BX:Z:AACTT +SRR2057595.10004161_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.1000475_TCTGG 0 chr19 13129234 255 28M * 0 0 * * XA:i:1 MD:Z:22T5 NM:i:1 UG:i:383 BX:Z:TCTGG +SRR2057595.10004845_GGATA 0 chr19 24394180 255 29M * 0 0 * * XA:i:2 MD:Z:11C15T1 NM:i:2 UG:i:481 BX:Z:GGATA +SRR2057595.10005100_TTAAA 16 chr19 4785463 255 21M * 0 0 * * XA:i:0 MD:Z:21 NM:i:0 UG:i:196 BX:Z:TTAAA +SRR2057595.10005761_CTAGA 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:43 BX:Z:CTAGA +SRR2057595.10005808_CAGTA 16 chr19 7545872 255 24M * 0 0 * * XA:i:0 MD:Z:24 NM:i:0 UG:i:276 BX:Z:CAGTA +SRR2057595.10005955_GACAT 16 chr19 4078302 255 67M * 0 0 * * XA:i:2 MD:Z:39A21A5 NM:i:2 UG:i:82 BX:Z:GACAT +SRR2057595.10006011_TCTCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:1 MD:Z:43A6 NM:i:1 UG:i:36 BX:Z:TCTCA +SRR2057595.10006546_GCTAC 16 chr19 4078446 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:97 BX:Z:GCTAC +SRR2057595.1000663_ACGCA 16 chr19 51217013 255 67M * 0 0 * * XA:i:2 MD:Z:3A46A16 NM:i:2 UG:i:647 BX:Z:ACGCA +SRR2057595.10006870_GTCCA 16 chr19 4078431 255 67M * 0 0 * * XA:i:1 MD:Z:11A55 NM:i:1 UG:i:107 BX:Z:GTCCA +SRR2057595.10006926_GGATA 0 chr19 24394180 255 29M * 0 0 * * XA:i:2 MD:Z:11C15T1 NM:i:2 UG:i:481 BX:Z:GGATA +SRR2057595.10006986_GCGGA 16 chr19 4078298 255 48M * 0 0 * * XA:i:1 MD:Z:43A4 NM:i:1 UG:i:19 BX:Z:GCGGA +SRR2057595.10007554_TATTG 0 chr19 13129322 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:404 BX:Z:TATTG +SRR2057595.10008217_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10008433_AGAGT 16 chr19 24180225 255 55M * 0 0 * * XA:i:1 MD:Z:41C13 NM:i:1 UG:i:474 BX:Z:AGAGT +SRR2057595.10008484_GGATC 16 chr19 4078298 255 61M * 0 0 * * XA:i:2 MD:Z:42G0A17 NM:i:2 UG:i:61 BX:Z:GGATC +SRR2057595.10008523_TATGG 0 chr19 13129322 255 56M * 0 0 * * XA:i:2 MD:Z:33C5A16 NM:i:2 UG:i:405 BX:Z:TATGG +SRR2057595.10008761_AGCTG 16 chr19 4806926 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:201 BX:Z:AGCTG +SRR2057595.10009130_CGGCC 16 chr19 52641725 255 20M * 0 0 * * XA:i:2 MD:Z:1C13A4 NM:i:2 UG:i:648 BX:Z:CGGCC +SRR2057595.10009213_CCCCG 16 chr19 4078298 255 38M * 0 0 * * XA:i:1 MD:Z:35C2 NM:i:1 UG:i:184 BX:Z:CCCCG +SRR2057595.10009394_TAGAT 16 chr19 5797239 255 48M * 0 0 * * XA:i:0 MD:Z:48 NM:i:0 UG:i:229 BX:Z:TAGAT +SRR2057595.10009607_TTGCA 0 chr19 41640185 255 61M * 0 0 * * XA:i:0 MD:Z:61 NM:i:0 UG:i:581 BX:Z:TTGCA +SRR2057595.1001008_TATCA 16 chr19 60837460 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:691 BX:Z:TATCA +SRR2057595.10010444_AGTGC 0 chr19 8798525 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:296 BX:Z:AGTGC +SRR2057595.10011106_CGGTC 16 chr19 28083193 255 60M * 0 0 * * XA:i:0 MD:Z:60 NM:i:0 UG:i:497 BX:Z:CGGTC +SRR2057595.10011603_ATGAA 0 chr19 56177103 255 20M * 0 0 * * XA:i:2 MD:Z:12G0T6 NM:i:2 UG:i:667 BX:Z:ATGAA +SRR2057595.10011612_GTTTC 0 chr19 58015463 255 56M * 0 0 * * XA:i:0 MD:Z:56 NM:i:0 UG:i:683 BX:Z:GTTTC +SRR2057595.10011924_CCCCG 16 chr19 4078298 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:184 BX:Z:CCCCG +SRR2057595.10012143_ATGAA 0 chr19 56007800 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:666 BX:Z:ATGAA +SRR2057595.10012487_TGGGT 16 chr19 4917801 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:202 BX:Z:TGGGT +SRR2057595.1001256_TGATG 16 chr19 3544146 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:1 BX:Z:TGATG +SRR2057595.10013018_GGGGC 16 chr19 4078446 255 67M * 0 0 * * XA:i:1 MD:Z:66T0 NM:i:1 UG:i:98 BX:Z:GGGGC +SRR2057595.10013156_GTCAA 16 chr19 4078298 255 57M * 0 0 * * XA:i:1 MD:Z:43A13 NM:i:1 UG:i:52 BX:Z:GTCAA +SRR2057595.10013642_AATTC 0 chr19 13129322 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:403 BX:Z:AATTC +SRR2057595.10013703_ATGTA 0 chr19 39042685 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:558 BX:Z:ATGTA +SRR2057595.10014121_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.1001419_CAGTA 16 chr19 7545872 255 24M * 0 0 * * XA:i:0 MD:Z:24 NM:i:0 UG:i:276 BX:Z:CAGTA +SRR2057595.10014588_AAGGT 0 chr19 8798514 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:294 BX:Z:AAGGT +SRR2057595.10014904_GCGAC 0 chr19 13129553 255 26M * 0 0 * * XA:i:2 MD:Z:3A0T21 NM:i:2 UG:i:399 BX:Z:GCGAC +SRR2057595.10015153_ATTAT 16 chr19 7546213 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:278 BX:Z:ATTAT +SRR2057595.10015867_ATCCG 16 chr19 37369602 255 35M * 0 0 * * XA:i:0 MD:Z:35 NM:i:0 UG:i:544 BX:Z:ATCCG +SRR2057595.10015949_CTTTC 16 chr19 4078298 255 53M * 0 0 * * XA:i:2 MD:Z:42G0A9 NM:i:2 UG:i:47 BX:Z:CTTTC +SRR2057595.10016080_TATCA 16 chr19 60837460 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:691 BX:Z:TATCA +SRR2057595.10016534_ACGCA 16 chr19 51217013 255 67M * 0 0 * * XA:i:2 MD:Z:3A46A16 NM:i:2 UG:i:647 BX:Z:ACGCA +SRR2057595.10016692_TGGAC 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:127 BX:Z:TGGAC +SRR2057595.10016737_GTTGT 16 chr19 12564367 255 21M * 0 0 * * XA:i:2 MD:Z:0A12A7 NM:i:2 UG:i:380 BX:Z:GTTGT +SRR2057595.10016874_TTGCA 0 chr19 44434953 255 43M * 0 0 * * XA:i:0 MD:Z:43 NM:i:0 UG:i:602 BX:Z:TTGCA +SRR2057595.10017358_AGATC 16 chr19 24352786 255 67M * 0 0 * * XA:i:1 MD:Z:25C41 NM:i:1 UG:i:479 BX:Z:AGATC +SRR2057595.10017433_GGCAG 0 chr19 13129322 255 22M * 0 0 * * XA:i:0 MD:Z:22 NM:i:0 UG:i:409 BX:Z:GGCAG +SRR2057595.10017742_TAGCG 0 chr19 25429456 255 20M * 0 0 * * XA:i:2 MD:Z:5A9A4 NM:i:2 UG:i:485 BX:Z:TAGCG +SRR2057595.10018091_TAACA 0 chr19 8799023 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:299 BX:Z:TAACA +SRR2057595.10018143_CCATC 0 chr19 5008676 255 30M * 0 0 * * XA:i:2 MD:Z:1A23G4 NM:i:2 UG:i:203 BX:Z:CCATC +SRR2057595.10018237_GCTTG 16 chr19 5797471 255 50M * 0 0 * * XA:i:0 MD:Z:50 NM:i:0 UG:i:226 BX:Z:GCTTG +SRR2057595.10018343_AATTC 0 chr19 13129322 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:403 BX:Z:AATTC +SRR2057595.10018510_TTCGT 0 chr19 47167097 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:623 BX:Z:TTCGT +SRR2057595.10018668_GAGAG 16 chr19 4078298 255 50M * 0 0 * * XA:i:1 MD:Z:43A6 NM:i:1 UG:i:35 BX:Z:GAGAG +SRR2057595.10018780_CTAGA 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:43 BX:Z:CTAGA +SRR2057595.10018930_TGATT 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:41 BX:Z:TGATT +SRR2057595.10019300_TTATA 16 chr19 4078298 255 49M * 0 0 * * XA:i:2 MD:Z:42G0A5 NM:i:2 UG:i:22 BX:Z:TTATA +SRR2057595.10019327_TTATA 16 chr19 4078298 255 49M * 0 0 * * XA:i:2 MD:Z:42G0A5 NM:i:2 UG:i:22 BX:Z:TTATA +SRR2057595.10019367_AGTAT 0 chr19 10600420 255 29M * 0 0 * * XA:i:0 MD:Z:29 NM:i:0 UG:i:328 BX:Z:AGTAT +SRR2057595.10019428_CTAAT 0 chr19 40194243 255 44M * 0 0 * * XA:i:2 MD:Z:5T20A17 NM:i:2 UG:i:562 BX:Z:CTAAT +SRR2057595.10019544_CCCCG 16 chr19 4078298 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:184 BX:Z:CCCCG +SRR2057595.100195_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.10019604_TCCAG 0 chr19 29662201 255 26M * 0 0 * * XA:i:1 MD:Z:0C25 NM:i:1 UG:i:509 BX:Z:TCCAG +SRR2057595.10019935_ACCGT 16 chr19 15980484 255 29M * 0 0 * * XA:i:1 MD:Z:12C16 NM:i:1 UG:i:431 BX:Z:ACCGT +SRR2057595.10020101_CTAGA 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:43 BX:Z:CTAGA +SRR2057595.10020134_CTTTC 16 chr19 4078298 255 45M * 0 0 * * XA:i:1 MD:Z:43A1 NM:i:1 UG:i:15 BX:Z:CTTTC +SRR2057595.10020463_GTCCT 16 chr19 59694502 255 21M * 0 0 * * XA:i:2 MD:Z:8T0A11 NM:i:2 UG:i:688 BX:Z:GTCCT +SRR2057595.10020843_TTAAA 16 chr19 4785463 255 21M * 0 0 * * XA:i:0 MD:Z:21 NM:i:0 UG:i:196 BX:Z:TTAAA +SRR2057595.10020890_AGTAT 0 chr19 10600420 255 29M * 0 0 * * XA:i:0 MD:Z:29 NM:i:0 UG:i:328 BX:Z:AGTAT +SRR2057595.10021356_CAGTA 16 chr19 61240266 255 26M * 0 0 * * XA:i:1 MD:Z:12C13 NM:i:1 UG:i:697 BX:Z:CAGTA +SRR2057595.10021429_CACTG 16 chr19 5798673 255 42M * 0 0 * * XA:i:0 MD:Z:42 NM:i:0 UG:i:249 BX:Z:CACTG +SRR2057595.10021850_AGTCT 0 chr19 13129558 255 37M * 0 0 * * XA:i:1 MD:Z:21A15 NM:i:1 UG:i:398 BX:Z:AGTCT +SRR2057595.10022647_TTGTT 16 chr19 61240268 255 24M * 0 0 * * XA:i:1 MD:Z:10C13 NM:i:1 UG:i:696 BX:Z:TTGTT +SRR2057595.10022675_GTGGA 16 chr19 35902240 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:535 BX:Z:GTGGA +SRR2057595.10023366_AGCTT 16 chr19 5799262 255 20M * 0 0 * * XA:i:0 MD:Z:20 NM:i:0 UG:i:245 BX:Z:AGCTT +SRR2057595.10023594_ATGAA 0 chr19 56007800 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:666 BX:Z:ATGAA +SRR2057595.10023941_TACGG 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:1T43A21 NM:i:2 UG:i:71 BX:Z:TACGG +SRR2057595.10023975_TCGTC 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:704 BX:Z:TCGTC +SRR2057595.10024016_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10024824_ACCGG 16 chr19 10924297 255 49M * 0 0 * * XA:i:0 MD:Z:49 NM:i:0 UG:i:336 BX:Z:ACCGG +SRR2057595.10025477_AGTAT 0 chr19 10600420 255 29M * 0 0 * * XA:i:0 MD:Z:29 NM:i:0 UG:i:328 BX:Z:AGTAT +SRR2057595.10025641_CGAAC 0 chr19 13129273 255 36M * 0 0 * * XA:i:2 MD:Z:25C1T8 NM:i:2 UG:i:386 BX:Z:CGAAC +SRR2057595.10025783_GCATA 16 chr19 4078298 255 57M * 0 0 * * XA:i:2 MD:Z:42G0A13 NM:i:2 UG:i:53 BX:Z:GCATA +SRR2057595.10025842_GATCT 16 chr19 4078433 255 67M * 0 0 * * XA:i:1 MD:Z:9A57 NM:i:1 UG:i:90 BX:Z:GATCT +SRR2057595.10026676_GGATC 16 chr19 4078298 255 61M * 0 0 * * XA:i:2 MD:Z:42G0A17 NM:i:2 UG:i:61 BX:Z:GGATC +SRR2057595.1002669_TATTC 16 chr19 5038295 255 24M * 0 0 * * XA:i:0 MD:Z:24 NM:i:0 UG:i:206 BX:Z:TATTC +SRR2057595.10027010_TGTTT 16 chr19 50526830 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:643 BX:Z:TGTTT +SRR2057595.10027016_GACGT 0 chr19 40194275 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:565 BX:Z:GACGT +SRR2057595.10027043_TATGG 0 chr19 13129322 255 56M * 0 0 * * XA:i:2 MD:Z:33C5A16 NM:i:2 UG:i:405 BX:Z:TATGG +SRR2057595.10027109_GGGGC 16 chr19 4078446 255 67M * 0 0 * * XA:i:1 MD:Z:66T0 NM:i:1 UG:i:98 BX:Z:GGGGC +SRR2057595.10027286_ATGAA 0 chr19 56177103 255 20M * 0 0 * * XA:i:2 MD:Z:12G0T6 NM:i:2 UG:i:667 BX:Z:ATGAA +SRR2057595.10027736_CCGAT 16 chr19 4078300 255 49M * 0 0 * * XA:i:2 MD:Z:40G0A7 NM:i:2 UG:i:44 BX:Z:CCGAT +SRR2057595.10027811_AAAGG 16 chr19 23347776 255 34M * 0 0 * * XA:i:1 MD:Z:9T24 NM:i:1 UG:i:467 BX:Z:AAAGG +SRR2057595.10028205_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10028311_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10029121_AGCCG 0 chr19 8799400 255 33M * 0 0 * * XA:i:0 MD:Z:33 NM:i:0 UG:i:302 BX:Z:AGCCG +SRR2057595.10029169_GGTTC 16 chr19 56896659 255 23M * 0 0 * * XA:i:2 MD:Z:0C0A21 NM:i:2 UG:i:674 BX:Z:GGTTC +SRR2057595.10029364_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10029373_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10029475_ACAAA 16 chr19 5688438 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:222 BX:Z:ACAAA +SRR2057595.10029622_TCAGG 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:407 BX:Z:TCAGG +SRR2057595.10029643_TATGG 0 chr19 13129322 255 56M * 0 0 * * XA:i:2 MD:Z:33C5A16 NM:i:2 UG:i:405 BX:Z:TATGG +SRR2057595.10030278_TATTG 0 chr19 13129322 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:404 BX:Z:TATTG +SRR2057595.10030560_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.10030639_GCGAC 0 chr19 13129553 255 26M * 0 0 * * XA:i:2 MD:Z:3A0T21 NM:i:2 UG:i:399 BX:Z:GCGAC +SRR2057595.10031002_AAGGT 0 chr19 8798514 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:294 BX:Z:AAGGT +SRR2057595.10031358_TAGTA 16 chr19 61240268 255 27M * 0 0 * * XA:i:1 MD:Z:10C16 NM:i:1 UG:i:702 BX:Z:TAGTA +SRR2057595.10031368_TAACA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:32 BX:Z:TAACA +SRR2057595.10032378_CAGCA 16 chr19 5801892 255 47M * 0 0 * * XA:i:0 MD:Z:47 NM:i:0 UG:i:258 BX:Z:CAGCA +SRR2057595.10032451_ACCGT 16 chr19 15980484 255 29M * 0 0 * * XA:i:1 MD:Z:12C16 NM:i:1 UG:i:431 BX:Z:ACCGT +SRR2057595.10032626_CATAA 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:401 BX:Z:CATAA +SRR2057595.10033730_AACTG 16 chr19 40362914 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:571 BX:Z:AACTG +SRR2057595.10034125_CTCTT 16 chr19 10761035 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:332 BX:Z:CTCTT +SRR2057595.1003424_GGATA 0 chr19 24394180 255 29M * 0 0 * * XA:i:2 MD:Z:11C15T1 NM:i:2 UG:i:481 BX:Z:GGATA +SRR2057595.10034498_ATGAA 0 chr19 56177103 255 20M * 0 0 * * XA:i:2 MD:Z:12G0T6 NM:i:2 UG:i:667 BX:Z:ATGAA +SRR2057595.10034756_ATGAA 0 chr19 56007800 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:666 BX:Z:ATGAA +SRR2057595.10034981_TTCGT 0 chr19 47167097 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:623 BX:Z:TTCGT +SRR2057595.10035130_GGATG 0 chr19 13129322 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:402 BX:Z:GGATG +SRR2057595.10035284_TACGT 0 chr19 40194308 255 34M * 0 0 * * XA:i:1 MD:Z:3A30 NM:i:1 UG:i:563 BX:Z:TACGT +SRR2057595.1003559_CCATA 0 chr19 13129322 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:408 BX:Z:CCATA +SRR2057595.10035727_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.1003581_TAACA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:32 BX:Z:TAACA +SRR2057595.10036346_AGTCG 0 chr19 17767763 255 67M * 0 0 * * XA:i:2 MD:Z:10A4A51 NM:i:2 UG:i:440 BX:Z:AGTCG +SRR2057595.10036388_CCAAG 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:37 BX:Z:CCAAG +SRR2057595.10037180_GGTAG 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:177 BX:Z:GGTAG +SRR2057595.10037421_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10037481_CATAA 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:401 BX:Z:CATAA +SRR2057595.10037635_TTCGT 0 chr19 47167097 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:623 BX:Z:TTCGT +SRR2057595.10037650_GATGA 0 chr19 48194137 255 59M * 0 0 * * XA:i:0 MD:Z:59 NM:i:0 UG:i:635 BX:Z:GATGA +SRR2057595.10037690_TCCCG 16 chr19 56004832 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:663 BX:Z:TCCCG +SRR2057595.10038156_AAACA 16 chr19 4078299 255 49M * 0 0 * * XA:i:1 MD:Z:42A6 NM:i:1 UG:i:33 BX:Z:AAACA +SRR2057595.10038276_TAGTA 0 chr19 10274608 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:324 BX:Z:TAGTA +SRR2057595.10038572_TGTTT 0 chr19 41595518 255 38M * 0 0 * * XA:i:1 MD:Z:25G12 NM:i:1 UG:i:579 BX:Z:TGTTT +SRR2057595.10038660_CAGAT 0 chr19 13129299 255 32M * 0 0 * * XA:i:2 MD:Z:1T18G11 NM:i:2 UG:i:396 BX:Z:CAGAT +SRR2057595.1003901_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:173 BX:Z:AGCCG +SRR2057595.1003905_CCTCA 0 chr19 8798505 255 48M * 0 0 * * XA:i:0 MD:Z:48 NM:i:0 UG:i:292 BX:Z:CCTCA +SRR2057595.10039161_CTTAC 16 chr19 4078298 255 54M * 0 0 * * XA:i:2 MD:Z:42G0A10 NM:i:2 UG:i:49 BX:Z:CTTAC +SRR2057595.10039249_TCTGT 16 chr19 45885126 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:608 BX:Z:TCTGT +SRR2057595.10039623_CGCCG 16 chr19 4078299 255 59M * 0 0 * * XA:i:1 MD:Z:42A16 NM:i:1 UG:i:58 BX:Z:CGCCG +SRR2057595.1004011_GTCGC 0 chr19 12080548 255 51M * 0 0 * * XA:i:2 MD:Z:28C20G1 NM:i:2 UG:i:345 BX:Z:GTCGC +SRR2057595.1004011_GTCGC 16 chr19 12084189 255 51M * 0 0 * * XA:i:2 MD:Z:1C20G28 NM:i:2 UG:i:352 BX:Z:GTCGC +SRR2057595.10040530_CCTTG 16 chr19 4078298 255 49M * 0 0 * * XA:i:1 MD:Z:43A5 NM:i:1 UG:i:24 BX:Z:CCTTG +SRR2057595.10041125_GACGC 16 chr19 32419974 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:514 BX:Z:GACGC +SRR2057595.10041133_ACAGA 16 chr19 4078299 255 22M * 0 0 * * XA:i:2 MD:Z:0C0C20 NM:i:2 UG:i:136 BX:Z:ACAGA +SRR2057595.10041495_TGTTT 0 chr19 41595518 255 38M * 0 0 * * XA:i:1 MD:Z:25G12 NM:i:1 UG:i:579 BX:Z:TGTTT +SRR2057595.10041520_TGACA 16 chr19 51217013 255 43M * 0 0 * * XA:i:2 MD:Z:3A1C37 NM:i:2 UG:i:646 BX:Z:TGACA +SRR2057595.10041577_GTCAA 16 chr19 4078298 255 57M * 0 0 * * XA:i:1 MD:Z:43A13 NM:i:1 UG:i:52 BX:Z:GTCAA +SRR2057595.10041629_GATGA 0 chr19 48194137 255 59M * 0 0 * * XA:i:0 MD:Z:59 NM:i:0 UG:i:635 BX:Z:GATGA +SRR2057595.10041662_GGGGC 16 chr19 4078446 255 67M * 0 0 * * XA:i:1 MD:Z:66T0 NM:i:1 UG:i:98 BX:Z:GGGGC +SRR2057595.10042115_GACAT 16 chr19 4078302 255 67M * 0 0 * * XA:i:2 MD:Z:39A21A5 NM:i:2 UG:i:82 BX:Z:GACAT +SRR2057595.10042422_TCTGG 0 chr19 13129234 255 28M * 0 0 * * XA:i:1 MD:Z:22T5 NM:i:1 UG:i:383 BX:Z:TCTGG +SRR2057595.1004285_CCAAC 0 chr19 8820047 255 33M * 0 0 * * XA:i:0 MD:Z:33 NM:i:0 UG:i:315 BX:Z:CCAAC +SRR2057595.10042887_AAAGG 16 chr19 23347776 255 34M * 0 0 * * XA:i:1 MD:Z:9T24 NM:i:1 UG:i:467 BX:Z:AAAGG +SRR2057595.10043446_CAGAC 0 chr19 5493782 255 41M * 0 0 * * XA:i:0 MD:Z:41 NM:i:0 UG:i:218 BX:Z:CAGAC +SRR2057595.10043768_GGCGT 16 chr19 5797126 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:233 BX:Z:GGCGT +SRR2057595.10043877_TGTTC 16 chr19 44470840 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:605 BX:Z:TGTTC +SRR2057595.10043959_GACGT 0 chr19 40194275 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:565 BX:Z:GACGT +SRR2057595.10044281_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10044428_TCTAC 0 chr19 8798544 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:288 BX:Z:TCTAC +SRR2057595.10044882_TTAAA 16 chr19 4785463 255 21M * 0 0 * * XA:i:0 MD:Z:21 NM:i:0 UG:i:196 BX:Z:TTAAA +SRR2057595.10044970_GTGGA 16 chr19 35902240 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:535 BX:Z:GTGGA +SRR2057595.10045235_CGGCC 16 chr19 52641725 255 20M * 0 0 * * XA:i:2 MD:Z:1C13A4 NM:i:2 UG:i:648 BX:Z:CGGCC +SRR2057595.10045410_TGGTC 16 chr19 4078419 255 67M * 0 0 * * XA:i:2 MD:Z:3T19A43 NM:i:2 UG:i:119 BX:Z:TGGTC +SRR2057595.10045875_TTGGC 16 chr19 36085526 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:536 BX:Z:TTGGC +SRR2057595.10045996_TAACA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:32 BX:Z:TAACA +SRR2057595.10046273_TCCAG 0 chr19 29662201 255 26M * 0 0 * * XA:i:1 MD:Z:0C25 NM:i:1 UG:i:509 BX:Z:TCCAG +SRR2057595.10047157_GTGGA 16 chr19 35902240 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:535 BX:Z:GTGGA +SRR2057595.10047180_CCTCA 0 chr19 8798505 255 48M * 0 0 * * XA:i:0 MD:Z:48 NM:i:0 UG:i:292 BX:Z:CCTCA +SRR2057595.10047439_AATTA 16 chr19 50733890 255 61M * 0 0 * * XA:i:2 MD:Z:55C0C4 NM:i:2 UG:i:645 BX:Z:AATTA +SRR2057595.10047481_ACGCA 16 chr19 5796074 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:230 BX:Z:ACGCA +SRR2057595.10047803_TTCGT 16 chr19 40658723 255 67M * 0 0 * * XA:i:1 MD:Z:53T13 NM:i:1 UG:i:576 BX:Z:TTCGT +SRR2057595.10047915_TCACA 0 chr19 13129273 255 45M * 0 0 * * XA:i:2 MD:Z:25C1T17 NM:i:2 UG:i:388 BX:Z:TCACA +SRR2057595.10048372_GTCCC 0 chr19 12347657 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:372 BX:Z:GTCCC +SRR2057595.10048372_GTCCC 0 chr19 12348043 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:374 BX:Z:GTCCC +SRR2057595.10048372_GTCCC 16 chr19 15035071 255 32M * 0 0 * * XA:i:1 MD:Z:29T2 NM:i:1 UG:i:428 BX:Z:GTCCC +SRR2057595.10048721_CTAAT 0 chr19 40194243 255 44M * 0 0 * * XA:i:2 MD:Z:5T20A17 NM:i:2 UG:i:562 BX:Z:CTAAT +SRR2057595.10048861_CTAGA 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:43 BX:Z:CTAGA +SRR2057595.10048933_GTCGC 0 chr19 12080548 255 51M * 0 0 * * XA:i:2 MD:Z:28C20G1 NM:i:2 UG:i:345 BX:Z:GTCGC +SRR2057595.10048933_GTCGC 16 chr19 12084189 255 51M * 0 0 * * XA:i:2 MD:Z:1C20G28 NM:i:2 UG:i:352 BX:Z:GTCGC +SRR2057595.10049852_AAAGG 16 chr19 23347776 255 34M * 0 0 * * XA:i:1 MD:Z:9T24 NM:i:1 UG:i:467 BX:Z:AAAGG +SRR2057595.10049896_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.10050290_AAAGA 0 chr19 42921257 255 44M * 0 0 * * XA:i:2 MD:Z:4G9G29 NM:i:2 UG:i:590 BX:Z:AAAGA +SRR2057595.10050342_TGAAC 16 chr19 5798677 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:250 BX:Z:TGAAC +SRR2057595.10050663_GACGT 0 chr19 40194275 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:565 BX:Z:GACGT +SRR2057595.10051008_TTCGT 16 chr19 40658723 255 67M * 0 0 * * XA:i:1 MD:Z:53T13 NM:i:1 UG:i:576 BX:Z:TTCGT +SRR2057595.10051043_GCCGA 16 chr19 4078298 255 64M * 0 0 * * XA:i:1 MD:Z:43A20 NM:i:1 UG:i:66 BX:Z:GCCGA +SRR2057595.1005223_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.10052510_ATACG 0 chr19 13129421 255 31M * 0 0 * * XA:i:1 MD:Z:11A19 NM:i:1 UG:i:393 BX:Z:ATACG +SRR2057595.10052666_ATCCG 16 chr19 37369602 255 35M * 0 0 * * XA:i:0 MD:Z:35 NM:i:0 UG:i:544 BX:Z:ATCCG +SRR2057595.10053188_CTTTC 16 chr19 4078298 255 45M * 0 0 * * XA:i:1 MD:Z:43A1 NM:i:1 UG:i:15 BX:Z:CTTTC +SRR2057595.10053341_TTGCA 0 chr19 41640185 255 61M * 0 0 * * XA:i:0 MD:Z:61 NM:i:0 UG:i:581 BX:Z:TTGCA +SRR2057595.10053375_ATGAA 0 chr19 56007800 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:666 BX:Z:ATGAA +SRR2057595.10053481_CTCCG 0 chr19 11195966 255 44M * 0 0 * * XA:i:1 MD:Z:0T43 NM:i:1 UG:i:338 BX:Z:CTCCG +SRR2057595.10053518_AATCG 16 chr19 4078303 255 67M * 0 0 * * XA:i:2 MD:Z:38A21A6 NM:i:2 UG:i:88 BX:Z:AATCG +SRR2057595.1005354_GAGAG 16 chr19 4078298 255 50M * 0 0 * * XA:i:1 MD:Z:43A6 NM:i:1 UG:i:35 BX:Z:GAGAG +SRR2057595.1005379_CCAAC 0 chr19 8820047 255 33M * 0 0 * * XA:i:0 MD:Z:33 NM:i:0 UG:i:315 BX:Z:CCAAC +SRR2057595.10053918_AGTAT 0 chr19 10600420 255 29M * 0 0 * * XA:i:0 MD:Z:29 NM:i:0 UG:i:328 BX:Z:AGTAT +SRR2057595.10053970_TATGG 0 chr19 13129322 255 56M * 0 0 * * XA:i:2 MD:Z:33C5A16 NM:i:2 UG:i:405 BX:Z:TATGG +SRR2057595.10054046_ATAAA 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:175 BX:Z:ATAAA +SRR2057595.10054497_AAAGA 0 chr19 42921257 255 44M * 0 0 * * XA:i:2 MD:Z:4G9G29 NM:i:2 UG:i:590 BX:Z:AAAGA +SRR2057595.10055409_TAGAC 0 chr19 23806079 255 67M * 0 0 * * XA:i:1 MD:Z:41T25 NM:i:1 UG:i:468 BX:Z:TAGAC +SRR2057595.10055838_CAGAC 0 chr19 5493782 255 41M * 0 0 * * XA:i:0 MD:Z:41 NM:i:0 UG:i:218 BX:Z:CAGAC +SRR2057595.10055912_TCTCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:1 MD:Z:43A6 NM:i:1 UG:i:36 BX:Z:TCTCA +SRR2057595.10055955_AGGAT 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:70 BX:Z:AGGAT +SRR2057595.10056203_TGCAA 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:703 BX:Z:TGCAA +SRR2057595.10056249_TCAGG 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:407 BX:Z:TCAGG +SRR2057595.10057756_CGAAC 0 chr19 13129273 255 36M * 0 0 * * XA:i:2 MD:Z:25C1T8 NM:i:2 UG:i:386 BX:Z:CGAAC +SRR2057595.10058241_TTAAA 16 chr19 4785463 255 21M * 0 0 * * XA:i:0 MD:Z:21 NM:i:0 UG:i:196 BX:Z:TTAAA +SRR2057595.10058572_TCTTA 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:705 BX:Z:TCTTA +SRR2057595.10058802_TAGTA 0 chr19 10274608 255 34M * 0 0 * * XA:i:1 MD:Z:9C24 NM:i:1 UG:i:324 BX:Z:TAGTA +SRR2057595.10058946_AGCCA 0 chr19 43786861 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:594 BX:Z:AGCCA +SRR2057595.10059254_GGTTC 16 chr19 56896659 255 23M * 0 0 * * XA:i:2 MD:Z:0C0A21 NM:i:2 UG:i:674 BX:Z:GGTTC +SRR2057595.10059642_CTAGC 16 chr19 4078298 255 62M * 0 0 * * XA:i:1 MD:Z:43A18 NM:i:1 UG:i:63 BX:Z:CTAGC +SRR2057595.10060193_ATAGA 0 chr19 38795805 255 62M * 0 0 * * XA:i:2 MD:Z:44C1A15 NM:i:2 UG:i:553 BX:Z:ATAGA +SRR2057595.10060270_ATGGT 16 chr19 12086074 255 45M * 0 0 * * XA:i:1 MD:Z:19T25 NM:i:1 UG:i:361 BX:Z:ATGGT +SRR2057595.10060437_TATAT 0 chr19 40194298 255 25M * 0 0 * * XA:i:1 MD:Z:13A11 NM:i:1 UG:i:567 BX:Z:TATAT +SRR2057595.10060447_CTTAC 16 chr19 4078298 255 54M * 0 0 * * XA:i:2 MD:Z:42G0A10 NM:i:2 UG:i:49 BX:Z:CTTAC +SRR2057595.10060568_GACGC 16 chr19 32419974 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:514 BX:Z:GACGC +SRR2057595.10060757_AGCCG 16 chr19 4078299 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:11 BX:Z:AGCCG +SRR2057595.10060776_TCTTG 16 chr19 4078423 255 67M * 0 0 * * XA:i:1 MD:Z:19A47 NM:i:1 UG:i:111 BX:Z:TCTTG +SRR2057595.10061005_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:173 BX:Z:AGCCG +SRR2057595.10061018_CACTG 16 chr19 5798673 255 42M * 0 0 * * XA:i:0 MD:Z:42 NM:i:0 UG:i:249 BX:Z:CACTG +SRR2057595.10061192_ATCCG 16 chr19 37369602 255 35M * 0 0 * * XA:i:0 MD:Z:35 NM:i:0 UG:i:544 BX:Z:ATCCG +SRR2057595.10061371_GCCTG 16 chr19 4078301 255 67M * 0 0 * * XA:i:2 MD:Z:40A21A4 NM:i:2 UG:i:75 BX:Z:GCCTG +SRR2057595.10061403_AAGTA 16 chr19 31326011 255 57M * 0 0 * * XA:i:1 MD:Z:23A33 NM:i:1 UG:i:512 BX:Z:AAGTA +SRR2057595.10061424_GCTAC 16 chr19 4078446 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:97 BX:Z:GCTAC +SRR2057595.10061535_ATCAA 0 chr19 20547357 255 23M * 0 0 * * XA:i:0 MD:Z:23 NM:i:0 UG:i:452 BX:Z:ATCAA +SRR2057595.10061661_ACAAA 16 chr19 5688438 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:222 BX:Z:ACAAA +SRR2057595.10061745_ATATA 0 chr19 13129330 255 22M * 0 0 * * XA:i:1 MD:Z:0T21 NM:i:1 UG:i:422 BX:Z:ATATA +SRR2057595.10062673_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.10062948_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10063079_GTCTA 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:124 BX:Z:GTCTA +SRR2057595.10063521_AGATG 16 chr19 61240273 255 22M * 0 0 * * XA:i:1 MD:Z:5C16 NM:i:1 UG:i:707 BX:Z:AGATG +SRR2057595.10064036_CAGCA 16 chr19 5801892 255 47M * 0 0 * * XA:i:0 MD:Z:47 NM:i:0 UG:i:258 BX:Z:CAGCA +SRR2057595.1006417_GATCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:34 BX:Z:GATCA +SRR2057595.10065256_AGGGG 0 chr19 13129331 255 37M * 0 0 * * XA:i:2 MD:Z:24C5A6 NM:i:2 UG:i:423 BX:Z:AGGGG +SRR2057595.10065460_ATGGT 0 chr19 8798504 255 50M * 0 0 * * XA:i:0 MD:Z:50 NM:i:0 UG:i:291 BX:Z:ATGGT +SRR2057595.10065656_GATGA 0 chr19 48194137 255 59M * 0 0 * * XA:i:0 MD:Z:59 NM:i:0 UG:i:635 BX:Z:GATGA +SRR2057595.10065671_ATGAT 16 chr19 61240306 255 36M * 0 0 * * XA:i:2 MD:Z:16A4T14 NM:i:2 UG:i:716 BX:Z:ATGAT +SRR2057595.10065713_ATGTC 0 chr19 37052811 255 64M * 0 0 * * XA:i:1 MD:Z:11C52 NM:i:1 UG:i:542 BX:Z:ATGTC +SRR2057595.10065896_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10066051_AACTT 16 chr19 4078420 255 64M * 0 0 * * XA:i:2 MD:Z:2T19A41 NM:i:2 UG:i:55 BX:Z:AACTT +SRR2057595.10066334_AGCCA 0 chr19 47324044 255 46M * 0 0 * * XA:i:0 MD:Z:46 NM:i:0 UG:i:627 BX:Z:AGCCA +SRR2057595.10066646_GGTTC 16 chr19 56896659 255 23M * 0 0 * * XA:i:2 MD:Z:0C0A21 NM:i:2 UG:i:674 BX:Z:GGTTC +SRR2057595.10066775_GTCAA 16 chr19 4078298 255 57M * 0 0 * * XA:i:1 MD:Z:43A13 NM:i:1 UG:i:52 BX:Z:GTCAA +SRR2057595.10066849_GGTCG 0 chr19 53686815 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:651 BX:Z:GGTCG +SRR2057595.10066970_TAACA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:32 BX:Z:TAACA +SRR2057595.10067118_TGATT 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:41 BX:Z:TGATT +SRR2057595.10067277_TAGAC 0 chr19 23806079 255 67M * 0 0 * * XA:i:1 MD:Z:41T25 NM:i:1 UG:i:468 BX:Z:TAGAC +SRR2057595.10067458_TACGG 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:1T43A21 NM:i:2 UG:i:71 BX:Z:TACGG +SRR2057595.10067669_GCTTG 16 chr19 5797471 255 50M * 0 0 * * XA:i:0 MD:Z:50 NM:i:0 UG:i:226 BX:Z:GCTTG +SRR2057595.10067834_GACGT 0 chr19 40194275 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:565 BX:Z:GACGT +SRR2057595.1006784_GTCCA 16 chr19 4078431 255 67M * 0 0 * * XA:i:1 MD:Z:11A55 NM:i:1 UG:i:107 BX:Z:GTCCA +SRR2057595.10067875_GCCGA 16 chr19 4078298 255 64M * 0 0 * * XA:i:1 MD:Z:43A20 NM:i:1 UG:i:66 BX:Z:GCCGA +SRR2057595.1006816_GATCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:34 BX:Z:GATCA +SRR2057595.10068282_TAACA 0 chr19 8799023 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:299 BX:Z:TAACA +SRR2057595.1006840_TTAAA 16 chr19 4785463 255 21M * 0 0 * * XA:i:0 MD:Z:21 NM:i:0 UG:i:196 BX:Z:TTAAA +SRR2057595.10068743_TTGGC 16 chr19 36085526 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:536 BX:Z:TTGGC +SRR2057595.10068873_CTTTT 0 chr19 56656380 255 54M * 0 0 * * XA:i:0 MD:Z:54 NM:i:0 UG:i:669 BX:Z:CTTTT +SRR2057595.10068999_AAGTA 16 chr19 31326011 255 57M * 0 0 * * XA:i:1 MD:Z:23A33 NM:i:1 UG:i:512 BX:Z:AAGTA +SRR2057595.10069428_AACGC 16 chr19 4078298 255 51M * 0 0 * * XA:i:1 MD:Z:43A7 NM:i:1 UG:i:42 BX:Z:AACGC +SRR2057595.10070094_TCTTA 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:705 BX:Z:TCTTA +SRR2057595.10070448_GGATC 16 chr19 4078298 255 61M * 0 0 * * XA:i:2 MD:Z:42G0A17 NM:i:2 UG:i:61 BX:Z:GGATC +SRR2057595.1007050_TATAC 0 chr19 8799026 255 60M * 0 0 * * XA:i:0 MD:Z:60 NM:i:0 UG:i:301 BX:Z:TATAC +SRR2057595.10070871_GAAGT 16 chr19 27988575 255 52M * 0 0 * * XA:i:0 MD:Z:52 NM:i:0 UG:i:494 BX:Z:GAAGT +SRR2057595.10071095_TGCAA 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:703 BX:Z:TGCAA +SRR2057595.10071228_CTAGC 16 chr19 4078298 255 62M * 0 0 * * XA:i:1 MD:Z:43A18 NM:i:1 UG:i:63 BX:Z:CTAGC +SRR2057595.10071424_CATTT 0 chr19 49059158 255 49M * 0 0 * * XA:i:2 MD:Z:30A7A10 NM:i:2 UG:i:639 BX:Z:CATTT +SRR2057595.10071487_CACTG 16 chr19 5798673 255 42M * 0 0 * * XA:i:0 MD:Z:42 NM:i:0 UG:i:249 BX:Z:CACTG +SRR2057595.10072063_CAGCA 16 chr19 5801892 255 47M * 0 0 * * XA:i:0 MD:Z:47 NM:i:0 UG:i:258 BX:Z:CAGCA +SRR2057595.10072068_GCCGA 16 chr19 4078298 255 64M * 0 0 * * XA:i:1 MD:Z:43A20 NM:i:1 UG:i:66 BX:Z:GCCGA +SRR2057595.10072172_TCTGC 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:178 BX:Z:TCTGC +SRR2057595.10072478_AGCTT 16 chr19 5799262 255 20M * 0 0 * * XA:i:0 MD:Z:20 NM:i:0 UG:i:245 BX:Z:AGCTT +SRR2057595.10072601_TTGGC 16 chr19 4078302 255 67M * 0 0 * * XA:i:2 MD:Z:39A21A5 NM:i:2 UG:i:79 BX:Z:TTGGC +SRR2057595.10072627_AATAC 0 chr19 12046308 255 33M * 0 0 * * XA:i:1 MD:Z:23A9 NM:i:1 UG:i:343 BX:Z:AATAC +SRR2057595.10072909_TCCAG 0 chr19 29662201 255 26M * 0 0 * * XA:i:1 MD:Z:0C25 NM:i:1 UG:i:509 BX:Z:TCCAG +SRR2057595.10072922_TGATT 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:41 BX:Z:TGATT +SRR2057595.10072936_TAGTA 16 chr19 61240268 255 27M * 0 0 * * XA:i:1 MD:Z:10C16 NM:i:1 UG:i:702 BX:Z:TAGTA +SRR2057595.10073020_TGTTT 0 chr19 41595518 255 38M * 0 0 * * XA:i:1 MD:Z:25G12 NM:i:1 UG:i:579 BX:Z:TGTTT +SRR2057595.10073067_CGTCT 0 chr19 13129329 255 35M * 0 0 * * XA:i:2 MD:Z:26C5A2 NM:i:2 UG:i:421 BX:Z:CGTCT +SRR2057595.10073092_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.1007313_ACATT 0 chr19 8799370 255 64M * 0 0 * * XA:i:0 MD:Z:64 NM:i:0 UG:i:297 BX:Z:ACATT +SRR2057595.10073348_CCTTG 16 chr19 4078298 255 49M * 0 0 * * XA:i:1 MD:Z:43A5 NM:i:1 UG:i:24 BX:Z:CCTTG +SRR2057595.10073515_GGATC 16 chr19 4078298 255 61M * 0 0 * * XA:i:2 MD:Z:42G0A17 NM:i:2 UG:i:61 BX:Z:GGATC +SRR2057595.10073784_TGATG 16 chr19 3544146 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:1 BX:Z:TGATG +SRR2057595.10073896_TTCGT 0 chr19 36922761 255 57M * 0 0 * * XA:i:1 MD:Z:5C51 NM:i:1 UG:i:538 BX:Z:TTCGT +SRR2057595.10073986_ATCTC 0 chr19 21910875 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:458 BX:Z:ATCTC +SRR2057595.10074218_AGTGC 0 chr19 8798525 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:296 BX:Z:AGTGC +SRR2057595.10074241_AAACA 16 chr19 4078299 255 49M * 0 0 * * XA:i:1 MD:Z:42A6 NM:i:1 UG:i:33 BX:Z:AAACA +SRR2057595.10074267_GCTAC 0 chr19 13129281 255 31M * 0 0 * * XA:i:2 MD:Z:17C1T11 NM:i:2 UG:i:391 BX:Z:GCTAC +SRR2057595.10074447_ATTAT 16 chr19 7546213 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:278 BX:Z:ATTAT +SRR2057595.10074845_GTTGT 16 chr19 12564367 255 21M * 0 0 * * XA:i:2 MD:Z:0A12A7 NM:i:2 UG:i:380 BX:Z:GTTGT +SRR2057595.10075349_ATGAT 16 chr19 61240306 255 36M * 0 0 * * XA:i:2 MD:Z:16A4T14 NM:i:2 UG:i:716 BX:Z:ATGAT +SRR2057595.10075433_GGGGC 16 chr19 4078446 255 67M * 0 0 * * XA:i:1 MD:Z:66T0 NM:i:1 UG:i:98 BX:Z:GGGGC +SRR2057595.10075536_GCATA 16 chr19 4078298 255 57M * 0 0 * * XA:i:2 MD:Z:42G0A13 NM:i:2 UG:i:53 BX:Z:GCATA +SRR2057595.10075745_GTATT 0 chr19 47323072 255 47M * 0 0 * * XA:i:1 MD:Z:0C46 NM:i:1 UG:i:628 BX:Z:GTATT +SRR2057595.10075769_ATTAT 16 chr19 7546213 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:278 BX:Z:ATTAT +SRR2057595.10075776_GCTAC 16 chr19 4078446 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:97 BX:Z:GCTAC +SRR2057595.10075922_AATAC 0 chr19 12046308 255 33M * 0 0 * * XA:i:1 MD:Z:23A9 NM:i:1 UG:i:343 BX:Z:AATAC +SRR2057595.10075965_GCATA 16 chr19 4078298 255 57M * 0 0 * * XA:i:2 MD:Z:42G0A13 NM:i:2 UG:i:53 BX:Z:GCATA +SRR2057595.10076063_AAATC 16 chr19 47891298 255 21M * 0 0 * * XA:i:2 MD:Z:4C11C4 NM:i:2 UG:i:633 BX:Z:AAATC +SRR2057595.10076706_TTGCA 0 chr19 44434953 255 43M * 0 0 * * XA:i:0 MD:Z:43 NM:i:0 UG:i:602 BX:Z:TTGCA +SRR2057595.1007711_TCCAG 0 chr19 29662201 255 26M * 0 0 * * XA:i:1 MD:Z:0C25 NM:i:1 UG:i:509 BX:Z:TCCAG +SRR2057595.10077136_ATGTC 0 chr19 37052811 255 64M * 0 0 * * XA:i:1 MD:Z:11C52 NM:i:1 UG:i:542 BX:Z:ATGTC +SRR2057595.10077390_GATCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:34 BX:Z:GATCA +SRR2057595.10077550_GATCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:34 BX:Z:GATCA +SRR2057595.10078118_AGCGC 16 chr19 27464368 255 52M * 0 0 * * XA:i:0 MD:Z:52 NM:i:0 UG:i:493 BX:Z:AGCGC +SRR2057595.10078219_GCATA 16 chr19 4078298 255 57M * 0 0 * * XA:i:2 MD:Z:42G0A13 NM:i:2 UG:i:53 BX:Z:GCATA +SRR2057595.10078231_ATGAT 16 chr19 61240306 255 36M * 0 0 * * XA:i:2 MD:Z:16A4T14 NM:i:2 UG:i:716 BX:Z:ATGAT +SRR2057595.10078601_AATGC 0 chr19 8798773 255 42M * 0 0 * * XA:i:0 MD:Z:42 NM:i:0 UG:i:306 BX:Z:AATGC +SRR2057595.10079068_TTGGC 16 chr19 4078302 255 67M * 0 0 * * XA:i:2 MD:Z:39A21A5 NM:i:2 UG:i:79 BX:Z:TTGGC +SRR2057595.10079368_GAGTG 0 chr19 8798502 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:289 BX:Z:GAGTG +SRR2057595.10079546_CGGCC 16 chr19 52641725 255 20M * 0 0 * * XA:i:2 MD:Z:1C13A4 NM:i:2 UG:i:648 BX:Z:CGGCC +SRR2057595.10079549_TCTTA 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:705 BX:Z:TCTTA +SRR2057595.10079658_CGAAC 0 chr19 13129273 255 36M * 0 0 * * XA:i:2 MD:Z:25C1T8 NM:i:2 UG:i:386 BX:Z:CGAAC +SRR2057595.10080019_TATGG 0 chr19 13129322 255 56M * 0 0 * * XA:i:2 MD:Z:33C5A16 NM:i:2 UG:i:405 BX:Z:TATGG +SRR2057595.10080347_GCAGA 0 chr19 8799844 255 53M * 0 0 * * XA:i:0 MD:Z:53 NM:i:0 UG:i:309 BX:Z:GCAGA +SRR2057595.10080699_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10080920_GCTAT 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:706 BX:Z:GCTAT +SRR2057595.10080946_TATAT 0 chr19 40194298 255 25M * 0 0 * * XA:i:1 MD:Z:13A11 NM:i:1 UG:i:567 BX:Z:TATAT +SRR2057595.10081070_CGCCG 16 chr19 4078299 255 59M * 0 0 * * XA:i:1 MD:Z:42A16 NM:i:1 UG:i:58 BX:Z:CGCCG +SRR2057595.10081274_AACGC 16 chr19 4078298 255 51M * 0 0 * * XA:i:1 MD:Z:43A7 NM:i:1 UG:i:42 BX:Z:AACGC +SRR2057595.10081914_AATAC 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:125 BX:Z:AATAC +SRR2057595.10082001_TCTTG 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:705 BX:Z:TCTTA +SRR2057595.1008200_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.1008230_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10082591_ACAGA 16 chr19 4078299 255 22M * 0 0 * * XA:i:2 MD:Z:0C0C20 NM:i:2 UG:i:136 BX:Z:ACAGA +SRR2057595.10082733_CCTTT 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:180 BX:Z:CCTTT +SRR2057595.1008274_ATCAA 0 chr19 20547357 255 23M * 0 0 * * XA:i:0 MD:Z:23 NM:i:0 UG:i:452 BX:Z:ATCAA +SRR2057595.100827_ACAGA 16 chr19 4078299 255 22M * 0 0 * * XA:i:2 MD:Z:0C0C20 NM:i:2 UG:i:136 BX:Z:ACAGA +SRR2057595.10082802_GCCGT 0 chr19 30861987 255 20M * 0 0 * * XA:i:2 MD:Z:3C2T13 NM:i:2 UG:i:511 BX:Z:GCCGT +SRR2057595.10082815_CCTCA 0 chr19 8798505 255 48M * 0 0 * * XA:i:0 MD:Z:48 NM:i:0 UG:i:292 BX:Z:CCTCA +SRR2057595.1008318_CTATT 16 chr19 56874049 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:672 BX:Z:CTATT +SRR2057595.10084177_ATCAA 0 chr19 20547357 255 23M * 0 0 * * XA:i:0 MD:Z:23 NM:i:0 UG:i:452 BX:Z:ATCAA +SRR2057595.1008424_GATAA 16 chr19 47057431 255 22M * 0 0 * * XA:i:2 MD:Z:6T10T4 NM:i:2 UG:i:622 BX:Z:GATAA +SRR2057595.10084297_CTTCA 0 chr19 13129483 255 30M * 0 0 * * XA:i:2 MD:Z:6G16A6 NM:i:2 UG:i:382 BX:Z:CTTCA +SRR2057595.10084300_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:173 BX:Z:AGCCG +SRR2057595.10084403_GACGT 0 chr19 40194275 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:565 BX:Z:GACGT +SRR2057595.10084641_ACGCA 16 chr19 51217013 255 67M * 0 0 * * XA:i:2 MD:Z:3A46A16 NM:i:2 UG:i:647 BX:Z:ACGCA +SRR2057595.10084680_CCCCG 16 chr19 4078298 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:184 BX:Z:CCCCG +SRR2057595.1008532_CCGCG 0 chr19 8991126 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:318 BX:Z:CCGCG +SRR2057595.1008533_AAGGT 0 chr19 8798514 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:294 BX:Z:AAGGT +SRR2057595.10085421_CGAGA 0 chr19 6375345 255 21M * 0 0 * * XA:i:2 MD:Z:10C5G4 NM:i:2 UG:i:265 BX:Z:CGAGA +SRR2057595.10085567_TATTC 16 chr19 5038295 255 24M * 0 0 * * XA:i:0 MD:Z:24 NM:i:0 UG:i:206 BX:Z:TATTC +SRR2057595.10085614_AATTA 16 chr19 50733889 255 62M * 0 0 * * XA:i:2 MD:Z:56C0C4 NM:i:2 UG:i:645 BX:Z:AATTA +SRR2057595.10085673_GAAGT 16 chr19 8360277 255 67M * 0 0 * * XA:i:2 MD:Z:10A52A3 NM:i:2 UG:i:285 BX:Z:GAAGT +SRR2057595.10086547_TTCGT 0 chr19 47167097 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:623 BX:Z:TTCGT +SRR2057595.10086746_TCTTC 16 chr19 4078299 255 36M * 0 0 * * XA:i:1 MD:Z:27A8 NM:i:1 UG:i:178 BX:Z:TCTGC +SRR2057595.1008697_TTGCA 0 chr19 44434953 255 43M * 0 0 * * XA:i:0 MD:Z:43 NM:i:0 UG:i:602 BX:Z:TTGCA +SRR2057595.10087261_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10087496_CACTG 16 chr19 5798673 255 42M * 0 0 * * XA:i:0 MD:Z:42 NM:i:0 UG:i:249 BX:Z:CACTG +SRR2057595.10087607_ATACG 0 chr19 13129421 255 31M * 0 0 * * XA:i:1 MD:Z:11A19 NM:i:1 UG:i:393 BX:Z:ATACG +SRR2057595.10087628_GGCGT 16 chr19 5797126 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:233 BX:Z:GGCGT +SRR2057595.10087735_TCTGC 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:178 BX:Z:TCTGC +SRR2057595.10088373_AAAGG 16 chr19 23347776 255 34M * 0 0 * * XA:i:1 MD:Z:9T24 NM:i:1 UG:i:467 BX:Z:AAAGG +SRR2057595.1008837_GCGGA 16 chr19 4078298 255 48M * 0 0 * * XA:i:1 MD:Z:43A4 NM:i:1 UG:i:19 BX:Z:GCGGA +SRR2057595.10088477_AGTCG 0 chr19 17767763 255 67M * 0 0 * * XA:i:2 MD:Z:10A4A51 NM:i:2 UG:i:440 BX:Z:AGTCG +SRR2057595.1008854_TTCGT 0 chr19 47167097 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:623 BX:Z:TTCGT +SRR2057595.10088899_CATAA 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:401 BX:Z:CATAA +SRR2057595.10088953_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.10088992_TATTC 16 chr19 5038295 255 24M * 0 0 * * XA:i:0 MD:Z:24 NM:i:0 UG:i:206 BX:Z:TATTC +SRR2057595.10089226_AAACA 16 chr19 4078299 255 49M * 0 0 * * XA:i:1 MD:Z:42A6 NM:i:1 UG:i:33 BX:Z:AAACA +SRR2057595.10089360_ACCTG 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:126 BX:Z:ACCTG +SRR2057595.10089372_CAAAC 0 chr19 34304620 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:531 BX:Z:CAAAC +SRR2057595.10089645_CATTT 0 chr19 49059158 255 49M * 0 0 * * XA:i:2 MD:Z:30A7A10 NM:i:2 UG:i:639 BX:Z:CATTT +SRR2057595.10089847_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10090471_ATTCT 0 chr19 3576314 255 23M * 0 0 * * XA:i:2 MD:Z:16T5T0 NM:i:2 UG:i:5 BX:Z:ATTCT +SRR2057595.10090859_TTGCA 0 chr19 44434953 255 43M * 0 0 * * XA:i:0 MD:Z:43 NM:i:0 UG:i:602 BX:Z:TTGCA +SRR2057595.10090930_CCTCG 0 chr19 13129273 255 46M * 0 0 * * XA:i:2 MD:Z:25C1T18 NM:i:2 UG:i:387 BX:Z:CCTCG +SRR2057595.10090939_GATGG 0 chr19 44417644 255 23M * 0 0 * * XA:i:2 MD:Z:15T4C2 NM:i:2 UG:i:601 BX:Z:GATGG +SRR2057595.10091088_TCGCT 16 chr19 4078298 255 64M * 0 0 * * XA:i:2 MD:Z:42G0A20 NM:i:2 UG:i:67 BX:Z:TCGCT +SRR2057595.10091189_AGCAT 0 chr19 36976220 255 56M * 0 0 * * XA:i:0 MD:Z:56 NM:i:0 UG:i:539 BX:Z:AGCAT +SRR2057595.10091303_GTCCA 16 chr19 4078431 255 67M * 0 0 * * XA:i:1 MD:Z:11A55 NM:i:1 UG:i:107 BX:Z:GTCCA +SRR2057595.10091425_AGCTT 16 chr19 5799262 255 20M * 0 0 * * XA:i:0 MD:Z:20 NM:i:0 UG:i:245 BX:Z:AGCTT +SRR2057595.10091775_TAGTA 0 chr19 10274608 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:324 BX:Z:TAGTA +SRR2057595.10091818_TCTGC 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:178 BX:Z:TCTGC +SRR2057595.10092607_CATAA 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:401 BX:Z:CATAA +SRR2057595.10092724_CTGAT 0 chr19 42626015 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:587 BX:Z:CTGAT +SRR2057595.10092926_GCTAC 16 chr19 4078446 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:97 BX:Z:GCTAC +SRR2057595.1009294_CCAAC 0 chr19 8820047 255 33M * 0 0 * * XA:i:0 MD:Z:33 NM:i:0 UG:i:315 BX:Z:CCAAC +SRR2057595.10093528_ATACG 0 chr19 13129421 255 31M * 0 0 * * XA:i:1 MD:Z:11A19 NM:i:1 UG:i:393 BX:Z:ATACG +SRR2057595.1009371_GCGGA 16 chr19 4078298 255 48M * 0 0 * * XA:i:1 MD:Z:43A4 NM:i:1 UG:i:19 BX:Z:GCGGA +SRR2057595.10094012_GTGGA 16 chr19 35902240 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:535 BX:Z:GTGGA +SRR2057595.1009415_TTAAA 16 chr19 4785463 255 21M * 0 0 * * XA:i:0 MD:Z:21 NM:i:0 UG:i:196 BX:Z:TTAAA +SRR2057595.10094858_ATTAT 0 chr19 4755116 255 67M * 0 0 * * XA:i:2 MD:Z:42C11T12 NM:i:2 UG:i:195 BX:Z:ATTAT +SRR2057595.10094990_TCCAG 0 chr19 29662201 255 26M * 0 0 * * XA:i:1 MD:Z:0C25 NM:i:1 UG:i:509 BX:Z:TCCAG +SRR2057595.10095705_AACAG 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:174 BX:Z:AACAG +SRR2057595.10095797_CGTCT 0 chr19 13129329 255 35M * 0 0 * * XA:i:2 MD:Z:26C5A2 NM:i:2 UG:i:421 BX:Z:CGTCT +SRR2057595.10095837_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10096002_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:173 BX:Z:AGCCG +SRR2057595.10096480_TTCGT 16 chr19 40658723 255 67M * 0 0 * * XA:i:1 MD:Z:53T13 NM:i:1 UG:i:576 BX:Z:TTCGT +SRR2057595.10096536_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10096738_CTAGA 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:43 BX:Z:CTAGA +SRR2057595.10096905_AATTA 16 chr19 50733890 255 61M * 0 0 * * XA:i:2 MD:Z:55C0C4 NM:i:2 UG:i:645 BX:Z:AATTA +SRR2057595.10096962_TGTTT 0 chr19 41595518 255 38M * 0 0 * * XA:i:1 MD:Z:25G12 NM:i:1 UG:i:579 BX:Z:TGTTT +SRR2057595.10096976_AGATC 16 chr19 24352786 255 67M * 0 0 * * XA:i:1 MD:Z:25C41 NM:i:1 UG:i:479 BX:Z:AGATC +SRR2057595.10097042_CACTG 16 chr19 5798673 255 42M * 0 0 * * XA:i:0 MD:Z:42 NM:i:0 UG:i:249 BX:Z:CACTG +SRR2057595.10097390_CCATC 0 chr19 5008676 255 30M * 0 0 * * XA:i:2 MD:Z:1A23G4 NM:i:2 UG:i:203 BX:Z:CCATC +SRR2057595.10097856_TATGG 0 chr19 13129322 255 56M * 0 0 * * XA:i:2 MD:Z:33C5A16 NM:i:2 UG:i:405 BX:Z:TATGG +SRR2057595.10097927_GCATA 16 chr19 4078298 255 57M * 0 0 * * XA:i:2 MD:Z:42G0A13 NM:i:2 UG:i:53 BX:Z:GCATA +SRR2057595.10098678_CTCCG 0 chr19 11195966 255 43M * 0 0 * * XA:i:1 MD:Z:0T42 NM:i:1 UG:i:338 BX:Z:CTCCG +SRR2057595.10098703_TAGTA 0 chr19 10274608 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:324 BX:Z:TAGTA +SRR2057595.10098714_GCAGG 0 chr19 47903921 255 21M * 0 0 * * XA:i:2 MD:Z:2T0T17 NM:i:2 UG:i:634 BX:Z:GCAGG +SRR2057595.10098837_AATTA 16 chr19 50733888 255 63M * 0 0 * * XA:i:2 MD:Z:57C0C4 NM:i:2 UG:i:645 BX:Z:AATTA +SRR2057595.10098952_AAATC 16 chr19 47891298 255 21M * 0 0 * * XA:i:2 MD:Z:4C11C4 NM:i:2 UG:i:633 BX:Z:AAATC +SRR2057595.10098_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10099454_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10099862_TGATT 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:41 BX:Z:TGATT +SRR2057595.10099876_AGATC 16 chr19 24352786 255 67M * 0 0 * * XA:i:1 MD:Z:25C41 NM:i:1 UG:i:479 BX:Z:AGATC +SRR2057595.10100225_GAAGT 16 chr19 8360277 255 67M * 0 0 * * XA:i:2 MD:Z:10A52A3 NM:i:2 UG:i:285 BX:Z:GAAGT +SRR2057595.10100273_CCTCG 0 chr19 13129273 255 46M * 0 0 * * XA:i:2 MD:Z:25C1T18 NM:i:2 UG:i:387 BX:Z:CCTCG +SRR2057595.10100361_AGAAT 16 chr19 24246487 255 44M * 0 0 * * XA:i:0 MD:Z:44 NM:i:0 UG:i:476 BX:Z:AGAAT +SRR2057595.10100419_TCTTG 16 chr19 4078423 255 67M * 0 0 * * XA:i:1 MD:Z:19A47 NM:i:1 UG:i:111 BX:Z:TCTTG +SRR2057595.10100887_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.10101063_ATGAA 0 chr19 56007800 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:666 BX:Z:ATGAA +SRR2057595.10101443_TATAT 0 chr19 40194298 255 25M * 0 0 * * XA:i:1 MD:Z:13A11 NM:i:1 UG:i:567 BX:Z:TATAT +SRR2057595.10101541_TCTAC 0 chr19 8798544 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:288 BX:Z:TCTAC +SRR2057595.10101549_AAAGA 0 chr19 42921257 255 44M * 0 0 * * XA:i:2 MD:Z:4G9G29 NM:i:2 UG:i:590 BX:Z:AAAGA +SRR2057595.10102419_GGGGC 16 chr19 4078446 255 67M * 0 0 * * XA:i:1 MD:Z:66T0 NM:i:1 UG:i:98 BX:Z:GGGGC +SRR2057595.10102630_CTATT 16 chr19 56874049 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:672 BX:Z:CTATT +SRR2057595.10102792_GCTAA 16 chr19 32596026 255 56M * 0 0 * * XA:i:0 MD:Z:56 NM:i:0 UG:i:517 BX:Z:GCTAA +SRR2057595.10102901_AACTT 16 chr19 4078420 255 64M * 0 0 * * XA:i:2 MD:Z:2T19A41 NM:i:2 UG:i:55 BX:Z:AACTT +SRR2057595.10102918_ATTAT 16 chr19 7546213 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:278 BX:Z:ATTAT +SRR2057595.10103451_GATCA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:34 BX:Z:GATCA +SRR2057595.10103512_TCTTG 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:705 BX:Z:TCTTA +SRR2057595.10103749_ATCAA 0 chr19 20547357 255 23M * 0 0 * * XA:i:0 MD:Z:23 NM:i:0 UG:i:452 BX:Z:ATCAA +SRR2057595.10103791_ACGAC 16 chr19 4078298 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:186 BX:Z:ACGAC +SRR2057595.10104373_ACAGA 16 chr19 4078299 255 22M * 0 0 * * XA:i:2 MD:Z:0C0C20 NM:i:2 UG:i:136 BX:Z:ACAGA +SRR2057595.10104410_GCAGG 16 chr19 8963035 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:317 BX:Z:GCAGG +SRR2057595.10104500_GTATT 0 chr19 47323072 255 47M * 0 0 * * XA:i:1 MD:Z:0C46 NM:i:1 UG:i:628 BX:Z:GTATT +SRR2057595.10105258_CGCCG 16 chr19 4078299 255 59M * 0 0 * * XA:i:1 MD:Z:42A16 NM:i:1 UG:i:58 BX:Z:CGCCG +SRR2057595.10105817_GCGGA 16 chr19 4078298 255 48M * 0 0 * * XA:i:1 MD:Z:43A4 NM:i:1 UG:i:19 BX:Z:GCGGA +SRR2057595.10105818_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:1 MD:Z:0C35 NM:i:1 UG:i:173 BX:Z:AGCCG +SRR2057595.10106002_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.10106174_CTGCC 16 chr19 3957282 255 22M * 0 0 * * XA:i:2 MD:Z:17C1G2 NM:i:2 UG:i:9 BX:Z:CTGCC +SRR2057595.10106174_CTGCC 16 chr19 3970897 255 22M * 0 0 * * XA:i:2 MD:Z:17C1G2 NM:i:2 UG:i:10 BX:Z:CTGCC +SRR2057595.10106212_AGCGC 16 chr19 27464368 255 52M * 0 0 * * XA:i:0 MD:Z:52 NM:i:0 UG:i:493 BX:Z:AGCGC +SRR2057595.10106528_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.1010661_ATTAT 16 chr19 7546213 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:278 BX:Z:ATTAT +SRR2057595.10106731_CCTCA 0 chr19 8798505 255 48M * 0 0 * * XA:i:0 MD:Z:48 NM:i:0 UG:i:292 BX:Z:CCTCA +SRR2057595.10106761_CCTCG 0 chr19 13129273 255 46M * 0 0 * * XA:i:2 MD:Z:25C1T18 NM:i:2 UG:i:387 BX:Z:CCTCG +SRR2057595.10106812_CCAAG 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:37 BX:Z:CCAAG +SRR2057595.10106892_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:1 MD:Z:25T10 NM:i:1 UG:i:173 BX:Z:AGCCG +SRR2057595.10107056_GGCAG 0 chr19 13129322 255 22M * 0 0 * * XA:i:0 MD:Z:22 NM:i:0 UG:i:409 BX:Z:GGCAG +SRR2057595.1010759_TCTAA 0 chr19 12086447 255 46M * 0 0 * * XA:i:0 MD:Z:46 NM:i:0 UG:i:364 BX:Z:TCTAA +SRR2057595.1010759_TCTAA 0 chr19 12086782 255 46M * 0 0 * * XA:i:0 MD:Z:46 NM:i:0 UG:i:369 BX:Z:TCTAA +SRR2057595.10108058_TCTCG 16 chr19 6446235 255 26M * 0 0 * * XA:i:0 MD:Z:26 NM:i:0 UG:i:266 BX:Z:TCTCG +SRR2057595.10108108_TGTTT 0 chr19 41595519 255 37M * 0 0 * * XA:i:1 MD:Z:24G12 NM:i:1 UG:i:580 BX:Z:TGTTT +SRR2057595.10108134_TCTTG 16 chr19 4078423 255 67M * 0 0 * * XA:i:1 MD:Z:19A47 NM:i:1 UG:i:111 BX:Z:TCTTG +SRR2057595.10108322_TATTG 0 chr19 13129322 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:404 BX:Z:TATTG +SRR2057595.10108355_TAACA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:32 BX:Z:TAACA +SRR2057595.10108853_GGCGT 16 chr19 5797126 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:233 BX:Z:GGCGT +SRR2057595.10108990_TTGGC 16 chr19 36085526 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:536 BX:Z:TTGGC +SRR2057595.10109290_TATGA 16 chr19 4797429 255 49M * 0 0 * * XA:i:0 MD:Z:49 NM:i:0 UG:i:197 BX:Z:TATGA +SRR2057595.10109512_CTCTT 16 chr19 10761035 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:332 BX:Z:CTCTT +SRR2057595.10109901_ATGAA 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:175 BX:Z:ATAAA +SRR2057595.10110057_TTGTT 16 chr19 61240268 255 24M * 0 0 * * XA:i:1 MD:Z:10C13 NM:i:1 UG:i:696 BX:Z:TTGTT +SRR2057595.10110159_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10110319_TAGTA 16 chr19 61240268 255 27M * 0 0 * * XA:i:1 MD:Z:10C16 NM:i:1 UG:i:702 BX:Z:TAGTA +SRR2057595.10110664_TATTG 0 chr19 13129322 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:404 BX:Z:TATTG +SRR2057595.10110928_ATCCG 16 chr19 37369602 255 35M * 0 0 * * XA:i:0 MD:Z:35 NM:i:0 UG:i:544 BX:Z:ATCCG +SRR2057595.10111198_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.10111819_CCTTT 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:180 BX:Z:CCTTT +SRR2057595.1011210_GCTAC 16 chr19 4078446 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:97 BX:Z:GCTAC +SRR2057595.10112142_GATCT 16 chr19 4078433 255 67M * 0 0 * * XA:i:1 MD:Z:9A57 NM:i:1 UG:i:90 BX:Z:GATCT +SRR2057595.10112493_TCTGT 16 chr19 45885126 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:608 BX:Z:TCTGT +SRR2057595.10112719_GATTG 0 chr19 8799830 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:314 BX:Z:GATTG +SRR2057595.10113412_TGTTT 0 chr19 41595518 255 38M * 0 0 * * XA:i:1 MD:Z:25G12 NM:i:1 UG:i:579 BX:Z:TGTTT +SRR2057595.1011350_CATGC 16 chr19 42000773 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:584 BX:Z:CATGC +SRR2057595.10113554_TGACC 0 chr19 12086441 255 52M * 0 0 * * XA:i:0 MD:Z:52 NM:i:0 UG:i:362 BX:Z:TGACC +SRR2057595.10113554_TGACC 0 chr19 12086776 255 52M * 0 0 * * XA:i:0 MD:Z:52 NM:i:0 UG:i:367 BX:Z:TGACC +SRR2057595.10113832_CTAGA 0 chr19 7609388 255 42M * 0 0 * * XA:i:0 MD:Z:42 NM:i:0 UG:i:279 BX:Z:CTAGA +SRR2057595.10113868_ATGAA 0 chr19 56007800 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:666 BX:Z:ATGAA +SRR2057595.10113897_CATGC 16 chr19 42000773 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:584 BX:Z:CATGC +SRR2057595.10114178_TTAAA 16 chr19 4785463 255 21M * 0 0 * * XA:i:0 MD:Z:21 NM:i:0 UG:i:196 BX:Z:TTAAA +SRR2057595.10114345_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10114548_GATGA 16 chr19 55763541 255 22M * 0 0 * * XA:i:2 MD:Z:2C14G4 NM:i:2 UG:i:662 BX:Z:GATGA +SRR2057595.10114848_GTTGT 16 chr19 12564367 255 21M * 0 0 * * XA:i:2 MD:Z:0A12A7 NM:i:2 UG:i:380 BX:Z:GTTGT +SRR2057595.10115148_CTATT 16 chr19 56874049 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:672 BX:Z:CTATT +SRR2057595.10115906_TGGAC 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:127 BX:Z:TGGAC +SRR2057595.10116129_GACTC 16 chr19 25016726 255 31M * 0 0 * * XA:i:0 MD:Z:31 NM:i:0 UG:i:483 BX:Z:GACTC +SRR2057595.10116129_GACTC 16 chr19 60850233 255 31M * 0 0 * * XA:i:0 MD:Z:31 NM:i:0 UG:i:693 BX:Z:GACTC +SRR2057595.10116527_CTTAC 16 chr19 4078298 255 54M * 0 0 * * XA:i:2 MD:Z:42G0A10 NM:i:2 UG:i:49 BX:Z:CTTAC +SRR2057595.10116529_GTCCT 16 chr19 59694502 255 21M * 0 0 * * XA:i:2 MD:Z:8T0A11 NM:i:2 UG:i:688 BX:Z:GTCCT +SRR2057595.1011687_TAACA 16 chr19 4078298 255 50M * 0 0 * * XA:i:2 MD:Z:42G0A6 NM:i:2 UG:i:32 BX:Z:TAACA +SRR2057595.10116904_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10117105_CGCTT 0 chr19 13129322 255 26M * 0 0 * * XA:i:0 MD:Z:26 NM:i:0 UG:i:406 BX:Z:CGCTT +SRR2057595.10117768_TATCA 16 chr19 60837460 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:691 BX:Z:TATCA +SRR2057595.1011778_TGTTT 0 chr19 41595518 255 38M * 0 0 * * XA:i:1 MD:Z:25G12 NM:i:1 UG:i:579 BX:Z:TGTTT +SRR2057595.10117819_CTAGC 16 chr19 4078298 255 62M * 0 0 * * XA:i:2 MD:Z:21C21A18 NM:i:2 UG:i:63 BX:Z:CTAGC +SRR2057595.10117947_AATCG 16 chr19 4078303 255 67M * 0 0 * * XA:i:2 MD:Z:38A21A6 NM:i:2 UG:i:88 BX:Z:AATCG +SRR2057595.10118173_TTCGT 0 chr19 36922761 255 57M * 0 0 * * XA:i:1 MD:Z:5C51 NM:i:1 UG:i:538 BX:Z:TTCGT +SRR2057595.10118916_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:173 BX:Z:AGCCG +SRR2057595.10119454_TGTTT 0 chr19 41595518 255 38M * 0 0 * * XA:i:1 MD:Z:25G12 NM:i:1 UG:i:579 BX:Z:TGTTT +SRR2057595.10119776_GTCAA 16 chr19 4078298 255 57M * 0 0 * * XA:i:1 MD:Z:43A13 NM:i:1 UG:i:52 BX:Z:GTCAA +SRR2057595.10119820_CAGAC 0 chr19 5493782 255 41M * 0 0 * * XA:i:0 MD:Z:41 NM:i:0 UG:i:218 BX:Z:CAGAC +SRR2057595.10119919_GAATG 16 chr19 33191616 255 28M * 0 0 * * XA:i:2 MD:Z:1G17G8 NM:i:2 UG:i:526 BX:Z:GAATG +SRR2057595.10120218_TGGTA 16 chr19 28085954 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:498 BX:Z:TGGTA +SRR2057595.10120257_TCAGG 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:407 BX:Z:TCAGG +SRR2057595.10120271_GTCCC 0 chr19 12347657 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:372 BX:Z:GTCCC +SRR2057595.10120271_GTCCC 0 chr19 12348043 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:374 BX:Z:GTCCC +SRR2057595.10120271_GTCCC 16 chr19 15035071 255 32M * 0 0 * * XA:i:1 MD:Z:29T2 NM:i:1 UG:i:428 BX:Z:GTCCC +SRR2057595.10120377_TATTC 16 chr19 5038295 255 24M * 0 0 * * XA:i:0 MD:Z:24 NM:i:0 UG:i:206 BX:Z:TATTC +SRR2057595.10120569_GCATA 16 chr19 4078298 255 57M * 0 0 * * XA:i:2 MD:Z:42G0A13 NM:i:2 UG:i:53 BX:Z:GCATA +SRR2057595.10120734_AGTAA 16 chr19 4078423 255 67M * 0 0 * * XA:i:2 MD:Z:19A24T22 NM:i:2 UG:i:110 BX:Z:AGTAA +SRR2057595.10120758_AACGC 16 chr19 4078298 255 51M * 0 0 * * XA:i:1 MD:Z:43A7 NM:i:1 UG:i:42 BX:Z:AACGC +SRR2057595.10120905_ACCTG 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:126 BX:Z:ACCTG +SRR2057595.10121117_TGTTC 0 chr19 61241963 255 27M * 0 0 * * XA:i:2 MD:Z:20A5G0 NM:i:2 UG:i:719 BX:Z:TGTTC +SRR2057595.10121638_AATTA 16 chr19 50733888 255 63M * 0 0 * * XA:i:2 MD:Z:57C0C4 NM:i:2 UG:i:645 BX:Z:AATTA +SRR2057595.10122214_TATTG 0 chr19 13129322 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:404 BX:Z:TATTG +SRR2057595.10122328_TAAAG 16 chr19 61240268 255 28M * 0 0 * * XA:i:1 MD:Z:10C17 NM:i:1 UG:i:710 BX:Z:TAAAG +SRR2057595.10122369_AGCTT 16 chr19 5799262 255 20M * 0 0 * * XA:i:0 MD:Z:20 NM:i:0 UG:i:245 BX:Z:AGCTT +SRR2057595.10122580_AAGGT 0 chr19 8798514 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:294 BX:Z:AAGGT +SRR2057595.101226_AACTT 16 chr19 4078420 255 64M * 0 0 * * XA:i:2 MD:Z:2T19A41 NM:i:2 UG:i:55 BX:Z:AACTT +SRR2057595.10123331_ATTAT 16 chr19 7546213 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:278 BX:Z:ATTAT +SRR2057595.10123377_CTAAT 0 chr19 40194243 255 44M * 0 0 * * XA:i:2 MD:Z:5T20A17 NM:i:2 UG:i:562 BX:Z:CTAAT +SRR2057595.1012355_CAACT 16 chr19 4078298 255 49M * 0 0 * * XA:i:2 MD:Z:42G0A5 NM:i:2 UG:i:25 BX:Z:CAACT +SRR2057595.10123626_TCTAC 0 chr19 8798544 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:288 BX:Z:TCTAC +SRR2057595.10124059_GATGA 16 chr19 55763541 255 22M * 0 0 * * XA:i:2 MD:Z:2C14G4 NM:i:2 UG:i:662 BX:Z:GATGA +SRR2057595.1012426_CATTT 0 chr19 49059158 255 49M * 0 0 * * XA:i:2 MD:Z:30A7A10 NM:i:2 UG:i:639 BX:Z:CATTT +SRR2057595.10125107_TATGG 0 chr19 13129322 255 56M * 0 0 * * XA:i:2 MD:Z:33C5A16 NM:i:2 UG:i:405 BX:Z:TATGG +SRR2057595.1012543_GGATC 16 chr19 4078298 255 61M * 0 0 * * XA:i:2 MD:Z:42G0A17 NM:i:2 UG:i:61 BX:Z:GGATC +SRR2057595.10125580_TGCTA 0 chr19 10830415 255 21M * 0 0 * * XA:i:2 MD:Z:0A10C9 NM:i:2 UG:i:334 BX:Z:TGCTA +SRR2057595.1012559_TGATT 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:41 BX:Z:TGATT +SRR2057595.10125877_CTATT 16 chr19 56874049 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:672 BX:Z:CTATT +SRR2057595.10125930_ATCGG 0 chr19 8799853 255 44M * 0 0 * * XA:i:0 MD:Z:44 NM:i:0 UG:i:310 BX:Z:ATCGG +SRR2057595.10125966_GATGA 16 chr19 55763541 255 22M * 0 0 * * XA:i:2 MD:Z:2C14G4 NM:i:2 UG:i:662 BX:Z:GATGA +SRR2057595.10126475_TCAGG 0 chr19 38669715 255 54M * 0 0 * * XA:i:0 MD:Z:54 NM:i:0 UG:i:550 BX:Z:TCAGG +SRR2057595.10126817_CAGAC 0 chr19 5493782 255 41M * 0 0 * * XA:i:0 MD:Z:41 NM:i:0 UG:i:218 BX:Z:CAGAC +SRR2057595.1012699_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10127386_GCAGA 0 chr19 8799844 255 53M * 0 0 * * XA:i:0 MD:Z:53 NM:i:0 UG:i:309 BX:Z:GCAGA +SRR2057595.10127482_AAGGT 0 chr19 8798514 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:294 BX:Z:AAGGT +SRR2057595.10127633_GGACC 16 chr19 4078298 255 28M * 0 0 * * XA:i:1 MD:Z:18C9 NM:i:1 UG:i:143 BX:Z:GGACC +SRR2057595.10127800_CCTTT 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:180 BX:Z:CCTTT +SRR2057595.10127806_TTTTC 16 chr19 4078296 255 67M * 0 0 * * XA:i:2 MD:Z:0T44A21 NM:i:2 UG:i:69 BX:Z:TTTTC +SRR2057595.10127925_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10127969_TCAGC 16 chr19 4078298 255 49M * 0 0 * * XA:i:1 MD:Z:43A5 NM:i:1 UG:i:23 BX:Z:TCAGC +SRR2057595.10128074_TATTC 16 chr19 5038295 255 24M * 0 0 * * XA:i:0 MD:Z:24 NM:i:0 UG:i:206 BX:Z:TATTC +SRR2057595.10128154_CTAGG 16 chr19 5798663 255 67M * 0 0 * * XA:i:2 MD:Z:0G0G65 NM:i:2 UG:i:242 BX:Z:CTAGG +SRR2057595.10128297_ATAAA 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:175 BX:Z:ATAAA +SRR2057595.10128465_TTGCA 0 chr19 44434953 255 43M * 0 0 * * XA:i:0 MD:Z:43 NM:i:0 UG:i:602 BX:Z:TTGCA +SRR2057595.10128496_AGCCG 16 chr19 4078298 255 39M * 0 0 * * XA:i:0 MD:Z:39 NM:i:0 UG:i:11 BX:Z:AGCCG +SRR2057595.1012868_AAACA 16 chr19 4078299 255 49M * 0 0 * * XA:i:1 MD:Z:42A6 NM:i:1 UG:i:33 BX:Z:AAACA +SRR2057595.1012871_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.10128720_CCTCA 0 chr19 8798505 255 48M * 0 0 * * XA:i:0 MD:Z:48 NM:i:0 UG:i:292 BX:Z:CCTCA +SRR2057595.10128783_TTGCA 0 chr19 44434953 255 43M * 0 0 * * XA:i:0 MD:Z:43 NM:i:0 UG:i:602 BX:Z:TTGCA +SRR2057595.10128887_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.10128931_AGCCG 16 chr19 4078299 255 36M * 0 0 * * XA:i:0 MD:Z:36 NM:i:0 UG:i:173 BX:Z:AGCCG +SRR2057595.10129033_ACAGA 16 chr19 4078299 255 22M * 0 0 * * XA:i:2 MD:Z:0C0C20 NM:i:2 UG:i:136 BX:Z:ACAGA +SRR2057595.10129046_TCAAA 16 chr19 4078298 255 38M * 0 0 * * XA:i:0 MD:Z:38 NM:i:0 UG:i:185 BX:Z:TCAAA +SRR2057595.10129405_ATGAA 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:175 BX:Z:ATAAA +SRR2057595.1012950_GGGAT 0 chr19 13129424 255 25M * 0 0 * * XA:i:1 MD:Z:8A16 NM:i:1 UG:i:394 BX:Z:GGGAT +SRR2057595.10129521_CGCCG 16 chr19 4078299 255 59M * 0 0 * * XA:i:1 MD:Z:42A16 NM:i:1 UG:i:58 BX:Z:CGCCG +SRR2057595.10129734_TATTG 0 chr19 13129322 255 28M * 0 0 * * XA:i:0 MD:Z:28 NM:i:0 UG:i:404 BX:Z:TATTG +SRR2057595.10130199_GGCAG 0 chr19 13129322 255 22M * 0 0 * * XA:i:0 MD:Z:22 NM:i:0 UG:i:409 BX:Z:GGCAG +SRR2057595.10130333_GTCGC 0 chr19 12080548 255 51M * 0 0 * * XA:i:2 MD:Z:28C20G1 NM:i:2 UG:i:345 BX:Z:GTCGC +SRR2057595.10130333_GTCGC 16 chr19 12084189 255 51M * 0 0 * * XA:i:2 MD:Z:1C20G28 NM:i:2 UG:i:352 BX:Z:GTCGC +SRR2057595.10130767_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.10130878_TTAGT 16 chr19 4078302 255 67M * 0 0 * * XA:i:2 MD:Z:39A21A5 NM:i:2 UG:i:81 BX:Z:TTAGT +SRR2057595.10131031_AGCTT 16 chr19 5799262 255 20M * 0 0 * * XA:i:0 MD:Z:20 NM:i:0 UG:i:245 BX:Z:AGCTT +SRR2057595.10131148_CAGCA 16 chr19 5801892 255 47M * 0 0 * * XA:i:0 MD:Z:47 NM:i:0 UG:i:258 BX:Z:CAGCA +SRR2057595.10131269_CATGC 0 chr19 13129570 255 22M * 0 0 * * XA:i:1 MD:Z:9A12 NM:i:1 UG:i:392 BX:Z:CATGC +SRR2057595.10131482_CTTTC 16 chr19 4078298 255 53M * 0 0 * * XA:i:2 MD:Z:42G0A9 NM:i:2 UG:i:47 BX:Z:CTTTC +SRR2057595.10132002_CCCGT 16 chr19 47349292 255 19M * 0 0 * * XA:i:2 MD:Z:8A3G6 NM:i:2 UG:i:629 BX:Z:CCCGT +SRR2057595.10132034_ATACC 16 chr19 35673152 255 21M * 0 0 * * XA:i:2 MD:Z:7T2A10 NM:i:2 UG:i:534 BX:Z:ATACC +SRR2057595.10132108_TCTAC 0 chr19 8798544 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:288 BX:Z:TCTAC +SRR2057595.10132494_CGAAC 0 chr19 13129273 255 36M * 0 0 * * XA:i:2 MD:Z:25C1T8 NM:i:2 UG:i:386 BX:Z:CGAAC +SRR2057595.10133086_TCTTA 16 chr19 61240268 255 27M * 0 0 * * XA:i:1 MD:Z:10C16 NM:i:1 UG:i:705 BX:Z:TCTTA +SRR2057595.10133357_GTCGC 0 chr19 12080548 255 51M * 0 0 * * XA:i:2 MD:Z:28C20G1 NM:i:2 UG:i:345 BX:Z:GTCGC +SRR2057595.10133357_GTCGC 16 chr19 12084189 255 51M * 0 0 * * XA:i:2 MD:Z:1C20G28 NM:i:2 UG:i:352 BX:Z:GTCGC +SRR2057595.1013336_TTGTT 16 chr19 61240268 255 24M * 0 0 * * XA:i:1 MD:Z:10C13 NM:i:1 UG:i:696 BX:Z:TTGTT +SRR2057595.10133381_AGAGT 16 chr19 24180225 255 55M * 0 0 * * XA:i:1 MD:Z:41C13 NM:i:1 UG:i:474 BX:Z:AGAGT +SRR2057595.10133440_CTGCG 16 chr19 61240305 255 36M * 0 0 * * XA:i:2 MD:Z:17A4T13 NM:i:2 UG:i:712 BX:Z:CTGCG +SRR2057595.10134016_GACAC 0 chr19 12348541 255 26M * 0 0 * * XA:i:0 MD:Z:26 NM:i:0 UG:i:378 BX:Z:GACAC +SRR2057595.10134016_GACAC 16 chr19 15034579 255 26M * 0 0 * * XA:i:0 MD:Z:26 NM:i:0 UG:i:429 BX:Z:GACAC +SRR2057595.10134072_GCAGG 16 chr19 8963035 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:317 BX:Z:GCAGG +SRR2057595.1013452_ATATA 0 chr19 13129330 255 22M * 0 0 * * XA:i:1 MD:Z:0T21 NM:i:1 UG:i:422 BX:Z:ATATA +SRR2057595.10134979_AGTAT 0 chr19 10600420 255 29M * 0 0 * * XA:i:0 MD:Z:29 NM:i:0 UG:i:328 BX:Z:AGTAT +SRR2057595.10135304_GTCCA 16 chr19 4078431 255 67M * 0 0 * * XA:i:1 MD:Z:11A55 NM:i:1 UG:i:107 BX:Z:GTCCA +SRR2057595.1013548_TGTTT 16 chr19 50526830 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:643 BX:Z:TGTTT +SRR2057595.10135758_TTGCA 0 chr19 44434953 255 43M * 0 0 * * XA:i:0 MD:Z:43 NM:i:0 UG:i:602 BX:Z:TTGCA +SRR2057595.10135807_CCCGT 16 chr19 47349292 255 19M * 0 0 * * XA:i:2 MD:Z:8A3G6 NM:i:2 UG:i:629 BX:Z:CCCGT +SRR2057595.1013613_GATTG 0 chr19 55501486 255 31M * 0 0 * * XA:i:1 MD:Z:8A22 NM:i:1 UG:i:659 BX:Z:GATTG +SRR2057595.1013616_ATAGA 0 chr19 38795805 255 62M * 0 0 * * XA:i:1 MD:Z:46A15 NM:i:1 UG:i:553 BX:Z:ATAGA +SRR2057595.10136267_TTCGC 16 chr19 4078428 255 67M * 0 0 * * XA:i:2 MD:Z:14A24T27 NM:i:2 UG:i:74 BX:Z:TTCGC +SRR2057595.10136714_ATCCG 16 chr19 37369602 255 35M * 0 0 * * XA:i:0 MD:Z:35 NM:i:0 UG:i:544 BX:Z:ATCCG +SRR2057595.10136765_GCTAT 0 chr19 57701750 255 67M * 0 0 * * XA:i:2 MD:Z:34T1T30 NM:i:2 UG:i:677 BX:Z:GCTAT +SRR2057595.10137151_ATTCT 0 chr19 3576314 255 23M * 0 0 * * XA:i:2 MD:Z:16T5T0 NM:i:2 UG:i:5 BX:Z:ATTCT +SRR2057595.10137262_ACAAA 16 chr19 4078298 255 49M * 0 0 * * XA:i:1 MD:Z:43A5 NM:i:1 UG:i:27 BX:Z:ACAAA +SRR2057595.10137317_GATTG 0 chr19 8799830 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:314 BX:Z:GATTG +SRR2057595.10137603_GCTAT 16 chr19 61240267 255 28M * 0 0 * * XA:i:1 MD:Z:11C16 NM:i:1 UG:i:706 BX:Z:GCTAT +SRR2057595.10137788_TAGTA 0 chr19 10274608 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:324 BX:Z:TAGTA +SRR2057595.1013793_GCGCA 0 chr19 8798523 255 30M * 0 0 * * XA:i:0 MD:Z:30 NM:i:0 UG:i:295 BX:Z:GCGCA +SRR2057595.10138029_GTTGT 16 chr19 12564367 255 21M * 0 0 * * XA:i:2 MD:Z:0A12A7 NM:i:2 UG:i:380 BX:Z:GTTGT +SRR2057595.10138656_CTAGC 16 chr19 4078298 255 62M * 0 0 * * XA:i:1 MD:Z:43A18 NM:i:1 UG:i:63 BX:Z:CTAGC +SRR2057595.10138697_GCCGT 0 chr19 30861987 255 20M * 0 0 * * XA:i:2 MD:Z:3C2T13 NM:i:2 UG:i:511 BX:Z:GCCGT +SRR2057595.10138764_GATGG 0 chr19 44417644 255 23M * 0 0 * * XA:i:2 MD:Z:15T4C2 NM:i:2 UG:i:601 BX:Z:GATGG +SRR2057595.10138897_TGATT 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:41 BX:Z:TGATT +SRR2057595.10139239_GAGTG 0 chr19 8798502 255 51M * 0 0 * * XA:i:0 MD:Z:51 NM:i:0 UG:i:289 BX:Z:GAGTG +SRR2057595.10139499_AGTCT 0 chr19 13129558 255 37M * 0 0 * * XA:i:1 MD:Z:21A15 NM:i:1 UG:i:398 BX:Z:AGTCT +SRR2057595.10139596_GCTAC 16 chr19 4078446 255 67M * 0 0 * * XA:i:0 MD:Z:67 NM:i:0 UG:i:97 BX:Z:GCTAC +SRR2057595.10139652_GGATG 0 chr19 13129322 255 27M * 0 0 * * XA:i:0 MD:Z:27 NM:i:0 UG:i:402 BX:Z:GGATG +SRR2057595.10139674_ATGAA 16 chr19 4078298 255 37M * 0 0 * * XA:i:0 MD:Z:37 NM:i:0 UG:i:175 BX:Z:ATAAA +SRR2057595.10139952_TAGTA 16 chr19 61240268 255 27M * 0 0 * * XA:i:1 MD:Z:10C16 NM:i:1 UG:i:702 BX:Z:TAGTA +SRR2057595.10140058_ATATA 0 chr19 13129330 255 22M * 0 0 * * XA:i:1 MD:Z:0T21 NM:i:1 UG:i:422 BX:Z:ATATA +SRR2057595.10141122_GTCCC 0 chr19 12347657 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:372 BX:Z:GTCCC +SRR2057595.10141122_GTCCC 0 chr19 12348043 255 32M * 0 0 * * XA:i:0 MD:Z:32 NM:i:0 UG:i:374 BX:Z:GTCCC +SRR2057595.10141122_GTCCC 16 chr19 15035071 255 32M * 0 0 * * XA:i:1 MD:Z:29T2 NM:i:1 UG:i:428 BX:Z:GTCCC +SRR2057595.10141295_TAGTA 0 chr19 10274608 255 34M * 0 0 * * XA:i:0 MD:Z:34 NM:i:0 UG:i:324 BX:Z:TAGTA +SRR2057595.10142441_GTCTA 16 chr19 4078420 255 67M * 0 0 * * XA:i:2 MD:Z:2T19A44 NM:i:2 UG:i:124 BX:Z:GTCTA +SRR2057595.10142482_GCCGA 16 chr19 7389994 255 48M * 0 0 * * XA:i:1 MD:Z:47C0 NM:i:1 UG:i:274 BX:Z:GCCGA +SRR2057595.10143102_AAGTA 16 chr19 31326011 255 57M * 0 0 * * XA:i:1 MD:Z:23A33 NM:i:1 UG:i:512 BX:Z:AAGTA +SRR2057595.10143387_GCATA 16 chr19 4078298 255 57M * 0 0 * * XA:i:2 MD:Z:42G0A13 NM:i:2 UG:i:53 BX:Z:GCATA +SRR2057595.10143631_ATTTT 0 chr19 45095257 255 57M * 0 0 * * XA:i:0 MD:Z:57 NM:i:0 UG:i:606 BX:Z:ATTTT +SRR2057595.10143687_CGCTT 0 chr19 13129322 255 26M * 0 0 * * XA:i:0 MD:Z:26 NM:i:0 UG:i:406 BX:Z:CGCTT +SRR2057595.10143772_TGATT 16 chr19 4078299 255 50M * 0 0 * * XA:i:2 MD:Z:41G0A7 NM:i:2 UG:i:41 BX:Z:TGATT +SRR2057595.10143909_CAGAC 0 chr19 5493782 255 41M * 0 0 * * XA:i:0 MD:Z:41 NM:i:0 UG:i:218 BX:Z:CAGAC +SRR2057595.10143972_CTTAC 16 chr19 4078298 255 54M * 0 0 * * XA:i:2 MD:Z:42G0A10 NM:i:2 UG:i:49 BX:Z:CTTAC +SRR2057595.10144216_CATAA 0 chr19 13129322 255 35M * 0 0 * * XA:i:1 MD:Z:33C1 NM:i:1 UG:i:401 BX:Z:CATAA +SRR2057595.10144367_ATGAT 16 chr19 61240306 255 36M * 0 0 * * XA:i:2 MD:Z:16A4T14 NM:i:2 UG:i:716 BX:Z:ATGAT diff -r 000000000000 -r a6477bafd522 test-data/group_in2.bam Binary file test-data/group_in2.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_in3.bam Binary file test-data/group_in3.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_in4.bam Binary file test-data/group_in4.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_in5.bam Binary file test-data/group_in5.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_in6.bam Binary file test-data/group_in6.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_out2.bam Binary file test-data/group_out2.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_out3.bam Binary file test-data/group_out3.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_out3.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/group_out3.tab Wed Jan 10 19:09:42 2018 -0500 @@ -0,0 +1,498 @@ +read_id contig position gene umi umi_count final_umi final_umi_count unique_id +chrM_81_583_3:0:0_2:0:0_f4/1_TT chrM 80 NA TT 1 TT 1 0 +chrM_110_732_3:0:0_2:0:0_160/1_CC chrM 109 NA CC 1 CC 1 1 +chrM_118_613_6:0:0_4:0:0_169/1_AG chrM 117 NA AG 1 AG 1 2 +chrM_149_684_2:0:0_2:0:0_6e/1_TA chrM 148 NA TA 1 TA 1 3 +chrM_152_616_6:0:0_2:0:0_84/1_GC chrM 151 NA GC 1 GC 1 4 +chrM_247_748_8:0:0_1:0:0_1b7/1_GC chrM 246 NA GC 1 GC 1 5 +chrM_280_772_1:0:0_9:0:0_31/1_CG chrM 279 NA CG 1 CG 1 6 +chrM_292_819_2:0:0_1:0:0_18d/1_CA chrM 291 NA CA 1 CA 1 7 +chrM_390_890_3:0:0_6:0:0_bf/1_CA chrM 389 NA CA 1 CA 1 8 +chrM_447_921_2:0:0_3:0:0_1cc/1_AC chrM 446 NA AC 1 AC 1 9 +chrM_469_983_2:0:0_1:0:0_121/1_AC chrM 468 NA AC 1 AC 1 10 +chrM_541_1074_3:0:0_8:1:0_22/1_AC chrM 540 NA AC 1 AC 1 11 +chrM_8_556_3:0:0_2:0:0_1b1/1_AC chrM 556 NA AC 1 AC 1 12 +chrM_112_577_2:0:0_4:0:0_17b/1_CG chrM 577 NA CG 1 CG 1 13 +chrM_627_1063_6:0:0_4:1:0_12a/1_AA chrM 626 NA AA 1 AA 1 14 +chrM_164_650_5:0:0_3:0:0_164/1_AC chrM 650 NA AC 1 AC 1 15 +chrM_200_695_6:0:0_2:0:0_5a/1_TA chrM 695 NA TA 3 TA 3 16 +chrM_168_695_5:0:0_2:0:0_1af/1_TA chrM 695 NA TA 3 TA 3 16 +chrM_170_695_1:0:0_6:0:0_1e3/1_TA chrM 695 NA TA 3 TA 3 16 +chrM_209_705_3:0:0_2:0:0_1b3/1_TA chrM 705 NA TA 1 TA 1 17 +chrM_288_807_3:0:0_5:0:0_8e/1_AG chrM 807 NA AG 1 AG 1 18 +chrM_818_1274_1:1:0_4:0:0_124/1_TC chrM 817 NA TC 1 TC 1 19 +chrM_442_872_4:0:0_6:0:0_146/1_AT chrM 872 NA AT 1 AT 1 20 +chrM_460_919_3:0:0_1:0:0_0/1_GA chrM 919 NA GA 1 GA 1 21 +chrM_384_950_5:0:0_2:0:0_34/1_GG chrM 950 NA GG 1 GG 1 22 +chrM_476_980_5:0:0_2:0:0_133/1_TA chrM 980 NA TA 1 TA 1 23 +chrM_552_999_4:0:0_5:0:0_bb/1_AT chrM 999 NA AT 1 AT 1 24 +chrM_512_1030_5:0:0_7:1:0_a5/1_AC chrM 1030 NA AC 1 AC 1 25 +chrM_1103_1634_5:0:0_3:0:0_36/1_TA chrM 1102 NA TA 1 TA 1 26 +chrM_648_1148_3:0:0_3:1:0_1d1/1_TT chrM 1148 NA TT 1 TT 1 27 +chrM_1167_1797_8:0:0_2:1:0_14e/1_GT chrM 1166 NA GT 1 GT 1 28 +chrM_619_1169_3:0:0_5:1:0_33/1_CC chrM 1169 NA CC 1 CC 1 29 +chrM_668_1248_0:0:0_4:0:0_fa/1_TT chrM 1248 NA TT 1 TT 1 30 +chrM_1288_1759_4:0:0_2:1:0_a1/1_TA chrM 1287 NA TA 1 TA 1 31 +chrM_1327_1786_4:0:0_3:1:0_1ad/1_CA chrM 1326 NA CA 1 CA 1 32 +chrM_866_1351_5:1:0_6:0:0_1c9/1_TA chrM 1351 NA TA 1 TA 1 33 +chrM_853_1366_6:1:0_4:0:0_13b/1_TC chrM 1366 NA TC 1 TC 1 34 +chrM_1399_1851_2:0:0_11:1:0_85/1_AT chrM 1398 NA AT 1 AT 1 35 +chrM_946_1444_4:1:0_6:0:0_113/1_TT chrM 1444 NA TT 1 TT 1 36 +chrM_943_1485_2:1:0_5:0:0_1e2/1_AT chrM 1485 NA AT 1 AT 1 37 +chrM_1022_1501_4:0:0_5:0:0_132/1_GA chrM 1501 NA GA 1 GA 1 38 +chrM_1504_2013_4:1:0_3:0:0_10b/1_AT chrM 1503 NA AT 1 AT 1 39 +chrM_1505_1934_1:1:0_5:0:0_157/1_TC chrM 1504 NA TC 1 TC 1 40 +chrM_997_1511_2:1:0_1:0:0_1d8/1_AT chrM 1511 NA AT 1 AT 1 41 +chrM_1521_2070_5:1:0_3:0:0_1a/1_TA chrM 1520 NA TA 1 TA 1 42 +chrM_1575_2126_5:1:0_2:1:0_18c/1_TA chrM 1574 NA TA 1 TA 1 43 +chrM_1063_1598_5:0:0_5:0:0_f5/1_TG chrM 1598 NA TG 1 TG 1 44 +chrM_1605_2128_3:1:0_4:1:0_1ea/1_TT chrM 1604 NA TT 1 TT 1 45 +chrM_1065_1609_3:0:0_5:0:0_18e/1_TA chrM 1609 NA TA 1 TA 1 46 +chrM_1140_1619_1:0:0_4:0:0_74/1_TT chrM 1619 NA TT 1 TT 1 47 +chrM_1111_1626_4:0:0_2:0:0_186/1_AC chrM 1626 NA AC 1 AC 1 48 +chrM_1664_2135_0:1:0_3:1:0_179/1_TT chrM 1663 NA TT 1 TT 1 49 +chrM_1699_2147_1:0:0_7:0:0_10c/1_AA chrM 1698 NA AA 1 AA 1 50 +chrM_1706_2240_3:0:0_6:0:0_99/1_GT chrM 1705 NA GT 1 GT 1 51 +chrM_1756_2309_0:0:0_3:0:0_c3/1_AT chrM 1755 NA AT 1 AT 1 52 +chrM_1344_1758_4:0:0_3:1:0_75/1_TT chrM 1758 NA TT 1 TT 1 53 +chrM_1223_1777_5:0:0_6:1:0_42/1_TG chrM 1777 NA TG 1 TG 1 54 +chrM_1790_2351_4:0:0_4:0:0_f3/1_CC chrM 1789 NA CC 1 CC 1 55 +chrM_1308_1807_4:0:0_4:1:0_45/1_TA chrM 1807 NA TA 1 TA 1 56 +chrM_1814_2315_6:0:0_3:0:0_12d/1_AA chrM 1813 NA AA 1 AA 1 57 +chrM_1862_2304_4:0:0_1:0:0_19d/1_AA chrM 1861 NA AA 1 AA 1 58 +chrM_1363_1869_5:0:0_3:1:0_aa/1_TA chrM 1869 NA TA 1 TA 1 59 +chrM_1363_1887_2:0:0_3:0:0_d3/1_CA chrM 1887 NA CA 1 CA 1 60 +chrM_1428_1904_7:0:0_1:0:0_1b0/1_AC chrM 1904 NA AC 1 AC 1 61 +chrM_1480_1975_4:0:0_6:0:0_cd/1_CT chrM 1975 NA CT 1 CT 1 62 +chrM_1998_2552_3:0:0_2:0:0_35/1_AC chrM 1997 NA AC 1 AC 1 63 +chrM_1427_2035_6:0:0_1:0:0_7d/1_TA chrM 2035 NA TA 1 TA 1 64 +chrM_2065_2529_4:1:0_2:0:0_69/1_CG chrM 2064 NA CG 1 CG 1 65 +chrM_2079_2546_5:1:0_2:0:0_37/1_TA chrM 2078 NA TA 1 TA 1 66 +chrM_1576_2102_2:1:0_1:0:0_3/1_CT chrM 2102 NA CT 1 CT 1 67 +chrM_2118_2553_2:0:0_6:0:0_e6/1_CC chrM 2117 NA CC 1 CC 1 68 +chrM_1707_2180_3:0:0_2:0:0_152/1_TT chrM 2180 NA TT 1 TT 1 69 +chrM_1950_2421_3:0:0_3:0:0_13/1_TC chrM 2421 NA TC 1 TC 1 70 +chrM_2426_2953_5:0:0_5:0:0_83/1_AT chrM 2425 NA AT 1 AT 1 71 +chrM_2497_3082_1:0:0_6:0:0_1a4/1_AG chrM 2496 NA AG 1 AG 1 72 +chrM_2021_2503_3:1:0_4:0:0_4d/1_AC chrM 2503 NA AC 1 AC 1 73 +chrM_2510_3002_4:0:0_1:0:0_1d/1_AA chrM 2509 NA AA 1 AA 1 74 +chrM_2528_3123_6:1:0_4:0:0_122/1_AT chrM 2527 NA AT 1 AT 1 75 +chrM_2646_3108_4:1:0_3:0:0_68/1_GA chrM 2645 NA GA 1 GA 1 76 +chrM_2063_2659_1:1:0_8:0:0_93/1_GG chrM 2659 NA GG 1 GG 1 77 +chrM_2160_2679_2:0:0_4:0:0_177/1_TC chrM 2679 NA TC 1 TC 1 78 +chrM_2823_3302_7:0:0_9:0:0_de/1_CC chrM 2822 NA CC 1 CC 1 79 +chrM_2441_2926_5:0:0_3:0:0_25/1_AA chrM 2926 NA AA 1 AA 1 80 +chrM_2958_3490_6:0:0_5:1:0_192/1_TT chrM 2957 NA TT 1 TT 1 81 +chrM_2988_3473_2:0:0_4:1:0_c6/1_GT chrM 2987 NA GT 1 GT 1 82 +chrM_3032_3426_3:0:0_3:1:0_17c/1_AG chrM 3031 NA AG 1 AG 1 83 +chrM_2563_3061_1:1:0_3:0:0_70/1_GG chrM 3061 NA GG 1 GG 1 84 +chrM_2692_3152_1:1:0_4:0:0_5d/1_TC chrM 3152 NA TC 1 TC 1 85 +chrM_2612_3153_2:1:0_1:0:0_1df/1_AT chrM 3153 NA AT 1 AT 1 86 +chrM_3202_3692_4:0:0_2:0:0_95/1_AA chrM 3201 NA AA 1 AA 1 87 +chrM_2733_3206_4:0:0_1:0:0_4e/1_GT chrM 3206 NA GT 1 GT 1 88 +chrM_2706_3256_4:0:0_9:0:0_16e/1_GA chrM 3256 NA GA 1 GA 1 89 +chrM_2864_3323_4:0:0_5:0:0_fb/1_CT chrM 3323 NA CT 1 CT 1 90 +chrM_2856_3324_2:0:0_6:0:0_4/1_GC chrM 3324 NA GC 1 GC 1 91 +chrM_2902_3333_5:0:0_3:0:0_4a/1_CC chrM 3333 NA CC 1 CC 1 92 +chrM_2890_3410_3:0:0_6:1:0_1b5/1_TA chrM 3410 NA TA 1 TA 1 93 +chrM_2875_3422_5:0:0_7:1:0_c5/1_GG chrM 3422 NA GG 1 GG 1 94 +chrM_2890_3430_4:0:0_8:1:0_19a/1_TT chrM 3430 NA TT 1 TT 1 95 +chrM_2949_3443_10:0:0_6:1:0_139/1_CT chrM 3443 NA CT 1 CT 1 96 +chrM_2872_3457_8:0:0_3:1:0_134/1_TG chrM 3457 NA TG 1 TG 1 97 +chrM_2948_3460_5:0:0_7:1:0_1c1/1_TG chrM 3460 NA TG 1 TG 1 98 +chrM_3015_3496_4:0:0_5:1:0_137/1_AT chrM 3496 NA AT 1 AT 1 99 +chrM_2965_3521_2:0:0_2:1:0_eb/1_TA chrM 3521 NA TA 1 TA 1 100 +chrM_3561_4058_2:0:0_2:0:0_14d/1_TG chrM 3560 NA TG 1 TG 1 101 +chrM_3590_4068_6:0:0_6:1:0_b8/1_GG chrM 3589 NA GG 1 GG 1 102 +chrM_3123_3605_5:0:0_4:0:0_14/1_AG chrM 3605 NA AG 1 AG 1 103 +chrM_3158_3635_7:0:0_0:0:0_1c/1_AG chrM 3635 NA AG 1 AG 1 104 +chrM_3253_3640_6:1:0_8:0:0_10a/1_GT chrM 3640 NA GT 1 GT 1 105 +chrM_3663_4221_4:0:0_1:0:0_1a9/1_TA chrM 3662 NA TA 1 TA 1 106 +chrM_3679_4214_3:1:0_4:0:0_54/1_AG chrM 3678 NA AG 1 AG 1 107 +chrM_3688_4175_3:1:0_8:0:0_fd/1_AA chrM 3687 NA AA 1 AA 1 108 +chrM_3216_3691_4:1:0_4:0:0_110/1_TT chrM 3691 NA TT 1 TT 1 109 +chrM_3701_4211_4:0:0_3:0:0_e0/1_TA chrM 3700 NA TA 1 TA 1 110 +chrM_3213_3703_4:1:0_6:0:0_199/1_AG chrM 3703 NA AG 1 AG 1 111 +chrM_3738_4170_4:1:0_3:0:0_3b/1_CA chrM 3737 NA CA 1 CA 1 112 +chrM_3752_4208_9:1:0_5:0:0_11f/1_CA chrM 3751 NA CA 1 CA 1 113 +chrM_3171_3761_4:0:0_6:0:0_38/1_AT chrM 3761 NA AT 1 AT 1 114 +chrM_3213_3785_3:1:0_10:0:0_ed/1_GT chrM 3785 NA GT 1 GT 1 115 +chrM_3793_4298_6:0:0_5:0:0_51/1_TC chrM 3792 NA TC 1 TC 1 116 +chrM_3805_4310_4:1:0_3:0:0_161/1_TA chrM 3804 NA TA 1 TA 1 117 +chrM_3817_4309_3:1:0_3:0:0_86/1_AT chrM 3816 NA AT 1 AT 1 118 +chrM_3863_4305_5:0:0_5:0:0_e9/1_AC chrM 3862 NA AC 1 AC 1 119 +chrM_3399_3957_5:1:0_2:1:0_43/1_AT chrM 3957 NA AT 1 AT 1 120 +chrM_3449_3959_6:0:0_3:1:0_16a/1_GA chrM 3959 NA GA 1 GA 1 121 +chrM_3971_4563_4:0:0_5:0:0_150/1_GG chrM 3970 NA GG 1 GG 1 122 +chrM_3541_3989_4:0:0_2:1:0_90/1_GA chrM 3989 NA GA 1 GA 1 123 +chrM_3597_3996_3:0:0_7:0:0_62/1_GT chrM 3996 NA GT 1 GT 1 124 +chrM_3450_3998_6:0:0_2:0:0_1ba/1_GC chrM 3998 NA GC 1 GC 1 125 +chrM_3637_4061_5:0:0_2:1:0_111/1_AT chrM 4061 NA AT 1 AT 1 126 +chrM_4068_4556_3:0:0_3:0:0_6a/1_CA chrM 4067 NA CA 1 CA 1 127 +chrM_4074_4533_2:0:0_2:0:0_9e/1_TA chrM 4073 NA TA 1 TA 1 128 +chrM_4082_4621_7:0:0_3:0:0_198/1_AA chrM 4081 NA AA 1 AA 1 129 +chrM_4160_4692_5:0:0_4:1:0_1f0/1_AT chrM 4159 NA AT 1 AT 1 130 +chrM_4199_4751_2:0:0_1:1:0_15c/1_GC chrM 4198 NA GC 1 GC 1 131 +chrM_4229_4769_4:0:0_6:0:0_cc/1_GA chrM 4228 NA GA 1 GA 1 132 +chrM_3770_4283_1:0:0_8:0:0_13f/1_AA chrM 4283 NA AA 1 AA 1 133 +chrM_4374_4823_2:0:0_5:0:0_dc/1_TA chrM 4373 NA TA 1 TA 1 134 +chrM_4382_4881_2:0:0_5:0:0_195/1_AT chrM 4381 NA AT 1 AT 1 135 +chrM_4384_4922_5:0:0_4:1:0_1e1/1_TA chrM 4383 NA TA 1 TA 1 136 +chrM_4056_4522_5:0:0_4:0:0_e3/1_AT chrM 4522 NA AT 1 AT 1 137 +chrM_4553_5058_6:0:0_4:0:0_a3/1_AT chrM 4552 NA AT 1 AT 1 138 +chrM_4565_5190_3:0:0_4:0:0_19c/1_AC chrM 4564 NA AC 1 AC 1 139 +chrM_4106_4574_6:0:0_2:0:0_64/1_TC chrM 4574 NA TC 1 TC 1 140 +chrM_4705_5072_2:0:0_2:1:0_b7/1_TT chrM 4704 NA TT 1 TT 1 141 +chrM_4724_5279_5:0:0_8:0:0_2e/1_AA chrM 4723 NA AA 1 AA 1 142 +chrM_4840_5344_4:1:0_4:0:0_3c/1_GT chrM 4839 NA GT 1 GT 1 143 +chrM_4509_4974_5:0:0_3:0:0_ea/1_CT chrM 4974 NA CT 1 CT 1 144 +chrM_4444_5010_5:1:0_5:1:0_119/1_TG chrM 5010 NA TG 1 TG 1 145 +chrM_4491_5013_5:1:0_6:1:0_1b/1_GA chrM 5013 NA GA 1 GA 1 146 +chrM_4490_5017_8:0:0_4:0:0_8d/1_AG chrM 5017 NA AG 1 AG 1 147 +chrM_5058_5566_6:0:0_5:0:0_1c6/1_AG chrM 5057 NA AG 1 AG 1 148 +chrM_5062_5522_7:0:0_5:0:0_e8/1_TC chrM 5061 NA TC 1 TC 1 149 +chrM_4566_5086_5:1:0_6:1:0_19/1_AT chrM 5086 NA AT 1 AT 1 150 +chrM_5119_5706_5:0:0_4:0:0_140/1_TG chrM 5118 NA TG 1 TG 1 151 +chrM_5121_5635_3:0:0_3:0:0_165/1_AG chrM 5120 NA AG 1 AG 1 152 +chrM_5166_5697_4:0:0_4:0:0_ec/1_TC chrM 5165 NA TC 1 TC 1 153 +chrM_4583_5175_2:0:0_4:0:0_10d/1_TT chrM 5175 NA TT 1 TT 1 154 +chrM_5184_5616_6:0:0_4:0:0_1ee/1_AA chrM 5183 NA AA 1 AA 1 155 +chrM_4713_5198_5:0:0_1:0:0_f0/1_TA chrM 5198 NA TA 1 TA 1 156 +chrM_4833_5239_3:1:0_3:0:0_56/1_AG chrM 5239 NA AG 1 AG 1 157 +chrM_5292_5797_3:0:0_3:0:0_130/1_CC chrM 5291 NA CC 1 CC 1 158 +chrM_5315_5855_2:0:0_5:0:0_136/1_TA chrM 5314 NA TA 1 TA 1 159 +chrM_5342_5909_5:0:0_5:0:0_5/1_GC chrM 5341 NA GC 1 GC 1 160 +chrM_4886_5348_5:0:0_3:0:0_d7/1_TC chrM 5348 NA TC 1 TC 1 161 +chrM_5355_5797_0:0:0_5:0:0_1d3/1_TA chrM 5354 NA TA 1 TA 1 162 +chrM_5362_5890_4:0:0_3:0:0_1dd/1_AT chrM 5361 NA AT 1 AT 1 163 +chrM_4974_5463_3:0:0_4:0:0_8/1_CT chrM 5463 NA CT 1 CT 1 164 +chrM_5125_5611_2:0:0_2:0:0_11b/1_AT chrM 5611 NA AT 1 AT 1 165 +chrM_5163_5634_8:0:0_3:0:0_1bf/1_TT chrM 5634 NA TT 1 TT 1 166 +chrM_5185_5662_3:0:0_4:0:0_a4/1_AT chrM 5662 NA AT 1 AT 1 167 +chrM_5683_6115_5:0:0_5:0:0_b5/1_CT chrM 5682 NA CT 1 CT 1 168 +chrM_5306_5717_1:0:0_4:0:0_1ab/1_TT chrM 5717 NA TT 1 TT 1 169 +chrM_5181_5752_4:0:0_4:0:0_128/1_GT chrM 5752 NA GT 1 GT 1 170 +chrM_5828_6286_5:0:0_4:0:0_17f/1_TT chrM 5827 NA TT 1 TT 1 171 +chrM_5310_5842_3:0:0_6:0:0_1d0/1_TA chrM 5842 NA TA 1 TA 1 172 +chrM_5845_6332_2:0:0_3:0:0_db/1_AA chrM 5844 NA AA 1 AA 1 173 +chrM_5856_6361_4:0:0_4:0:0_1d4/1_AA chrM 5855 NA AA 1 AA 1 174 +chrM_5410_5865_9:0:0_5:0:0_57/1_TA chrM 5865 NA TA 1 TA 1 175 +chrM_5411_5901_3:0:0_2:0:0_197/1_TT chrM 5901 NA TT 1 TT 1 176 +chrM_5423_5960_2:0:0_6:0:0_1ca/1_TA chrM 5960 NA TA 1 TA 1 177 +chrM_5461_5962_4:0:0_8:0:0_1a8/1_AT chrM 5962 NA AT 1 AT 1 178 +chrM_5990_6598_5:0:0_2:0:0_16f/1_TC chrM 5989 NA TC 1 TC 1 179 +chrM_6047_6502_2:0:0_6:0:0_1f3/1_CA chrM 6046 NA CA 1 CA 1 180 +chrM_6078_6610_7:0:0_6:0:0_29/1_GG chrM 6077 NA GG 1 GG 1 181 +chrM_5565_6084_3:0:0_3:0:0_194/1_TC chrM 6084 NA TC 1 TC 1 182 +chrM_6087_6574_4:0:0_2:0:0_12e/1_CC chrM 6086 NA CC 1 CC 1 183 +chrM_5710_6256_2:0:0_6:0:0_12/1_CT chrM 6256 NA CT 2 CT 2 184 +chrM_5852_6256_2:0:0_4:0:0_58/1_CT chrM 6256 NA CT 2 CT 2 184 +chrM_5770_6275_1:0:0_3:0:0_a6/1_AG chrM 6275 NA AG 1 AG 1 185 +chrM_5839_6288_5:0:0_5:0:0_1c0/1_GG chrM 6288 NA GG 1 GG 1 186 +chrM_5829_6293_3:0:0_1:0:0_bd/1_GC chrM 6293 NA GC 1 GC 1 187 +chrM_5793_6362_4:0:0_5:0:0_1a5/1_CC chrM 6362 NA CC 1 CC 1 188 +chrM_5942_6425_5:0:0_1:0:0_16d/1_AC chrM 6425 NA AC 1 AC 1 189 +chrM_5989_6468_5:0:0_1:0:0_16c/1_TG chrM 6468 NA TG 1 TG 1 190 +chrM_5994_6492_3:0:0_3:0:0_153/1_TA chrM 6492 NA TA 1 TA 1 191 +chrM_6011_6514_3:0:0_7:0:0_1ef/1_TG chrM 6514 NA TG 1 TG 1 192 +chrM_6134_6610_4:0:0_2:0:0_d6/1_GG chrM 6610 NA GG 1 GG 1 193 +chrM_6073_6613_4:0:0_6:0:0_e/1_CT chrM 6613 NA CT 1 CT 1 194 +chrM_6655_7146_3:0:0_6:0:0_55/1_TT chrM 6654 NA TT 1 TT 1 195 +chrM_6155_6684_4:0:0_6:0:0_156/1_GG chrM 6684 NA GG 1 GG 1 196 +chrM_6698_7198_5:0:1_4:0:0_19e/1_AT chrM 6697 NA AT 1 AT 1 197 +chrM_6222_6744_5:0:0_2:0:0_116/1_CA chrM 6744 NA CA 1 CA 1 198 +chrM_6258_6750_3:0:0_2:0:0_b0/1_TT chrM 6750 NA TT 1 TT 1 199 +chrM_6331_6811_6:0:0_5:0:0_b6/1_GT chrM 6811 NA GT 1 GT 1 200 +chrM_6370_6852_5:0:0_4:0:0_5c/1_AT chrM 6852 NA AT 1 AT 1 201 +chrM_6507_6920_3:0:0_5:0:1_1d2/1_AA chrM 6920 NA AA 1 AA 1 202 +chrM_6479_6981_2:0:0_5:0:1_18f/1_TG chrM 6981 NA TG 1 TG 1 203 +chrM_6530_7006_2:0:0_2:0:1_101/1_TA chrM 7006 NA TA 1 TA 1 204 +chrM_7048_7515_9:0:0_4:0:0_7e/1_AT chrM 7047 NA AT 1 AT 1 205 +chrM_7081_7575_3:0:0_0:0:0_10/1_CC chrM 7080 NA CC 1 CC 1 206 +chrM_7119_7629_5:0:0_7:0:0_1cf/1_GC chrM 7118 NA GC 1 GC 1 207 +chrM_7270_7705_10:0:0_4:0:0_b4/1_GC chrM 7269 NA GC 1 GC 1 208 +chrM_7272_7803_5:0:0_4:0:0_fe/1_TA chrM 7271 NA TA 1 TA 1 209 +chrM_6787_7349_2:0:1_1:0:0_141/1_TG chrM 7349 NA TG 1 TG 1 210 +chrM_7365_7794_3:0:0_5:0:0_1c7/1_GA chrM 7364 NA GA 1 GA 1 211 +chrM_7413_7973_6:0:0_5:0:0_158/1_TC chrM 7412 NA TC 1 TC 1 212 +chrM_6786_7414_0:0:1_6:0:0_1b6/1_AT chrM 7414 NA AT 1 AT 1 213 +chrM_7472_7886_4:0:0_4:0:0_8b/1_GG chrM 7471 NA GG 1 GG 1 214 +chrM_7505_7965_4:0:0_2:0:0_3a/1_CT chrM 7504 NA CT 1 CT 1 215 +chrM_7095_7547_2:0:0_3:0:0_149/1_GG chrM 7547 NA GG 1 GG 1 216 +chrM_7564_8021_3:0:0_2:0:0_1/1_AA chrM 7563 NA AA 1 AA 1 217 +chrM_7067_7564_4:0:0_5:0:0_107/1_AA chrM 7564 NA AA 1 AA 1 218 +chrM_7016_7567_4:0:0_5:0:0_27/1_TT chrM 7567 NA TT 1 TT 1 219 +chrM_7200_7687_2:0:0_6:0:0_c/1_TA chrM 7687 NA TA 1 TA 1 220 +chrM_7267_7793_1:0:0_4:0:0_41/1_AC chrM 7793 NA AC 1 AC 1 221 +chrM_7797_8230_4:0:0_5:0:0_f6/1_GC chrM 7796 NA GC 1 GC 1 222 +chrM_7821_8285_3:0:0_3:0:0_1f1/1_TT chrM 7820 NA TT 1 TT 1 223 +chrM_7859_8392_2:0:0_3:0:0_114/1_TT chrM 7858 NA TT 1 TT 1 224 +chrM_7881_8395_3:0:0_2:0:0_173/1_CT chrM 7880 NA CT 1 CT 1 225 +chrM_7339_7893_4:0:0_3:0:0_19b/1_TA chrM 7893 NA TA 1 TA 1 226 +chrM_7424_7912_5:0:0_4:0:0_ba/1_GG chrM 7912 NA GG 1 GG 1 227 +chrM_7467_7949_6:0:0_9:0:0_7f/1_AT chrM 7949 NA AT 1 AT 1 228 +chrM_7956_8469_8:0:0_5:1:0_14a/1_AA chrM 7955 NA AA 1 AA 1 229 +chrM_7963_8508_6:0:0_5:1:0_73/1_AC chrM 7962 NA AC 1 AC 1 230 +chrM_8018_8547_2:0:0_2:1:0_39/1_TA chrM 8017 NA TA 1 TA 1 231 +chrM_8039_8593_3:0:0_3:1:0_71/1_AG chrM 8038 NA AG 1 AG 1 232 +chrM_7594_8044_3:0:0_5:0:0_e1/1_TT chrM 8044 NA TT 1 TT 1 233 +chrM_7570_8062_4:0:0_1:0:0_1ed/1_TC chrM 8062 NA TC 1 TC 1 234 +chrM_7576_8099_6:0:0_4:0:0_1bc/1_GT chrM 8099 NA GT 1 GT 1 235 +chrM_8114_8625_4:0:0_3:1:0_17a/1_AT chrM 8113 NA AT 1 AT 1 236 +chrM_7655_8131_5:0:0_7:0:0_19f/1_CG chrM 8131 NA CG 1 CG 1 237 +chrM_8139_8619_2:0:0_5:1:0_166/1_CC chrM 8138 NA CC 1 CC 1 238 +chrM_7825_8191_3:0:0_2:0:0_1ec/1_CT chrM 8191 NA CT 1 CT 1 239 +chrM_7641_8206_4:0:0_5:0:0_13a/1_GT chrM 8206 NA GT 1 GT 1 240 +chrM_7613_8208_2:0:0_5:0:0_52/1_AT chrM 8208 NA AT 1 AT 1 241 +chrM_7752_8261_10:0:0_5:0:0_142/1_GA chrM 8261 NA GA 1 GA 1 242 +chrM_7739_8294_7:0:0_2:0:0_144/1_CC chrM 8294 NA CC 1 CC 1 243 +chrM_8303_8710_5:1:0_6:0:0_87/1_CA chrM 8302 NA CA 1 CA 1 244 +chrM_8309_8842_3:1:0_5:0:0_115/1_GA chrM 8308 NA GA 1 GA 1 245 +chrM_7965_8408_4:0:0_0:0:0_143/1_CC chrM 8408 NA CC 1 CC 1 246 +chrM_8031_8436_7:0:0_5:0:0_c0/1_TC chrM 8436 NA TC 1 TC 1 247 +chrM_8438_8973_3:1:0_2:0:0_1d5/1_AA chrM 8437 NA AA 1 AA 1 248 +chrM_8511_8920_4:0:0_3:0:0_ca/1_AT chrM 8510 NA AT 1 AT 1 249 +chrM_8532_8954_3:0:0_4:0:0_cb/1_TT chrM 8531 NA TT 1 TT 1 250 +chrM_8565_8973_6:0:0_1:0:0_66/1_GA chrM 8564 NA GA 1 GA 1 251 +chrM_8575_9003_7:0:0_2:0:0_c7/1_TA chrM 8574 NA TA 1 TA 1 252 +chrM_8125_8594_2:0:0_2:1:0_187/1_TA chrM 8594 NA TA 1 TA 1 253 +chrM_8144_8650_4:0:0_3:1:0_14b/1_GG chrM 8650 NA GG 1 GG 1 254 +chrM_8117_8669_3:0:0_5:0:0_16b/1_AG chrM 8669 NA AG 1 AG 1 255 +chrM_8296_8802_4:1:0_9:0:0_1db/1_GC chrM 8802 NA GC 1 GC 1 256 +chrM_8807_9277_3:0:0_4:0:1_108/1_TG chrM 8806 NA TG 1 TG 1 257 +chrM_8825_9275_3:0:0_5:0:1_1da/1_AA chrM 8824 NA AA 1 AA 1 258 +chrM_8868_9400_4:0:0_3:0:1_1e8/1_AC chrM 8867 NA AC 1 AC 1 259 +chrM_8885_9340_4:0:0_5:0:1_c8/1_TT chrM 8884 NA TT 1 TT 1 260 +chrM_8938_9405_8:0:0_2:0:1_1c3/1_CG chrM 8937 NA CG 1 CG 1 261 +chrM_8396_8947_8:1:0_2:0:0_155/1_AC chrM 8947 NA AC 1 AC 1 262 +chrM_8461_8951_7:0:0_1:0:0_1cb/1_TG chrM 8951 NA TG 1 TG 1 263 +chrM_8987_9581_3:0:0_3:0:0_1cd/1_TA chrM 8986 NA TA 1 TA 1 264 +chrM_9007_9456_2:0:0_4:0:0_23/1_AA chrM 9006 NA AA 1 AA 1 265 +chrM_8520_9009_2:0:0_7:0:0_193/1_TA chrM 9009 NA TA 1 TA 1 266 +chrM_8547_9016_5:0:0_4:0:0_32/1_TG chrM 9016 NA TG 1 TG 1 267 +chrM_9075_9551_5:0:1_3:0:0_4b/1_CC chrM 9074 NA CC 1 CC 1 268 +chrM_9102_9561_3:0:1_3:0:0_2b/1_TA chrM 9101 NA TA 1 TA 1 269 +chrM_9113_9686_4:0:1_1:0:0_78/1_CA chrM 9112 NA CA 1 CA 1 270 +chrM_8621_9133_8:0:0_5:0:0_145/1_GC chrM 9133 NA GC 1 GC 1 271 +chrM_9166_9651_6:0:1_4:0:0_15b/1_TT chrM 9165 NA TT 1 TT 1 272 +chrM_9168_9638_4:0:1_5:0:0_175/1_TC chrM 9167 NA TC 1 TC 1 273 +chrM_9172_9645_0:0:1_5:0:0_e7/1_CA chrM 9171 NA CA 1 CA 1 274 +chrM_9233_9739_3:0:1_5:0:0_91/1_CT chrM 9232 NA CT 1 CT 1 275 +chrM_9249_9749_6:0:1_2:0:0_12c/1_TA chrM 9247 NA TA 1 TA 1 276 +chrM_9338_9899_1:0:0_2:0:0_a8/1_TT chrM 9337 NA TT 1 TT 1 277 +chrM_8914_9385_2:0:0_3:0:1_135/1_GC chrM 9385 NA GC 1 GC 1 278 +chrM_8829_9414_6:0:0_2:0:1_1bd/1_GT chrM 9414 NA GT 1 GT 1 279 +chrM_8916_9457_5:0:0_5:0:0_167/1_AC chrM 9457 NA AC 1 AC 1 280 +chrM_8970_9471_5:0:0_6:0:0_1ac/1_CT chrM 9471 NA CT 1 CT 1 281 +chrM_8928_9477_3:0:0_6:0:0_48/1_GG chrM 9477 NA GG 2 GG 2 282 +chrM_8929_9477_7:0:0_5:0:0_f9/1_GG chrM 9477 NA GG 2 GG 2 282 +chrM_9536_10046_5:0:0_1:0:0_d5/1_TA chrM 9535 NA TA 1 TA 1 283 +chrM_9579_10013_2:0:0_3:0:0_9/1_AC chrM 9578 NA AC 1 AC 1 284 +chrM_9672_10188_5:0:0_9:0:0_11/1_CG chrM 9671 NA CG 1 CG 1 285 +chrM_9689_10220_2:0:0_2:0:0_151/1_AA chrM 9688 NA AA 1 AA 1 286 +chrM_9159_9700_2:0:1_4:0:0_2/1_GA chrM 9700 NA GA 1 GA 1 287 +chrM_9103_9709_4:0:1_5:0:0_24/1_AT chrM 9709 NA AT 1 AT 1 288 +chrM_9768_10236_5:0:0_4:0:0_1d9/1_GT chrM 9767 NA GT 1 GT 1 289 +chrM_9271_9836_3:0:0_2:0:0_77/1_CT chrM 9836 NA CT 1 CT 1 290 +chrM_9847_10370_4:0:0_5:0:0_bc/1_CA chrM 9846 NA CA 1 CA 1 291 +chrM_9856_10428_4:0:0_4:1:0_120/1_GC chrM 9855 NA GC 1 GC 1 292 +chrM_9866_10380_3:0:0_5:0:0_8f/1_CC chrM 9865 NA CC 1 CC 1 293 +chrM_9895_10412_1:0:0_7:1:0_7a/1_AT chrM 9894 NA AT 1 AT 1 294 +chrM_9904_10380_7:0:0_3:0:0_1e0/1_TG chrM 9903 NA TG 1 TG 1 295 +chrM_9917_10459_5:0:0_6:0:0_9c/1_CG chrM 9916 NA CG 1 CG 1 296 +chrM_9489_9919_2:0:0_2:0:0_b9/1_GG chrM 9919 NA GG 1 GG 1 297 +chrM_9452_9927_4:0:0_5:0:0_190/1_GA chrM 9927 NA GA 1 GA 1 298 +chrM_9932_10484_3:0:0_5:0:0_148/1_CT chrM 9931 NA CT 1 CT 1 299 +chrM_9983_10519_5:0:0_2:0:0_dd/1_AC chrM 9982 NA AC 1 AC 1 300 +chrM_9402_10032_4:0:0_5:0:0_a2/1_GG chrM 10032 NA GG 1 GG 1 301 +chrM_10074_10591_1:0:0_0:0:0_14c/1_TA chrM 10073 NA TA 1 TA 1 302 +chrM_10103_10594_8:0:0_6:1:0_1e9/1_AT chrM 10102 NA AT 1 AT 1 303 +chrM_10141_10683_4:0:0_4:0:0_6f/1_TA chrM 10140 NA TA 1 TA 1 304 +chrM_10226_10692_5:1:0_4:0:0_103/1_CT chrM 10225 NA CT 1 CT 1 305 +chrM_9669_10227_5:0:0_4:0:0_4f/1_TG chrM 10227 NA TG 1 TG 1 306 +chrM_9735_10257_0:0:0_7:0:0_b2/1_TA chrM 10257 NA TA 1 TA 1 307 +chrM_10326_10838_5:0:0_3:0:0_f2/1_TT chrM 10325 NA TT 1 TT 1 308 +chrM_9825_10334_2:0:0_6:0:0_123/1_AC chrM 10334 NA AC 1 AC 1 309 +chrM_10009_10476_0:0:0_5:0:0_176/1_CA chrM 10476 NA CA 1 CA 1 310 +chrM_9995_10574_6:0:0_4:1:0_7c/1_GG chrM 10574 NA GG 1 GG 1 311 +chrM_10582_11133_4:0:0_3:0:0_21/1_GC chrM 10581 NA GC 1 GC 1 312 +chrM_10590_11016_2:0:0_6:0:0_81/1_AA chrM 10589 NA AA 1 AA 1 313 +chrM_10607_11150_5:0:0_2:0:0_8a/1_CC chrM 10606 NA CC 1 CC 1 314 +chrM_10041_10610_6:0:0_4:0:0_ef/1_GG chrM 10610 NA GG 1 GG 1 315 +chrM_10679_11090_2:0:0_6:0:0_1b4/1_GA chrM 10678 NA GA 1 GA 1 316 +chrM_10690_11211_4:0:0_3:0:0_184/1_AC chrM 10689 NA AC 1 AC 1 317 +chrM_10692_11227_2:0:0_4:0:0_49/1_TG chrM 10691 NA TG 1 TG 1 318 +chrM_10325_10814_5:1:0_6:0:0_129/1_GC chrM 10814 NA GC 1 GC 1 319 +chrM_10834_11259_4:0:0_3:0:0_47/1_GA chrM 10833 NA GA 1 GA 1 320 +chrM_10846_11368_1:0:0_6:0:0_fc/1_TG chrM 10845 NA TG 1 TG 1 321 +chrM_10850_11370_7:0:0_4:0:0_15e/1_TT chrM 10849 NA TT 1 TT 1 322 +chrM_10491_10996_4:0:0_6:0:0_1a6/1_AG chrM 10996 NA AG 1 AG 1 323 +chrM_11001_11542_4:0:0_9:0:0_15/1_GT chrM 11000 NA GT 1 GT 1 324 +chrM_11067_11640_2:0:0_7:0:0_d2/1_TT chrM 11066 NA TT 1 TT 1 325 +chrM_10600_11115_3:0:0_3:0:0_10f/1_TT chrM 11115 NA TT 1 TT 1 326 +chrM_11179_11645_6:0:1_1:0:0_181/1_TC chrM 11178 NA TC 1 TC 1 327 +chrM_10716_11197_6:0:0_5:0:0_67/1_CA chrM 11197 NA CA 1 CA 1 328 +chrM_10649_11218_5:0:0_8:0:0_50/1_GT chrM 11218 NA GT 1 GT 1 329 +chrM_10624_11233_2:0:0_3:0:0_2f/1_CT chrM 11233 NA CT 1 CT 1 330 +chrM_11239_11740_6:0:0_1:0:0_3e/1_CT chrM 11238 NA CT 1 CT 1 331 +chrM_10825_11310_5:0:0_6:0:0_1ce/1_TA chrM 11310 NA TA 1 TA 1 332 +chrM_11423_11976_7:0:0_3:0:0_189/1_GA chrM 11422 NA GA 1 GA 1 333 +chrM_11522_12121_5:0:0_2:0:0_13d/1_GG chrM 11521 NA GG 1 GG 1 334 +chrM_11080_11521_6:0:0_7:0:0_ff/1_GT chrM 11521 NA GT 1 GT 1 335 +chrM_11076_11558_1:0:0_3:0:0_92/1_AT chrM 11558 NA AT 1 AT 1 336 +chrM_11197_11674_1:0:0_3:0:0_c1/1_TT chrM 11674 NA TT 1 TT 1 337 +chrM_11227_11680_5:0:1_3:0:0_76/1_GT chrM 11680 NA GT 1 GT 1 338 +chrM_11759_12257_5:0:0_8:0:0_9b/1_TT chrM 11758 NA TT 1 TT 1 339 +chrM_11236_11767_1:0:0_2:0:0_e2/1_TT chrM 11767 NA TT 1 TT 1 340 +chrM_11780_12326_4:0:0_7:0:0_10e/1_GG chrM 11779 NA GG 1 GG 1 341 +chrM_11789_12313_6:0:0_5:0:0_1f2/1_CG chrM 11788 NA CG 1 CG 1 342 +chrM_11815_12366_6:0:0_1:0:0_1be/1_AA chrM 11814 NA AA 1 AA 1 343 +chrM_11829_12278_5:0:0_8:0:0_11a/1_CA chrM 11828 NA CA 1 CA 1 344 +chrM_11834_12212_5:0:0_4:0:0_131/1_CC chrM 11833 NA CC 1 CC 1 345 +chrM_11305_11834_4:0:0_5:0:0_5e/1_TA chrM 11834 NA TA 1 TA 1 346 +chrM_11315_11842_2:0:0_8:0:0_126/1_TA chrM 11842 NA TA 1 TA 1 347 +chrM_11890_12397_8:0:0_2:0:0_30/1_CC chrM 11889 NA CC 1 CC 1 348 +chrM_11898_12351_4:0:0_1:0:0_a9/1_AC chrM 11897 NA AC 1 AC 1 349 +chrM_11902_12409_3:0:0_4:0:0_ad/1_AA chrM 11901 NA AA 1 AA 1 350 +chrM_11932_12361_4:0:0_6:0:0_44/1_TC chrM 11931 NA TC 1 TC 1 351 +chrM_11938_12488_7:0:0_6:0:0_ee/1_GG chrM 11937 NA GG 1 GG 1 352 +chrM_11949_12385_7:0:0_5:0:0_118/1_TT chrM 11948 NA TT 1 TT 1 353 +chrM_11551_11956_7:0:0_5:0:0_d0/1_CT chrM 11956 NA CT 1 CT 1 354 +chrM_11468_11961_6:0:0_4:0:0_c9/1_TA chrM 11961 NA TA 1 TA 1 355 +chrM_11495_11994_2:0:0_5:0:0_100/1_AG chrM 11994 NA AG 1 AG 1 356 +chrM_11413_12020_5:0:0_2:0:0_4c/1_AG chrM 12020 NA AG 1 AG 1 357 +chrM_11598_12060_3:0:0_5:0:0_89/1_AT chrM 12060 NA AT 1 AT 1 358 +chrM_11616_12104_4:0:0_7:0:0_6/1_TG chrM 12104 NA TG 1 TG 1 359 +chrM_12150_12654_2:0:0_3:0:0_12b/1_CA chrM 12149 NA CA 1 CA 1 360 +chrM_12160_12669_8:0:0_2:0:0_1ae/1_CC chrM 12159 NA CC 1 CC 1 361 +chrM_11701_12198_5:0:0_2:0:0_96/1_AA chrM 12198 NA AA 1 AA 1 362 +chrM_11745_12199_2:0:0_3:0:0_7/1_CA chrM 12199 NA CA 1 CA 1 363 +chrM_11818_12273_4:0:0_4:0:0_196/1_TA chrM 12273 NA TA 1 TA 1 364 +chrM_11791_12301_5:0:0_6:0:0_17d/1_GA chrM 12301 NA GA 1 GA 1 365 +chrM_12330_12776_4:0:0_4:0:0_63/1_AT chrM 12329 NA AT 1 AT 1 366 +chrM_12336_12829_5:0:0_5:0:0_f/1_AC chrM 12335 NA AC 1 AC 1 367 +chrM_11746_12340_4:0:0_2:0:0_6c/1_GG chrM 12340 NA GG 1 GG 1 368 +chrM_12356_12830_4:0:0_5:0:0_104/1_AA chrM 12355 NA AA 1 AA 1 369 +chrM_12371_12852_6:0:0_5:0:0_28/1_CT chrM 12370 NA CT 1 CT 1 370 +chrM_12405_12863_2:0:0_3:0:0_154/1_AC chrM 12404 NA AC 1 AC 1 371 +chrM_12413_12965_5:0:0_6:0:0_174/1_AT chrM 12412 NA AT 1 AT 1 372 +chrM_12006_12431_4:0:0_4:0:0_15a/1_AT chrM 12431 NA AT 1 AT 1 373 +chrM_11969_12486_5:0:0_4:0:0_2a/1_AC chrM 12486 NA AC 1 AC 1 374 +chrM_12495_12966_1:0:0_3:0:0_18a/1_TT chrM 12494 NA TT 1 TT 1 375 +chrM_12089_12585_1:0:0_7:0:0_1e4/1_TT chrM 12585 NA TT 1 TT 1 376 +chrM_12621_13150_2:0:0_4:0:0_9a/1_CA chrM 12620 NA CA 1 CA 1 377 +chrM_12209_12659_3:0:0_3:0:0_1a7/1_GT chrM 12659 NA GT 1 GT 1 378 +chrM_12702_13164_4:0:0_3:0:0_88/1_AC chrM 12701 NA AC 1 AC 1 379 +chrM_12712_13263_5:0:0_6:0:0_1a3/1_AT chrM 12711 NA AT 1 AT 1 380 +chrM_12284_12729_4:0:0_4:0:0_61/1_TA chrM 12729 NA TA 1 TA 1 381 +chrM_12259_12743_4:0:0_2:0:0_1a2/1_TT chrM 12743 NA TT 1 TT 1 382 +chrM_12260_12784_5:0:0_2:0:0_6d/1_TG chrM 12784 NA TG 1 TG 1 383 +chrM_12316_12791_2:0:0_3:0:0_1d6/1_TT chrM 12791 NA TT 1 TT 1 384 +chrM_12262_12819_3:0:0_5:0:0_ce/1_AT chrM 12819 NA AT 1 AT 1 385 +chrM_12329_12832_4:0:0_4:0:0_1c5/1_TC chrM 12832 NA TC 1 TC 1 386 +chrM_12323_12847_3:0:0_3:0:0_163/1_GA chrM 12847 NA GA 1 GA 1 387 +chrM_12350_12855_4:0:0_3:0:0_3f/1_AG chrM 12855 NA AG 1 AG 1 388 +chrM_12495_12939_2:0:0_4:0:0_117/1_TC chrM 12939 NA TC 1 TC 1 389 +chrM_12949_13437_4:0:0_3:0:0_1a0/1_GG chrM 12948 NA GG 1 GG 1 390 +chrM_12449_12988_5:0:0_4:0:0_79/1_AA chrM 12988 NA AA 1 AA 1 391 +chrM_13032_13560_5:0:0_4:0:0_b1/1_GC chrM 13031 NA GC 1 GC 1 392 +chrM_12548_13045_4:0:0_3:0:0_be/1_TA chrM 13045 NA TA 1 TA 1 393 +chrM_12698_13119_5:0:0_3:0:0_a/1_GC chrM 13119 NA GC 1 GC 1 394 +chrM_12657_13157_3:0:0_2:0:0_a0/1_TG chrM 13157 NA TG 1 TG 1 395 +chrM_12629_13160_5:0:0_5:0:0_125/1_GG chrM 13160 NA GG 1 GG 1 396 +chrM_12642_13191_2:0:0_5:0:0_1dc/1_AT chrM 13191 NA AT 1 AT 1 397 +chrM_13236_13695_2:0:0_9:0:0_1c8/1_TT chrM 13235 NA TT 1 TT 1 398 +chrM_13271_13696_3:0:0_5:0:0_188/1_AC chrM 13270 NA AC 1 AC 1 399 +chrM_12763_13284_4:0:0_1:0:0_d1/1_AG chrM 13284 NA AG 1 AG 1 400 +chrM_13300_13825_6:0:0_3:0:0_ac/1_AA chrM 13299 NA AA 2 AA 2 401 +chrM_13300_13872_3:0:0_5:1:0_c2/1_AA chrM 13299 NA AA 2 AA 2 401 +chrM_12845_13364_2:0:0_5:0:0_109/1_AG chrM 13364 NA AG 1 AG 1 402 +chrM_13502_13977_3:0:0_5:0:0_ae/1_TT chrM 13501 NA TT 1 TT 1 403 +chrM_13627_14137_3:0:0_3:0:0_162/1_AC chrM 13626 NA AC 1 AC 1 404 +chrM_13055_13651_7:0:0_3:0:0_8c/1_TA chrM 13651 NA TA 1 TA 1 405 +chrM_13676_14200_0:1:0_4:0:0_2d/1_CA chrM 13675 NA CA 1 CA 1 406 +chrM_13138_13709_5:0:0_7:0:0_72/1_TT chrM 13709 NA TT 1 TT 1 407 +chrM_13242_13718_6:0:0_4:0:0_11d/1_AT chrM 13718 NA AT 1 AT 1 408 +chrM_13287_13740_4:0:0_5:1:0_1b8/1_GT chrM 13740 NA GT 1 GT 1 409 +chrM_13767_14345_8:0:0_5:0:0_f8/1_GC chrM 13766 NA GC 1 GC 1 410 +chrM_13776_14348_5:0:0_6:0:0_12f/1_GG chrM 13775 NA GG 1 GG 1 411 +chrM_13243_13775_6:0:0_5:1:0_1e/1_TG chrM 13775 NA TG 1 TG 1 412 +chrM_13792_14355_3:0:0_5:0:0_d/1_CA chrM 13791 NA CA 1 CA 1 413 +chrM_13799_14307_7:0:0_6:0:0_82/1_TT chrM 13798 NA TT 1 TT 1 414 +chrM_13306_13843_4:0:0_6:0:0_b/1_AT chrM 13843 NA AT 1 AT 1 415 +chrM_13398_13896_7:0:0_5:1:0_1e5/1_GC chrM 13896 NA GC 1 GC 1 416 +chrM_13409_13926_6:0:0_2:0:0_13c/1_AG chrM 13926 NA AG 1 AG 1 417 +chrM_13997_14575_4:0:0_2:0:0_97/1_AC chrM 13996 NA AC 1 AC 1 418 +chrM_13416_13997_2:0:0_2:0:0_d9/1_TA chrM 13997 NA TA 1 TA 1 419 +chrM_13540_14114_4:1:0_4:0:0_1b2/1_TC chrM 14114 NA TC 1 TC 1 420 +chrM_14120_14578_4:0:0_4:0:0_26/1_CG chrM 14119 NA CG 1 CG 1 421 +chrM_14134_14579_3:0:0_5:0:0_da/1_CA chrM 14133 NA CA 1 CA 1 422 +chrM_13711_14180_4:0:0_8:0:0_13e/1_GG chrM 14180 NA GG 1 GG 1 423 +chrM_14205_14774_5:0:0_5:0:0_105/1_CT chrM 14204 NA CT 1 CT 1 424 +chrM_13714_14222_2:0:0_1:0:0_182/1_GA chrM 14222 NA GA 1 GA 1 425 +chrM_13745_14261_7:0:0_4:0:0_b3/1_TA chrM 14261 NA TA 1 TA 1 426 +chrM_13794_14264_5:0:0_5:0:0_102/1_TA chrM 14264 NA TA 1 TA 1 427 +chrM_13867_14400_3:0:0_1:0:0_14f/1_AA chrM 14400 NA AA 1 AA 1 428 +chrM_13938_14449_2:0:0_6:0:0_1eb/1_GT chrM 14449 NA GT 1 GT 1 429 +chrM_13943_14452_3:0:0_4:0:0_11c/1_GT chrM 14452 NA GT 1 GT 1 430 +chrM_14006_14481_5:0:0_3:0:0_df/1_TA chrM 14481 NA TA 1 TA 1 431 +chrM_14135_14581_5:0:0_5:0:0_1de/1_AT chrM 14581 NA AT 1 AT 1 432 +chrM_14612_15169_2:0:0_5:0:0_e4/1_TT chrM 14611 NA TT 1 TT 1 433 +chrM_14668_15109_2:0:0_3:0:0_3d/1_TG chrM 14667 NA TG 1 TG 1 434 +chrM_14696_15126_2:0:0_2:0:0_a7/1_AA chrM 14695 NA AA 1 AA 1 435 +chrM_14727_15156_4:0:0_4:0:0_178/1_GA chrM 14726 NA GA 1 GA 1 436 +chrM_14245_14783_3:0:0_1:0:0_9d/1_GG chrM 14783 NA GG 1 GG 1 437 +chrM_14314_14802_5:0:0_3:0:0_191/1_GG chrM 14802 NA GG 1 GG 1 438 +chrM_14817_15261_5:0:0_4:0:0_65/1_AT chrM 14816 NA AT 1 AT 1 439 +chrM_14863_15320_2:0:0_4:0:0_d8/1_AA chrM 14862 NA AA 1 AA 1 440 +chrM_14869_15387_2:0:0_7:0:0_1f/1_AC chrM 14868 NA AC 1 AC 1 441 +chrM_14888_15340_3:0:0_3:0:0_112/1_GG chrM 14887 NA GG 1 GG 1 442 +chrM_14368_14894_6:0:0_5:0:0_40/1_TC chrM 14894 NA TC 1 TC 1 443 +chrM_14410_14941_6:0:0_3:0:0_e5/1_TA chrM 14941 NA TA 1 TA 1 444 +chrM_14962_15543_5:0:0_8:0:0_46/1_CA chrM 14961 NA CA 1 CA 1 445 +chrM_14456_14987_3:0:0_5:0:0_ab/1_GT chrM 14987 NA GT 1 GT 1 446 +chrM_15070_15569_5:0:0_5:0:0_cf/1_CA chrM 15069 NA CA 1 CA 1 447 +chrM_15140_15686_2:0:0_5:0:0_185/1_CC chrM 15139 NA CC 1 CC 1 448 +chrM_15142_15661_7:0:0_3:0:0_11e/1_CT chrM 15141 NA CT 1 CT 1 449 +chrM_15192_15694_5:0:0_2:0:0_f7/1_CA chrM 15191 NA CA 1 CA 1 450 +chrM_15211_15685_5:0:0_2:0:0_1d7/1_AC chrM 15210 NA AC 1 AC 1 451 +chrM_15225_15786_3:0:0_6:0:0_17e/1_TT chrM 15224 NA TT 1 TT 1 452 +chrM_15258_15810_4:0:0_6:0:0_5f/1_TT chrM 15257 NA TT 1 TT 1 453 +chrM_14817_15317_5:0:0_2:0:0_59/1_GC chrM 15317 NA GC 1 GC 1 454 +chrM_15324_15836_4:0:0_3:0:0_94/1_AC chrM 15323 NA AC 1 AC 1 455 +chrM_15365_15880_4:1:0_3:0:0_80/1_CA chrM 15364 NA CA 1 CA 1 456 +chrM_15408_15863_2:1:0_3:0:0_1e6/1_AG chrM 15407 NA AG 1 AG 1 457 +chrM_15439_15924_0:0:0_4:0:0_172/1_TG chrM 15438 NA TG 1 TG 1 458 +chrM_14931_15457_2:0:0_4:1:0_1a1/1_AA chrM 15457 NA AA 1 AA 1 459 +chrM_15547_16054_5:0:0_1:0:0_af/1_GA chrM 15546 NA GA 1 GA 1 460 +chrM_15134_15560_3:0:0_4:0:0_1c4/1_GG chrM 15560 NA GG 1 GG 1 461 +chrM_15575_16118_5:0:0_4:0:0_1aa/1_GC chrM 15574 NA GC 1 GC 1 462 +chrM_15052_15629_7:0:0_9:0:0_15f/1_GA chrM 15629 NA GA 1 GA 1 463 +chrM_15698_16224_5:0:0_6:0:0_138/1_AC chrM 15697 NA AC 1 AC 1 464 +chrM_15247_15721_2:1:0_6:0:0_17/1_TC chrM 15721 NA TC 1 TC 1 465 +chrM_15218_15763_2:0:0_2:0:0_171/1_AG chrM 15763 NA AG 1 AG 1 466 +chrM_15860_16419_5:0:0_8:0:0_53/1_TA chrM 15859 NA TA 1 TA 1 467 +chrM_15864_16329_3:0:0_3:0:0_18b/1_CG chrM 15863 NA CG 1 CG 1 468 +chrM_15507_15962_3:0:0_0:0:0_5b/1_GT chrM 15962 NA GT 1 GT 1 469 +chrM_15430_15985_5:0:0_3:0:0_2c/1_AT chrM 15985 NA AT 1 AT 1 470 +chrM_15706_16238_6:0:0_7:0:0_7b/1_AG chrM 16238 NA AG 1 AG 1 471 +chrM_16252_16701_2:1:0_5:0:0_147/1_TA chrM 16251 NA TA 1 TA 1 472 +chrM_15753_16280_5:0:0_6:0:0_180/1_TG chrM 16280 NA TG 1 TG 1 473 +chrM_15777_16347_2:0:0_6:0:0_16/1_AT chrM 16347 NA AT 1 AT 1 474 +chrM_16366_16911_6:1:0_6:1:1_168/1_CC chrM 16365 NA CC 1 CC 1 475 +chrM_16370_16833_6:1:0_6:0:1_1bb/1_AC chrM 16369 NA AC 1 AC 1 476 +chrM_16402_16876_2:0:0_1:0:1_98/1_CA chrM 16401 NA CA 1 CA 1 477 +chrM_16426_16953_3:0:0_2:2:0_60/1_AA chrM 16425 NA AA 1 AA 1 478 +chrM_15986_16442_3:0:0_7:0:0_15d/1_TG chrM 16442 NA TG 1 TG 1 479 +chrM_16030_16460_6:0:0_4:1:0_20/1_AT chrM 16460 NA AT 1 AT 1 480 +chrM_15909_16463_3:0:0_7:1:0_106/1_GT chrM 16463 NA GT 1 GT 1 481 +chrM_15984_16500_4:0:0_1:0:0_6b/1_GT chrM 16500 NA GT 1 GT 1 482 +chrM_16070_16513_3:0:0_7:1:0_1e7/1_TA chrM 16513 NA TA 1 TA 1 483 +chrM_16064_16572_4:0:0_6:1:0_f1/1_TG chrM 16572 NA TG 1 TG 1 484 +chrM_16245_16740_4:1:0_2:0:1_127/1_TA chrM 16742 NA TA 1 TA 1 485 +chrM_16315_16834_5:0:0_3:0:1_18/1_TT chrM 16834 NA TT 1 TT 1 486 +chrM_16229_16843_2:1:0_2:0:1_c4/1_GC chrM 16843 NA GC 1 GC 1 487 +chrM_16420_16905_5:0:0_2:0:1_159/1_GT chrM 16905 NA GT 2 GT 2 488 +chrM_16420_16905_5:0:0_2:0:1_159/1_GT chrM 16905 NA GT 2 GT 2 488 +chrM_16419_16923_3:0:0_3:0:1_183/1_CG chrM 16923 NA CG 1 CG 1 489 +chrM_16585_16993_5:0:1_5:1:0_1b9/1_GA chrM 16993 NA GA 1 GA 1 490 diff -r 000000000000 -r a6477bafd522 test-data/group_out4.bam Binary file test-data/group_out4.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_out4.tab diff -r 000000000000 -r a6477bafd522 test-data/group_out5.bam Binary file test-data/group_out5.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/group_out6.bam Binary file test-data/group_out6.bam has changed diff -r 000000000000 -r a6477bafd522 test-data/out_R1.fastq.gz Binary file test-data/out_R1.fastq.gz has changed diff -r 000000000000 -r a6477bafd522 test-data/out_R2.fastq.gz Binary file test-data/out_R2.fastq.gz has changed diff -r 000000000000 -r a6477bafd522 test-data/out_SE.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/out_SE.fastq Wed Jan 10 19:09:42 2018 -0500 @@ -0,0 +1,288 @@ +@HISEQ:105:C2UE1ACXX:3:1101:11160:2245_TTA 1:N:0:CAGATC +AAAAGTAGTTAATATATTAGATTTGTTTGATAGTGGTAGTATATATTTTTTATTTTAGTATTTAGGAGGTAGAGGTAGATGAATTTTTGAGTTTAAAG ++ +BBBFFFFBFFFFFIIIIIIBIIIIBIFFBFIIBBFFFFFFFFIIIIIIIIIIIIIIIBFIIIIIB7BBBFBBFFFF77 +# stdin : <_io.TextIOWrapper name='input_read1.gz' encoding='ascii'> +# stdlog : <_io.TextIOWrapper name='/tmp/tmpAvWtMd/files/000/dataset_20.dat' mode='a' encoding='UTF-8'> +# stdout : <_io.TextIOWrapper name='out1.gz' encoding='ascii'> +# timeit_file : None +# timeit_header : None +# timeit_name : all +# whitelist : None +2018-01-03 17:06:39,365 INFO Starting barcode extraction +2018-01-03 17:06:39,377 INFO Input Reads: 100 +2018-01-03 17:06:39,377 INFO Reads output: 100 +# job finished in 0 seconds at Wed Jan 3 17:06:39 2018 -- 0.60 0.05 0.00 0.00 -- 9eecaab6-9562-4b29-8a2e-ae6c1034a14c diff -r 000000000000 -r a6477bafd522 test-data/out_single.log --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/out_single.log Wed Jan 10 19:09:42 2018 -0500 @@ -0,0 +1,38 @@ +# output generated by extract --bc-pattern=XXXNNN --stdin=/tmp/tmpAvWtMd/files/000/dataset_13.dat --stdout /tmp/tmpAvWtMd/files/000/dataset_14.dat --3prime --quality-filter-threshold 10 --quality-encoding phred33 --log=/tmp/tmpAvWtMd/files/000/dataset_15.dat +# job started at Wed Jan 3 17:06:20 2018 on tzk-ThinkPad-T450s -- 8897b3b4-c6e2-4ee3-b6c0-71044a63b7e1 +# pid: 22794, system: Linux 4.10.0-42-generic #46~16.04.1-Ubuntu SMP Mon Dec 4 15:57:59 UTC 2017 x86_64 +# blacklist : None +# compresslevel : 6 +# error_correct_cell : False +# extract_method : string +# filter_cell_barcode : None +# filter_cell_barcodes : False +# log2stderr : False +# loglevel : 1 +# pattern : XXXNNN +# pattern2 : None +# prime3 : True +# quality_encoding : phred33 +# quality_filter_mask : None +# quality_filter_threshold : 10 +# random_seed : None +# read2_in : None +# read2_out : False +# read2_stdout : False +# reads_subset : None +# reconcile : False +# retain_umi : None +# short_help : None +# stderr : <_io.TextIOWrapper name='' mode='w' encoding='UTF-8'> +# stdin : <_io.TextIOWrapper name='/tmp/tmpAvWtMd/files/000/dataset_13.dat' mode='r' encoding='UTF-8'> +# stdlog : <_io.TextIOWrapper name='/tmp/tmpAvWtMd/files/000/dataset_15.dat' mode='a' encoding='UTF-8'> +# stdout : <_io.TextIOWrapper name='/tmp/tmpAvWtMd/files/000/dataset_14.dat' mode='w' encoding='UTF-8'> +# timeit_file : None +# timeit_header : None +# timeit_name : all +# whitelist : None +2018-01-03 17:06:20,360 INFO Starting barcode extraction +2018-01-03 17:06:20,362 INFO Input Reads: 100 +2018-01-03 17:06:20,362 INFO Reads output: 72 +2018-01-03 17:06:20,362 INFO filtered: umi quality: 28 +# job finished in 0 seconds at Wed Jan 3 17:06:20 2018 -- 0.58 0.06 0.00 0.00 -- 8897b3b4-c6e2-4ee3-b6c0-71044a63b7e1 diff -r 000000000000 -r a6477bafd522 test-data/t_R1.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/t_R1.fastq Wed Jan 10 19:09:42 2018 -0500 @@ -0,0 +1,400 @@ +@HISEQ:105:C2UE1ACXX:3:1101:11160:2245 1:N:0:CAGATC +AAAAGTAGTTAATATATTAGATTTGTTTGATAGTGGTAGTATATATTTTTTATTTTAGTATTTAGGAGGTAGAGGTAGATGAATTTTTGAGTTTAAAGTTA ++ +BBBFFFFBFFFFFIIIIIIBIIIIBIFFBFIIBBFFFFFFFFIIIIIIIIIIIIIIIBFIIIIIB7BBBFBBFFFF77 + Extract UMI from fastq files + + macros.xml + + + samtools + + 0: + --gene-tag '$gene_tag' + #end if + #if $input.is_of_type("sam"): + --in-sam + #end if + -I '$input_file' -S deduped.bam && + samtools sort deduped.bam -@ \${GALAXY_SLOTS:-1} -o '$output' -O BAM + ]]> + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ", +replacing with e.g ":". + +Alternatively, if your UMIs are encoded in a tag, you can specify this +by setting the option --extract-umi-method=tag and set the tag name +with the --umi-tag option. For example, if your UMIs are encoded in +the 'UM' tag, provide the following options: +"--extract-umi-method=tag --umi-tag=UM" + +The start postion of a read is considered to be the start of its alignment +minus any soft clipped bases. A read aligned at position 500 with +cigar 2S98M will be assumed to start at postion 498. + + +Methods +------- + +dedup can be run with multiple methods to identify groups of reads with +the same (or similar) UMI(s). All methods start by identifying the +reads with the same mapping position. + +The simpliest method, "unique", groups reads with the exact same +UMI. The network-based methods, "cluster", "adjacency" and +"directional", build networks where nodes are UMIs and edges connect +UMIs with an edit distance <= threshold (usually 1). The groups of +reads are then defined from the network in a method-specific manner. + + "unique" + Reads group share the exact same UMI + + "percentile" + Reads group share the exact same UMI. UMIs with counts < 1% of the + median counts for UMIs at the same position are ignored. + + "cluster" + Identify clusters of connected UMIs (based on hamming distance + threshold). Each network is a read group + + "adjacency" + Cluster UMIs as above. For each cluster, select the node(UMI) + with the highest counts. Visit all nodes one edge away. If all + nodes have been visted, stop. Otherise, repeat with remaining + nodes until all nodes have been visted. Each step + defines a read group. + + "directional" (default) + Identify clusters of connected UMIs (based on hamming distance + threshold) and umi A counts >= (2* umi B counts) - 1. Each + network is a read group. + +Options +------- + +--extract-umi-method (choice) + How are the UMIs encoded in the read? + + Options are: + + - "read_id" (default) + UMIs contained at the end of the read separated as + specified with --umi-separator option + + - "tag" + UMIs contained in a tag, see --umi-tag option + +--umi-separator (string) + Separator between read id and UMI. See --extract-umi-method above + +--umi-tag (string) + Tag which contains UMI. See --extract-umi-method above + +--edit-distance-threshold (int) + For the adjacency and cluster methods the threshold for the + edit distance to connect two UMIs in the network can be + increased. The default value of 1 works best unless the UMI is + very long (>14bp) + +--paired + BAM is paired end - output both read pairs. This will also + force the use of the template length to determine reads with + the same mapping coordinates. + +--spliced-is-unique + Causes two reads that start in the same position on the same + strand and having the same UMI to be considered unique if one is + spliced and the other is not. (Uses the 'N' cigar operation to test + for splicing) + +--soft-clip-threshold (int) + Mappers that soft clip, will sometimes do so rather than mapping a + spliced read if there is only a small overhang over the exon + junction. By setting this option, you can treat reads with at least + this many bases soft-clipped at the 3' end as spliced. + +--multimapping-detection-method (string, choice) + If the sam/bam contains tags to identify multimapping reads, you can + specify for use when selecting the best read at a given loci. + Supported tags are "NH", "X0" and "XT". If not specified, the read + with the highest mapping quality will be selected + +--read-length + Use the read length as as a criteria when deduping, for e.g sRNA-Seq + +--whole-contig + Consider all alignments to a single contig together. This is useful if + you have aligned to a transcriptome multi-fasta + +--subset (float, [0-1]) + Only consider a fraction of the reads, chosen at random. This is useful + for doing saturation analyses. + +--chrom + Only consider a single chromosome. This is useful for debugging purposes + +--per-contig (string) + Deduplicate per contig (field 3 in BAM; RNAME). + All reads with the same contig will be + considered to have the same alignment position. This is useful + if your library prep generates PCR duplicates with non identical + alignment positions such as CEL-Seq. In this case, you would + align to a reference transcriptome with one transcript per gene + +--per-gene (string) + Deduplicate per gene. As above except with this option you can + align to a reference transcriptome with more than one transcript + per gene. You need to also provide --gene-transcript-map option. + This will also add a metacontig ('MC') tag to the reads if used + in conjunction with --output-bam + +--gene-transcript-map (string) + File mapping genes to transripts (tab separated), e.g: + + gene1 transcript1 + gene1 transcript2 + gene2 transcript3 + +--gene-tag (string) + Deduplicate per gene. As per --per-gene except here the gene + information is encoded in the bam read tag specified so you do + not need to supply --gene-transcript-map + +--output-bam (string, filename) + Output a tagged bam file to stdout or -S + +-i, --in-sam/-o, --out-sam + By default, inputs are assumed to be in BAM format and output are output + in BAM format. Use these options to specify the use of SAM format for + inputs or outputs. + +-I (string, filename) input file name + The input file must be sorted and indexed. + +-S (string, filename) output file name + +-L (string, filename) log file name + +Usage +----- + umi_tools dedup -I infile.bam -S grouped.bam -- + + ]]> + +