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,2720AG99.2% 438.2 / ‑4.9 116V236A (GTT→GCT) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (47/68);  total (48/68)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-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.

Only 100 of 116 total aligned reads displayed.

TTCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAACTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGATTTCG  >  REL606/3328224‑3328319
                                                |                                               
ttCAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCt                                               >  1:1953388/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:6009936/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             <  1:2476202/51‑1 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             <  1:1525370/51‑1 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             <  1:5709653/51‑1 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:2013035/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:3075983/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:8984963/1‑51 (MQ=255)
  cAGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTc                                             >  1:1259089/1‑51 (MQ=255)
   aGGGTAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCg                                            <  1:10362507/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)
       tAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAttt                                        <  1:2260914/51‑1 (MQ=255)
       tAAAACCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGAttt                                        >  1:864548/1‑51 (MQ=255)
        aaaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTc                                       <  1:9197355/51‑1 (MQ=255)
         aaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCa                                      <  1:4630237/51‑1 (MQ=255)
         aaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCa                                      <  1:8484096/51‑1 (MQ=255)
         aaaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCa                                      <  1:6082886/51‑1 (MQ=255)
          aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCAc                                     >  1:10296766/1‑51 (MQ=255)
          aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCAc                                     <  1:8329316/51‑1 (MQ=255)
          aaCCGCGTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCAc                                     <  1:4906729/51‑1 (MQ=255)
               cgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACt                                >  1:9757042/1‑51 (MQ=255)
               cgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACt                                <  1:1318058/51‑1 (MQ=255)
               cgTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACt                                >  1:7006102/1‑51 (MQ=255)
                gTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACtt                               <  1:3385737/51‑1 (MQ=255)
                gTGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACtt                               <  1:1055500/51‑1 (MQ=255)
                 tGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTc                              <  1:8231883/51‑1 (MQ=255)
                 tGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTc                              >  1:7844287/1‑51 (MQ=255)
                 tGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTc                              <  1:3411963/51‑1 (MQ=255)
                 tGGAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGACTTCACGGACTTc                              <  1:5981704/51‑1 (MQ=38)
                  ggAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCa                             >  1:10118145/1‑51 (MQ=255)
                   gAACTCCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCAt                            <  1:9867159/51‑1 (MQ=255)
                   gAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCAt                            >  1:10176979/1‑51 (MQ=255)
                   gAACACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCAt                            <  1:1300764/51‑1 (MQ=255)
                    aaCACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATc                           <  1:9493168/51‑1 (MQ=255)
                    aaCACCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATc                           <  1:5346489/51‑1 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          >  1:17784/1‑51 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          >  1:2776712/1‑51 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          <  1:310603/51‑1 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          >  1:7458143/1‑51 (MQ=255)
                     acacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          <  1:6747698/51‑1 (MQ=255)
                      cacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCg                          <  1:4168385/50‑1 (MQ=255)
                      cacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGc                         <  1:748511/51‑1 (MQ=255)
                      cacCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGc                         >  1:9005993/1‑51 (MQ=255)
                       acCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGcc                        >  1:9889688/1‑51 (MQ=255)
                       acCTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGcc                        >  1:1886220/1‑51 (MQ=255)
                        ccTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGccc                       <  1:6704048/51‑1 (MQ=255)
                        ccTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGccc                       <  1:8561055/51‑1 (MQ=255)
                        ccTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGccc                       <  1:7474227/51‑1 (MQ=255)
                         cTTCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCg                      <  1:2622942/51‑1 (MQ=255)
                          ttCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCgg                     >  1:7137118/1‑51 (MQ=255)
                           tCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGa                    <  1:8780830/51‑1 (MQ=255)
                           tCTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCAGGa                    >  1:3645736/1‑51 (MQ=255)
                            cTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGAt                   <  1:218977/51‑1 (MQ=255)
                            cTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGAt                   >  1:10487565/1‑51 (MQ=255)
                            cTGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGAt                   <  1:480822/51‑1 (MQ=255)
                             tGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATa                  <  1:8610887/51‑1 (MQ=255)
                             tGCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATa                  <  1:10445837/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:3171981/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 >  1:6619885/1‑51 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:1578318/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:2914495/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:6861725/51‑1 (MQ=255)
                              gCCAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATaa                 <  1:6813345/51‑1 (MQ=255)
                               ccAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAg                >  1:3534678/1‑51 (MQ=255)
                               ccAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAg                >  1:9603782/1‑51 (MQ=255)
                               ccAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAg                <  1:9328414/51‑1 (MQ=255)
                               ccAGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAg                >  1:410603/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:4570946/51‑1 (MQ=255)
                                 aGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGcc              >  1:9846187/1‑51 (MQ=255)
                                 aGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGcc              <  1:8218006/51‑1 (MQ=255)
                                 aGGTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGcc              >  1:464239/1‑51 (MQ=255)
                                   gTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGa            <  1:9305470/51‑1 (MQ=255)
                                   gTTACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGa            <  1:7360934/51‑1 (MQ=255)
                                    ttACGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGaa           >  1:5242467/1‑51 (MQ=255)
                                      aCGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGa            <  1:8606195/48‑1 (MQ=255)
                                      aCGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAAcc         >  1:193718/1‑51 (MQ=255)
                                       cGGCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCg        <  1:3555483/51‑1 (MQ=255)
                                        ggCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGa       >  1:6700627/1‑51 (MQ=255)
                                        ggCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGa       <  1:8632742/51‑1 (MQ=255)
                                        ggCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGa       >  1:9872844/1‑51 (MQ=255)
                                        ggCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGa       >  1:528154/1‑51 (MQ=255)
                                        ggCCACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGa       <  1:7923984/51‑1 (MQ=255)
                                          ccACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAtt     <  1:6187159/51‑1 (MQ=255)
                                          ccACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAtt     <  1:372550/51‑1 (MQ=255)
                                           cACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAttt    <  1:8617226/51‑1 (MQ=255)
                                           cACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAttt    <  1:5768121/51‑1 (MQ=255)
                                           cACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAttt    <  1:252211/51‑1 (MQ=255)
                                           cACGAGCTTCGATTTCACGGACTTCATCGCCCGGATAAGCCGAACCGAt      >  1:10858182/1‑49 (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 < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: