Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,075,492 C→G 7.3% T70T (ACC→ACG)  mreB → cell wall structural actin‑like protein MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,075,4920CG 7.3% 160.6 / 2.5 54T70T (ACC→ACG) mreBcell wall structural actin‑like protein MreB
Reads supporting (aligned to +/- strand):  ref base C (22/28);  new base G (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

CCCAGGGTCACTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTG  >  NC_002947/1075347‑1075630
                                                                                                                                                 |                                                                                                                                          
cccAGGGTCACTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCAt                                                                                                                                         <  1:133890/149‑1 (MQ=255)
 ccAGGGTCACTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATg                                                                                                                                        <  6:141856/149‑1 (MQ=255)
 ccAGGGTCACTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATg                                                                                                                                        <  5:125002/149‑1 (MQ=255)
   aGGGTCACTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATgg                                                                                                                                       >  4:301353/1‑148 (MQ=255)
         aCTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATGGGTACGGATGGCAAcaac                                                                                                                                <  7:164890/149‑4 (MQ=255)
          cTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCTTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAg                                                                                                                                >  4:171304/1‑148 (MQ=255)
              cccATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGaaaa                                                                                                                            <  1:115066/148‑1 (MQ=255)
               ccATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGc                                                                                                                          <  6:281312/149‑1 (MQ=255)
                 aTGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGcgt                                                                                                                        >  6:158284/1‑149 (MQ=255)
                       aaGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGc                                                                                                                  <  3:301353/149‑1 (MQ=255)
                       aaGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGc                                                                                                                  >  7:253336/1‑149 (MQ=255)
                       aaGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGc                                                                                                                  >  1:5458/1‑149 (MQ=255)
                                    gCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGa                                                                                                      <  5:70865/148‑1 (MQ=255)
                                      aTGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGc                                                                                                   >  1:154766/1‑149 (MQ=255)
                                                      tttCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGc                                                                                   >  7:74234/1‑149 (MQ=255)
                                                        tCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTg                                                                                      <  5:234525/144‑1 (MQ=255)
                                                            tCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACg                                                                             >  2:40539/1‑149 (MQ=255)
                                                                   gggTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGc                                                                       >  6:21091/1‑148 (MQ=255)
                                                                    ggTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGccc                                                                          <  7:20711/144‑1 (MQ=255)
                                                                    ggTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCCGGG‑CGTACGgcccagcat                                                                    <  5:23742/149‑7 (MQ=255)
                                                                      tACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACa                                                                   <  3:171304/149‑1 (MQ=255)
                                                                        cTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACAtt                                                                 <  4:162938/149‑1 (MQ=255)
                                                                            cAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAg                                                             <  8:74234/149‑1 (MQ=255)
                                                                             aaCACCCTTATTTACGTCCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGc                                                            >  5:211025/1‑149 (MQ=255)
                                                                             aaCACACTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAAc                                                                                              <  2:40353/115‑1 (MQ=255)
                                                                             aaCACACTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAAc                                                                                              >  1:40353/1‑115 (MQ=255)
                                                                                  ccTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAg                                                             >  4:197923/1‑143 (MQ=255)
                                                                                  ccTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAg                                                             <  3:197923/143‑1 (MQ=255)
                                                                                  ccTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTc                                                       >  5:186186/1‑149 (MQ=255)
                                                                                  ccTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTc                                                       >  3:286944/1‑149 (MQ=255)
                                                                                    ttATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATg                                                                                                                                        >  7:114497/1‑66 (MQ=255)
                                                                                    ttATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATg                                                                                                                                        <  8:114497/66‑1 (MQ=255)
                                                                                      aTTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTc                                                    <  3:288815/148‑1 (MQ=255)
                                                                                          aCGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGAt                                                <  8:77762/148‑1 (MQ=255)
                                                                                          aCGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGAt                                                <  2:5458/148‑1 (MQ=255)
                                                                                            gtgcgtgAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAAc                                                                                              <  3:150472/100‑1 (MQ=255)
                                                                                            gtgcgtgAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAAc                                                                                              >  4:150472/1‑100 (MQ=255)
                                                                                            gtgcgtgAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGaa                                             <  5:158284/149‑1 (MQ=255)
                                                                                                    gcgcGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCg                                     >  1:160405/1‑149 (MQ=255)
                                                                                                        ggTATCGTCCAGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCAt                                 >  6:103267/1‑149 (MQ=255)
                                                                                                         gTATCGTCCTGAATGAGCCCTCGGTTGTTTCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCAtt                                <  8:253336/149‑1 (MQ=255)
                                                                                                         gTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCAtt                                <  5:21091/149‑1 (MQ=255)
                                                                                                               tCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGAc                          >  8:200343/1‑149 (MQ=255)
                                                                                                                     aTGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGc                    <  4:217870/149‑1 (MQ=255)
                                                                                                                          cccTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTg               <  2:154766/149‑1 (MQ=255)
                                                                                                                             tCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGACCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGt                                    <  7:47468/125‑1 (MQ=255)
                                                                                                                             tCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAACGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAATGCCGGCGt                                    >  8:47468/1‑125 (MQ=255)
                                                                                                                               tctggttgCCATGGGTACGGATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCTATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGc              <  4:206089/141‑1 (MQ=255)
                                                                                                                               tctggttgCCATGGGTACGGATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCAATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGc              >  3:206089/5‑145 (MQ=255)
                                                                                                                                   gttgCCATGGGTACGGATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAAt      >  5:207728/1‑149 (MQ=255)
                                                                                                                                       ccATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGcc                             >  7:185933/1‑122 (MQ=255)
                                                                                                                                       ccATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGcc                             <  8:185933/122‑1 (MQ=255)
                                                                                                                                       ccATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTg  >  5:115343/1‑149 (MQ=255)
                                                                                                                                       ccATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTg  <  2:160405/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
CCCAGGGTCACTTTCCCATGTTCAAGAAACTGCGTGGCATGTTTTCCAGCGATCTTTCTATCGATCTGGGTACTGCCAACACCCTTATTTACGTGCGTGAGCGCGGTATCGTCCTGAATGAGCCCTCGGTTGTTGCCATCCGTACCCATGGCAACCAGAAAAGCGTCGTCGCCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTG  >  NC_002947/1075347‑1075630

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