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

GCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTT  >  NC_000913/1478179‑1478473
                                                                                                                                                 |                                                                                                                                                     
gCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGt                                                                                                                                                  <  2:367204/151‑1 (MQ=255)
 cGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGtt                                                                                                                                                 <  2:617130/151‑1 (MQ=255)
             gTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGGGGTTTCATTATGTa                                                                                                                                       >  1:109349/1‑149 (MQ=255)
               cGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGGAAAGGGGGGGGGGTTACATTACGTAtt                                                                                                                                     >  1:679147/1‑149 (MQ=255)
                  aGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAAGGGGGGGGGGTACATTACGGAtttg                                                                                                                                   >  1:630802/1‑146 (MQ=255)
                    tgtgCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGGGAGAAAAGGGGGGGAGGTTACATTACGTATTTGTACCt                                                                                                                              >  1:202664/1‑151 (MQ=255)
                         ttcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCATCACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAATCCGTTTTCTGCTATGGCGGTTGCCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtg                                                                                                                          <  1:42867/150‑1 (MQ=255)
                          tcttATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCtggtgg                                                                                                                         <  2:727094/150‑1 (MQ=255)
                             tATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTggg                                                                                                                        <  2:114224/148‑1 (MQ=255)
                                gAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATc                                                                                                                   <  1:308300/150‑1 (MQ=255)
                                   cgcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGt                                                                                                                 <  2:617484/149‑1 (MQ=255)
                                    gcTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGAGAGAAAAAGGGGGGAGGTTACATTACGTATTCGTACCTGGTGGGGGATCGTCAt                                                                                                              >  1:867029/1‑151 (MQ=255)
                                      tGAAATATAATGCTGGCAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGGTTTCTGCTATGGCGGTTGTCTGAAGGAAAAAGGGGAGGGGTGGCAGTGCGAATTTTGACCTAGTTGGGCATCGtgctt                                                                                                             >  2:120637/1‑146 (MQ=255)
                                         aaTATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGGGAGGAAAGGGGGGGGGGGTACATTTAATATTCGTTTCTCGGGGGGccacggattttgtt                                                                                                         >  2:268619/1‑138 (MQ=255)
                                               aTGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCaa                                                                                                   <  1:676932/151‑1 (MQ=255)
                                                 gCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATg                                                                                                 <  1:895868/151‑1 (MQ=255)
                                                      aaaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGGGGGAAAAGGGGGGGAGGGTACATTACGTGTTTGTACCTGGTGGGGGGTCGTCATTTGTTTTTTTATGGTg                                                                                              >  2:773444/1‑149 (MQ=255)
                                                       aaaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACTGAGAGAAAAAGGGGGGGGGGTACATTACTTATTTGTTCCCGGTGGGGCATCCTCATTTTTTTTTTAAttggggtt                                                                                           >  2:671974/1‑144 (MQ=255)
                                                        aaaTCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACCGCGAGAAAAGGGGGGGGGGTTACATTACGTAATTGTACCTTGTGGGGGATCGGCATTTGTTTTTTAATGGTGGttt                                                                                          >  2:382319/1‑151 (MQ=255)
                                                                ggCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCa                                                                                  <  1:396315/151‑1 (MQ=255)
                                                                 gCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGGAAAGGGGGGGGGGGTTCATTACGTATTCGTACCTGGTGGGGGATCGTCATTTTTTTTTTAATGGTGGTTTTAGGTCCAc                                                                                 >  2:429650/1‑151 (MQ=255)
                                                                  cAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTAACGGGGGGAAAAGGGGGGGGGGGTACATTACGTATTTGTACCTGGTGGGGGATCGTCATTTTTTTTTTAATGGTGATTTTAGGTc                                                                                    >  1:761657/1‑147 (MQ=255)
                                                                       cTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGAGAGAAAAGGGGGGGAGGGTTCATTACGTATTTGTACCTGGTGGGGGATCGTCCTTTTTTTTTTTATGGTGATTTTTGGTCCCCGc                                                                               >  2:537497/1‑147 (MQ=255)
                                                                          ttCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGGGAGGGTACATTACGTATTCGTACCTGGTGGGGGATCGTCATTTGTTTTTTTATGGGGGTTTCAGGGCCCCGc                                                                               >  1:887104/1‑144 (MQ=255)
                                                                             gcgTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCCGAGAGGAAAGGGGGGGAGGTTACATTACGTATTCGTACCTGGTGGGGGATCGTCCTTTGTTTTTTAATGGTGATTTTAGGTCCACGCTGGGtgacgt                                                                     >  2:900415/1‑151 (MQ=255)
                                                                                ttgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGCGAGAAAAGGGGGGGGGGTTACATTTCTCATTCGTACCTGGTGGGGGATCGGCATTTGTTTTTTTATGGTGGTTTTTGGGCCACCCTGGGtggcgt                                                                     >  2:358074/1‑148 (MQ=255)
                                                                                 tgttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTATCGGGGAGAAAAGGGGGGGGGGTTAACTTATTGATTTGTACCTGGGGGGGGATCGGCATTTTTTTTTTAATGGGGGTTTTTGGTCCACGCTGGAtgacgt                                                                     >  2:59560/1‑147 (MQ=255)
                                                                                  gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAGGGGGGGAGGGTATATTTCGTATTTGCACCCGGTGGGGCGACGTCCTTTGTTTTTTAATGGGGGATTTAGGTTCACGGTGGGtgacggtg                                                                   >  2:462043/1‑148 (MQ=255)
                                                                                  gttAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTTACGGGGAGAAAAGGGGGGGGGGGTACATTTAGTATTTGTAACTGGGGGGGGATCGTCCTTTTTTTTTTTATTGTGATTTTAGGGCCCCGCCGGGtgacgttgtcg                                                                >  2:62175/1‑151 (MQ=255)
                                                                                           aGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTgcg                                                       <  1:862214/151‑1 (MQ=255)
                                                                                              aaCGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTg                                                    <  2:343318/151‑1 (MQ=255)
                                                                                                 gCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATa                                                 <  1:428159/151‑1 (MQ=255)
                                                                                                 gCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATa                                                 <  1:306236/151‑1 (MQ=255)
                                                                                                    gATAAACCGTTTTCTGCTATGGCGGTTTCCGGGGAGAAAAGGGGGGGGGGGTACATTACGTATTTGTACCTGGGGGGGGATCGTCATTTTTTTTTTAATGGGGATTTTAGGGCCCCGCTGGATGACGTTGGCGTTTTTTTCGCTGGTAAg                                               >  1:748198/1‑150 (MQ=255)
                                                                                                                   gCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGaa                                <  1:174103/150‑1 (MQ=255)
                                                                                                                      aTGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGaa                            <  2:94261/151‑1 (MQ=255)
                                                                                                                      aTGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGa                             <  2:126146/150‑1 (MQ=255)
                                                                                                                           ggTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTg                       <  2:219569/151‑1 (MQ=255)
                                                                                                                              tACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTa                     <  1:646873/150‑1 (MQ=255)
                                                                                                                                   gagagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCtt                >  2:480481/1‑150 (MQ=255)
                                                                                                                                     gagaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTtata             <  1:692681/151‑1 (MQ=255)
                                                                                                                                      agaAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTatat            >  2:95079/1‑151 (MQ=255)
                                                                                                                                         aaaGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTg         >  2:614323/1‑151 (MQ=255)
                                                                                                                                                gTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACAttt  <  2:146155/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCGATTAATCGTTGTCGGAGTGTGCTTCTTATGAGCGCTGAAATATAATGCTGGAAAAATCTCTGGCAACACTTTTCGCGTTGTTAATTTTAGCAACGCTGATAAACCGTTTTCTGCTATGGCGGTTACCGGAGAGAAAAGGGGGTGAGGTTACATTACGTATTCGTACCTGGTGGGGCATCGTCATTTGTTTTTCAATGGTGATTTCAGGTCCACGCTGGATGACGTTGACGTTTTTTGCGCTGATAAGTTTTCTGGCATTGAAAGAATATTGTACGCTTATATCTGTACATTT  >  NC_000913/1478179‑1478473

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