Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,592,2900TG25.0% 115.2 / 11.4 74F212V (TTC→GTC) tsrmethyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base T (19/35);  new base G (17/1);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.54e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCG  >  NC_000913/4592055‑4592537
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            
caATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCATGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCTTCAGCGATAACAATGCCTCCGACAGCCAGGCGATGTGGATTCTGGTGTGCTTGATGATCTTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCt                                                                                                                                                                                                                                               <  1:510334/251‑1 (MQ=255)
caATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCt                                                                                                                                                                                                                                               <  1:351071/251‑1 (MQ=255)
caATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCt                                                                                                                                                                                                                                               <  1:437208/251‑1 (MQ=255)
 aaTGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTTAAGCCTc                                                                                                                                                                                                                                              >  1:129046/1‑251 (MQ=255)
        tGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGCCACAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTCCTGGCGGTTACCTTCTACGCCTAGGTCGGTATTaaaaagtctcttgt                                                                                                                                                                                                                                        >  2:282931/1‑239 (MQ=255)
         ggCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCTTCTTCGCCGTCTGGTTCGGTATTTAAACCTCGCTGGtgg                                                                                                                                                                                                                                      >  2:520821/1‑249 (MQ=255)
            ggAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAAAGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATATTCGCCGGCCGGGTCGGGATTAAAGCCCCGCTGGTAgggc                                                                                                                                                                                                                                   >  1:299932/1‑248 (MQ=255)
               gCTGATCCAACTGTTAGGTGCAGGCAAAATCACCGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTACATTTTCCCGTCATGTTTGGGTATTAAAGctcgccgggggacgga                                                                                                                                                                                                                                  >  2:444850/1‑235 (MQ=255)
               gCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCTCTTTGCCCGCTGGGTCGGGATTAAACCCCCGCTGGtcaagcca                                                                                                                                                                                                                                 >  1:254570/1‑243 (MQ=255)
               gCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTgg                                                                                                                                                                                                                                         <  1:444850/242‑1 (MQ=255)
                 tGATCCAACTGTTAGGTGCAGGCAAAATCAACGCGTTCTTTGATCAGCCGACCCAGGGATATCAGGAGGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGGTCGGCTCCATGATATCGCCGTCAGCGCTCACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGAGTGATGATCGTCGTACTGGCGGGCCCTTTCTCCGTCGTGTTCGGTATTAAAGACTccccatcatcgccgca                                                                                                                                                                                                                               >  1:531102/1‑235 (MQ=255)
                  gATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGGTGGGAGCGCCAATg                                                                                                                                                                                                                             >  1:138744/1‑250 (MQ=255)
                        aCTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGGTCGGGATTAAAACCTCGCTGGGAGCGCCAATGAATCGcc                                                                                                                                                                                                                      >  1:337229/1‑251 (MQ=255)
                         cTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGTTCGGGATTTAAGCCTCGCTGGTAGCGcccctgaattccct                                                                                                                                                                                                                     >  1:232151/1‑239 (MQ=255)
                         cTGTTAGGAGCAGGCAAAATCAACGAGTTATTTGATCAACCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGACGTACTGGCGGACATTTTTGCCGTCTGGGTCGGGGTTAAAGACTCGCTGTTAGCGCCAAAGaaaagcct                                                                                                                                                                                                                     >  2:344624/1‑245 (MQ=255)
                           gTTAGGTGCAGGCAAAATCATCGAGTTCGTTGATCATCCGACCCAGGGATATCAGGACGGGTTCGAGAAGCAGTATGTTGCTTACATTGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATTTGGATTCTTGTTGGGGTTATGATCTTCGTACTGGCGATCATCTTCGCCGTCTGGGTCGGCATTAAATCCTCGCTGGTAGCGCCGATGAATGGCCTg                                                                                                                                                                                                                    <  2:308941/250‑1 (MQ=255)
                           gTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGTGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTg                                                                                                                                                                                                                    <  1:498315/250‑1 (MQ=255)
                              aGGTGCAGGCAAATTCAACGAGTTCTTAGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATAGCTCCGTTAGCGATAACAATGCCGCCTACAGCCAGGAGATGTGGATTCGGTTGGGCGTGAAGATCGGCGTACTGGGGGTCATATTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAAgcgcgtgggtg                                                                                                                                                                                                                <  2:156558/251‑8 (MQ=255)
                                      gCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTTGGTATTAAAACCTCGCCGGGAGCGCCAATTGATCGCCTGATTGACACCAtt                                                                                                                                                                                                       >  1:197266/1‑251 (MQ=255)
                                       cAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAAGGAATCGCCTGATTGACGGCATTc                                                                                                                                                                                                      <  1:140130/251‑1 (MQ=255)
                                        aaaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTCGGTATTAAAGCCTCGCTGGTAGCGCCCATGAATCGCCTGATTGACAGCATTCg                                                                                                                                                                                                     >  1:218146/1‑251 (MQ=255)
                                                 gAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGc                                                                                                                                                                                            <  1:520821/251‑1 (MQ=255)
                                                   gTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGTTCGGTATTTAAACCTCGCTGGGAGCGCCCATGCATCGCCTGGTTGGCAACCTTTGTCATATTGc                                                                                                                                                                                            >  1:411163/1‑249 (MQ=255)
                                                         tGATCAGCCGACCCAGCGATATCAGGACGGTTTCGCGAAGCAGTATGTGGCTTACATGGAGCACAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTTAAGACTCGCTGGTCGCGCCCATTAATCGCCTGATTGACAGCATTCGTCAGATTGCAggcgg                                                                                                                                                                                      >  2:30678/1‑249 (MQ=255)
                                                             cAGCCGACTCAGGGATATCTGGACGGTGTCGATATGCAGTATGTGGCTGACATGGAGCAAGACGAGCGGTTCCATGATATCTCCGTCAGCGATAACATTGCCTCCTACAGCTAGGCGTTGTGTACTCTGGTGTGCGTGATGATCGTCGGACGGGTGGTCTTGTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGTTAGCGCCATTGAATCTCCTGATTAACAGCATTCTTCTGATTGCAGTCGGCGTTCt                                                                                                                                                                               <  2:349154/251‑1 (MQ=255)
                                                                 cGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTg                                                                                                                                                                           <  1:117682/251‑1 (MQ=255)
                                                                  gACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGa                                                                                                                                                                          <  1:196688/251‑1 (MQ=255)
                                                                   aCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGaa                                                                                                                                                                         <  2:476478/251‑1 (MQ=255)
                                                                          gATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAAGCGCCTGATTGACAGCATTCGTCAGATTGCAGGCGGCGATCTGGTGAAACCGAtt                                                                                                                                                                  <  2:74373/251‑1 (MQ=255)
                                                                             atccggACGGGTTCGAGAAGCAGTATGTTGCTTTCATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCTGTAATAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGCTCTGGTGAATCCGATTGa                                                                                                                                                                <  2:378691/246‑1 (MQ=255)
                                                                               cAGGACGGTTTCGATAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAATCCGATTGAGGt                                                                                                                                                             <  2:138744/251‑1 (MQ=255)
                                                                               cAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCACAACGTTCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCAAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTCCTGCCGGTCATCTTCGCCCTCTAGGTCGGTTCTAAAGCCTCCCTGGTAGCGCCCAAGACTCCACCGATTTCCAGCAATCGTACTATTGCAGgcggcggattggggcaaacgattttggg                                                                                                                                                              >  1:228375/1‑229 (MQ=255)
                                                                                 ggACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCTGGTTCGGGATTAAAGCCCCGCTGGTAGCGCCCATGAATCGCCCGATTTACCGCATTCGTCATATTTCCGGCGGCGATCTGGTggaaacggttggggggg                                                                                                                                                           >  2:228251/1‑235 (MQ=255)
                                                                                           gagaAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCTGGGTTGGTATTTAAGCCCCCCTGGTAGGGCCACAGAATCGCCTGGTTGACAGCATTTGTCATATTGCAGGCGGCGGTCCGGGGAAACCGATTGAGGGGGATGGCTCTa                                                                                                                                                  >  1:317023/1‑250 (MQ=255)
                                                                                           gagaAGCAGTATGTGGCTTACATAGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTTGGTATTAAAAACCCGCTGGGAGCGCCAATGAATCCCCTGATTGACAGCATTCGTCCTATTTCAAGGGGCGATCTGGGGAAACCGGTTGAGGGGGGTGGGTCTa                                                                                                                                                  >  1:246504/1‑250 (MQ=255)
                                                                                              aaGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATga                                                                                                                                              <  2:424977/251‑1 (MQ=255)
                                                                                                            ttACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGCCTGGGTCGGGATTAAAGACCCGCTGGGAGCGCACAAGAATCGCCCTATTTGCACCATTCGGCCTATTTACGGCGGCGGTCTGGGGAAAACGGTTTAGGGGGGTAGGTCTAATGAGAGGGGGGAACTgg                                                                                                                               >  2:97164/1‑251 (MQ=255)
                                                                                                             tACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTgg                                                                                                                               <  1:211712/250‑1 (MQ=255)
                                                                                                                         aaCGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAACTTCGCCGTCCGGTTTGGTATTAAAAACCCCATGGTAGGGCCAATGACTCGCCTGATTTGCAACAATTGTCATATTTCCGGCGGCGGTTTTGGGGAAACGATTGTGGGGGATGGCTCTAATGAGATGGGGCCACAGGAAGAGAATTTgtg                                                                                                                  >  1:332552/1‑249 (MQ=255)
                                                                                                                            gATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCCTTTTTGCCGTCTGGGTTGGTATTTAAGCCTCGCTGGTAGCGGCAAAGAATTGCCCGGTTGACAGGATTCGTCTTATTTCAGGGGGGGATCTGGTGGAAACGAATGAGGGGGATGGGTCTTATGAGAAGGGGCAAATGGCAGAGAGGTTTCGcc                                                                                                                >  1:437237/1‑250 (MQ=255)
                                                                                                                            gATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCCTCCGGTTTGGGATTTAAGCCCCGCCGGGAGCGCCCATGAAACACCCGGTTGGCAGGATTTGTCATTTTGCAGGGGGCGGACTTGGGGAACCCATTGAGGGGGGTGGGTCTAATGAAATgggggcacgtggagatagatttggccc                                                                                                                >  2:189562/1‑228 (MQ=255)
                                                                                                                            gATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTTGGTATTAAAGACTCGCTGGTCGCGCCCATGAATAGCCTGAATTACAGCATTCGTCATATTTCAGGCGGCGGTCTGGGGAAAACGATTGAGGTGGATGGGTCTAATGAGGTGGGGCAACTGGCcgagaggtttcggcc                                                                                                               >  1:326281/1‑242 (MQ=255)
                                                                                                                             aTCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCCCCGTCTGGGTCGGTATTAAAACCTCCCTGGTAGAGCCCATGAATCGCCTGATTGACAGCATTCCTCATATTTCAGGCGGCGGTCTGGGGGAACCGGTTGAGGGGGGTGGCTCTAATGAGATGGGGCAACCGGCAGAGAGTTTGGGccct                                                                                                              >  1:452065/1‑249 (MQ=255)
                                                                                                                                 gCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGGGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGGCATATTGCAGGCGGCGAGCTGGTGAAACCGATTGAGGTTGATGGCTTTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGc                                                                                                          <  1:506675/251‑1 (MQ=255)
                                                                                                                                    ccATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCa                                                                                                         <  2:218146/249‑1 (MQ=255)
                                                                                                                                        gATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTTGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGgag                                                                                                    <  1:216904/250‑1 (MQ=255)
                                                                                                                                             cGCCGTCAGCGATAACAATGCCGCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGa                                                                                              <  2:129046/251‑1 (MQ=255)
                                                                                                                                               ccGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTCGGGATTAAAACCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTTCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGGTGGCTCTAATGAGATGGGGGCAATGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATg                                                                                            >  1:166400/1‑251 (MQ=255)
                                                                                                                                                    aGCGATACCAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCTCCGTCTGGGTCGGGATTAAAGCCTCGCTGGTAACGCCCATGAAACGCCTGGATTACAGCCTTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGGTTGCGCCATATGCAGGGAGAGCTGATGCGTAc                                                                                       >  1:183723/1‑251 (MQ=255)
                                                                                                                                                     gCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTAcc                                                                                      <  1:113827/251‑1 (MQ=255)
                                                                                                                                                         tAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAACCTCGCTGGTAGCGCCCATGAATCGCCTGATTGACAGCATTCGTCCTATTGCAGGCGGCGATCTGGTGAAAACGAATGTGGTGGATGGCTCTAATGAGAAGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCggcg                                                                                  >  1:424694/1‑248 (MQ=255)
                                                                                                                                                                  ctcctACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGACATCTTCGCCGTCTGGGTCGGGATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCCGATTGACAGCATTCGTCATATTTCAGGCGGCGATCTGGGGGAACCGCATGGGGGGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATTCCGGGGGGGCTGATGCGTACCGTCGGGGAtgtgg                                                                         >  1:511633/1‑250 (MQ=255)
                                                                                                                                                                    cctACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGAAGTgcgc                                                                       <  1:147456/251‑1 (MQ=255)
                                                                                                                                                                    cctACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTgc                                                                         <  2:232151/249‑1 (MQ=255)
                                                                                                                                                                       aCAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAAc                                                                   <  1:535032/251‑1 (MQ=255)
                                                                                                                                                                            cAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGCCATCTTCGCCGTCTGGGTCGGTATTACAGCCTCGCTGGTGACGCCAATGAAACGCCTGATTGACAGCATTCGTCAAATTGCAGGCGGGGATCTGGGGAAAACGCTTGAGGGGGAGGGGTCTAATTAGATGGGGCAAATGGCAGCGAGTTTGTGCCCTATACCGGGAGAGATGGTGAGTACCGGCGGGGATGTGCGCGACGGGGc                                                              >  2:493656/1‑251 (MQ=255)
                                                                                                                                                                                      ggATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATc                                                    <  1:228251/251‑1 (MQ=255)
                                                                                                                                                                                      ggATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGACGTCTGGCTCGGTATTAAAACCTCGCTGGTTGCGGCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGGTTGAGGGGGATGGCTCTAATGAGATGGGGCAACTGGCCGAAAGTTTGCGCCATATGCAGGGAGAGCTGATGCGGACCGTCGGGGGAGTGCGCCACGGGGCCAAAGCCAtt                                                    >  1:439525/1‑250 (MQ=255)
                                                                                                                                                                                          tCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCAGTCTAGCTCGGTATTAAAGCCCCGCTCGTAGCGCCAAAGAAAAGCCTGATTGACAGAATGCGTCATCTTTCACGCGGCGATCTTGTGAAAACGAATTAGGGGGATGGCCATAATGAGATTGGGCAAACGGCAGAAAGTTTTCGCCATATGCAGGGAGAGCAGGAGCGTAACCTCGGGGGTGTGCGCAAACGGGGcaaagccctctat                                                  >  1:352363/1‑238 (MQ=255)
                                                                                                                                                                                                         gatCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTTGGTATTAAAACCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCCTATTGCAGGCGGCGATCTGGTGAAACCGATTGGGGTGGGTGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGGGAGCTGATGGGTACCGTCGGGGATGTGCGCAACGGGGCCAATGCCCTCTATAACGGTGCCAGCgga                                  >  1:93969/1‑248 (MQ=255)
                                                                                                                                                                                                          atCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAAt                                <  2:8662/251‑1 (MQ=255)
                                                                                                                                                                                                                   cTGGCGGTCCTCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTAGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCaggtc                       >  2:160937/1‑247 (MQ=255)
                                                                                                                                                                                                                   cTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGcc                            <  1:160937/246‑1 (MQ=255)
                                                                                                                                                                                                                       cGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGTCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAAt                   <  1:30678/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggTCATCTTCGACGTCTGGTTCGGTATTAAAGCCTCGCTGGCAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGACGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAAt                   >  2:221999/1‑250 (MQ=255)
                                                                                                                                                                                                                             aTCTTCGACGTATGGTTCGGTCTTAACGCCTCGCTGGCAGCGCCACTGACCCGCCCGATTGACAGCATTCGACATCTTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGACAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGGGATGTGCGCAACCGGGACAATGCCATCTATAGCGGTGCCAGCGAACTAGCCACCGGCAAtcccg               >  2:224796/1‑246 (MQ=255)
                                                                                                                                                                                                                                ttCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGctct           >  2:377102/1‑251 (MQ=255)
                                                                                                                                                                                                                                 tCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGATAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAAGCGCCACCGGCAATAACGAt             <  2:218988/248‑1 (MQ=255)
                                                                                                                                                                                                                                 tCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGGCGGTGATGTGCTCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGAtct           <  2:246504/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGTCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGAtctctc        <  1:344624/251‑1 (MQ=255)
                                                                                                                                                                                                                                      gTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTc     <  1:97164/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGAAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTcg    >  1:462011/1‑251 (MQ=255)
                                                                                                                                                                                                                                           ggTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTcgcg  <  1:282931/250‑1 (MQ=255)
                                                                                                                                                                                                                                           ggTTCGGGATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCAGTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCATGGAGAGTTGATGCGTACCGTCGGTGATGTCCGCATCGTGGCCAATGCCATCTATAGCGTTGCCAGCGAAATCGCCACCGGCAAGAACGATCTCTCTTcgc   <  2:503532/249‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                            
CAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGGATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCGCCATATGCAGGGAGAGCTGATGCGTACCGTCGGTGATGTGCGCAACGGGGCCAATGCCATCTATAGCGGTGCCAGCGAAATCGCCACCGGCAATAACGATCTCTCTTCGCG  >  NC_000913/4592055‑4592537

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