Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,128,408 T→C 33.4% D155G (GAC→GGC)  yghO ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,128,4080TC33.4% 46.1 / 40.8 46D155G (GAC→GGC) yghOpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (13/17);  new base C (10/5);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTA  >  W3110S.gb/3128342‑3128473
                                                                  |                                                                 
cccAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCaaa                                                               <  1:1776679/71‑1 (MQ=255)
   aGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCTTCATCGCACAGGGTGGTGTGTCAAAAcc                                                            >  1:623287/1‑71 (MQ=255)
   aGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAAcc                                                            >  1:3539717/1‑71 (MQ=255)
   aGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAAcc                                                            >  1:799994/1‑71 (MQ=255)
      tgttCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCtt                                                          >  1:1723879/1‑70 (MQ=255)
      tgttCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCtt                                                          >  1:3813304/1‑70 (MQ=255)
      tgttCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCtt                                                          >  1:320988/1‑70 (MQ=255)
             tatGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAg                                                   >  1:3181967/1‑70 (MQ=255)
                 cgcgGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGgtgtgt                                                                   >  1:3999564/1‑50 (MQ=255)
                 cgcgGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGgtgtgt                                                                   >  1:383929/1‑50 (MQ=255)
                 cgcgGCATACCACGGTTTGCCGTGTGCCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATc                                              <  1:2762090/71‑1 (MQ=255)
                      cATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCtt                                         <  1:541845/71‑1 (MQ=255)
                          ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTcc                                      <  1:304191/70‑1 (MQ=255)
                          ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTcc                                      <  1:2994320/70‑1 (MQ=255)
                          ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                     >  1:2480924/1‑71 (MQ=255)
                          ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCt                                     >  1:3938280/1‑71 (MQ=255)
                           cACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCtg                                    <  1:828963/71‑1 (MQ=255)
                            aCGGTTTGCCGTGTGGCATCATCGCACAGGGTGgtgtgt                                                                   >  1:3526874/1‑39 (MQ=255)
                            aCGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTaa                                                >  1:1940294/1‑58 (MQ=255)
                            aCGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCtga                                   >  1:1646654/1‑71 (MQ=255)
                            aCGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCtga                                   >  1:1384907/1‑71 (MQ=255)
                               gTTTGCCGTGTGGCATCATCG‑‑CAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCtgatt                                 >  1:1053158/1‑68 (MQ=255)
                                   gCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGata                             <  1:3427147/70‑1 (MQ=255)
                                       tgtgGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAg                        <  1:1215144/71‑1 (MQ=255)
                                           gCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCt                     >  1:1596923/1‑70 (MQ=255)
                                             atcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGaa                  <  1:2921785/71‑1 (MQ=255)
                                              tcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                     <  1:2483839/51‑1 (MQ=255)
                                              tcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                     <  1:3980456/51‑1 (MQ=255)
                                              tcatcGCACAGGGTGGTGTGGCAAAACCTTCAATCAGTAATCCGCTTTCCt                                     <  1:1387962/51‑1 (MQ=255)
                                                atcGCAGAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg                >  1:523828/1‑70 (MQ=255)
                                                atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgat                                  >  1:2123484/1‑52 (MQ=255)
                                                atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtga                                   >  1:884721/1‑51 (MQ=255)
                                                atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg                >  1:1863046/1‑70 (MQ=255)
                                                atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg                >  1:3688340/1‑70 (MQ=255)
                                                atcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATc                                              >  1:3857043/1‑40 (MQ=255)
                                                atcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg                >  1:3383146/1‑70 (MQ=255)
                                                atcGCACAGGGTGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtga                                   >  1:3511569/1‑50 (MQ=39)
                                                  cGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGAc             >  1:1855888/1‑71 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCttt                                        <  1:4072750/40‑1 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCttt                                        <  1:405881/40‑1 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCttt                                        <  1:4044453/40‑1 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCttt                                        <  1:2165596/40‑1 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCttt                                        <  1:2506359/40‑1 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCtt                                         <  1:956813/39‑1 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCtt                                         <  1:2312630/39‑1 (MQ=255)
                                                      caGGGTGGTGTGTCAAAACCTTCAATCAGGAATCCGCttt                                        <  1:328031/40‑1 (MQ=255)
                                                       aGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTg        <  1:3302415/71‑1 (MQ=255)
                                                         ggtggtGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                     <  1:2614750/40‑1 (MQ=255)
                                                           tggtgTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATAtt                           <  1:868447/48‑1 (MQ=255)
                                                             gtgtgtCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTa  <  1:381596/71‑1 (MQ=255)
                                                             gtgtgcCAAAACCTTCAATCAGTAATCCGCTTTCCtgattga                               >  1:2131900/1‑42 (MQ=255)
                                                                  tCAAAACCTTCAATCAGTAATACGCTTTCCTGATTGatat                            >  1:2688339/1‑40 (MQ=255)
                                                                  |                                                                 
CCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTA  >  W3110S.gb/3128342‑3128473

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: