Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,065,6050AT23.1% 267.3 / 69.5 159intergenic (+40/‑23)ECDH1_RS05245/ECDH1_RS05250tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (72/48);  new base T (25/11);  total (99/60)
Fisher's exact test for biased strand distribution p-value = 3.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

CAAAAAATCTGCAGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATT  >  DH1/1065457‑1065753
                                                                                                                                                       |                                                                                                                                                      
gaataaaTCTGGAGTACAGTAAGTTTCCCGATGCATACGTAGCTCATCTTGATAGAGTAATCGGCTAGGAAGCGTTGGG‑TCTGAGGTTGTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATcagc                                                                                                                                                      <  3:526187/147‑1 (MQ=18)
  aaaacTCTGCTGTAACGTAAGTATCGCGCTGAAGACGTAGCTCAGCGGTATAGAGGACTCGGCTACGACGGGTAGGG‑TCTCAGCTTGTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATcagcag                                                                                                                                                    <  7:385405/151‑1 (MQ=9)
       tCTGCAGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTaaa                                                                                                                                                >  5:753463/1‑150 (MQ=255)
        cTGCAGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATAATACTAAATATGCATAAATGAACGGTAAG‑‑AATATCAACAGTAatat                                                                                                                                              >  4:521381/1‑151 (MQ=11)
          gCAGTAAAGTCAGTTTCGCGATGCATGCGTAGGTGAGCTGGATAGAGTACTCGGCTACGATCCGTTGGG‑TCGGCGGGTGTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAaaata                                                                                                                                             <  7:599361/150‑1 (MQ=18)
          gCAGTAAAGTAAGATTCCCTAGGCATCCGTAGCTCAGCTGGATAGACTACTCGGCTACGAACCGTGCGG‑TCGTAGGTTGTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAaaataa                                                                                                                                            <  7:574143/151‑1 (MQ=21)
               aaagtaaGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑c                                                                                                                                     >  5:770832/1‑151 (MQ=255)
               aaagtaaGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATCTTACTAAATATGGCGAAAGTAGGGTTAGC‑‑AAGATGAGCAGTaataaaacgttg                                                                                                                                       >  4:279929/1‑141 (MQ=14)
                     agtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATg                                                                                                                               <  6:752389/151‑1 (MQ=255)
                     agtTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATg                                                                                                                               <  6:757433/151‑1 (MQ=255)
                      gatgtcctttGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGCGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGc                                                                                                                              <  5:358823/142‑1 (MQ=255)
                       ggtCCCGATGCCTCCGTTGGTAATCTGGATAGAGTACTCGGCCCCGAAGCGAGCGG‑TTGGAGGTTCGAATCCTGCGGGCTGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCa                                                                                                                             <  5:91558/149‑1 (MQ=21)
                        ccgactttgacacagTAGATCAGATTGTTTGAGTACTCGTCTATGTACCGATCTG‑GCTTTGCGTCTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCAt                                                                                                                            <  6:956792/137‑1 (MQ=9)
                         tCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACGTGATATGGCTGTAGTAGCGGTATC‑‑AAGATCAGCAGTAAAATAA‑ATTT‑CCCGATGCAt                                                                                                                            >  3:711668/1‑150 (MQ=255)
                           cccctcgccgctgaGCTTATCTGGCTAGACGACTCTGCTCCGATCTGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATTCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCg                                                                                                                         <  4:939903/138‑1 (MQ=18)
                            cGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATTGATCTATTAAGGTGAGG‑‑AATATCAGCCGGAAAATAA‑AGAT‑CCCGCTGGCTCCGt                                                                                                                        >  5:553561/1‑151 (MQ=14)
                             gATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTATTACTTGATCTTGCTTTATGATCGGTATA‑‑AATATCAGCAGTAAAATAA‑ATAA‑CCCCATGcgtaagt                                                                                                                        >  8:799865/1‑144 (MQ=17)
                                     cGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCCATGCCTACGTAGCGCaagt                                                                                                               >  1:271710/1‑148 (MQ=255)
                                            aaGCTGGATAGACTACTCGGCGACGATCGCTCCGC‑TCGGAGCGTCTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATaga                                                                                                        <  7:934979/150‑1 (MQ=17)
                                               cTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGATCTTAGAATCCTCCAGGATGCACAATCTCTTAATTGATATGGATTTAGTAGGGGTATC‑‑AATATCATAAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGt                                                                                                      >  8:21078/1‑150 (MQ=18)
                                               cTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGAACCATCTCTTAGGTGCTATGGGGTTAGGAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCAGATGCATCCGTAGCTCAGCTGGATAGAGtc                                                                                                     >  7:904320/1‑150 (MQ=21)
                                                   aTAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTATTACTTGATATTGCTTTAGTAGGGGTTTC‑‑AATATCTGGAGTAATATAA‑ATTT‑CCCGATGCATAGGTTGATCAGCTGGATTGAGTAATCg                                                                                                 >  4:333739/1‑151 (MQ=21)
                                                   aTAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATAATACTAGATATGCCTTTATGAAGGTTAAC‑‑AATATCAACAGAAAAATAA‑ATTA‑CACGATGCATACCTAGAAGACCTCGATAGAGAACTCg                                                                                                 >  8:392275/1‑151 (MQ=9)
                                                          aCTCGGCTACGAGCGCTGCGG‑TCTGAGGGTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGa                                                                                          <  3:312816/151‑1 (MQ=17)
                                                            gcGGCTACG‑AGCGAGCGGCTAGGAGAGTCTAATCCACGCTTATTCACGATCTCTTACTTGATTTGTCTTTAGTAGCGTTATA‑‑AATATCTGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAAc                                                                                        <  2:435024/150‑1 (MQ=9)
                                                                        gcgcggg‑TCGGAGCGTCTAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGa                                                                            <  7:4852/146‑1 (MQ=255)
                                                                             gg‑TCGGAGGTTCGAATCCTCCCGGATGCACCATATCGTCCTAGATATGGCGTAATTAGCGGAAAC‑‑AAGATCAACAGTAAATAAA‑ATTT‑CCCGATGCATCAGAAACAGAGCGGTATAGAGTACGCGGCTACGAACCCAGCGCTCGGAggaa                                                                        >  4:200942/1‑148 (MQ=255)
                                                                                 cGGAGTTTCGAATCCTCCCGTATTCACCATCTCTTAGTTGATATGGCTTTAGTAGCGTTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGa                                                                     <  8:62929/150‑1 (MQ=255)
                                                                                    aGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑ACCGATGCATACGTAGCTCAGCAGGATAGACTACTACGCAACGAACCCAGCGGTCGGAGGGTCGAAta                                                                  >  6:397869/1‑149 (MQ=14)
                                                                                     ggTTCGAATCCTCCCGGATGCACCATCTTTTACTTGATATGGCTTTAGTAGCGTTATA‑‑AATATCAGAAGTAAAATAA‑ATTT‑TCCTATGCATTCTTAGCTCATCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAAtcct                                                                <  8:903298/151‑1 (MQ=14)
                                                                                     ggTTCGAATCCTCCCGGATGCACCATCTGTTACTTGATATGGCTATAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGCTGCGTCCGTAGCGCAGCGGGTTAGAGTACTCGGTTCGGAGCCGCGCGGGCGGAGGTTCGATtccg                                                                >  4:963383/1‑150 (MQ=14)
                                                                                         cGAATCCTCCCGGATGCACCATCTCTTGATAGATATGGCGTTATGAACGGTAAC‑‑AATATCATCAGGAATTTAA‑ATAT‑CCCTATTCGGCGGTAGATGATAAGGATAGAGTGCTCGGCTGCGCGCCGCGCGGGGGGGGGTTGGatgcggcggg                                                            >  6:882690/1‑142 (MQ=2)
                                                                                           aaTCCTCCCGGATGCACCATCTATTACTTGATATGTCTTTAGTAGCTGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTccc                                                             <  1:753565/148‑1 (MQ=255)
                                                                                           aaTCCGCCCGGCCGCACCCTCTCTAACTAGACCGGGCTTTACTAGCTCTATA‑‑AATATCTGCAGTAAAATAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGa                                                          <  8:1015815/151‑1 (MQ=11)
                                                                                             acctccCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATg                                                        <  3:358894/150‑1 (MQ=39)
                                                                                                ttttttatttaataTCTCTTACATGAAAAGCTTTTATTTGCTTTATA‑‑ATTATTACCAATAATATTA‑TTTT‑CCGTTTTCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATg                                                        <  4:1055152/135‑1 (MQ=0)
                                                                                                  ccGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATG‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGGTGGAATCGTCCCGGCTGCAccc                                                   >  8:468523/1‑150 (MQ=255)
                                                                                                    ggATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑ACCGATGCATCCGTAGCTAAGCAGGATAGAGTACTAGGCTACGAAACGAGCCGTCGGAGGTTCGAAAACTACCGGATGCACCat                                                  >  8:427787/1‑150 (MQ=14)
                                                                                                        gCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATtc                                              >  6:514307/1‑150 (MQ=40)
                                                                                                            cctctcTTACTGGCTACCGCCTTCGGCGCGGGAGC‑‑AATATCAGCAGTAGATTAATTTTT‑CCCGCTCCATCCGTAGCTCAGCTGGTCCGCCTTGACTGCTGCGAAGAGCTGGGGCGTAGCTTCGAATCCTCCCGGATGCACCATATTCTCCg                                          <  3:236778/149‑1 (MQ=2)
                                                                                                               ctctTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATcagca                                                                                                                                                     >  6:1007573/1‑42 (MQ=35)
                                                                                                                tctTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTa                                                                                                                                                  >  8:552104/1‑44 (MQ=35)
                                                                                                                  ttACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCtt                                   >  7:737375/1‑151 (MQ=255)
                                                                                                                   tACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTc                                  >  7:24954/1‑151 (MQ=255)
                                                                                                                   tACTAGAGATGGTTTTTAGAGCGTTAAC‑‑AAAATCAGCAGTAAATTAA‑TTTT‑CCCTCTGCATACGTAGCTCATCTGTATAGACTACTCGGCTACGAACCGATCGGTCTGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATc                                     <  7:7096/148‑1 (MQ=2)
                                                                                                                    aCTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTc                                  >  2:2535/1‑150 (MQ=255)
                                                                                                                       tGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGc                               >  7:683524/1‑150 (MQ=255)
                                                                                                                       tGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGc                               >  7:562007/1‑150 (MQ=255)
                                                                                                                        gATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑GCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACGGAGCGGTCCCCGGTTCAATTCCTGCCCGACGGACCATAGtctcggcacgcgtgtgca                              >  3:138486/1‑136 (MQ=11)
                                                                                                                         atatGGCTTTAGTAGCGGTATC‑‑AATAAAAGCAGTAAAATAA‑ATTT‑CCCGATGAATCAGTAGCGCAGGTTGATAGAGTACTCGGCTACGAACCGAGAAGTAGGAGGTTCGAATCCTCCCGGATGCCCGATATTCTCCGAAATCTTCGGcg                              >  8:841414/1‑148 (MQ=14)
                                                                                                                          tatGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATg                           >  4:120931/1‑151 (MQ=255)
                                                                                                                          tatGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATCCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTAGGAGGTTCGAATCCTCCCGGATGCACAATATTCTCAGAAATCTTCAGCAATg                           >  8:49385/1‑151 (MQ=255)
                                                                                                                          tacggCTTTAGTAGCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATAATACTAAATATTCAgaa                              >  6:94224/4‑148 (MQ=14)
                                                                                                                           atGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATg                           >  5:763629/1‑150 (MQ=255)
                                                                                                                           acggCTTTAGTAGCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATACTCCGAGATATTCATCAAt                            >  3:936967/3‑151 (MQ=14)
                                                                                                                            tGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTACAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGa                          >  7:154015/1‑150 (MQ=255)
                                                                                                                            cggCTTTAGTAGCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGGGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATACGAAATATTCatatatt                           >  3:820386/2‑145 (MQ=14)
                                                                                                                              gCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAgg                       >  7:834855/1‑151 (MQ=255)
                                                                                                                               cTTTAGTAGCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCCCCATATGCTCCGAGATATTCCGCAAGGacg                        >  4:584853/1‑149 (MQ=14)
                                                                                                                               cTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCCAAATCTTCAGCAATGAAGGt                      >  6:464420/1‑151 (MQ=255)
                                                                                                                                tttAGTAGCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTATACGAAATATTCAGCAATGaacg                       >  2:873836/1‑149 (MQ=11)
                                                                                                                                 ttactcgCGATATC‑‑AAAATCTCCAGTAAGATTA‑TTTT‑CCAGAGGCCTCCGTAGCGCAGCTTTATATACGAATCGGCTACTACCCGTCGGGGCTCCGCTTGTAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTa                     <  1:677354/144‑1 (MQ=2)
                                                                                                                                 ttAGTAGCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGAAATGAAGGt                      >  4:802633/1‑151 (MQ=11)
                                                                                                                                 ttAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTAt                    >  1:587862/1‑151 (MQ=255)
                                                                                                                                 ttAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTAt                    >  1:1024766/1‑151 (MQ=255)
                                                                                                                                  tagtagCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTACGAGATATTGATCaagtatcggt                     >  1:917502/1‑143 (MQ=14)
                                                                                                                                   agtagCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCTACATATTGATGAATGAAGGttt                    >  1:341769/1‑149 (MQ=14)
                                                                                                                                   agtagCGGTATCAAAAAATCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTACGATATATGCAGCAATGAAGGTAt                    >  1:381569/1‑151 (MQ=11)
                                                                                                                                   agtagCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTTACCGGATAGA‑GCGTAGGACACCTCGAGTGACAACTCGTCACAGCGCAGAGCGGTCGGAGGTTCGAATCGTCGCGGATGCACCATATTCTCCGAAGTCTTCCGcgaggccggttggg                  >  7:163688/1‑138 (MQ=2)
                                                                                                                                   agtagCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCg                  >  5:958182/1‑151 (MQ=255)
                                                                                                                                   agtagCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAAACTTCAGCAATGAAGGTATCg                  >  1:476330/1‑151 (MQ=255)
                                                                                                                                   agtagCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑ACCGATTAGTGCGTAGCTCAGCAGGAAAGAGTACAAGGCTCACAACCGAGAGGTCGGAGGTTCGAATCCTCCCGGATGCAACATATTGTCCGAAATCTTCAGCAATGAAGGTATCg                  >  7:163685/1‑151 (MQ=14)
                                                                                                                                   aggcgggatata‑‑AATATCCGCAGTACACTAA‑TTTA‑CCCTATGCCGCCGGAGCTCAGCTGGCGAGAGTACTCGGCCACGACCCGATCGGGCGGAGGGTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATc                   <  1:184494/142‑1 (MQ=2)
                                                                                                                                    tttggtttata‑‑AAAATCAGCAGTAAACTAA‑AATT‑CGCGGTGCTTCCTTAGCTCCTCTGGTTTGAGCACGCGCCCTCGACGCGCCGGCGCTGCGGTTGTAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCg                  <  3:1065631/143‑1 (MQ=9)
                                                                                                                                     taggtttata‑‑AATATCTGCAGTAAACTAA‑TTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGa                 <  8:523738/144‑1 (MQ=11)
                                                                                                                                      agCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACGATATTCTCCGAAATCTTCAGCAATGAAGGTATCGa                 >  4:356243/1‑149 (MQ=255)
                                                                                                                                        gtctata‑‑AAAATCTGCAGTAAATTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAagag             <  2:341600/148‑1 (MQ=11)
                                                                                                                                        cggatcaaaaaaatCTGCAGTAAAGTAA‑CTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGaa                <  3:454266/138‑1 (MQ=14)
                                                                                                                                        cGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑AACGATGCATCCGTAGCTCAGATGGATAGAGTACTCGGCTACGAAGCGAGCGGTCGGAGGTTCGAATCCTCACGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGatg               >  2:995096/1‑147 (MQ=14)
                                                                                                                                         tttata‑‑AATATCAGCAGTAAATTAA‑TTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAagag             <  4:667583/148‑1 (MQ=255)
                                                                                                                                          gtgtctaaaTATCAGCAGTAAACTAA‑CTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAagag             <  3:407001/144‑1 (MQ=255)
                                                                                                                                          gTATC‑‑AAGCTCACCACTAAAATAA‑ATTT‑CCCGATTCCTCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGt           <  7:398073/151‑1 (MQ=14)
                                                                                                                                           tataaaaaaatCAGCAGTAAACTAA‑CTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAgg                       <  1:746446/131‑1 (MQ=14)
                                                                                                                                            aTC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGt           >  2:698316/1‑149 (MQ=255)
                                                                                                                                             tC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGc        >  1:1072416/1‑151 (MQ=255)
                                                                                                                                             tC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CACGATGCATCCGTAGCTCAGATGGATAGAGAACTCGGCTAAGAACCGAGCGGTCGTAGCTAAGAATCCACCAGGATGCACCATATTCTACGAAATCTTCAGAAATGAAGGTATAGAATAGGTAGc        >  8:754334/1‑151 (MQ=255)
                                                                                                                                             cc‑‑AAAATCAGCAGTAAACTAA‑TTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTCATGGTCTACGCAGCTTTCGTTCGGTGTTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGc        <  3:99780/150‑1 (MQ=14)
                                                                                                                                              c‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCg                                                                                      >  3:689947/1‑72 (MQ=255)
                                                                                                                                              c‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATtct                                             >  4:666997/1‑113 (MQ=40)
                                                                                                                                              c‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCg       >  7:105401/1‑151 (MQ=255)
                                                                                                                                               taaaTATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGa                                                                                     <  6:77643/72‑1 (MQ=255)
                                                                                                                                               aaaaaatCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTCCGAAATATTCAGCAATGAAGGTATCCaa                >  4:454008/6‑143 (MQ=14)
                                                                                                                                               aaaaaatCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATAGTCCTAGATATTGAGGAAGTAACGGATCGAagatctacc        >  6:466356/6‑145 (MQ=11)
                                                                                                                                                  atatCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGtaagc                                     >  1:842604/1‑117 (MQ=255)
                                                                                                                                                   tatCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGa     >  5:346560/1‑150 (MQ=255)
                                                                                                                                                   aatCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTGTCCGAAATATTCAGCAAGTAAGGTTTGGAAGATCTAGAGGa     >  5:595000/2‑150 (MQ=11)
                                                                                                                                                    atCTGCAGTAAAGTAA‑GTTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTACGAGATATTCAGCAATTAAGGTATAGAATAGGTAGAGTAAt   >  1:919000/1‑151 (MQ=14)
                                                                                                                                                      cagcagTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTc                                                                                                  >  6:834262/1‑54 (MQ=255)
                                                                                                                                                      cagcagTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAAtt  >  8:249141/1‑150 (MQ=255)
                                                                                                                                                       agcagTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGc                                                                                   >  3:14560/1‑68 (MQ=40)
                                                                                                                                                       |                                                                                                                                                      
CAAAAAATCTGCAGTAAAGTAAGTTTCCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGG‑TCGGAGGTTCGAATCCTCCCGGATGCACCATCTCTTACTTGATATGGCTTTAGTAGCGGTATC‑‑AATATCAGCAGTAAAATAA‑ATTT‑CCCGATGCATCCGTAGCTCAGCTGGATAGAGTACTCGGCTACGAACCGAGCGGTCGGAGGTTCGAATCCTCCCGGATGCACCATATTCTCCGAAATCTTCAGCAATGAAGGTATCGAAGAGGTAGCGGAATT  >  DH1/1065457‑1065753

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑