Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 3,328,657 G→A 11.6% P108S (CCG→TCG)  mreB ← cell wall structural complex MreBCD, actin‑like component MreB

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6063,328,6570G→A11.6% 6.7 ‑1.0 ‑0.143322mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (21/17):  total (25/18)
Fisher's exact test strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.37e-01

GGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTT  >  REL606/3328622‑3328692
                                   |                                   
ggTCGCGCCAACCGGCACACAAACCAGAACGCGCgg                                     >  1:1795307/1‑36 (MQ=255)
   cgcgCCAACCGGCACACAAACCAGAACGCGCGGGCt                                  >  1:5470941/1‑36 (MQ=255)
    gcgcCAACCGGCACACAAACCAGAACGCGCGGGCtt                                 <  1:5725481/36‑1 (MQ=255)
     cgcCAACCGGCACACAAACCAGAACGCGCGGGCTTg                                >  1:2119324/1‑36 (MQ=255)
      gcCAAACGGCACACAAACCAGAACGCGCGGGCTTgg                               >  1:6082736/1‑36 (MQ=255)
       ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGa                              <  1:4990215/36‑1 (MQ=255)
       ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGa                              >  1:403876/1‑36 (MQ=255)
       ccAACCGGCACACAAACCAGAACGCGCGGGCTTGGa                              >  1:3542606/1‑36 (MQ=255)
         aaCCGGCACACAAACCAGAACGCGCGGGCTTGGACg                            >  1:2489126/1‑36 (MQ=255)
          aCCGGCACACAAACCAGAACGCGCGGGCTTGGACGc                           <  1:5130182/36‑1 (MQ=255)
          aCCGGCACACAAACCAGAACGCGCGAGCTTGGACGc                           >  1:211998/1‑36 (MQ=255)
           ccGGCACACAAACCAGAACGCGCGGGCTTGGACGCa                          <  1:3772938/36‑1 (MQ=255)
           ccGGCACACAAACCAGAACGCGCGGGCTTGGACGCa                          <  1:1022897/36‑1 (MQ=255)
           acGGCACACAAACCAGAACGCGCGGGCTTGGACGCa                          >  1:7439000/2‑36 (MQ=255)
            cGGCACACAAACCAGAACGCGCGAGCTTGGACGCAt                         >  1:1458314/1‑36 (MQ=255)
             ggCACACAAACCAGAACGCGCGGGCTTGGACGCATa                        >  1:6951080/1‑36 (MQ=255)
             ggCACACAAACCAGAACGCGCGAGCTTGGACGCATa                        >  1:3355141/1‑36 (MQ=255)
             ggCACAAAAACCAGAACGCGCGGGCTTGGACGCATa                        >  1:6166301/1‑36 (MQ=255)
               cacacaAACCAGAACGCGCGGGCTTGGACGCATaaa                      >  1:3611465/1‑36 (MQ=255)
               aacacaAACCAGAACGCGCGGGCTTGGACGCATaaa                      <  1:5544331/35‑1 (MQ=255)
                acacaAACCAGAACGCGCGGGCTTGGACGCATAAAg                     >  1:4210173/1‑36 (MQ=255)
                acacaAACCAGAACGCGCGGGCTTGGACGCATAAAg                     >  1:7445058/1‑36 (MQ=255)
                 cacaAACCAGAACGCGCGGGCTTGGACGCATAAAGc                    <  1:5690676/36‑1 (MQ=255)
                  acaAACCAGAACGCGCGGGCTTGGACGCATAAAGCt                   <  1:1198604/36‑1 (MQ=255)
                  acaAACCAGAACGCGCGGGCTTGGACGCATAAAGCt                   <  1:5255795/36‑1 (MQ=255)
                  acaAACCAGAACGCGCGGGCTTGGACGCATAAAGCt                   >  1:5809425/1‑36 (MQ=255)
                   caAACCAGAACGCGCGGGCTTGGACGCATAAAGCTg                  >  1:678080/1‑36 (MQ=255)
                    aaaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGt                 <  1:3160869/36‑1 (MQ=255)
                    aaaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGt                 >  1:1266324/1‑36 (MQ=255)
                     aaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGtt                >  1:1266661/1‑36 (MQ=255)
                     aaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGtt                >  1:3955541/1‑36 (MQ=255)
                     aaCCAGAACGCGCGGGCTTGGACGCATAAAGCTGtt                <  1:1918724/36‑1 (MQ=255)
                     aaCCAGAACGCGCGAGCTTGGACGCATAAAGCTGtt                <  1:6675898/36‑1 (MQ=255)
                        cAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCt             <  1:1192542/36‑1 (MQ=255)
                          gAACGCGCGGGCTTGGACGCATAAAGCTGTTGCtgt           >  1:2340392/1‑36 (MQ=255)
                           aaCGCGCGGGCTTGGACGCATAAAGCTGTTGCtgtg          <  1:4987362/36‑1 (MQ=255)
                            aCGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGc         <  1:441703/36‑1 (MQ=255)
                            aCGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGc         <  1:6892817/36‑1 (MQ=255)
                             cgcgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCa        >  1:528843/1‑36 (MQ=255)
                              gcgcgAGCTTGGACGCATAAAGCTGTTGCTGTGCAc       >  1:2864785/1‑36 (MQ=255)
                               cgcgGGCTTGGACGCATAAAGCTGTTTCTGTGCCCt      <  1:5496209/36‑1 (MQ=37)
                               cgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACt      >  1:3855699/1‑36 (MQ=255)
                               cgcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACt      <  1:2136337/36‑1 (MQ=255)
                                gcgGGCTTGGACGCATAAAGCTGTTGCTGTGCACtt     >  1:5524428/1‑36 (MQ=255)
                                  gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt   >  1:5107378/1‑36 (MQ=255)
                                  gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt   >  1:5028427/1‑36 (MQ=255)
                                  gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt   <  1:6134219/36‑1 (MQ=255)
                                  gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt   <  1:1544599/36‑1 (MQ=255)
                                  gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt   >  1:293831/1‑36 (MQ=255)
                                  gggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGt   <  1:2063131/36‑1 (MQ=255)
                                   ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt  <  1:631009/36‑1 (MQ=255)
                                   ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt  <  1:3728441/36‑1 (MQ=255)
                                   ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt  <  1:2795548/36‑1 (MQ=255)
                                   ggCTTGGACGCATAAAGCTGTTGCTGTGCACTTGtt  <  1:2640500/36‑1 (MQ=255)
                                   |                                   
GGTCGCGCCAACCGGCACACAAACCAGAACGCGCGGGCTTGGACGCATAAAGCTGTTGCTGTGCACTTGTT  >  REL606/3328622‑3328692

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG