Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,400,385 A→G 6.5% S186S (TCT→TCC 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
*W3110S.gb3,400,3850AG 6.5% 166.3 / 2.7 61S186S (TCT→TCCmreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACC  >  W3110S.gb/3400317‑3400444
                                                                    |                                                           
agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACag                                                            <  1:778980/70‑1 (MQ=255)
agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACag                                                            <  1:2324519/70‑1 (MQ=255)
agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACag                                                            <  1:3769060/70‑1 (MQ=255)
agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACag                                                            <  1:2918028/70‑1 (MQ=255)
agagAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTAGAAACGGTCACCACCAATGCGCACag                                                            <  1:356117/70‑1 (MQ=255)
 gagaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACagaa                                                          <  1:1055331/71‑1 (MQ=255)
 gagaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACagaa                                                          <  1:1388948/71‑1 (MQ=255)
 gagaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACagaa                                                          <  1:2750565/71‑1 (MQ=255)
  agaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACag                                                            >  1:2754919/1‑68 (MQ=255)
   gaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACagaa                                                          >  1:2295291/1‑69 (MQ=255)
   gaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACagaa                                                          >  1:700980/1‑69 (MQ=255)
   gaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACagaa                                                          >  1:1706087/1‑69 (MQ=255)
   gaACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACagaa                                                          >  1:3175564/1‑69 (MQ=255)
     aCCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAagag                                                       <  1:1145505/70‑1 (MQ=255)
         tAATTACGTCGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTa                                                  <  1:681361/71‑1 (MQ=255)
         tAATTACGACGCACATAGTTGGTGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTa                                                  <  1:2240528/71‑1 (MQ=255)
         tAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTa                                                  <  1:3811499/71‑1 (MQ=255)
          aaTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTa                                                  <  1:2249604/70‑1 (MQ=255)
           aTTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTaa                                                 <  1:3845711/70‑1 (MQ=255)
                   gCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCacac                                          >  1:2416108/1‑69 (MQ=255)
                      caTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                     <  1:3184040/71‑1 (MQ=255)
                      caTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                     <  1:3748602/71‑1 (MQ=255)
                             gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAATCCACACCGTTCAAg                                  >  1:235632/1‑67 (MQ=255)
                             gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACcaca                                           >  1:1660073/1‑58 (MQ=255)
                             gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAg                                  >  1:1632030/1‑67 (MQ=255)
                             gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAg                                  >  1:2528259/1‑67 (MQ=255)
                             gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAg                                  >  1:1235310/1‑67 (MQ=255)
                             gatgatAGCTTCGTCGAAACGGTCAACACCAATGCGCACAGAAGAGGAGTAAACcaca                                           >  1:1578757/1‑58 (MQ=255)
                             aatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACcaca                                           >  1:2351955/2‑58 (MQ=255)
                              atgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGAt                             <  1:1047714/71‑1 (MQ=255)
                                 atAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATaa                           <  1:3417272/70‑1 (MQ=255)
                                 atAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATaa                           <  1:2685324/70‑1 (MQ=255)
                                   aGCTTCGTCGAAACGGTCACCACCAATGCGCACagaaga                                                        >  1:2161630/1‑39 (MQ=255)
                                   aGCTTCGTCGAAACGGTCACCACCAATGCGCACagaaga                                                        >  1:1381745/1‑39 (MQ=255)
                                   aGCTTCGTCGAAACGGTCACCACCAATGCGCACagaaga                                                        >  1:1535147/1‑39 (MQ=255)
                                   aGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAg                        >  1:1541922/1‑71 (MQ=255)
                                   aGCTTCGTCGAAACGGTCACCAACAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAg                                  >  1:1600098/1‑61 (MQ=255)
                                   aGCTTAGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCaa                                   >  1:3233377/1‑60 (MQ=255)
                                    gCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGc                       <  1:3250257/71‑1 (MQ=255)
                                     cTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGga                                <  1:166992/61‑1 (MQ=255)
                                       tcgtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATaa                           >  1:1506060/1‑64 (MQ=255)
                                       tcgtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATaa                           >  1:1984272/1‑64 (MQ=255)
                                        cgtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTaaa                                                >  1:1569531/1‑42 (MQ=255)
                                        cgtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTaaa                                                >  1:1792773/1‑42 (MQ=255)
                                         gtctAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                  <  1:2285290/71‑1 (MQ=255)
                                         gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGt                                                   >  1:372682/1‑38 (MQ=255)
                                         gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGt                                                   >  1:3591270/1‑38 (MQ=255)
                                         gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                  <  1:442043/71‑1 (MQ=255)
                                         gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                  <  1:1547087/71‑1 (MQ=255)
                                         gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                  <  1:3663666/71‑1 (MQ=255)
                                         gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                  <  1:1893145/71‑1 (MQ=255)
                                          tcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTc                 <  1:526648/71‑1 (MQ=255)
                                            gAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTc                 <  1:2846724/69‑1 (MQ=255)
                                             aaaCGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                     <  1:3135176/48‑1 (MQ=255)
                                             aaaCGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                     <  1:3154212/48‑1 (MQ=255)
                                             aaaCGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGt              <  1:798059/71‑1 (MQ=255)
                                               aCGGTCACCACCAATGCGCACGGAAGAGGAGTAAACCACGCCGTTCAAGGAGATAACAGCAACTTCAg               >  1:1046731/1‑68 (MQ=255)
                                               aCGGTCACCACCAATGCGCACGGAAGAGGAGTAAACCACGCCGTTCAAGGAGATAACAGCAACTTCAg               >  1:138800/1‑68 (MQ=255)
                                                cGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTgg            >  1:2768360/1‑70 (MQ=255)
                                                  gTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACa                         >  1:3560734/1‑55 (MQ=255)
                                                    caccacCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTacca       >  1:3600277/1‑71 (MQ=255)
                                                       caccaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac         <  1:298636/66‑1 (MQ=255)
                                                        accaATGCGCACAGAAGAGGAGTAAACCACACCGTTCaa                                   >  1:940372/1‑39 (MQ=255)
                                                        accaATGCGCACAGAAGAGGAGTAAACCACACCGTTCaa                                   >  1:1894539/1‑39 (MQ=255)
                                                        accaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAAc                          <  1:1710088/48‑1 (MQ=255)
                                                        accaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAAc                          <  1:2369749/48‑1 (MQ=255)
                                                        accaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAAc                          <  1:1978737/48‑1 (MQ=255)
                                                         ccaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTgg            >  1:2672866/1‑61 (MQ=255)
                                                         ccaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTgg            >  1:2698983/1‑61 (MQ=255)
                                                         ccaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccaccacc  >  1:2039054/1‑71 (MQ=255)
                                                         ccaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccaccacc  <  1:2391604/71‑1 (MQ=255)
                                                         ccaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccaccacc  <  1:2158993/71‑1 (MQ=255)
                                                          caATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCaa                     >  1:237449/1‑51 (MQ=255)
                                                          caATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCaa                     >  1:2890036/1‑51 (MQ=255)
                                                            aTGCGCACGGAAGAGGAGTAAACCACGCCGTTCAAGGAGATa                            <  1:2834357/42‑1 (MQ=38)
                                                            aTGCGCACGGAAGAGGAGTAAACCACGCCGTTCAAGGAGATa                            <  1:1720824/42‑1 (MQ=38)
                                                             tGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGgag                               <  1:3654716/38‑1 (MQ=255)
                                                                    agaagaGGAGTAAACCACACCGTTCAAGGAGATAACAg                        >  1:3046775/1‑38 (MQ=255)
                                                                    agaagaGGAGTAAACCACACCGTTCAAGGAGATAACAg                        >  1:780887/1‑38 (MQ=255)
                                                                    agaagaGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAg               <  1:3729598/47‑1 (MQ=255)
                                                                    |                                                           
AGAGAACCGTAATTACGACGCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACC  >  W3110S.gb/3400317‑3400444

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: