Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,328,272 A→G 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.7% 282.4 / ‑4.7 76V236A (GTT→GCT) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (29/46);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

TTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCG  >  REL606/3328224‑3328319
                                                |                                               
ttCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCt                                               >  1:1953388/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:1259089/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:6009936/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:3075983/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             <  1:1525370/51‑1 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:3197774/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:2013035/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             <  1:5709653/51‑1 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             <  1:2476202/51‑1 (MQ=255)
    gggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGa                                           <  1:4255198/51‑1 (MQ=255)
     ggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAt                                          <  1:3899928/51‑1 (MQ=255)
     ggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAt                                          >  1:3066418/1‑51 (MQ=255)
     ggTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGAt                                          >  1:5801350/1‑51 (MQ=255)
      gTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAtt                                         <  1:4564830/51‑1 (MQ=255)
      gTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAtt                                         <  1:2174307/51‑1 (MQ=255)
       tAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAttt                                        <  1:2260914/51‑1 (MQ=255)
       tAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAttt                                        >  1:864548/1‑51 (MQ=255)
       tAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAttt                                        >  1:163421/1‑51 (MQ=255)
         aaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCa                                      <  1:4630237/51‑1 (MQ=255)
         aaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCa                                      <  1:6082886/51‑1 (MQ=255)
          aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCAc                                     <  1:4906729/51‑1 (MQ=255)
           aCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACg                                    >  1:2778980/1‑51 (MQ=255)
              gcgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGAc                                 >  1:6629404/1‑51 (MQ=255)
               cgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACt                                <  1:1318058/51‑1 (MQ=255)
               cgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACt                                >  1:7006102/1‑51 (MQ=255)
                gTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACtt                               <  1:1055500/51‑1 (MQ=255)
                gTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACtt                               >  1:4627436/1‑51 (MQ=255)
                gTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACtt                               <  1:3385737/51‑1 (MQ=255)
                 tGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTc                              <  1:3411963/51‑1 (MQ=255)
                 tGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGACTTCACGGACTTc                              <  1:5981704/51‑1 (MQ=38)
                   gAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCAt                            <  1:1300764/51‑1 (MQ=255)
                    aaCACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATc                           <  1:5346489/51‑1 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          >  1:17784/1‑51 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          <  1:6747698/51‑1 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          <  1:310603/51‑1 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          >  1:2776712/1‑51 (MQ=255)
                      cacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          <  1:4168385/50‑1 (MQ=255)
                      cacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGc                         <  1:748511/51‑1 (MQ=255)
                       acCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGcc                        >  1:1886220/1‑51 (MQ=255)
                        ccTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGccc                       <  1:6704048/51‑1 (MQ=255)
                        ccTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGccc                       <  1:3886169/51‑1 (MQ=255)
                         cTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCg                      <  1:2622942/51‑1 (MQ=255)
                          ttCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCgg                     >  1:7137118/1‑51 (MQ=255)
                           tCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGa                    <  1:5905553/51‑1 (MQ=255)
                           tCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCAGGa                    >  1:3645736/1‑51 (MQ=255)
                            cTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGAt                   <  1:218977/51‑1 (MQ=255)
                            cTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGAt                   <  1:480822/51‑1 (MQ=255)
                             tGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATa                  >  1:70957/1‑51 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:3171981/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 >  1:6619885/1‑51 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:2914495/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:6813345/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:1578318/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:6861725/51‑1 (MQ=255)
                               ccAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAg                >  1:410603/1‑51 (MQ=255)
                               ccAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAg                >  1:3534678/1‑51 (MQ=255)
                                cAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGc               <  1:5653157/51‑1 (MQ=255)
                                cAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGc               <  1:500761/51‑1 (MQ=255)
                                 aGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATCAGcc              >  1:2869863/1‑51 (MQ=255)
                                 aGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGcc              >  1:464239/1‑51 (MQ=255)
                                 aGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGcc              <  1:4570946/51‑1 (MQ=255)
                                   gTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGa            >  1:5131376/1‑51 (MQ=255)
                                    ttACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGaa           >  1:5242467/1‑51 (MQ=255)
                                      aCGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAAcc         >  1:193718/1‑51 (MQ=255)
                                       cGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCg        <  1:3555483/51‑1 (MQ=255)
                                        ggCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGa       >  1:528154/1‑51 (MQ=255)
                                        ggCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGa       >  1:6700627/1‑51 (MQ=255)
                                          ccACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAtt     <  1:6187159/51‑1 (MQ=255)
                                          ccACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAtt     <  1:372550/51‑1 (MQ=255)
                                           cACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAttt    <  1:5768121/51‑1 (MQ=255)
                                           cACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAttt    <  1:252211/51‑1 (MQ=255)
                                            aCGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTc   <  1:6442046/51‑1 (MQ=255)
                                            aCGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTc   <  1:721505/51‑1 (MQ=255)
                                            aCGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTc   <  1:498963/51‑1 (MQ=255)
                                            aCGAGCTTCGATTTCACGGACTTCATCGCCCAGATAAGCCGAACCGATTTc   <  1:3071952/51‑1 (MQ=255)
                                             cGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCg  <  1:2040300/51‑1 (MQ=255)
                                                |                                               
TTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCG  >  REL606/3328224‑3328319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: