Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,508,6620TA46.8% 54.5 / 15.8 47E384D (GAA→GAT) yhfWphosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (7/18);  new base A (22/0);  total (29/18)
Fisher's exact test for biased strand distribution p-value = 1.07e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508515‑3508809
                                                                                                                                                     |                                                                                                                                                   
cGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGCCCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTca                                                                                                                                                  <  2:77339/151‑1 (MQ=255)
    gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTcacaca                                                                                                                                              <  2:844632/151‑1 (MQ=255)
        ttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGGGC‑GCGGGAAAAATTACAAACgg                                                                                                                                           >  1:265636/1‑150 (MQ=255)
          aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATCACAAACggga                                                                                                                                         >  2:16784/1‑148 (MQ=255)
          aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTa                                                                                                                                         <  2:826175/150‑1 (MQ=255)
              ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGGAAAAATCACACACCGGAGcccc                                                                                                                                    >  1:111746/1‑151 (MQ=255)
                       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:845213/151‑1 (MQ=255)
                                                                                                                                                  aaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATc  <  1:94050/151‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGG‑TGGCGC‑GCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCATC  >  NC_000913/3508515‑3508809

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