Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,478,3240TG41.8% 78.6 / 23.6 56G33G (GGT→GGG) ynbBputative CDP‑diglyceride synthase
Reads supporting (aligned to +/- strand):  ref base T (8/24);  new base G (23/0);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 4.74e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACA  >  NC_000913/1478183‑1478470
                                                                                                                                              |                                                                                                                                                  
ttAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGATTTA‑GGCGGGAAAACGGGGGGAgg                                                                                                                                                >  1:832479/1‑146 (MQ=255)
  aaTCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCC‑CGGGAGAAAAAGGGGGGGGGGTACAtt                                                                                                                                         >  2:551278/1‑151 (MQ=255)
        tGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTa                                                                                                                                    <  1:839016/150‑1 (MQ=255)
            ggAGTGTGCTTCTTATGAGCGCTGAAATTTAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGGAAAAGGGGGGGGGGTACATTACGTATTTGt                                                                                                                               >  2:984080/1‑151 (MQ=255)
             gAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTa                                                                                                                              <  2:201900/151‑1 (MQ=255)
                 gtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtg                                                                                                                          <  1:654784/151‑1 (MQ=255)
                  tgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAC‑CGGGAGAAAAGGGGGGGAGGGTAAATTTCGTATTCGTAccctg                                                                                                                         >  1:790077/1‑148 (MQ=255)
                      tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGATACC‑GGGGAGGAAAGGGGGGGAGGGTTCATTATGTATTCTTGCCCgggggg                                                                                                                     >  1:560433/1‑151 (MQ=255)
                            gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCC‑GGAGGGGAAAAGGGGGGGGGGTAAATTTCATAATCGCACCCGGTggggg                                                                                                                   >  2:93327/1‑146 (MQ=255)
                            gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATc                                                                                                                <  2:151654/150‑1 (MQ=255)
                            gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCg                                                                                                               <  1:848615/151‑1 (MQ=255)
                                 cTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACT‑CGGGGGAAAAGGGGGGGGGGGTACATTTCGTATTCGGACCTGGGGGGGGGTCGtgct                                                                                                           >  1:251780/1‑147 (MQ=255)
                                  tGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAC‑GGAGAGAAAAGGGGGGGAGGGTACATTACGTATTTGTACCTGGGGGGGGATCGTCAttt                                                                                                         >  2:603022/1‑151 (MQ=255)
                                    aaaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCC‑GGGGAGAAAAAGGGGGGGGGTTACATTACGTATTCGTACCTGGTGGGGCGTCGTCATTTGt                                                                                                       >  2:105876/1‑151 (MQ=255)
                                    aaaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCAtt                                                                                                          <  1:564753/148‑1 (MQ=255)
                                        ataATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAAC‑CGGGGGAAAAAGGGGGGGGGGTATATTACGTATTCGTACCTGGGGGGGGATTGTCATtttttttt                                                                                                   >  2:935130/1‑151 (MQ=255)
                                         taATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGCGGGAAAAGGGGGGGGGGGTATATTACGTGTTCGTACCTGGTGGGGGATCGTCCTTTGtttttt                                                                                                  >  2:852613/1‑150 (MQ=255)
                                          aaTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAAGGGGGGGGGGTAAACTACATTTTCTTACCTGGGGGGGGGGCGGCATTTGttttt                                                                                                   >  1:853653/1‑149 (MQ=255)
                                           aTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGGGAAAAGGGGGGGGGGGTACATTACGTATTTGTACCCGGTGGGGGAGCGTCATTTGtttttt                                                                                                  >  2:131011/1‑148 (MQ=255)
                                                       tctcTGGCAACACTTTTCTCGTTGTTAATTTTAGCAACGCTTTTAAACCGTTTTCTGCTATGGCGGTTACCTGGTGTGAAAAGGGGG‑GAGGGTGGCTTATGTTGTTGTAACTTTTGGGGTATTGTCACTTGTTTTTTTATTGTGGttttt                                                                                     >  1:213140/1‑148 (MQ=255)
                                                            ggCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGt                                                                                  <  1:321349/148‑1 (MQ=255)
                                                               aaCACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTGACC‑GGAGAGGAAAAGGGGGGGGGGTGCATTACTCATTATGACCTGGTGGGGCGTCCTCCTttttttttttat                                                                                               >  2:848012/1‑128 (MQ=255)
                                                                 cacTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGGAAAAGGGGGGAGGGGATATTATGTGTATGTACCTTGGGGGGGAACGTCTTTTGTTTTTTACTGGTGGTTTTCGGTTCAc                                                                              >  2:413472/1‑147 (MQ=255)
                                                                     tttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGGCCCCGCTGGa                                                                        <  1:857540/149‑1 (MQ=255)
                                                                       tCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtga                                                                     <  2:347567/150‑1 (MQ=255)
                                                                        cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGTTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgac                                                                    <  2:27295/150‑1 (MQ=255)
                                                                        cgcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGAtgacg                                                                   <  1:333464/151‑1 (MQ=255)
                                                                         gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTAC‑GGGGGGGAAAGGGGGGGGGGGTACATTTAGTATTCGTAACTGGGGGGGGATCGGCCTTTTTTTTTTAATGGGGGTTTTAGGGCCACGCTGGGtgaag                                                                   >  2:711902/1‑148 (MQ=255)
                                                                          cgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGAGAAAAGGGGGGGGGGGTACAATACGGATTCGTACCTGGTGGGGGATCGTCATTTTTTTTTTAATTGGGATTTTAGGTCCACGCAGGGtgacg                                                                   >  2:702606/1‑149 (MQ=255)
                                                                            ttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGAGAAAAAGGGGGGGGGGTACATTACGTATTCGTACCTGGGGGGGGATCGTCATTTGTTTTTTAATGGGGGTTTTAGGGCCACGCTGGGtggcggt                                                                 >  2:334039/1‑147 (MQ=255)
                                                                                     ttAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGGGGAAAGGGGGGGAGGGTACATTACGGATTTGTACCTGGTGGGGGATCGTCATTTGTTTTTTCATGGTGGTTTCAGGGCCACGCTGGGTGACGTTGAGGtttttt                                                       >  2:146898/1‑150 (MQ=255)
                                                                                      tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGAGAAAAAGGGGGGGGGGTACATTAAGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTAATGGTGATTTTAGGGGCACGCTGGGTGACGTTGGCGtttttt                                                       >  2:840027/1‑149 (MQ=255)
                                                                                      tAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGGGGGGGTACATTACGTATTCGTACCTGGGGGGGGATCGTCATTTTTTTTTTAATTGGGGTTTTTGGTCCCCCCTGGGTTACGTTTGCGgtttttt                                                      >  2:793819/1‑149 (MQ=255)
                                                                                         cAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGGGAGAAAAAGGGGGGGGGGTACATTACGTATTTGTACCTGGTGGGGCATCGTCATTTGTTTTTTTATGGTGGTTTTAGGTCCACGCTGGGTGACGTTGACGtttttt                                                       >  2:27637/1‑146 (MQ=255)
                                                                                          aaCGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTg                                                 <  2:611252/151‑1 (MQ=255)
                                                                                                 aTAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGtt                                          <  1:264499/151‑1 (MQ=255)
                                                                                                      ccGTTTTCTGCTATGGCGGTTACC‑GGGGAGGAAAGGGGGGGGGGGTACATTACTTATTCGTACCTGGGGGGGGATAGTCATTTGTTTTTTAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTTTGCTGATAAGtttttt                                      >  1:76397/1‑148 (MQ=255)
                                                                                                        gTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGc                                   <  2:845005/151‑1 (MQ=255)
                                                                                                          tgtCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGc                                   <  1:268886/147‑1 (MQ=255)
                                                                                                                  aTGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGaa                         <  2:38282/151‑1 (MQ=255)
                                                                                                                    ggCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAata                       <  1:427038/151‑1 (MQ=255)
                                                                                                                    ggCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAat                        >  1:98858/1‑150 (MQ=255)
                                                                                                                       ggTTACC‑GGAGAGAAAAGGTGGTGATGTTACATTACGTATTCTTACCTGGTGGGGCATAGTCATTTGTTTTTCAATGGTGATTTCTGGTCCACGCTTGATGACGTTGACGTTTGTTGCGCTGATAAGTTTTCTGGCTTTCAAAGAATATTg                    <  1:601252/151‑1 (MQ=255)
                                                                                                                       ggTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTTCCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTg                    <  2:276544/151‑1 (MQ=255)
                                                                                                                          tACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTAc                 >  2:716142/1‑151 (MQ=255)
                                                                                                                           aCC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACg                >  2:790963/1‑151 (MQ=255)
                                                                                                                            cc‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGc               <  1:876370/151‑1 (MQ=255)
                                                                                                                             c‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCt              >  2:882333/1‑151 (MQ=255)
                                                                                                                                 agagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtat           >  2:272070/1‑151 (MQ=255)
                                                                                                                                  gagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtat           >  1:28049/1‑150 (MQ=255)
                                                                                                                                   agaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtata          <  2:142248/150‑1 (MQ=255)
                                                                                                                                      aaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTg      >  2:36170/1‑151 (MQ=255)
                                                                                                                                       aaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGTTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGt     <  1:847612/151‑1 (MQ=255)
                                                                                                                                       aaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGt     <  2:832877/151‑1 (MQ=255)
                                                                                                                                          ggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTAc   >  2:382954/1‑150 (MQ=255)
                                                                                                                                          ggggTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACa  <  2:82945/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACC‑GGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACA  >  NC_000913/1478183‑1478470

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