Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,080,607 A→G 14.8% H124R (CAC→CGC)  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,080,6070AG14.8% 61.9 / 5.2 27H124R (CAC→CGC) mreBcell wall structural actin‑like protein MreB
Reads supporting (aligned to +/- strand):  ref base A (9/14);  new base G (2/2);  total (11/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

CCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGGCGCGGGTGCCCGTGAAGTGTTCCTGATCGAAGAGCC  >  NC_002947/1080472‑1080742
                                                                                                                                       |                                                                                                                                       
ccGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCa                                                                                                                                         >  2:38438/1‑136 (MQ=255)
     ggTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACgaga                                                                                                                                    >  1:230705/1‑136 (MQ=255)
      gTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGaa                                                                                                                                   <  1:227650/136‑1 (MQ=255)
                          cTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGcc                                                                                                               <  1:226213/136‑1 (MQ=255)
                                  tACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTg                                                                                                       >  1:372468/1‑136 (MQ=255)
                                     gCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTg                                                                                                    >  1:60526/1‑136 (MQ=255)
                                       ccGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCATGCTGAt                                                                                                  >  2:385504/1‑136 (MQ=255)
                                          gCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCGCGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTg                                                                                               >  1:255343/1‑136 (MQ=255)
                                          gCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCGCGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTg                                                                                               >  1:257780/1‑136 (MQ=255)
                                          gCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTg                                                                                               <  2:444612/136‑1 (MQ=255)
                                           cAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGc                                                                                              >  2:216347/1‑136 (MQ=255)
                                           cAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGc                                                                                              <  1:442256/136‑1 (MQ=255)
                                            aaCATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGcg                                                                                             >  1:502939/1‑136 (MQ=255)
                                                           cGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGa                                                                              <  2:465699/136‑1 (MQ=255)
                                                            gTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGAc                                                                             >  1:247606/1‑136 (MQ=255)
                                                                gATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCGCGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAg                                                                         <  2:257780/136‑1 (MQ=255)
                                                                gATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCGCGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAg                                                                         <  2:255343/136‑1 (MQ=255)
                                                                gATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAg                                                                         <  2:248072/136‑1 (MQ=255)
                                                                  tGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGt                                                                       <  1:216347/136‑1 (MQ=255)
                                                                                      gACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCcgcg                                                   <  2:247606/136‑1 (MQ=255)
                                                                                           cAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCg                                              <  1:118899/136‑1 (MQ=255)
                                                                                           cAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCATGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCg                                              <  1:385504/136‑1 (MQ=255)
                                                                                                      aaaaaaTGTTGCAATATTTTATCAACTAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGgcgc                                   <  1:234044/136‑1 (MQ=255)
                                                                                                      aaaaaaTGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGgcgc                                   >  2:281791/1‑136 (MQ=255)
                                                                                                      aaaaaaTGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGgcgc                                   <  2:101017/136‑1 (MQ=255)
                                                                                                      aaaaaaTGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTAGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGgcgc                                   >  2:192832/1‑136 (MQ=255)
                                                                                                            tGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGGCGCGGGTGc                             <  1:103443/136‑1 (MQ=255)
                                                                                                                 cAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGGCGCGGGTGCCCGTg                        <  2:60526/136‑1 (MQ=255)
                                                                                                                                       aCGAGAAGAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGGCGCGGGTGCCCGTGAAGTGTTCCTGATCGAAGAGcc  <  2:316088/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                       
CCGTCGGTACCGAAGCCAAACGCATGCTGGGCCGTACGCCCGGCAACATTGCAGCCATTCGTCCGATGAAGGACGGCGTCATTGCCGACTTCAGCGTCTGCGAAAAAATGTTGCAATATTTTATCAACAAGGTTCACGAGAACAGCTTCCTGCAGCCGAGCCCGCGCGTGCTGATCTGCGTGCCGTGCAAGTCGACCCAGGTCGAGCGCCGCGCCATCCGCGAGTCGGCCCTCGGCGCGGGTGCCCGTGAAGTGTTCCTGATCGAAGAGCC  >  NC_002947/1080472‑1080742

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