Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,328,272 A→G 100% V236A (GTT→GCT)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,328,2720AG98.3% 96.7 / ‑6.4 30V236A (GTT→GCT) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/12);  total (17/13)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATG  >  REL606/3328135‑3328415
                                                                                                                                         |                                                                                                                                               
gTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCAc                                                                                                                                     >  2:473914/1‑150 (MQ=255)
   ccGGCGGGCACTGTTCCAGTGCGACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGGGCTTCGAGGATTTCATTTGAGTTAAGGGTAAAACCGCGTGGAAGACCTTCTGCCATGTTGCGGCCACGAGCTTCGACTTCACGGa                                                                                                                                  >  1:457321/1‑150 (MQ=255)
     ggcggGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACt                                                                                                                                >  1:109229/1‑150 (MQ=255)
          gCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATc                                                                                                                           >  2:444391/1‑150 (MQ=255)
           cACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                                                                                                                          >  2:254823/1‑150 (MQ=255)
           cACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                                                                                                                          >  2:662129/1‑150 (MQ=255)
             cTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGcc                                                                                                                        >  1:230539/1‑150 (MQ=255)
               gTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCg                                                                                                                      >  2:297118/1‑150 (MQ=255)
               gTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACCAGCTTCGATTTCACGGACTTCATCGCCCg                                                                                                                      >  1:96102/1‑150 (MQ=255)
               gTTCCAGTGCAACAATTACAGCGCACACAATACCGGTCAGCGGTTCCTGCAGTACTTCGAGGATTTCATTGGAGTTCAGGATAACACCGCGTGGAACACCTTCTGCCAGGTTACGGCAACGAGCTTCGATTTCCCGGACTTAATCGCCCg                                                                                                                      >  1:606225/1‑150 (MQ=255)
                           ccATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAAc                                                                                                          >  1:108028/1‑150 (MQ=255)
                                        cacaATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCtt                                                                                             >  2:651622/1‑150 (MQ=255)
                                            aTACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGGGCTTGATa                                                                                         >  1:70772/1‑150 (MQ=255)
                                                 ggTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTc                                                                                    <  2:93605/150‑1 (MQ=255)
                                                      gCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCgg                                                                               <  2:234196/150‑1 (MQ=255)
                                                                cAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCAcc                                                                     <  2:109551/150‑1 (MQ=255)
                                                                cAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCAcc                                                                     <  2:315114/150‑1 (MQ=255)
                                                                           cgcttgcATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTGACGGCCACGAACTTCGACTGCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAAc                                                          <  1:559239/144‑1 (MQ=255)
                                                                              tttCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGt                                                       >  2:343144/1‑150 (MQ=255)
                                                                                             gggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACAt                                        >  1:491162/1‑150 (MQ=255)
                                                                                             gggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGACCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACAt                                        >  1:600510/1‑150 (MQ=255)
                                                                                             gggTAAAACCGCGGGGAACACCTTCTGCCAGGTTACGGCCACGGGCTTCGATTTCACGCACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACAt                                        >  2:400525/1‑150 (MQ=255)
                                                                                                 aaaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGtt                                    <  2:150163/150‑1 (MQ=255)
                                                                                                   aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTTATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTtga                                  <  1:404267/150‑1 (MQ=255)
                                                                                                    aCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTtgat                                 >  1:444317/1‑150 (MQ=255)
                                                                                                       gcgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTgatgat                              <  1:472570/150‑1 (MQ=255)
                                                                                                       gcgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTgatgat                              <  2:96389/150‑1 (MQ=255)
                                                                                                       gcgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTGCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAgatgatgat                              <  1:239818/150‑1 (MQ=255)
                                                                                                       gcgTGGAACACCTTCTGCCAGGTTACGACCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTgatgat                              <  2:602277/150‑1 (MQ=255)
                                                                                                                           ggTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTcac          >  2:523661/1‑150 (MQ=255)
                                                                                                                               aCGGCCACGAGCGTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCaccacc      <  1:522160/150‑1 (MQ=255)
                                                                                                                                 ggCCCCGAGCTTCTCTTTCACGGACTTCATCGCCCAGATAAGCCGAAGCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGTTTCGTCGAAACGGTCACCCCCaa    <  1:400322/150‑1 (MQ=255)
                                                                                                                                   ccACGGGCTTCGATTTCACGGACTTCAGCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGGTTCGCCGGTCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCaccccctacg  <  1:40044/150‑5 (MQ=255)
                                                                                                                                   ccACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATg  <  2:492557/150‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                               
GTTCCGGCGGGCACTGTTCCAGTGCAACCATTACCGCGCTCACAATACCGGTCAGCGGTTCCTGCAGTGCTTCGAGGATTTCATTGGAGTTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCGTGCTTGATACGTTCTGCGGTGGCTTCACCGATCAGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATG  >  REL606/3328135‑3328415

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