Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,328,951 A→G S10P (TCC→CCC)  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,9510AG95.5% 260.2 / 1.5 67S10P (TCC→CCC) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (40/24);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

TTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCG  >  REL606/3328878‑3329021
                                                                         |                                                                      
ttCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa                                                                      <  2:866807/76‑1 (MQ=255)
ttCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa                                                                      >  1:1156165/1‑76 (MQ=255)
ttCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa                                                                      >  2:320427/1‑76 (MQ=255)
      aCGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGc                                                                >  2:820798/1‑76 (MQ=255)
      aCGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGc                                                                >  2:1612931/1‑76 (MQ=255)
      aCGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGc                                                                >  2:1168427/1‑76 (MQ=255)
       cGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGcc                                                               <  1:489397/76‑1 (MQ=255)
        gATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCa                                                              >  2:373138/1‑76 (MQ=255)
         aTGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCAc                                                             >  1:198797/1‑76 (MQ=255)
          tgtcTTGTTCTTTTAGTTAAATGAGGGTGTTCGCAGTTCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACg                                                            <  2:1387399/73‑1 (MQ=255)
            ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa                                                          >  1:564137/1‑76 (MQ=255)
            ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa                                                          <  1:411470/76‑1 (MQ=255)
            ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa                                                          >  1:969454/1‑76 (MQ=255)
            ccTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaa                                                          <  2:1203979/76‑1 (MQ=255)
             cTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaaa                                                         >  1:1331221/1‑76 (MQ=255)
                gTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAttt                                                      >  1:858630/1‑76 (MQ=255)
                 tCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAtttt                                                     >  1:1337760/1‑76 (MQ=255)
                  ccTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAACttttt                                                    >  1:293831/1‑76 (MQ=255)
                  ccTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAttttt                                                    >  1:293819/1‑76 (MQ=255)
                   cTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAtttttt                                                   >  1:902170/1‑76 (MQ=255)
                     tttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCa                                                 >  2:149792/1‑76 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               >  1:1609562/1‑76 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               >  1:1612778/1‑76 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               >  1:933182/1‑76 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               >  1:752102/1‑76 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               >  1:665169/1‑76 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               >  1:1041827/1‑76 (MQ=255)
                         cATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACac                                             >  2:433383/1‑75 (MQ=255)
                           tAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATAc                                           <  1:1064735/76‑1 (MQ=255)
                             aaTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTa                                         <  1:684000/76‑1 (MQ=255)
                              aTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTaa                                        <  1:178053/76‑1 (MQ=255)
                                gAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAgg                                      >  2:505743/1‑76 (MQ=255)
                                gAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAgg                                      <  1:1591985/76‑1 (MQ=255)
                                 aGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAggg                                     >  2:1101476/1‑76 (MQ=255)
                                  gggTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGa                                    <  2:1421194/76‑1 (MQ=255)
                                    gTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATa                                  <  1:550462/76‑1 (MQ=255)
                                     tATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATaa                                 >  2:1524256/1‑76 (MQ=255)
                                      aTTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAAt                                >  2:170315/1‑76 (MQ=255)
                                      aTTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAAt                                >  2:163756/1‑76 (MQ=255)
                                         cGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCt                             >  2:224528/1‑76 (MQ=255)
                                         cGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCt                             >  1:571906/1‑76 (MQ=255)
                                          gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg                            <  1:1130432/76‑1 (MQ=255)
                                          gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg                            >  1:1374797/1‑76 (MQ=255)
                                           taGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGa                           >  2:1025146/2‑76 (MQ=255)
                                            aGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaa                          <  1:135914/76‑1 (MQ=255)
                                            aGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaa                          <  2:1168382/76‑1 (MQ=255)
                                             gtgcccAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         <  2:1439796/73‑1 (MQ=255)
                                             gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         >  1:1374087/1‑76 (MQ=255)
                                             gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         <  2:1067699/76‑1 (MQ=255)
                                              tACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAg                        >  2:894834/1‑76 (MQ=255)
                                               aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc                       >  1:1474558/1‑76 (MQ=255)
                                               aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc                       >  1:699312/1‑76 (MQ=255)
                                               aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc                       >  2:177639/1‑76 (MQ=255)
                                                cccAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCt                      <  1:44379/76‑1 (MQ=255)
                                                   aGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg                   <  1:1287688/76‑1 (MQ=255)
                                                   aGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg                   <  1:1328818/76‑1 (MQ=255)
                                                   aGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTggg                   <  1:204216/76‑1 (MQ=255)
                                                      tCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCg                <  2:584661/76‑1 (MQ=255)
                                                      tCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCg                >  1:1558427/1‑76 (MQ=255)
                                                       cAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCgg               <  2:1180116/76‑1 (MQ=255)
                                                        aaTGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGga              <  2:19207/76‑1 (MQ=255)
                                                          tGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaa            <  1:189576/76‑1 (MQ=255)
                                                           tgACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaa           <  2:46695/75‑1 (MQ=255)
                                                           ggACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaa           >  1:497844/1‑76 (MQ=255)
                                                             aCAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaaga         <  1:1430807/76‑1 (MQ=255)
                                                              cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa        >  2:1626276/1‑76 (MQ=255)
                                                              cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa        >  2:1614445/1‑76 (MQ=255)
                                                              cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa        >  1:148916/1‑76 (MQ=255)
                                                              cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa        >  2:1260576/1‑76 (MQ=255)
                                                               aaGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaaa       >  2:158316/1‑76 (MQ=255)
                                                                 gTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAAt     <  1:1622236/76‑1 (MQ=255)
                                                                   cATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATcc   >  2:1360436/1‑76 (MQ=255)
                                                                    aTTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCg  <  2:1174774/76‑1 (MQ=255)
                                                                         |                                                                      
TTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCG  >  REL606/3328878‑3329021

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