Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,514,494 A→G 13.1% V181A (GTC→GCC)  ftsK ← DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,514,4940AG13.1% 91.7 / 4.1 37V181A (GTC→GCC) ftsKDNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base A (17/15);  new base G (2/3);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

CGTCCTCGCGCAGCTGTGCCTCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGC  >  NC_002947/4514348‑4514640
                                                                                                                                                  |                                                                                                                                                  
cGTCCTCGCGCAGCGGTGCCTCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACgg                                                                                                                                                 <  1:377483/150‑1 (MQ=255)
     tCGCGCAGCTGTGCCTCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAGCCACGACAGGTCGGTGAAGACGGTCAgg                                                                                                                                            <  1:429361/150‑1 (MQ=255)
         gcAGCTGTGCCTCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCga                                                                                                                                        >  2:448469/1‑150 (MQ=255)
                   cTCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGCTGAAGACGGTCAGGCCcaacaggaacag                                                                                                                              >  2:374966/1‑150 (MQ=255)
                    tCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCgaaca                                                                                                                                     <  1:731641/142‑1 (MQ=255)
                    tCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCgaaca                                                                                                                                     >  2:731641/1‑142 (MQ=255)
                                             cTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGGCGGGGGTCATGGCTATCGTGAACAGGGCAATGAACATCAGCGTGCTGCCCt                                                                                                    <  1:575271/150‑1 (MQ=255)
                                                               cAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTgc                                                                                  >  1:416031/1‑150 (MQ=255)
                                                                aGCCCCTTGCACCAGTTCGAGCAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCa                                                                                                                >  2:573237/1‑119 (MQ=255)
                                                                aGCCCCTTGCACCAGTTCGAGCAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCa                                                                                                                <  1:573237/119‑1 (MQ=255)
                                                                            cAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGGCGGTCAGGc                                                                                                                                           >  1:650694/1‑80 (MQ=255)
                                                                            cAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGGCGGTCAGGc                                                                                                                                           <  2:650694/80‑1 (MQ=255)
                                                                            cAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGGCGGTCAGGc                                                                                                                                           <  2:650655/80‑1 (MQ=255)
                                                                            cAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGGCGGTCAGGc                                                                                                                                           >  1:650655/1‑80 (MQ=255)
                                                                                  gAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTc                                                                                                                                               >  2:313868/1‑70 (MQ=255)
                                                                                  gAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTc                                                                                                                                               <  1:313868/70‑1 (MQ=255)
                                                                                  gAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTcagcag                                                                                       <  1:239277/126‑1 (MQ=255)
                                                                                  gAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTcagcag                                                                                       >  2:239277/1‑126 (MQ=255)
                                                                                   aaCAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTcag                                                                                          >  1:44428/1‑122 (MQ=255)
                                                                                   aaCAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTcag                                                                                          <  2:44428/122‑1 (MQ=255)
                                                                                      aGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCg                                                           <  2:416031/150‑1 (MQ=255)
                                                                                        gTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCa                                                                                                                                              <  2:749332/65‑1 (MQ=255)
                                                                                        gTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCa                                                                                                                                              >  1:749332/1‑65 (MQ=255)
                                                                                           gAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAg                                                      >  2:176002/1‑150 (MQ=255)
                                                                                             gggTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTg                                                    <  1:448469/150‑1 (MQ=255)
                                                                                                            tCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTca                                                                                           >  2:350976/1‑96 (MQ=255)
                                                                                                            tCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTca                                                                                           <  1:350976/96‑1 (MQ=255)
                                                                                                              aGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAAc                                   >  2:705086/1‑150 (MQ=255)
                                                                                                               ggTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTcc                                                >  2:518763/1‑136 (MQ=255)
                                                                                                               ggTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTcc                                                <  1:518763/136‑1 (MQ=255)
                                                                                                                gTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAggg                                                                                                                            <  2:245993/59‑1 (MQ=255)
                                                                                                                gTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAggg                                                                                                                            >  1:245993/1‑59 (MQ=255)
                                                                                                                gTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCa                                                                            <  2:744357/107‑1 (MQ=255)
                                                                                                                gTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAAACGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCa                                                                            >  1:744357/1‑107 (MQ=255)
                                                                                                                     tCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTc                                                 >  1:593955/1‑129 (MQ=255)
                                                                                                                     tCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTc                                                 <  2:593955/129‑1 (MQ=255)
                                                                                                                                               aaGACGGTCAGGCCGAACAGGAACAGGGCAATGTACGTCAGCGTGCTGCCCGTCGCGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATgc  >  1:447552/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTCCTCGCGCAGCTGTGCCTCAAGGCGCTTGCGCTCGTTGCGCGCTTCCCACCAGCGGGTGGCAGCCCCTTGCACCAGTTCGAACAGGTCGAGGGTGATCTTGCCGGTCAGGTCCATCACCTTGAACCACGACAGGTCGGTGAAGACGGTCAGGCCGAACAGGAACAGGGCAATGAACATCAGCGTGCTGCCCTGCACGTTCAGCAGGTTGCGCGCCAGGTCGCCCAGGCTTTCGCCCAGTGCTCCACCGGCGGAGAACGGCAGGCTGGCGGGCGGGTGGAAGTGGATATGC  >  NC_002947/4514348‑4514640

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