Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH13,448,5200AC26.0% 591.4 / 52.2 283intergenic (‑80/+76)malZ/proYalpha‑glycosidase/proline‑specific permease ProY
Reads supporting (aligned to +/- strand):  ref base A (143/65);  new base C (2/71);  total (145/138)
Fisher's exact test for biased strand distribution p-value = 1.33e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 301 total aligned reads displayed.

GTAATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGC  >  DH1/3448370‑3448670
                                                                                                                                                      |                                                                                                                                                      
gtAATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGgga                                                                                                                                                        <  5:968297/151‑1 (MQ=255)
   aTGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggagg                                                                                                                                                      >  8:484583/1‑150 (MQ=255)
     gAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagg                                                                                                                                                   <  5:742520/151‑1 (MQ=255)
      aGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGggaggagga                                                                                                                                                  <  7:929643/151‑1 (MQ=255)
            tGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACgggg                                                                                                                                             <  5:1273095/150‑1 (MQ=255)
             gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACggggg                                                                                                                                            <  4:972116/150‑1 (MQ=255)
             gATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                           >  1:274235/1‑151 (MQ=255)
                ctttgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTCTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCgg                                                                                                                                         >  3:217920/1‑150 (MQ=255)
                  ttgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGACGGGGGc                                                                                                                                           <  7:142616/146‑1 (MQ=255)
                   tgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCAGCTTGCCATAGCCTAATGCAGGATCCCTCATCAGCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGTTac                                                                                                                                     >  2:187072/1‑150 (MQ=255)
                   tgctttgTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATa                                                                                                                                      >  1:311593/1‑150 (MQ=255)
                       ttgTTTAACAAATGGGGGCACCGGCAGGTGCCATTCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTTTATTTGGGGGAGGAGGACGGGGGCGGATAAc                                                                                                                                    <  1:1064011/148‑1 (MQ=255)
                           ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGCTCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                              >  4:371129/1‑150 (MQ=255)
                           ttAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCgg                                                                                                                              >  6:947873/1‑150 (MQ=255)
                            tAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                             >  5:682551/1‑150 (MQ=255)
                            tAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGa                                                                                                                             >  3:1176157/1‑150 (MQ=255)
                             aaCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                           >  7:679899/1‑151 (MQ=255)
                              aCAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAAt                                                                                                                           >  8:1067453/1‑150 (MQ=255)
                                aaaTGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGc                                                                                                                        >  1:1135219/1‑151 (MQ=255)
                                   ggggggCACCGGCAGGTTCCATTAATTTAACATAATAATTCCCCTTTTATCAAAATCTGATCACCTTGCCAAAGCATATAGCAGTATCACTAATCATGTGAATGATTTTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                     <  5:1107130/150‑1 (MQ=255)
                                    gggggCACCGGCATGTGCCATTCTTTTAACATCATAATTCCCCTGTGATTAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTTAAGTATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCgg                                                                                                                     <  8:749686/150‑1 (MQ=255)
                                      ggccACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  <  5:233855/148‑1 (MQ=255)
                                       gccACCGGCAGGTGCCATGCATTAAACATAATAATTCCCCTCTGATCAAAATCTGATCACCTTGCCATAGCCTAATGCCGGATCACTCATCATCTGAAGTATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCa                                                                                                                 <  6:677628/149‑1 (MQ=255)
                                        gccccGGCAGGTGCCATTAATTTAACATATTAATTCCCCTGTGATTAAAATCTGATCACCTTGCCATACCCTAATGCAGGATCACTCATCTTCTTAAGTATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAg                                                                                                                <  2:344818/148‑1 (MQ=255)
                                        gCACCGGCAGGTGCCATTAATTTAACATCATAATTCCCCTTTTATGAAAATCTGATCCCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGTTTATTTGGGGGCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAg                                                                                                                <  4:223509/151‑1 (MQ=255)
                                        gCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  >  3:825507/1‑149 (MQ=255)
                                        gCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGCGGCGGATAACTTTCGGAATAGCCGGATc                                                                                                                  >  6:310020/1‑149 (MQ=255)
                                          cccGGCAGGTGCAATGCATTTAACATCATAATCCCCCTGTCATGAAAAGCTCAACCCCTGGCCACCGCCTAATGCAGCAGCACTCATCATCTCAAGCATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGa                                                                                                              <  8:1014901/150‑1 (MQ=255)
                                            cGGCAGGTGCCATGCATTAAACATAATAATTCCCCTGTGATGAAAATCTGATCACCTTCCCATACCCTAATGCAGGATCACTCATCATCTGAAGGATTTTTGGGGGCGGAGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGAtt                                                                                                            <  2:1217592/151‑1 (MQ=255)
                                             ggCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGAttt                                                                                                           >  5:512387/1‑151 (MQ=255)
                                                aGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcct                                                                                                        >  5:84922/1‑151 (MQ=255)
                                                 ggTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGACCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctc                                                                                                       >  8:1120819/1‑151 (MQ=255)
                                                   tGCCATTCATTAAACATCATAATTCCCCTTTTATTAAAATCTTATCACCTTGCCCTACCCTAATCCAGGATCACTCATCCTCTCAAGGATTTTTGGGGGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTtcctcc                                                                                                      <  4:151629/150‑1 (MQ=255)
                                                         tcATATAACATCAAAAATCCCCTGCGATCAAAATCTGATCACCTTCCCATCCCCTACTGCACGAACACTCATCATCTTACGGATTTTTGGGGGCGGAGGCGCGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                <  4:566067/149‑1 (MQ=255)
                                                         gCATATAACATCATAATTCCCGTGTGATGAAAATCTGAACACCTTGCCATAGCCTAATGCAGGATCACTCATAATCTGAAGGATATGTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAAcc                                                                                                <  7:195199/150‑1 (MQ=255)
                                                         gCATATAACATCATAATTCCCCTGTTATGAAAATCTGATCACCTTTCCATAGCCTAATCCAGGATCACTCATCATCTGAAGGATCTTTGGGGGAGGCGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCg                                                                                               <  2:483603/151‑1 (MQ=255)
                                                           aTTTAACACCATAATTCCGCTGTCATGAAAATCTGTTCACCTTGTCATAGCCTCTTCCAGCATCACTATTCTTGTGACGGATATTGGGGGGCGCCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGc                                                                                             <  6:973901/151‑1 (MQ=255)
                                                            tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGcg                                                                                            >  4:674666/1‑150 (MQ=255)
                                                            tttAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCt                                                                                            >  3:1021804/1‑151 (MQ=255)
                                                                aCATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                         >  3:1108996/1‑150 (MQ=255)
                                                                 catcatAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAc                                                                                         >  1:50155/1‑149 (MQ=255)
                                                                   acctAATGCCCCTGTGATCAAAATCGGCTCCCCTGCCCATACCCTCAGGCCCGACCACTCATCATCTGAAGGATATTGGGGGGCGGCGCCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGaa                                                                                     <  3:574877/148‑1 (MQ=255)
                                                                      tAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                  >  1:437491/1‑151 (MQ=255)
                                                                      tAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                  >  3:183681/1‑151 (MQ=255)
                                                                       aaTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGc                                                                                  >  1:1261161/1‑150 (MQ=255)
                                                                        aTTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTCATCCAGCATCACTCATCCTCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCaa                                                                                <  5:659035/151‑1 (MQ=255)
                                                                         ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                 >  1:1102614/1‑149 (MQ=255)
                                                                         ttCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCa                                                                                 >  5:132475/1‑149 (MQ=255)
                                                                          tCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTACTCCTGAACAGGCTTACAGAACGCAAgg                                                                              >  6:879003/1‑151 (MQ=255)
                                                                           ccccTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGACCCGGCTTACAGAACGCAAgg                                                                              >  5:1007894/1‑150 (MQ=255)
                                                                             ccTGTGATCAAAATCTGATCACCTTGCCATAGCCTAAGGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                           <  2:922934/151‑1 (MQ=255)
                                                                              cTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                           >  7:7899/1‑150 (MQ=255)
                                                                              cTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGAtt                                                                           >  4:478963/1‑150 (MQ=255)
                                                                                 tgatgaAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATACCTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGCACGCAAGGATTACTg                                                                       >  6:189112/1‑151 (MQ=255)
                                                                                     caaaccgcgcgcACCTGGCCATCGCCGAATGCACGATCCCTCCTCCTCGGTCGGATAGTGGGCGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGttt                                                                   <  4:798621/140‑1 (MQ=255)
                                                                                      aaaaTGTGAGCACCTTCCCATCGCCTTCTGCTGGATCACTCATCTTCTGAAGTATTTTTGGGGGCGGCGGAGGGGGGGGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                  <  2:760445/151‑1 (MQ=255)
                                                                                      aaaaTCTGCTCACCTTGCCATAGCCTAATGCACCATCACTCATCATCTCAAGGATTTTTGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGtttt                                                                  <  8:117905/151‑1 (MQ=255)
                                                                                       aaaTCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCCTATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTc                                                                 >  8:108532/1‑151 (MQ=255)
                                                                                             gctCACCTTCCCATCCCCTAATCCAGGATCACTCAGCATCTGAAGGTTATTGGGGGGAGGCGGCCGGCGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTcagcca                                                            <  5:754675/148‑1 (MQ=255)
                                                                                               tCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                         >  5:30886/1‑151 (MQ=255)
                                                                                               ccaccgacccATCGCCTACTGCAGGAGCACTCATCAGCTGAAGGTTATTTGGGGGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCt                                                         <  5:1041465/143‑1 (MQ=255)
                                                                                                 cccTAGCCACACCCTAATCCAGGAGCACTCATCAGCTGCAGGATCTTTGGGCGCGGCGGCCGCCCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGa                                                       <  5:603292/150‑1 (MQ=255)
                                                                                                 aCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTgc                                                       >  5:416851/1‑150 (MQ=255)
                                                                                                 aCCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTg                                                        >  3:467770/1‑150 (MQ=255)
                                                                                                  ccTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTAAAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGAc                                                      >  7:629299/1‑151 (MQ=255)
                                                                                                   catGCCATAGCCTAATGCAGGATCACTCATAATCTGAAGGATATTTGGAGGAGGAGGACGGGGGCGGATACCTTTCGGAATAGCCGCATCAGCATGTCCTCCTGAACCGGCTTAGAGAACGAAAGGATTAATGGTTTTgagccacgtcacg                                                     >  4:289485/3‑144 (MQ=255)
                                                                                                      gCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                  >  1:472065/1‑151 (MQ=255)
                                                                                                      cccATAGCCTAATGCAGGAGCACTCATCATCTGAAGGATATTTGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATc                                                  <  5:785845/150‑1 (MQ=255)
                                                                                                        cATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGt                                                >  5:1242686/1‑151 (MQ=255)
                                                                                                         aaccccTAATGCCGGACCACTCCTCATCTGAAGCATATTTGGGGGAGGCGGACGGCGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTg                                               <  8:325779/147‑1 (MQ=255)
                                                                                                            gCCTAATGCAGGATCACTCCTCATCTGCAGGATATGGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                            <  2:1129565/151‑1 (MQ=255)
                                                                                                            cccTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGCGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTggcg                                            <  4:723015/150‑1 (MQ=255)
                                                                                                             ccTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCTGCTCAGGATTTCCTCCTGAACCGGCTTACCGAACGCAAGGCTTACTGGTTTTCAGCCAGCTGAGGATCGTggcgg                                           >  2:612402/1‑151 (MQ=255)
                                                                                                               tAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                         >  7:907383/1‑151 (MQ=255)
                                                                                                               aactGCAGGAGCACTCATCATCTCACGGATATTGGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGgcggcg                                         <  7:530616/148‑1 (MQ=255)
                                                                                                                 aTGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGtt                                       >  2:527073/1‑151 (MQ=255)
                                                                                                                   gCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                      >  8:244814/1‑150 (MQ=255)
                                                                                                                   gCAGGAGCACTCATCATCTCAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGttt                                      <  7:1259763/150‑1 (MQ=255)
                                                                                                                    cagctccACTCATCATGTGAAGGATATTGGGGGGAGGAGGACCGCGGCGGATAACTTTCGGAATAGCGGGAGCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaa                                    <  3:773134/145‑1 (MQ=255)
                                                                                                                     aGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTaaa                                   <  8:411029/151‑1 (MQ=255)
                                                                                                                         acccTCATCATCTGAAGGATATTTGGGGGAGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATc                               <  3:967300/148‑1 (MQ=255)
                                                                                                                          cACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                              >  8:684141/1‑151 (MQ=255)
                                                                                                                          cACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                              >  2:600576/1‑151 (MQ=255)
                                                                                                                          cACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                              >  3:315812/1‑151 (MQ=255)
                                                                                                                           aCTCATCCTCTGAAGGATATTTGGGGGCGGCGGCCGCGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCa                             <  1:180051/151‑1 (MQ=255)
                                                                                                                            cTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATcc                              >  7:1036058/1‑149 (MQ=255)
                                                                                                                             cccaccgcccacggATATGGGGGGGAGGAGGCCGGGGCCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGc                           <  4:1199731/139‑1 (MQ=255)
                                                                                                                              cctcatcTGAAGGATATTGGGGGGAGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGcc                          <  1:329241/149‑1 (MQ=255)
                                                                                                                              catcatCTGACGGATAGTTGGGGGCGGCCCCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAg                            <  1:967636/149‑1 (MQ=255)
                                                                                                                                 catcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAAt                       >  6:1126054/1‑151 (MQ=255)
                                                                                                                                  atcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                      >  4:884894/1‑151 (MQ=255)
                                                                                                                                   tcTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATc                      >  8:178009/1‑150 (MQ=255)
                                                                                                                                           atatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTGTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcac              >  5:996293/1‑150 (MQ=255)
                                                                                                                                            tatTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGcaca             >  1:1190387/1‑150 (MQ=255)
                                                                                                                                             atgtGGGGGCGGCGGCCGGGCGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaac           <  2:273720/148‑1 (MQ=255)
                                                                                                                                               ttGGGGGCGGCGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCacaaca          <  4:1147745/150‑1 (MQ=255)
                                                                                                                                                   gggCGGAGGCCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCa     <  6:470750/151‑1 (MQ=255)
                                                                                                                                                    ggaggcggcCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc    <  1:1007364/151‑1 (MQ=255)
                                                                                                                                                    ggaggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAc    <  7:1179860/151‑1 (MQ=255)
                                                                                                                                                      aggaggaCGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCAcgc  >  1:804196/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GTAATGAGCAGTTGATCTTTGCTTTGTTTAACAAATGGGGGCACCGGCAGGTGCCATGCATTTAACATCATAATTCCCCTGTGATGAAAATCTGATCACCTTGCCATAGCCTAATGCAGGATCACTCATCATCTGAAGGATATTTGGGGGAGGAGGACGGGGGCGGATAACTTTCGGAATAGCCGGATCAGGATTTCCTCCTGAACCGGCTTACAGAACGCAAGGATTACTGGTTTTCAGCCAGCTGACGATCGTGGCGGCGTTTAAACATCCAGCCAATCAACAGCACAACAATCCACGC  >  DH1/3448370‑3448670

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: