Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,400,403 A→G 6.4% G180G (GGT→GGC 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,4030AG 6.4% 165.6 / 2.5 62G180G (GGT→GGCmreBcell wall structural complex MreBCD, actin‑like component MreB
Reads supporting (aligned to +/- strand):  ref base A (30/28);  new base G (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGCTT  >  W3110S.gb/3400336‑3400473
                                                                   |                                                                      
gCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCacac                                                                       >  1:2416108/1‑69 (MQ=255)
   caTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                                                  <  1:3748602/71‑1 (MQ=255)
   caTAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                                                  <  1:3184040/71‑1 (MQ=255)
          gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAATCCACACCGTTCAAg                                                               >  1:235632/1‑67 (MQ=255)
          gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACcaca                                                                        >  1:1660073/1‑58 (MQ=255)
          gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAg                                                               >  1:1235310/1‑67 (MQ=255)
          gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAg                                                               >  1:2528259/1‑67 (MQ=255)
          gatgatAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAg                                                               >  1:1632030/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:2685324/70‑1 (MQ=255)
              atAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATaa                                                        <  1:3417272/70‑1 (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)
                      gtctAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                                               <  1:2285290/71‑1 (MQ=255)
                      gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                                               <  1:1547087/71‑1 (MQ=255)
                      gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                                               <  1:1893145/71‑1 (MQ=255)
                      gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                                               <  1:3663666/71‑1 (MQ=255)
                      gtcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACtt                                               <  1:442043/71‑1 (MQ=255)
                       tcgAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTc                                              <  1:526648/71‑1 (MQ=255)
                         gAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTc                                              <  1:2846724/69‑1 (MQ=255)
                          aaaCGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                                                  <  1:3154212/48‑1 (MQ=255)
                          aaaCGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTc                                                                  <  1:3135176/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:2391604/71‑1 (MQ=255)
                                      ccaATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccaccacc                               >  1:2039054/1‑71 (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:1720824/42‑1 (MQ=38)
                                         aTGCGCACGGAAGAGGAGTAAACCACGCCGTTCAAGGAGATa                                                         <  1:2834357/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)
                                                     gaggagTAAACCACACCGTTCAAGGAGATAACAGCAACTTCa                                             >  1:42818/1‑42 (MQ=255)
                                                     gaggagTAAACCACACCGTTCAAGGAGATAACAGCAACTTCa                                             >  1:256517/1‑42 (MQ=255)
                                                     gaggagTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAg                                            >  1:1769597/1‑43 (MQ=255)
                                                      aggagTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAAccacca                 >  1:1848192/1‑69 (MQ=255)
                                                          gTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAAc           <  1:2291892/71‑1 (MQ=255)
                                                           tAAACCACACCGTTCAAGGAGATAACAGCAACTTCAg                                            <  1:3584584/37‑1 (MQ=255)
                                                           tAAACCACACCGTTCAAGGAGATAACAGCAACTTCAg                                            <  1:3327460/37‑1 (MQ=255)
                                                           tAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccacca                                 >  1:510245/1‑48 (MQ=255)
                                                           tAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccacca                                 >  1:1308319/1‑48 (MQ=255)
                                                           tAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac                                      >  1:2821670/1‑43 (MQ=255)
                                                           tAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTac                                      >  1:1687899/1‑43 (MQ=255)
                                                             aaCCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAACCgg        >  1:2237821/1‑71 (MQ=255)
                                                                cacaCCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCg     >  1:1499068/1‑71 (MQ=255)
                                                                 acacCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccac                                   >  1:243760/1‑40 (MQ=255)
                                                                 acacCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGc    >  1:3462503/1‑71 (MQ=255)
                                                                 acacCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGc    >  1:3189616/1‑71 (MQ=255)
                                                                  cacCGTTCAAGGAGATAACAGCAACTTCAGTGGTaccaccacc                               <  1:121812/43‑1 (MQ=255)
                                                                   acCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGCtt  <  1:2892819/71‑1 (MQ=255)
                                                                   |                                                                      
GCACATAGTTGATGATAGCTTCGTCGAAACGGTCACCACCAATGCGCACAGAAGAGGAGTAAACCACACCGTTCAAGGAGATAACAGCAACTTCAGTGGTACCACCACCGATATCAACCACCATAGAACCGGTCGCTT  >  W3110S.gb/3400336‑3400473

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: