Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,513,904 A→G 7.8% V363A (GTC→GCC)  ftsK ← DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,513,9040AG 7.8% 136.6 / 2.3 51V363A (GTC→GCC) ftsKDNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base A (25/22);  new base G (2/2);  total (27/24)
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.41e-01

TCCTTGGCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGCTTGCTCGGCT  >  NC_002947/4513760‑4514045
                                                                                                                                                |                                                                                                                                             
tCCTTGGCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGcc                                                                                                                                           >  4:57224/1‑149 (MQ=255)
  cTTGGCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGTGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCg                                                                                                                                          <  1:60310/148‑1 (MQ=255)
      gCCAGGTTGGCGATGCGGCTGACCTTCACGCTAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGcc                                                                                                                                       >  6:33131/1‑147 (MQ=255)
      gCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGGCGTCGGCCag                                                                                                                                     <  5:61633/149‑2 (MQ=255)
      gCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGcc                                                                                                                                       >  1:105723/1‑147 (MQ=255)
      gCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGcc                                                                                                                                       >  7:7369/1‑147 (MQ=255)
      gCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCaa                                                                                                                                     <  7:18936/149‑1 (MQ=255)
        cAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGCCCGACGCCGGCCAAg                                                                                                                                    >  4:62073/1‑148 (MQ=255)
              ggCGATGCGGCTGACCTTCACGCTAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCt                                                                                                                              <  5:33131/148‑1 (MQ=255)
                      ggCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATAc                                                                                                                     >  2:247808/1‑149 (MQ=255)
                           aCCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGAt                                                                                                                >  8:104832/1‑149 (MQ=255)
                                    ccAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTg                                                                                                       >  4:330904/1‑149 (MQ=255)
                                     cAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGt                                                                                                      >  2:345385/1‑149 (MQ=255)
                                      aGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGtt                                                                                                     <  3:330904/149‑1 (MQ=255)
                                           ccggcTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGCATCCCCCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGACGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGGTAGGcc                                                                                                >  8:250537/1‑149 (MQ=255)
                                             ggcTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCgg                                                                                              >  8:80210/1‑149 (MQ=255)
                                              gcTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCggg                                                                                             <  3:57224/149‑1 (MQ=255)
                                              gcTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCggg                                                                                             <  1:247808/149‑1 (MQ=255)
                                                  gATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTc                                                                                           >  7:51468/1‑147 (MQ=255)
                                                        gTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATg                                                                                    <  5:214704/148‑1 (MQ=255)
                                                           aCGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGtt                                                                                                     <  5:41229/128‑1 (MQ=255)
                                                           aCGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGtt                                                                                                     >  6:41229/1‑128 (MQ=255)
                                                                tAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCCCACCGAACTCCGTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAg                                                                           >  7:283153/1‑149 (MQ=255)
                                                                        gggccgggGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGt                                                                   >  5:73734/1‑149 (MQ=255)
                                                                          gtggcggTGGATTCGATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGc                                                                 >  3:40450/6‑149 (MQ=255)
                                                                          gccgggGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGc                                                                 <  6:73734/149‑1 (MQ=255)
                                                                              gggTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTc                                                             >  4:243326/1‑149 (MQ=255)
                                                                              cggTGGATTCGATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTc                                                             >  1:322770/2‑149 (MQ=255)
                                                                                     ttcgaTCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGgc                                                       <  4:40450/144‑1 (MQ=255)
                                                                                     tCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGgcg                                                      <  1:345385/149‑1 (MQ=255)
                                                                                        aaTCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTACCTTCCACGGCGCTg                                                   >  5:105899/1‑149 (MQ=255)
                                                                                          tCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTc                                                 <  2:322770/149‑1 (MQ=255)
                                                                                          tCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCCGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTc                                                 >  1:1097/1‑149 (MQ=255)
                                                                                             aCCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGa                                                                                      <  6:197423/109‑1 (MQ=255)
                                                                                             aCCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGa                                                                                      >  5:197423/1‑109 (MQ=255)
                                                                                             aCCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCCGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTcgac                                              <  2:1097/149‑1 (MQ=255)
                                                                                                        cacaCCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAgcggc                                    <  1:332237/148‑1 (MQ=255)
                                                                                                               aaCTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCtt                            <  7:250537/149‑1 (MQ=255)
                                                                                                                  tCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGgcgc                                  >  5:144327/1‑140 (MQ=255)
                                                                                                                  tCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGgcgc                                  <  6:144327/140‑1 (MQ=255)
                                                                                                                    cTTGAGCTTGATTTCCAGCAGCTGGTCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTcgacga                                            <  8:62901/128‑1 (MQ=255)
                                                                                                                    cTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTcc                        <  5:148877/148‑1 (MQ=255)
                                                                                                                    cTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTcc                        <  2:105723/148‑1 (MQ=255)
                                                                                                                    cTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTcc                        <  8:7369/148‑1 (MQ=255)
                                                                                                                    cTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCt                       >  5:126780/1‑149 (MQ=255)
                                                                                                                               tttCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTACCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGc            <  6:105899/149‑1 (MQ=255)
                                                                                                                                  ccAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGCTTg         >  8:14409/1‑149 (MQ=255)
                                                                                                                                        cttGGCCGGCGTCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGCTTGCTCGGc   >  8:24994/3‑149 (MQ=255)
                                                                                                                                         ttGGCCGGCGTCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGCTTGCTCGGc   >  4:152091/2‑148 (MQ=255)
                                                                                                                                         cTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGCTTGCTCGGCt  >  7:143049/1‑149 (MQ=255)
                                                                                                                                           ggCCGGCGTCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGCTTGCTCGGc   <  3:152091/146‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TCCTTGGCCAGGTTGGCGATGCGGCTGACCTTCACGCCAGCGGCCGGCTGGATTTCGTAACGGGTAATCACCGGGCCGGGGTGGATCGAATCCACCGCCACTTCCACACCGAACTCCTTGAGCTTGATTTCCAGCAGCTGGCCGACGCCGGCCAAGGATTCTGGCGAATACTCGATCTTCTTCTGTTCGGCCGGGTCGAGGATGGAAATCGAGGGCAAGGTGCCTTCCACGGCGCTGTCGACGAACAGCGGCGCCTGCTTTTCCTTCATCGCCCGCTTGCTCGGCT  >  NC_002947/4513760‑4514045

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