Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,386,8620GA60.8% 33.7 / 104.2 97intergenic (+31/‑72)yfbL/yfbMputative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (6/32);  new base A (59/0);  total (65/32)
Fisher's exact test for biased strand distribution p-value = 6.23e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.39e-12

ATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACC  >  NC_000913/2386722‑2387007
                                                                                                                                             |                                                                                                                                                           
aTCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGcac                                                                                                                                                 <  1:822879/151‑1 (MQ=255)
   gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAGAGACATATACatt                                                                                                                                             >  1:295700/1‑151 (MQ=255)
   gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGCACatt                                                                                                                                             >  2:1234577/1‑151 (MQ=255)
   gTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAAAGACAAGgacaat                                                                                                                                              >  2:262804/1‑149 (MQ=255)
        aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAACGACAAGGACaatttttt                                                                                                                                        >  1:79151/1‑151 (MQ=255)
        aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCACAAAAAAAGACaaggaaaatttttt                                                                                                                                           >  1:1433770/1‑138 (MQ=255)
        aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGAAATGCACattttttt                                                                                                                                        >  2:1328045/1‑151 (MQ=255)
        aaaCAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGCTATAACTTTGTTATACAAAAGTAAATAGTAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGGAAAGCACattttttt                                                                                                                                        >  2:1627731/1‑151 (MQ=255)
           cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAAAACACatttttttttt                                                                                                                                     >  2:654753/1‑150 (MQ=255)
           cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGAAAatttttttttt                                                                                                                                     >  1:860560/1‑150 (MQ=255)
           cAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGACAAGGAAatttttttttt                                                                                                                                     >  2:1135170/1‑150 (MQ=255)
              tACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCt                                                                                                                                  <  2:64071/151‑1 (MQ=255)
               aCCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCAAAAAAAACACAAGGAAATTTTTTTTTTCCt                                                                                                                                  >  2:776471/1‑150 (MQ=255)
                  aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATc                                                                                                                              <  2:295499/151‑1 (MQ=255)
                   cTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGGAAAAAAGACAAGGACAATTTTTTTTTCCTcct                                                                                                                               >  1:894630/1‑147 (MQ=255)
                       cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCCGAAAAAACAAAATACACAATTTTTTTTTCCCCCTCTGttt                                                                                                                       >  1:1469799/1‑151 (MQ=255)
                           gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGTTCCAGAAAAAAAGAAAAGCAAAATTTTTTTTTCCCCCTCTGtttttta                                                                                                                   >  2:264585/1‑148 (MQ=255)
                           gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACAAGCAACTTTTTTTTTTCCCTCTCTTTTTTTAg                                                                                                                   >  2:1084761/1‑151 (MQ=255)
                           gTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGAAAATTTTTTTTTTCCTCCTCCTTTTTTAg                                                                                                                   >  1:1265130/1‑151 (MQ=255)
                            tGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGa                                                                                                                  <  2:1191152/151‑1 (MQ=255)
                                          gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTccc                                                                                                                                   <  2:289493/122‑1 (MQ=255)
                                             tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTc                                                                                              <  2:1170839/151‑1 (MQ=255)
                                               aaTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCACAAAAAGAGACATGCACATTTTTTTTTTCTTCCTCTGTTTTTAGATGTGaacaaagtaactcaa                                                                                              >  1:1718461/1‑138 (MQ=255)
                                                aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAAc                                                                                          <  2:713613/151‑1 (MQ=255)
                                                   aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTGGAGGCGGGTCCAGAAAAAAAGAAAAAACAATTTTTTTTTTTCTCTTATGTTTTTATATGTGGAACATGAAATTTaaaagg                                                                                       >  1:1186573/1‑147 (MQ=255)
                                                   aTCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGACAAACAAACTTTTTTTTTTCTTTTCTTTTTTTAAATGCGAAAAAACAAAttttaacgg                                                                                        >  2:1643454/1‑144 (MQ=255)
                                                       gAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTTAGAAAAAAAGACATACACATTTTTTTTTTTCTCATTTGTTTTTAAATTTGAAAAAAGaaaatttaaaggggga                                                                                     >  1:14893/1‑138 (MQ=255)
                                                            tGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAaaata                                                                              <  2:748375/151‑1 (MQ=255)
                                                             ggCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGACAAAGACAATTTTTTTTTCCTCCTTTTTTTTTAGATGCGAAAAATGAAAATTTAAGGAGGGAaaaaaa                                                                             >  1:1425433/1‑151 (MQ=255)
                                                              gCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCCGAAAAAAAGACAATCACATTTTTTTTTTTCTCCTCTGTTTTTAAAAGCGAAAACAGAAAATCAACCGGGGAAaaataaa                                                                            >  1:667742/1‑151 (MQ=255)
                                                                tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATg                                                                          <  2:735753/151‑1 (MQ=255)
                                                                  aGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCGGAAAAAAGAGATGCACATTTTATTTTTTCTCCTCTGTTTTTAAATGCGGAAAATGAAAATTTAAGGGGGGAAACAAAAATgg                                                                         >  1:1272075/1‑151 (MQ=255)
                                                                   ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCAAAAAAAAGACATACACACTTTTTTTTTTCCTCCTTTTTTTTTAAAGCGTAACATTAAAATCCACGGAGGGAAAAAAAATtggg                                                                       >  1:1716724/1‑150 (MQ=255)
                                                                   ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAACGACATGCACATTTTTTTTTTTCTCCTCCGTTTTTAGATGCGGAAAATAAAAATTAACAGAGGAAAAAAAAATaggg                                                                       >  2:531434/1‑150 (MQ=255)
                                                                   ggTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCCAAAAAAACGACAAGCCCACTTTTTTTTTCCTCATTTTTTTTAAAAAGGGAAAAAAGAAAATTCAAAGGGGGAAAAAAAAATggg                                                                        >  2:591918/1‑151 (MQ=255)
                                                                    gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAAAAAAACACATTTTATTTTTCCTTTTCTTGTTTTAAAAGGGGAAAAAGAAAATTAAAGGGGGGAAAAAAAAATGGGa                                                                       >  2:551691/1‑151 (MQ=255)
                                                                    gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGGAAAAACGACAAGGCCATTTTTTTTTTCCTCCTCTTTTTTTAAAAGGGGAAAAAGAAAATTTACCGGGGGAAAAAAAAATGGGa                                                                       >  2:1023907/1‑151 (MQ=255)
                                                                    gTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGGCCCGCCAAATTTTTTTTTCCTCCTCCTTTTTTAAATAGGGAAAAAGAAAATTAAACGGGGGAAAAAAAAATGGGa                                                                       >  2:20213/1‑151 (MQ=255)
                                                                     tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGGAAAAAACACAAGAACAATTTTTTTTTCCTCATCTTGTTTTTAATTGGGAAAAAGAAAATTAAACGGGGGAAAAAAAAATGGGaa                                                                      >  1:336139/1‑151 (MQ=255)
                                                                        tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAt                                                                 <  1:827238/151‑1 (MQ=255)
                                                                        tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGAt                                                                 <  1:153822/151‑1 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGAAATGCAAATTTTTTTTTTCCCCATCTGTTTTTAAAAGGGAAAAAATAAAATTAACCGGGGGAAAAAAAATTGGAATGATTgggtt                                                           >  1:216953/1‑148 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAACACAAGAACATTTTTTTTTTCCTCCTCTGTTTTTAGATTCGGAACAAGAAATTTAAAGGAGGAAAAATAAAATGGGATGATTTGCtt                                                           >  1:1580206/1‑150 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAAACAAGCACATTTTATTTTTTCTCCTCTGTTTTTAGAAGAGAAACAATAAAATCTACAGAGGGAAAATAAAAGGGGATGaattggta                                                           >  1:209209/1‑144 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGACAAGCAAAATTTTTTTTTCCTCATCCGTTTTTAAAAGGGGAAAAAAAAAATTAAACGAGGAAAAAAAAATAGGGATGAATGGCTa                                                           >  2:789956/1‑151 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCTAAAAAACGACATGAACATTTTTTTTTTCCTCCTCTGTTTTAAAAAGCGAAAAATGAAAATCAAAGGAGGGAAAATAAAATGGAAAGATTGGGTa                                                           >  1:192004/1‑151 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCACATTTTATTTTTTCTCCTCTGTTTTTAGAAGCGGAAAATGAAAATCAAAGGAGGAAAAAAAAAAGGGAATTATTTGCTa                                                           >  1:813306/1‑151 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACCACACGCACCTTTTTTTTTTCCCCCTTTGTTTTTAGATTCGGAAAATCAAATTCTACGGGGGGAAAAAAAATAGGGATTATTTGCtt                                                           >  2:1187576/1‑150 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAAAAGAGAAACAAGAAATTTTAAGGGGGAAAAATAAATTGGGAAGATTGGGTa                                                           >  1:729713/1‑151 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt                                                            <  1:1515701/150‑1 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGGCAAGGACATTTTTTTTTTCCTCATCTTTTTTTAGAAGCGGAACAAAAAATTTAACGGGGGGAAAAAAAATGGGGATGATTTGCtt                                                           >  2:418670/1‑150 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAACGAAAAGCAAATTTTTTTTTTCCTCCTCTGTTTTAAGATGCGAAAAAAGAAATTTAAAGGAGGAAAAAAAAATTGGAATGTTTGGCtt                                                           >  2:649302/1‑150 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAAAAAAAACCAAATATTCATTTTATTTTTTCTCCTCTGGTTTTAGATAAGAAACATGAAAATTTACTGAGGAAAAAAAAAAGGGGAAGATTTGCTa                                                           >  1:1734079/1‑151 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGACAAGAACATTTTTTTTTTCCTCTTCTGTTTTAAGATGGGAAACATGAAAATTAACCGGGGAAAAATAAATTGGGATGAtttggtt                                                           >  1:486867/1‑146 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAACAAATGCACATTTTATTTTTTCTCATCTGTTTTAAGATGGGGAACATGAAAATCAAAGGCGGAAAAAAAAATAGGGATGATTGGcga                                                           >  1:127670/1‑149 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCAAAAAAACCACAAAGACATTTTTTTTTTCCTCCTCTGTTTTAAAATGCGAAAAAAGAAAATCAAACGGGGGAAAAAAAATTGGGATGAATTGGc                                                             >  2:1317772/1‑150 (MQ=255)
                                                                              gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCACAATTTTTTTTTTCTCATCTGTTTTTAGATGCGAAACATGAAAATCAAAGGAGGGAAAATAAATGGGGATGATTGGCTa                                                           >  1:1429311/1‑151 (MQ=255)
                                                                               tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAt                                                          <  1:583608/151‑1 (MQ=255)
                                                                               tgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTAt                                                          <  1:1030933/151‑1 (MQ=255)
                                                                                gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAATGCACATTTTTTTTTTCCTCATCTGTTTTTAAAAGCGAAACAATAAAATTAAAGGGGGAAAAAAAAAAGGGGATGATTGGGTAtt                                                         >  2:156130/1‑151 (MQ=255)
                                                                                gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGaaaaa                                                                                <  1:567956/129‑1 (MQ=255)
                                                                                gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAAAACATTTTTTTTTTCCTCATCTGTTTTTATAAACGAAAAAATAAAATTAACCGAGGGAAAATAAATTGGGATGATTGGGTaat                                                         >  1:367498/1‑149 (MQ=255)
                                                                                gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGAAAAGCAAATTTTTTTTTTCCTCATCTTTTTTTAAATGAGGAAAAATAAAATTAAAGGGGGGAAAAAAAATTGGGattatttgttttt                                                           >  2:1680940/1‑139 (MQ=255)
                                                                                gtTACAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCACATTTTTTTTTTCCTCATCTGTTTTTAGATGCGAAACATGAAAATTAAAGGAGGGAAAAAAAAAGGGGATTATTGGGTAtt                                                         >  2:210158/1‑151 (MQ=255)
                                                                                       cTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACAATTTTTTTTTCCTCATCTTTTTTTAAATGCGAAAAAAGAAAATCAAACGAGGAAAAAAAAAAGGGGAAGGTTTGCTAtttttccga                                                  >  2:254629/1‑146 (MQ=255)
                                                                                         ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGGCAAGCACATTTTTTTTTTTCTCCTCTTTTTTTAAATGCGAAAAAAGAAAATTAAAGGGGGGAAAAAAAATGGGGATGATTTGGTATTTTTTAGaaa                                                >  2:631882/1‑151 (MQ=255)
                                                                                         ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAACACAAAAAAAATTTTTTTTTTCTCCTCTTGTTTTAAAAAGGAAAAATTAAAATTTAACGAGGGAAAAAAAATGGGGATGGTTTGCTCTTTTTTAGaaa                                                >  1:476220/1‑151 (MQ=255)
                                                                                           gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGACATGCCAAATTTTTTTTTCCCCCTCTGTTTTTAAAAGCGAAACAAGAAAATTAAAGGCGGGAAAATAAAATGGGATTATTGGGtttttttttgaaaat                                              >  2:509750/1‑142 (MQ=255)
                                                                                           gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCACAAAAAACGGCCAGAACATTTTTTTTTTCCTCCACTTTTTTTAAATGCCAAACAAGCAAATTAACCGGGGGAAAACAAATGGGAATTGTTTGCTTTTTTTGCAGaatt                                               >  2:1626600/1‑149 (MQ=255)
                                                                                            ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCCCCTCTTTTTTTAAATGCGAAAAAAAAAAATTAAAGGGGGGAAAAAAAATGGGGATGATTTGGTATTTTTTAGaaaatt                                             >  2:167729/1‑148 (MQ=255)
                                                                                            ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTg                                             <  2:1524046/151‑1 (MQ=255)
                                                                                            ttATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCACATTTTATTTTTTCTCATCTGTTTTTAGAATCGAAACATGAAAATTAAAAGAGGGAAAAAAAATAGGGATGATTGGGTAtttttgtgaaaatt                                             >  1:371040/1‑141 (MQ=255)
                                                                                              ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc                                           <  1:156512/151‑1 (MQ=255)
                                                                                              ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGAc                                           <  2:1421919/151‑1 (MQ=255)
                                                                                              ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAAAGAACATTTTTTTTTTTCTTCTCTGTTTTTAAAAGGGGAACATTAAATTTTAAGGGGGAAAAAAAAAATGGAATGATTTGGTATTTTTCAGGAAATGaa                                           >  1:111093/1‑150 (MQ=255)
                                                                                                acaacaGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTc                                         <  2:146727/151‑1 (MQ=255)
                                                                                                 caacaGTAAATAATAAAAATTTCTCGAGGCGGGTTCAGAAAAAAAGACAAGTACATTTTATTTTTCCTCCTCTGTTTTTAGAAGAGAAAAAAGAAAATCAAACGGGGGAAAAAAAATGGGAAAGATTGGGCATTTTTTAGAAAATGAATCa                                        >  2:1048724/1‑151 (MQ=255)
                                                                                                   acaGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATGCACATTTTTTTTTTCCTCCTCTGTTTTTAAATGGGAAAAATTAAATTTAAAGGGGGGAAAAAAAAATGGGAATGTTTGGTTTTTTTTGAAAAAttaatcaaa                                       >  2:736035/1‑144 (MQ=255)
                                                                                                   acaGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAAGGACATTTTTTTTTTCCTCCTCTGTTTTTTGATGCGAAAAAAGAAAATTAAAGGAGGGAAAAAAAAAAGGGATGATTTGGTATTTTTTAGaaaattaatcaaa                                      >  1:1540800/1‑141 (MQ=255)
                                                                                                    caGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAACGAAATGAACATTTTTTTTTACCCCCTCTGTTTTTAAAATCGAAACAAGAAATTTTAAGGGGGGAAAATAAAATGGGAAGATTGGCACTTTTTTTGAAAATTAATCAGaa                                     >  2:319310/1‑151 (MQ=255)
                                                                                                       tAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAg                                  <  1:1634283/151‑1 (MQ=255)
                                                                                                         aataatAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAt                                <  2:1226390/151‑1 (MQ=255)
                                                                                                          ataataAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGAtt                               <  2:1215076/151‑1 (MQ=255)
                                                                                                                   aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTg                      <  2:530227/151‑1 (MQ=255)
                                                                                                                   aTTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTg                      <  1:332239/151‑1 (MQ=255)
                                                                                                                             gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc              <  1:962302/149‑1 (MQ=255)
                                                                                                                             gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTc              <  2:895144/149‑1 (MQ=255)
                                                                                                                             gCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCt             <  1:150458/150‑1 (MQ=255)
                                                                                                                                 gTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGa        >  1:112903/1‑151 (MQ=255)
                                                                                                                                    ccAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaa      <  1:1206328/151‑1 (MQ=255)
                                                                                                                                      aGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa    >  2:1285954/1‑151 (MQ=255)
                                                                                                                                       gAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGaaaaa    >  2:41577/1‑150 (MQ=255)
                                                                                                                                        aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAcc  <  1:1271614/151‑1 (MQ=255)
                                                                                                                                        aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc   <  1:1732570/150‑1 (MQ=255)
                                                                                                                                        aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc   <  1:771263/150‑1 (MQ=255)
                                                                                                                                        aaaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc   <  2:257753/150‑1 (MQ=255)
                                                                                                                                         aaaaGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc   >  2:448602/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                           
ATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCACATTTTATTTTCCCTCATCTGTTTTAAGATGCGAAACATGAAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACC  >  NC_000913/2386722‑2387007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: