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

CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508515‑3508804
                                                                                                                                                |                                                                                                                                                   
cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGCCCATTTTGTGGCGGTGGCGCGCGGAAAAATTca                                                                                                                                             <  2:77339/151‑1 (MQ=255)
    gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcacaca                                                                                                                                         <  2:844632/151‑1 (MQ=255)
     tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGGGGgaaaaataaaac                                                                                                                                          >  2:146581/1‑142 (MQ=255)
        ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACAAACgg                                                                                                                                      >  1:265636/1‑150 (MQ=255)
          aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACAAACggga                                                                                                                                    >  2:16784/1‑148 (MQ=255)
          aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTa                                                                                                                                    <  2:826175/150‑1 (MQ=255)
              ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGcccc                                                                                                                               >  1:111746/1‑151 (MQ=255)
                       tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAATTTACACACGGGAGCCCCCCCATc                                                                                                                         >  1:347461/1‑148 (MQ=255)
                       tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCAAACACCGGAGCCCCcac                                                                                                                            >  1:426138/1‑145 (MQ=255)
                            atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATTACACACCGGAGCCCCCACATCAGaaa                                                                                                                    >  1:842922/1‑148 (MQ=255)
                                     cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGCCCCCCCATCAGAAAGCgtgggggg                                                                                                          >  1:138310/1‑145 (MQ=255)
                                           ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                  <  1:133552/151‑1 (MQ=255)
                                              gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGGCGGAAAAATTAAAAACGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACAACCAGc                                                                                                >  1:407449/1‑150 (MQ=255)
                                                gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                             <  1:208270/151‑1 (MQ=255)
                                                  gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAAAAAACGGGGGGGGCACAcc                                                                                                    >  1:961170/1‑142 (MQ=255)
                                                         ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACAAACAGGAGCCCCCCCATCAAAAAACGGGGGGCACAccccccgccgc                                                                                             >  2:27420/1‑135 (MQ=255)
                                                          cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCTGAGGGAGGAAAGAAAAGCGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                                                                   <  2:609003/151‑1 (MQ=255)
                                                          cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGGAAAATAAAACACAGGAGCCCCCCCACCAGAAAGCGGGGTGCGCACACCGAGCTGCGTAGc                                                                                        >  2:528140/1‑146 (MQ=255)
                                                          cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACAATTTTGTGGCGGTGGCGCGCCGAAAAAAAAAAACAAGGTAACCCCCCCACTTGAAAGCCTGAGGCGCTCACCGAGCCGCGTAGCGaac                                                                                    >  2:660885/1‑148 (MQ=255)
                                                                  cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTAAACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATcccct                                                                             >  1:914071/1‑147 (MQ=255)
                                                                    aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAACAAACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACAACGCGCAGCGTAGGGGACATcccctaccg                                                                         >  2:793805/1‑145 (MQ=255)
                                                                         cagaGGGTGTCTCCGGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCGCATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                    <  2:1018377/148‑1 (MQ=255)
                                                                         cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                      <  2:365289/149‑1 (MQ=255)
                                                                          aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                   <  2:246551/151‑1 (MQ=255)
                                                                              ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAATCACACAACGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCCGCCGCGGAACGGTCATCCCCTTCTGGTaaa                                                                    >  1:274431/1‑146 (MQ=255)
                                                                                tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATAACACACGGCAGCCCCCACCTCACAAAGAGGGGTGGGCACACCGCGCTGCGTAGCGATCATCCCCTCCTTATAAAACaac                                                               >  2:632615/1‑147 (MQ=255)
                                                                                 gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                            <  1:770478/151‑1 (MQ=255)
                                                                                     ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCAAACACGGGAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCGACCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc                                                         >  1:884204/1‑150 (MQ=255)
                                                                                        gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                                     <  2:134507/151‑1 (MQ=255)
                                                                                        gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGGGCCAGGAAAAAATACCACAATGTAGCCCCCCAAATAGAAAACCTGGTGCGCACAACGCGCAGCGTACCGAACAACCCCTGCTGATAAACcaccacatggccc                                                      >  2:876029/1‑143 (MQ=255)
                                                                                            aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                 <  2:842826/151‑1 (MQ=255)
                                                                                                 agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                            >  1:201097/1‑151 (MQ=255)
                                                                                                   ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGTGGGCACACCCAGCTGCGTAGCGGGCATCCCTTGCTGATAAAACAGCACTGGCACCACCTCCCgg                                             >  1:695288/1‑148 (MQ=255)
                                                                                                     aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAGCCCCCCCATCAAAAAACGGGGGGCGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACCGGCACCCCTTCCCgggga                                          >  1:749313/1‑147 (MQ=255)
                                                                                                       gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACGGGAACCCCCACATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCCCCTCCCGGGGAtgg                                       >  1:34741/1‑150 (MQ=255)
                                                                                                         aagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTTACACACAGTAGCCCCCCCACCAGAAAACCAGGTGCGCCCCACCAGACGCGCGGCGCTCATCCCCCGCCGATAAACCAACCCCCGCGCCCCCTCGCTGGTGtggtgg                                    >  1:65681/1‑151 (MQ=255)
                                                                                                           gaaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAAAAACACACGGTAGCCCCCCCACCAAAAAAAGCGGGGGGGCCAACCAGCCGCGGTGCGATCATCCCCTGCGTATAAAACAACACAGGcaccccatctccgg                                            >  1:148663/1‑131 (MQ=255)
                                                                                                             aaagaGCGACCATTTTGTGGCGGTTGGGGCGGGAGAAAATTACACACGGTAGACCCCCCAACATAAAGAGAGGGGCGCACACCGAGGTGCGGAGCGATTAACACTCGCTGATAAACAAGCACCGGCCCCACtctgcg                                              >  2:536019/1‑132 (MQ=255)
                                                                                                              aagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCCCGCCGCGGAGCGATCATCCCTTTCTGATAAAACAGCACCGGCACCACCTCGCGGGGATGGGGGGtgtg                               >  2:880255/1‑151 (MQ=255)
                                                                                                                  gCGACCATTTTGTGGCGGTGGGGCGGGGGAAAAAATACACACGGTAACCCCCCCAAATGAAAGCGTGGGGCGCACACCCAGATGCGGAGCGAGCATACCTTGCTGGTAAACCCGAACCGGCACCACTTCGCCGGGATGGTGGCTGTGGcc                            >  2:741814/1‑150 (MQ=255)
                                                                                                                     aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                        <  1:957106/151‑1 (MQ=255)
                                                                                                                           ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                   <  1:457986/150‑1 (MQ=255)
                                                                                                                                      ggcgCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt        <  1:61173/151‑1 (MQ=255)
                                                                                                                                          cgcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc    >  1:781823/1‑151 (MQ=255)
                                                                                                                                             cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg  <  1:976218/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508515‑3508804

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