Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,111,5710AG53.2% 102.8 / inf 496C325R (TGT→CGT) hflKprotease modulator HflK
Reads supporting (aligned to +/- strand):  ref base A (117/114);  new base G (113/150);  total (231/265)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.46e-01

Only 100 of 498 total aligned reads displayed.

CTGGCACCGTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCACGCGGCTGAAC  >  DH1/4111423‑4111721
                                                                                                                                                    |                                                                                                                                                      
  ggCACCGTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCt                                                                                                                                                      <  7:585406/149‑1 (MQ=255)
    cACCGTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTcgcg                                                                                                                                                  <  3:1327168/151‑1 (MQ=255)
      ccGTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTcgcg                                                                                                                                                  >  1:1567368/1‑149 (MQ=255)
      ccGTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAg                                                                                                                                                <  2:595134/151‑1 (MQ=255)
       cGTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGt                                                                                                                                               <  3:554374/151‑1 (MQ=255)
        gTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTGAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTa                                                                                                                                              <  4:1594359/151‑1 (MQ=255)
          tATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTaa                                                                                                                                             <  2:1720025/150‑1 (MQ=255)
           aTCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCCGCACTTTGCGGGTGTTACCCCACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAAtt                                                                                                                                           >  8:199087/1‑151 (MQ=255)
           aTCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAAtt                                                                                                                                           >  1:1009542/1‑151 (MQ=255)
            tCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAAttt                                                                                                                                          >  5:231271/1‑151 (MQ=255)
            tCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAAttt                                                                                                                                          >  8:1067242/1‑151 (MQ=255)
             cGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAAttt                                                                                                                                          <  3:343413/150‑1 (MQ=255)
             cGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTc                                                                                                                                         <  3:1668282/151‑1 (MQ=255)
              gCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTcc                                                                                                                                        >  5:1589414/1‑151 (MQ=255)
              gCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTcc                                                                                                                                        <  7:1635902/151‑1 (MQ=255)
                   tCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGgcg                                                                                                                                    <  1:1150379/150‑1 (MQ=255)
                   tCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGgc                                                                                                                                     <  6:16581/149‑1 (MQ=255)
                       cgcAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGCCt                                                                                                                               >  7:924551/1‑151 (MQ=255)
                        gcAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCtt                                                                                                                              <  1:1469947/151‑1 (MQ=255)
                        gcAGGGCCGTAACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGATGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTAACAAACACTTTTTCAATCGTATCGAAATACAGACGCTCGCGAGTAATTTCCCGCGCGGCtt                                                                                                                              >  1:859803/1‑151 (MQ=255)
                          aGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCttt                                                                                                                             >  1:1828235/1‑150 (MQ=255)
                           ggggCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTtat                                                                                                                           <  1:447523/151‑1 (MQ=255)
                           ggggCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTtat                                                                                                                           <  8:257732/151‑1 (MQ=255)
                           ggggCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGATACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTtat                                                                                                                           <  6:1341920/151‑1 (MQ=255)
                             ggCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTatat                                                                                                                         >  5:1533929/1‑151 (MQ=255)
                                gTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTcc                                                                                                                      >  6:826927/1‑151 (MQ=255)
                                  tACCACCTTTCAGCATCTGGGCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCg                                                                                                                     <  8:382041/150‑1 (MQ=255)
                                   accaccTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACGCTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCgg                                                                                                                    <  5:455777/150‑1 (MQ=255)
                                     caccTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGc                                                                                                                   >  4:1547994/1‑149 (MQ=255)
                                     caccTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGc                                                                                                                   >  6:1289428/1‑149 (MQ=255)
                                       ccTTTCAGCATCTGGTCTAACTTAAGAACCATCACGTTGCCACCTTTATCTTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGaa                                                                                                               <  8:1556408/151‑1 (MQ=255)
                                       ccTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGaa                                                                                                               <  4:311484/151‑1 (MQ=255)
                                       ccTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGa                                                                                                                >  7:1000980/1‑150 (MQ=255)
                                           tCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGtt                                                                                                            <  2:1261745/150‑1 (MQ=255)
                                             aGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAg                                                                                                         <  7:1378513/151‑1 (MQ=255)
                                               cATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGc                                                                                                        >  2:858311/1‑150 (MQ=255)
                                                aTCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCa                                                                                                       >  5:167786/1‑150 (MQ=255)
                                                  cTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCaa                                                                                                      >  5:747786/1‑149 (MQ=255)
                                                      tCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAgcgcg                                                                                                >  6:1855875/1‑151 (MQ=255)
                                                       cTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAgcgcgc                                                                                               >  7:1850460/1‑151 (MQ=255)
                                                         aaCGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAgcgcgc                                                                                               <  6:1903463/149‑1 (MQ=255)
                                                         aaCGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAgcgcg                                                                                                <  5:315566/148‑1 (MQ=255)
                                                              cAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTc                                                                                         >  1:1600588/1‑150 (MQ=255)
                                                                gAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCAcc                                                                                      <  1:290184/151‑1 (MQ=255)
                                                                   ccATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTg                                                                                   <  5:1099413/151‑1 (MQ=255)
                                                                       cAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCt                                                                               <  6:1778158/151‑1 (MQ=255)
                                                                           ttGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTcc                                                                            >  8:1295872/1‑150 (MQ=255)
                                                                            tGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTcc                                                                            >  2:1021321/1‑149 (MQ=255)
                                                                            tGCCACCTTACGCGCTAACCAGCACTTGGCGGGTGATACCCAACACTTTTTCCATCGGCTCGAAATACAGACCCTCGCGAGTAATATCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCCCCCTGAGCTTCCa                                                                           <  2:700312/150‑1 (MQ=255)
                                                                             gCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAgg                                                                         >  2:540572/1‑151 (MQ=255)
                                                                              ccACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAgg                                                                         >  4:771104/1‑150 (MQ=255)
                                                                               cACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGAt                                                                       <  6:323286/151‑1 (MQ=255)
                                                                                aCCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATg                                                                      <  4:1152166/151‑1 (MQ=255)
                                                                                aCCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATg                                                                      <  7:1930608/151‑1 (MQ=255)
                                                                                 ccTTTATCGTTAACCAGCACTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGt                                                                    >  4:1393847/1‑151 (MQ=255)
                                                                                     tATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCt                                                                  <  4:54097/150‑1 (MQ=255)
                                                                                        ggTTAAGCAGCACTTTGCGGGTGTTTCCCAACACTTTTTCCATCGTCTCTATATACAGACACTCGCGATTAATTTCCGGCTCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCGACTTCTCCCTGAGCTTCGAGGATGGTCTGGGc                                                              <  6:642329/150‑1 (MQ=255)
                                                                                           tAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCtt                                                           >  7:1462520/1‑151 (MQ=255)
                                                                                           tAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCt                                                            <  2:334958/150‑1 (MQ=255)
                                                                                            aaCCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTg                                                          >  2:1399933/1‑151 (MQ=255)
                                                                                            aaCCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTg                                                          >  1:675100/1‑151 (MQ=255)
                                                                                                  cACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTacgc                                                     >  3:854253/1‑150 (MQ=255)
                                                                                                    cTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTacgcacg                                                  >  7:109425/1‑151 (MQ=255)
                                                                                                      ttGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCAcgcg                                                >  3:1050106/1‑151 (MQ=255)
                                                                                                       tGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCAcgcg                                                <  7:762296/150‑1 (MQ=255)
                                                                                                       tGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACgcgc                                               >  4:946125/1‑151 (MQ=255)
                                                                                                        gCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGcc                                              >  3:124197/1‑151 (MQ=255)
                                                                                                           ggTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCctct                                           <  6:1487201/151‑1 (MQ=255)
                                                                                                            gtgtTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCtt                                          <  8:1713983/151‑1 (MQ=255)
                                                                                                             tgtTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTc                                         >  1:1339392/1‑151 (MQ=255)
                                                                                                              gtTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCGGGCGGGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTc                                         <  3:464924/150‑1 (MQ=255)
                                                                                                              gtTACCCAACACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTc                                         <  4:1691017/150‑1 (MQ=255)
                                                                                                              gtTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTc                                         <  4:1275577/150‑1 (MQ=255)
                                                                                                                     aaCACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATAc                                 >  8:1181429/1‑151 (MQ=255)
                                                                                                                     aaCACTTTTTCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATAc                                 >  4:1009207/1‑151 (MQ=255)
                                                                                                                        acTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGtt                              >  2:1511033/1‑151 (MQ=255)
                                                                                                                           ttttCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTg                             <  4:719617/149‑1 (MQ=255)
                                                                                                                             ttCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTgcgc                          <  1:503266/150‑1 (MQ=255)
                                                                                                                              tCCATCGTCTCGATATACAGACGCTCGCGAGTAATTTACGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGcc                         <  7:378718/150‑1 (MQ=255)
                                                                                                                               ccATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTg                       <  5:875996/151‑1 (MQ=255)
                                                                                                                                cATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTgc                      >  2:113813/1‑150 (MQ=255)
                                                                                                                                cATCGTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTg                       >  3:1794743/1‑150 (MQ=255)
                                                                                                                                cATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTg                       >  6:704120/1‑150 (MQ=255)
                                                                                                                                cATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTg                       >  8:571786/1‑150 (MQ=255)
                                                                                                                                  tCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGAcc                    >  3:1122275/1‑151 (MQ=255)
                                                                                                                                   cGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCg                   <  3:1519677/151‑1 (MQ=255)
                                                                                                                                    gTCTCGATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGt                  <  6:670558/151‑1 (MQ=255)
                                                                                                                                     tctcGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGtt                 >  6:497272/1‑151 (MQ=255)
                                                                                                                                      ctagATATACAGACGCTCTCGAGTAATTTCAGGCGCGGCTTTATATTCCGGCAGAATTTTAGCAAAGCGATCCACATCACCCTGACCTTCCAGGAAGGTCACGGCCTTGTACTCACGCTCCACTACGAGGAAACGTTGCGCCTGACCGtt                 >  4:1856/4‑150 (MQ=255)
                                                                                                                                       tcGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCg               >  5:1392141/1‑151 (MQ=255)
                                                                                                                                        cGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCg               <  3:197105/150‑1 (MQ=255)
                                                                                                                                        cGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGc              >  8:768015/1‑151 (MQ=255)
                                                                                                                                         gATATACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCa             <  6:1146534/151‑1 (MQ=255)
                                                                                                                                         gATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGc              >  5:1219867/1‑150 (MQ=255)
                                                                                                                                          atatACAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCAc            <  1:1705095/151‑1 (MQ=255)
                                                                                                                                          atatACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCAc            <  5:649898/151‑1 (MQ=255)
                                                                                                                                           tataCAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCAAGCGCATCTTCGAGGATACGTTGCGCCTGACCGTTCGCAcg           >  6:27514/1‑151 (MQ=255)
                                                                                                                                              aCAGACGCTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCAcgcg         <  8:866080/150‑1 (MQ=255)
                                                                                                                                              aCAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCACGCgg        >  2:196612/1‑151 (MQ=255)
                                                                                                                                                  acacTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTGCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCACGCGGCTGa    <  7:823236/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CTGGCACCGTTATCGCTCTTCGCCGCAGGGGCGTTACCACCTTTCAGCATCTGGTCTAACGGCAGAACCATCAGGTTGCCACCTTTATCGTTAACCAGCACTTTGCGGGTGTTACCCAACACTTTTTCCATCGTCTCGATATACAGACACTCGCGAGTAATTTCCGGCGCGGCTTTATATTCCGGCAGAAGTTTAGCAAAGCGCGCCACTTCACCCTGAGCTTCCAGGATGGTCTGGGCCTTGTACGCACGCGCCTCTTCGAGGATACGTTGCGCCTGACCGTTCGCACGCGGCTGAAC  >  DH1/4111423‑4111721

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