Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6570AG42.9% 71.0 / 10.8 49F386S (TTC→TCC) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (7/21);  new base G (21/0);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 3.24e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508513‑3508806
                                                                                                                                                |                                                                                                                                                     
cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGgaaaa                                                                                                                                                   >  2:330805/1‑149 (MQ=255)
             aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGc                                                                                                                                    <  1:769461/151‑1 (MQ=255)
              aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGTAGcc                                                                                                                                   >  1:544768/1‑151 (MQ=255)
                 cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTAAACAAGGTAGccccc                                                                                                                                >  2:122709/1‑151 (MQ=255)
                 cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTAAAAACGGGAGccccc                                                                                                                                >  2:571053/1‑151 (MQ=255)
                 cGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGgcccc                                                                                                                                >  2:58297/1‑151 (MQ=255)
                    agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCcaca                                                                                                                             <  2:217748/151‑1 (MQ=255)
                    agagaTCAATATGCAGATGCTCATGGGCGAGGGCGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTCACACACGGGAGCCCCcac                                                                                                                              >  1:217735/1‑150 (MQ=255)
                          cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATCACACACGGGAGCCCCCacaacaaa                                                                                                                        >  2:674793/1‑148 (MQ=255)
                                gCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg                                                                                                                 <  2:448212/151‑1 (MQ=255)
                                  aGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtgg                                                                                                                <  1:487320/150‑1 (MQ=255)
                                       cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCACCTCAGAAAACgtgggg                                                                                                              >  1:531497/1‑145 (MQ=255)
                                          aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCacac                                                                                                       <  2:104837/151‑1 (MQ=255)
                                          aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGAGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTAGAGCGGCGAAAAATTCAAAAACCGCAACCCCCCCAACAAAAACCgaggggggc                                                                                                           >  1:298954/1‑142 (MQ=255)
                                              gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACGGGAGGCCCCCCCTCAGAAAACgcgggggggcccccc                                                                                                     >  2:39564/1‑138 (MQ=255)
                                               cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTAAAAACCGGAGCCCCCCCATCAGAAAGAGGGGGGCGGCCAcc                                                                                                      >  2:408484/1‑147 (MQ=255)
                                                 aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg                                                                                                <  1:185815/151‑1 (MQ=255)
                                                  gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTTGCGCCGGGGAAAAATCACAAACGGAAGCCCCCACATCACAAAACGCGGGGGGCCCAcc                                                                                                      >  1:61938/1‑144 (MQ=255)
                                                     tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCACACCGAgcggcg                                                                                              >  1:303220/1‑149 (MQ=255)
                                                     tgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCAAAAACCGGAGACCCCCCAACCCAAAGAGGGGGGCGcaca                                                                                                        >  2:296260/1‑139 (MQ=255)
                                                      gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                                                            <  2:762432/150‑1 (MQ=255)
                                                        aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                         <  2:311300/151‑1 (MQ=255)
                                                          cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                                                                          <  1:545189/148‑1 (MQ=255)
                                                               cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACAACGGAGGCCCCACATCAGAAAAAGTGGGGCGCACACCCAGCTGCGTAGCGATCATcc                                                                                  >  1:374033/1‑151 (MQ=255)
                                                                 gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGTGGGCGCCGAAAAAATCAAAAACAGTAACCCCCCCATCAGAAAGCCTGGGGCGGACCCCCGGCCACGGAGCGAACAccccc                                                                                >  1:630454/1‑150 (MQ=255)
                                                                  aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                               <  2:81994/151‑1 (MQ=255)
                                                                         cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCACATCAGAAAACGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGAt                                                                         >  1:627669/1‑150 (MQ=255)
                                                                           cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                      <  2:293395/151‑1 (MQ=255)
                                                                            aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACGGTAGCCCCCACATCAGAAAACGTGGTGGGCACAACCAGCAGCGTAGCGATCATCCCTTGCTGATaaa                                                                      >  2:492362/1‑150 (MQ=255)
                                                                              aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGGGGGGGGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAcc                                                                    >  1:240761/1‑150 (MQ=255)
                                                                              aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                   <  1:35678/151‑1 (MQ=255)
                                                                                    tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGTAGCCCCCCCCTCAAAAAGCGGGGGGCGCACACCCAGCTGCGGAGCGATCATCCCCTTCTGGTAAAACAGCa                                                                >  1:733052/1‑148 (MQ=255)
                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGAAAAAATTACACACAGTAGCCCCCACCTCACAAAGCGGGGGGCGCCCACCCAGCTGCGTAGCGAGCATCCCTTGCTGGTAAAACAACACTgg                                                            >  2:728190/1‑149 (MQ=255)
                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCACATCAGAAAACGTGGGGGGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCAc                                                               >  2:369223/1‑146 (MQ=255)
                                                                                       ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                            <  1:78024/149‑1 (MQ=255)
                                                                                             aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACGGGAGCCCCCCCATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGGTAAAACAACACTGGcacc                                                        >  1:130886/1‑147 (MQ=255)
                                                                                              aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                   <  1:645063/151‑1 (MQ=255)
                                                                                                  gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                               <  2:260204/151‑1 (MQ=255)
                                                                                                    gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGt                                             <  2:720782/151‑1 (MQ=255)
                                                                                                             gaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGTAGCCCCCACATCAGAAAGCGGGGGGGGCACACCCAGCCGCGTAGCGGTCATCCCTTGCTGGTAAACCAACACTGGCACCACCTCGCGGGGAtggtgg                                      >  1:237507/1‑149 (MQ=255)
                                                                                                                    gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCa                             <  2:672105/151‑1 (MQ=255)
                                                                                                                     cGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGGAGCCCCCACATCAGAAAGCGGGGGGGGGACACCGCGCTGCGTCGCGAACATCCCTTGCTGATAAAACAGCAATGGCACCCCTTCGCGGGGATGGTGGGTGTGAcc                              >  2:410994/1‑149 (MQ=255)
                                                                                                                             ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                     <  1:145909/150‑1 (MQ=255)
                                                                                                                               gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGAt                   <  1:523710/150‑1 (MQ=255)
                                                                                                                                          gcgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAt       <  2:698693/151‑1 (MQ=255)
                                                                                                                                           cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc      >  2:149361/1‑151 (MQ=255)
                                                                                                                                           cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc      >  2:715475/1‑151 (MQ=255)
                                                                                                                                             cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg    >  2:392224/1‑151 (MQ=255)
                                                                                                                                             cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg    <  1:107008/151‑1 (MQ=255)
                                                                                                                                               gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc  >  2:435023/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCC  >  NC_000913/3508513‑3508806

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