Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6810AC45.2% 42.1 / 12.0 43V378G (GTG→GGG) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (7/16);  new base C (19/0);  total (26/17)
Fisher's exact test for biased strand distribution p-value = 2.01e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCT  >  NC_000913/3508538‑3508831
                                                                                                                                                 |                                                                                                                                                      
tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGC‑GCGGGAAAATTTACACACGGGAGCCCCCCCATc                                                                                                                                                    >  1:347461/1‑148 (MQ=255)
tCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATCAAACACCGGAGCCCCcac                                                                                                                                                       >  1:426138/1‑145 (MQ=255)
     atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAAATTACACACCGGAGCCCCCACATCAGaaa                                                                                                                                               >  1:842922/1‑148 (MQ=255)
              cTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATTACACACCGTAGCCCCCCCATCAGAAAGCgtgggggg                                                                                                                                     >  1:138310/1‑145 (MQ=255)
                    ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                                             <  1:133552/151‑1 (MQ=255)
                       gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGCGGCGGAAAAATT‑AAAAACGGTAGCCCCCCCATCAGAAAACGGGGGGCGCACAACCAGc                                                                                                                           >  1:407449/1‑150 (MQ=255)
                         gggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGc                                                                                                                        <  1:208270/151‑1 (MQ=255)
                           gtgtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGC‑GCGGGAAAAATCACACACGGGAGCCCCCCCATCAAAAAACGGGGGGGGCACAcc                                                                                                                               >  1:961170/1‑142 (MQ=255)
                                  ccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATTACAAACAGGAGCCCCCCCATCAAAAAACGGGGGGCACAccccccgccgc                                                                                                                        >  2:27420/1‑135 (MQ=255)
                                   cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCTGAGGGAGGAAAGAAAAGCGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                                                                                              <  2:609003/151‑1 (MQ=255)
                                   cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGGAAAATAAAACACAGGAGCCCCCCCACCAGAAAGCGGGGTGCGCACACCGAGCTGCGTAGc                                                                                                                   >  2:528140/1‑146 (MQ=255)
                                   cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACAATTTTGTGGCGG‑TGGCGC‑GCCGAAAAAAAAAAACAAGGTAACCCCCCCACTTGAAAGCCTGAGGCGCTCACCGAGCCGCGTAGCGaac                                                                                                               >  2:660885/1‑148 (MQ=255)
                                           cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATTAAACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATcccct                                                                                                        >  1:914071/1‑147 (MQ=255)
                                             aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGC‑GCGGGAAAAATAACAAACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACAACGCGCAGCGTAGGGGACATcccctaccg                                                                                                    >  2:793805/1‑145 (MQ=255)
                                                  cagaGGGTGTCTCCGGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCGCATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaaa                                                                                               <  2:1018377/148‑1 (MQ=255)
                                                  cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                                                                                                 <  2:365289/149‑1 (MQ=255)
                                                   aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                                              <  2:246551/151‑1 (MQ=255)
                                                       ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGG‑GCGGGAAAAATCACACAACGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCCGCCGCGGAACGGTCATCCCCTTCTGGTaaa                                                                                               >  1:274431/1‑146 (MQ=255)
                                                         tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATAACACACGGCAGCCCCCACCTCACAAAGAGGGGTGGGCACACCGCGCTGCGTAGCGATCATCCCCTCCTTATAAAACaac                                                                                          >  2:632615/1‑147 (MQ=255)
                                                          gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                                                                                       <  1:770478/151‑1 (MQ=255)
                                                              ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCAAACACGGGAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCGACCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc                                                                                    >  1:884204/1‑150 (MQ=255)
                                                                 gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca                                                                                <  2:134507/151‑1 (MQ=255)
                                                                 gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑CGGGGC‑CAGGAAAAAATACCACAATGTAGCCCCCCAAATAGAAAACCTGGTGCGCACAACGCGCAGCGTACCGAACAACCCCTGCTGATAAACcaccacatggccc                                                                                 >  2:876029/1‑143 (MQ=255)
                                                                     aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                                                            <  2:842826/151‑1 (MQ=255)
                                                                          agggaggAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAATTCACACACGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                                                       >  1:201097/1‑151 (MQ=255)
                                                                            ggaggaAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGGGGTGGGCACACCCAGCTGCGTAGCGGGCATCCCTTGCTGATAAAACAGCACTGGCACCACCTCCCgg                                                                        >  1:695288/1‑148 (MQ=255)
                                                                              aggaAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATTACACACGGGAGCCCCCCCATCAAAAAACGGGGGGCGCACACCGAGCTGCGGAGCGATCATCCCCTGCTGATAAACCAGCACCGGCACCCCTTCCCgggga                                                                     >  1:749313/1‑147 (MQ=255)
                                                                                gaaagaaaAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATTACACACGGGAACCCCCACATCAGAAAACGTGGGGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCCCCTCCCGGGGAtgg                                                                  >  1:34741/1‑150 (MQ=255)
                                                                                  aagaaaagaGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTTACACACAGTAGCCCCCCCACCAGAAAACCAGGTGCGCCCCACCAGACGCGCGGCGCTCATCCCCCGCCGATAAACCAACCCCCGCGCCCCCTCGCTGGTGtggtgg                                                               >  1:65681/1‑151 (MQ=255)
                                                                                    gaaaagaGCGACCATTTTGTGGCGG‑TGGGGC‑GCGGGAAAAAAAACACACGGTAGCCCCCCCACCAAAAAAAGCGGGGGGGCCAACCAGCCGCGGTGCGATCATCCCCTGCGTATAAAACAACACAGGcaccccatctccgg                                                                       >  1:148663/1‑131 (MQ=255)
                                                                                      aaagaGCGACCATTTTGTGGCGGTTGGGGC‑G‑GGAGAAAATTACACACGGTAGACCCCCCAACATAAAGAGAGGGGCGCACACCGAGGTGCGGAGCGATTAACACTCGCTGATAAACAAGCACCGGCCCCACtctgcg                                                                         >  2:536019/1‑132 (MQ=255)
                                                                                       aagaGCGACCATTTTGTGGCGG‑TGGGGC‑GCGGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCCCGCCGCGGAGCGATCATCCCTTTCTGATAAAACAGCACCGGCACCACCTCGCGGGGATGGGGGGtgtg                                                          >  2:880255/1‑151 (MQ=255)
                                                                                           gCGACCATTTTGTGGCGG‑TGGGGC‑GGGGGAAAAAATACACACGGTAACCCCCCCAAATGAAAGCGTGGGGCGCACACCCAGATGCGGAGCGAGCATACCTTGCTGGTAAACCCGAACCGGCACCACTTCGCCGGGATGGTGGCTGTGGcc                                                       >  2:741814/1‑150 (MQ=255)
                                                                                              aCCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATg                                                   <  1:957106/151‑1 (MQ=255)
                                                                                                    ttGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                                              <  1:457986/150‑1 (MQ=255)
                                                                                                               ggcgC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGt                                   <  1:61173/151‑1 (MQ=255)
                                                                                                                   c‑gcgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc                               >  1:781823/1‑151 (MQ=255)
                                                                                                                      cgGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg                             <  1:976218/151‑1 (MQ=255)
                                                                                                                          aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt                         >  1:660019/1‑151 (MQ=255)
                                                                                                                           aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc                        <  1:94050/151‑1 (MQ=255)
                                                                                                                           aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc                        <  1:845213/151‑1 (MQ=255)
                                                                                                                                           gCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCGGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAatcat        >  2:970418/1‑151 (MQ=255)
                                                                                                                                            cccccACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAatcatc       <  2:529525/151‑1 (MQ=255)
                                                                                                                                             ccccACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCt      <  1:316014/151‑1 (MQ=255)
                                                                                                                                                 acaTCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCt  <  1:880355/151‑1 (MQ=255)
                                                                                                                                                 acaTCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGACACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGAGCCGCCATCACGACCAGGCAATCATCTGGCt  <  2:953138/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
TCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATCACGACCAGGCAATCATCTGGCT  >  NC_000913/3508538‑3508831

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: ‑