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.2% 232.5 / 1.8 62S10P (TCC→CCC) mreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (33/26);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCG  >  REL606/3328878‑3329021
                                                                         |                                                                      
ttCAATACGATGCCTTGTCCTTTTACATAAATGAGTGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa                                                                      >  1:1329226/1‑76 (MQ=255)
ttCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGgaa                                                                      <  2:424100/76‑1 (MQ=255)
    aTACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACAt                                                                  >  2:1050395/1‑76 (MQ=255)
    aTACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACAt                                                                  <  1:1813140/76‑1 (MQ=255)
     tACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTTATTGGGAAACagg                                                                 >  2:836558/1‑74 (MQ=255)
     tACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATg                                                                 <  2:501404/76‑1 (MQ=255)
      aCGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGc                                                                <  1:1503930/76‑1 (MQ=255)
         aTGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCAc                                                             <  1:974442/76‑1 (MQ=255)
         aTGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCAc                                                             >  1:423901/1‑76 (MQ=255)
           gCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGa                                                           >  2:1779596/1‑76 (MQ=255)
             cTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGaaa                                                         >  1:1058261/1‑76 (MQ=255)
               tGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAtt                                                       >  2:980977/1‑76 (MQ=255)
               tGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAtt                                                       >  1:767100/1‑76 (MQ=255)
                gTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAttt                                                      >  1:581802/1‑76 (MQ=255)
                gTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAAttt                                                      >  1:517079/1‑76 (MQ=255)
                 tCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAAtttt                                                     <  1:97889/76‑1 (MQ=255)
                    ttttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTc                                                  <  1:1042830/76‑1 (MQ=255)
                     tttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCa                                                 <  2:1545278/76‑1 (MQ=255)
                     tttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCa                                                 <  1:1602444/76‑1 (MQ=255)
                     tttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCa                                                 <  2:1580660/76‑1 (MQ=255)
                     tttACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCa                                                 <  1:634570/76‑1 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               >  2:845035/1‑76 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               <  1:1356793/76‑1 (MQ=255)
                       tACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAAc                                               <  2:672119/76‑1 (MQ=255)
                         cATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACAt                                             >  1:371214/1‑76 (MQ=255)
                          aTAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATa                                            <  2:940848/76‑1 (MQ=255)
                            aaaTGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACt                                          >  1:1729462/1‑76 (MQ=255)
                               tGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGCAATTTTTTCAACATACTAAg                                       >  2:1066527/1‑76 (MQ=255)
                                gAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAgg                                      <  1:265464/76‑1 (MQ=255)
                                gAGGGTATTCGCAGTACCCAGGTCATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAggg                                     >  2:1441551/1‑76 (MQ=255)
                                 aGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAggg                                     <  1:718404/76‑1 (MQ=255)
                                 aGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAggg                                     >  1:948300/1‑76 (MQ=255)
                                 aGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAggg                                     >  1:1563078/1‑76 (MQ=255)
                                 aGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAggg                                     >  2:1183565/1‑76 (MQ=255)
                                   ggTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGAt                                   <  2:51300/76‑1 (MQ=255)
                                    gTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATa                                  >  1:651746/1‑76 (MQ=255)
                                      aTTCGCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAAt                                >  1:906320/1‑76 (MQ=255)
                                          gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg                            <  1:1241772/76‑1 (MQ=255)
                                          gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg                            >  1:1288108/1‑76 (MQ=255)
                                          gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg                            >  2:713969/1‑76 (MQ=255)
                                          gCAGTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTg                            <  1:297959/76‑1 (MQ=255)
                                             gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         <  1:1353139/76‑1 (MQ=255)
                                             gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         >  1:553434/1‑76 (MQ=255)
                                             gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         <  1:1368475/76‑1 (MQ=255)
                                             gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         >  2:188349/1‑76 (MQ=255)
                                             gTACCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGaaa                         >  1:1259374/1‑76 (MQ=255)
                                               aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc                       >  1:110329/1‑76 (MQ=255)
                                               aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc                       >  1:1378535/1‑76 (MQ=255)
                                               aCCCAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGc                       >  1:827987/1‑76 (MQ=255)
                                                 ccAGGTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTg                     <  2:51884/76‑1 (MQ=255)
                                                    ggTCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTgggg                  <  2:1668093/76‑1 (MQ=255)
                                                      tCAATGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCg                >  2:687323/1‑76 (MQ=255)
                                                        aaTGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGga              >  1:410308/1‑76 (MQ=255)
                                                          tGGACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaa            <  2:3428/76‑1 (MQ=255)
                                                           tgACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCTgaaca           >  2:1829690/2‑74 (MQ=255)
                                                           tgACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaa           >  2:1348611/2‑76 (MQ=255)
                                                           ggACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACACACTAAGGGATAATCCTGTAATCTGGGGCGgaaaa           >  2:527321/1‑76 (MQ=255)
                                                            gACAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaag          <  2:1001903/76‑1 (MQ=255)
                                                             aaaaGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaaga         <  1:929804/74‑1 (MQ=255)
                                                             aCAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaaga         <  1:805810/76‑1 (MQ=255)
                                                             aCAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaaga         >  2:624692/1‑76 (MQ=255)
                                                              cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa        <  2:1666901/76‑1 (MQ=255)
                                                              cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa        <  1:327321/76‑1 (MQ=255)
                                                              cAAGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaa        <  1:1015719/76‑1 (MQ=255)
                                                               aaGTCATTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGgaaaagaaa       >  2:102659/1‑76 (MQ=255)
                                                                    aTTGGGAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCg  >  2:1428266/1‑76 (MQ=255)
                                                                         |                                                                      
TTCAATACGATGCCTTGTCCTTTTACATAAATGAGGGTATTCGCAGTACCCAGGTCAATGGACAAGTCATTGGAAAACATGCCACGAAATTTTTTCAACATACTAAGGGATAATCCTGAAAGCTGGGGCGGAAAAGAAAATCCG  >  REL606/3328878‑3329021

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