Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,512,756 G→A 9.3% G760G (GGC→GGT)  ftsK ← DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,512,7560GA 9.3% 114.9 / 2.9 43G760G (GGC→GGT) ftsKDNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base G (23/16);  new base A (2/2);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

ACTCGATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGTGCACCCGGATCGG  >  NC_002947/4512614‑4512902
                                                                                                                                              |                                                                                                                                                  
aCTCGATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCAccgccgccgc                                                                                                                                             >  2:737495/1‑150 (MQ=255)
    gATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCAccgccgccgccgcc                                                                                                                                         <  1:1061063/150‑1 (MQ=255)
         tGCGCGCGGCACGGTTGTAGCCGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCaa                                                                                                                                    >  2:74530/1‑150 (MQ=255)
                            gCCGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTc                                                                                                                 >  1:592640/1‑150 (MQ=255)
                            gCCGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTc                                                                                                                 <  1:143148/150‑1 (MQ=255)
                              cGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGa                                                                                                               >  2:1330130/1‑150 (MQ=255)
                                              ttGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATgtc                                                                                               <  1:737495/150‑1 (MQ=255)
                                                    ggCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGtt                                                                                         <  2:1240765/150‑1 (MQ=255)
                                                      cAGAAATGGAAGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGt                                                                                       >  2:881699/1‑150 (MQ=255)
                                                       aGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTa                                                                                      <  2:592640/150‑1 (MQ=255)
                                                        gAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAg                                                                                     >  2:492121/1‑150 (MQ=255)
                                                                   gcgACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCACCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGcc                                                                          >  1:1286480/1‑150 (MQ=255)
                                                                           tctcGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCt                                                                  >  1:112602/1‑150 (MQ=255)
                                                                             tcGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTc                                                                >  1:1387761/1‑150 (MQ=255)
                                                                             tcGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCACCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTc                                                                <  2:1286480/150‑1 (MQ=255)
                                                                                aGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCAc                                                             <  1:881699/150‑1 (MQ=255)
                                                                                 gCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACg                                                            >  1:1541196/1‑150 (MQ=255)
                                                                                       aTTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCACCGCCGCCGCCTTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCg                                                                                 <  1:1254864/123‑1 (MQ=255)
                                                                                       aTTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCACCGCCGCCGCCTTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCg                                                                                 >  2:1254864/1‑123 (MQ=255)
                                                                                           cACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTcc                                                                                  <  2:455508/118‑1 (MQ=255)
                                                                                           cACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTcc                                                                                  >  1:455508/1‑118 (MQ=255)
                                                                                            aCCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCAc                                                             <  2:977961/138‑1 (MQ=255)
                                                                                            aCCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCAc                                                             >  1:977961/1‑138 (MQ=255)
                                                                                            aCCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTg                                                 >  2:1451272/1‑150 (MQ=255)
                                                                                             ccGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGc                                                >  1:1111346/1‑150 (MQ=255)
                                                                                                     tcatACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACt                                        >  1:1467900/1‑150 (MQ=255)
                                                                                                     tcatACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCCACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACt                                        >  1:521100/1‑150 (MQ=255)
                                                                                                       atACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGcc                                                                          <  1:900345/114‑1 (MQ=255)
                                                                                                       atACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGcc                                                                          >  2:900345/1‑114 (MQ=255)
                                                                                                        tACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTtcg                                     <  2:1467900/150‑1 (MQ=255)
                                                                                                          cAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTtcgtc                                   <  2:112602/150‑1 (MQ=255)
                                                                                                                gTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGa                             >  2:279719/1‑150 (MQ=255)
                                                                                                                       tCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTa                                                                                      >  1:506915/1‑86 (MQ=255)
                                                                                                                       tCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTa                                                                                      <  2:506915/86‑1 (MQ=255)
                                                                                                                         ggCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATgtcg                                                                                              >  1:820875/1‑76 (MQ=255)
                                                                                                                         ggCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATgtcg                                                                                              <  2:820875/76‑1 (MQ=255)
                                                                                                                         ggCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATgc                    >  2:382406/1‑150 (MQ=255)
                                                                                                                            aTCGTCACCCTCGCCACCGCCGCCGCACCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGcc                 >  1:862999/1‑150 (MQ=255)
                                                                                                                              cGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCATCTTCGACCCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGt               <  2:862999/150‑1 (MQ=255)
                                                                                                                                 cACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGTGCa            <  1:1199382/150‑1 (MQ=255)
                                                                                                                                  aCCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTATCTGCTGGGGCGCCACGCAGCTCCCACACCTCGGCGCTGCTGTGCACTTCATCGTCGGACACGAATGCGCCGTGCAc           >  1:366214/1‑150 (MQ=255)
                                                                                                                                      tCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGC‑TTACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGTGCACCCGGa      >  1:921341/1‑150 (MQ=255)
                                                                                                                                         ccACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGTGCACCCGGATc    >  1:1443772/1‑150 (MQ=255)
                                                                                                                                           aCCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGTGCACCCGGATCgg  <  1:492121/150‑1 (MQ=255)
                                                                                                                                           aCCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGTGCACCCGGATCgg  <  2:564400/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
ACTCGATCATGCGCGCGGCACGGTTGTAGCCGATCTTCAGCTTGCGTTGCACGGCAGAAATGGAGGCGCGACGGCTCTCGAGCACGAATTGCACCGCCTCATCATACAGGGCGTCGGTTTCGGCATCGTCACCCTCGCCACCGCCGCCGCCGCCGTCAAAGCCGCTACCGGCCTCTTCGACGCCGTTGAGGATGTCGTCGTTGTAGTCCGGGGCGCCACGCAGCTTCCACGCCTCGACCGTGCGGTGCACTTCGTCGTCGGACACGAATGCGCCGTGCACCCGGATCGG  >  NC_002947/4512614‑4512902

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑