# HG changeset patch # User iuc # Date 1504042641 14400 # Node ID 860bc357b6787a2c58b27a4df85b2eb46baf462b planemo upload commit c8e46ecad0b1473097517e582ed6c43eb0635b36 diff -r 000000000000 -r 860bc357b678 macros.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/macros.xml Tue Aug 29 17:37:21 2017 -0400 @@ -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.0 + diff -r 000000000000 -r 860bc357b678 test-data/group_in1.sam --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/group_in1.sam Tue Aug 29 17:37:21 2017 -0400 @@ -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 860bc357b678 test-data/group_in2.bam Binary file test-data/group_in2.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_in3.bam Binary file test-data/group_in3.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_in4.bam Binary file test-data/group_in4.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_in5.bam Binary file test-data/group_in5.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_in6.bam Binary file test-data/group_in6.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_out1.bam Binary file test-data/group_out1.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_out2.bam Binary file test-data/group_out2.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_out3.bam Binary file test-data/group_out3.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_out3.tab --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/group_out3.tab Tue Aug 29 17:37:21 2017 -0400 @@ -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 860bc357b678 test-data/group_out4.bam Binary file test-data/group_out4.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_out4.tab diff -r 000000000000 -r 860bc357b678 test-data/group_out5.bam Binary file test-data/group_out5.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/group_out6.bam Binary file test-data/group_out6.bam has changed diff -r 000000000000 -r 860bc357b678 test-data/out_R1.fastq.gz Binary file test-data/out_R1.fastq.gz has changed diff -r 000000000000 -r 860bc357b678 test-data/out_R2.fastq.gz Binary file test-data/out_R2.fastq.gz has changed diff -r 000000000000 -r 860bc357b678 test-data/out_SE.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/out_SE.fastq Tue Aug 29 17:37:21 2017 -0400 @@ -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/tmpREDhtd/files/003/dataset_3773.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 +2017-08-29 14:20:44,451 INFO Starting barcode extraction +2017-08-29 14:20:44,461 INFO Input Reads: 100 +2017-08-29 14:20:44,461 INFO Reads output: 100 +# job finished in 0 seconds at Tue Aug 29 14:20:44 2017 -- 0.73 0.05 0.00 0.00 -- 5913ff4c-b1d1-4f65-abec-66ca6684bdbd diff -r 000000000000 -r 860bc357b678 test-data/out_single.log --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/out_single.log Tue Aug 29 17:37:21 2017 -0400 @@ -0,0 +1,37 @@ +# output generated by extract --bc-pattern=XXXNNN --stdin=/tmp/tmpREDhtd/files/003/dataset_3766.dat --stdout /tmp/tmpREDhtd/files/003/dataset_3767.dat --3prime --quality-filter-threshold 10 --quality-encoding phred33 --log=/tmp/tmpREDhtd/files/003/dataset_3768.dat +# job started at Tue Aug 29 14:20:26 2017 on gaius -- 1cfd46e7-5d74-4c25-8cdd-89cf8f41a393 +# pid: 21338, system: Linux 4.10.0-33-generic #37-Ubuntu SMP Fri Aug 11 10:55:28 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 +# short_help : None +# stderr : <_io.TextIOWrapper name='' mode='w' encoding='UTF-8'> +# stdin : <_io.TextIOWrapper name='/tmp/tmpREDhtd/files/003/dataset_3766.dat' mode='r' encoding='UTF-8'> +# stdlog : <_io.TextIOWrapper name='/tmp/tmpREDhtd/files/003/dataset_3768.dat' mode='a' encoding='UTF-8'> +# stdout : <_io.TextIOWrapper name='/tmp/tmpREDhtd/files/003/dataset_3767.dat' mode='w' encoding='UTF-8'> +# timeit_file : None +# timeit_header : None +# timeit_name : all +# whitelist : None +2017-08-29 14:20:26,493 INFO Starting barcode extraction +2017-08-29 14:20:26,495 INFO Input Reads: 100 +2017-08-29 14:20:26,495 INFO Reads output: 72 +2017-08-29 14:20:26,495 INFO filtered: umi quality: 28 +# job finished in 0 seconds at Tue Aug 29 14:20:26 2017 -- 0.75 0.02 0.00 0.00 -- 1cfd46e7-5d74-4c25-8cdd-89cf8f41a393 diff -r 000000000000 -r 860bc357b678 test-data/t_R1.fastq --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/test-data/t_R1.fastq Tue Aug 29 17:37:21 2017 -0400 @@ -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 $group_output: + --group-out '$group_out' + #end if + #if $input.is_of_type("sam"): + --in-sam + #end if + --output-bam + -I '$input_file' -S grouped.bam && + samtools sort grouped.bam -@ \${GALAXY_SLOTS:-1} -o '$output' -O BAM + ]]> + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + group_out + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + ", +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" + +By default, reads are considered identical if they have the same start +coordinate, are on the same strand, and have the same UMI. Optionally, +splicing status can be considered (see below). + +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 +------- + +group can be run with multiple methods to identify group 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. + +Note that the "percentile" method used with the dedup command is not +available with group. This is because this method does not group +similar UMIs as per the network methods. Instead it applies a +threshold for inclusion of the UMI in the output and excluded UMIs are +not assigned to a "true" UMI. + + "unique" + Reads group share the exact same UMI + + "cluster" + Identify clusters of connected UMIs (based on hamming distance + threshold). Each network is a read group + + "directional" + 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. + +The group command can be used to create two types of outfile: a tagged +BAM or a flatfile describing the read groups + +To generate the tagged-BAM file, use the option --output-bam and +provide a filename with the -S option. Alternatively, if you do not +provide a filename, the bam file will be outputted to the stdout. If +you have provided the --log/-L option to send the logging output +elsewhere, you can pipe the output from the group command directly to +e.g samtools sort like so: + +``umi_tools group -I inf.bam --group-out=grouped.tsv --output-bam --log=group.log --paired | samtools sort - -o grouped_sorted.bam`` + +The tagged-BAM file will have two tagged per read: + + - UG = Unique_id. + 0-indexed unique id number for each group of reads with the same genomic position and UMI or UMIs inferred to be from the same true UMI + errors + + - BX = Final UMI. + The inferred true UMI for the group + +To generate the flatfile describing the read groups, include the +--group-out= option. The columns of the read groups file are +below. The first five columns relate to the read. The final 3 columns +relate to the group. + + - read_id + read identifier + + - contig + alignment contig + + - position + Alignment position. Note that this position is not the start position of the read in the BAM file but the start of the read taking into account the read strand and cigar + + - umi + The read UMI + + - umi_count + The number of times this UMI is observed for reads at the same position + + - final_umi + The inferred true UMI for the group + + - final_umi_count + The total number of reads within the group + + - unique_id + The unique id for the 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 + +--method (choice, string) + Method used to identify PCR duplicates within reads. All methods + start by identifying the reads with the same mapping position + + Options are: + + - "unique" + Reads group share the exact same UMI + + - "cluster" + Identify clusters of connected UMIs (based on edit distance + threshold). Each network is a read group + + - "directional" + Identify clusters of connected UMIs (based on edit distance + threshold) and umi A counts >= (2* umi B counts) - 1. Each + network is a read group. + +--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 + +--group-out (string, filename) + Output a flatfile describing the read groups + +--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 group -I infile.bam --output-bam -S grouped.bam -L group.log -- + + ]]> + +