Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,128,408 T→C 25.7% 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,4080TC25.7% 122.1 / 49.5 74D155G (GAC→GGC) yghOpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (21/34);  new base C (5/14);  total (26/48)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.55e-01

TAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGC  >  W3110S.gb/3128339‑3128478
                                                                     |                                                                      
tAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTc                                                                       <  1:2546889/71‑1 (MQ=255)
  aCCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCaa                                                                     <  1:896080/71‑1 (MQ=255)
  aCCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCaa                                                                     <  1:1644942/71‑1 (MQ=255)
  aCCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCa                                                                      >  1:447401/1‑70 (MQ=255)
    ccAGTTATTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCaaaa                                                                   >  1:3140074/1‑71 (MQ=255)
          gttCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTc                                                              >  1:2716371/1‑70 (MQ=255)
          gttCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTc                                                              >  1:3042197/1‑70 (MQ=255)
          gttCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTc                                                              >  1:1511797/1‑70 (MQ=255)
            tCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           >  1:454333/1‑71 (MQ=255)
             cAATATGCGCGGCATTCCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           <  1:2396814/70‑1 (MQ=255)
             cAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           >  1:1410784/1‑70 (MQ=255)
             cAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           >  1:3026025/1‑70 (MQ=255)
             cAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           >  1:1010685/1‑70 (MQ=255)
             cAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           >  1:1170702/1‑70 (MQ=255)
               atatGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCaaa                                                                    <  1:3021170/59‑1 (MQ=255)
               atatGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAg                                                        <  1:1547825/71‑1 (MQ=255)
               atatGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAg                                                        <  1:3134353/71‑1 (MQ=255)
                  tGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTa                                                      >  1:2877669/1‑70 (MQ=255)
                     gcgGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATcc                                                  <  1:268364/71‑1 (MQ=255)
                       ggCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAACCCTTCAATCAGTAATCCGc                                                <  1:1847164/71‑1 (MQ=255)
                       ggCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGc                                                <  1:1423599/71‑1 (MQ=255)
                           tACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTc                                                              <  1:3416303/53‑1 (MQ=255)
                           tACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           <  1:2267246/56‑1 (MQ=255)
                           tACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           <  1:1253135/56‑1 (MQ=255)
                             ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                          <  1:4606467/71‑1 (MQ=255)
                             ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                          >  1:3801991/1‑71 (MQ=255)
                             ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                          >  1:1260114/1‑71 (MQ=255)
                             ccACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                          >  1:3928427/1‑71 (MQ=255)
                               aCGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           >  1:2590085/1‑52 (MQ=255)
                                 ggTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAg                                                        <  1:486226/53‑1 (MQ=255)
                                 ggTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAg                                                        <  1:605958/53‑1 (MQ=255)
                                 ggTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAg                                                        <  1:256894/53‑1 (MQ=255)
                                 ggTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgatt                                      >  1:316068/1‑71 (MQ=255)
                                 gctttGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAg                                                        <  1:2223768/51‑1 (MQ=255)
                                  gTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgattg                                     <  1:4591705/71‑1 (MQ=255)
                                  gTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgatt                                      >  1:4442749/1‑70 (MQ=255)
                                  gTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCt                                          <  1:3857442/66‑1 (MQ=255)
                                      gCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGata                                  <  1:3258283/70‑1 (MQ=255)
                                       ccGTGTGGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATAtt                                <  1:223352/71‑1 (MQ=255)
                                         gtgtGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgattga                                    <  1:536490/65‑1 (MQ=255)
                                          tgtgGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAAt                                                           <  1:2117457/41‑1 (MQ=255)
                                            tgGCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCtgattga                                    >  1:619582/1‑62 (MQ=255)
                                              gCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtga                                        >  1:2553946/1‑56 (MQ=255)
                                              gCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCt                          >  1:827862/1‑70 (MQ=255)
                                              gCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTg                         <  1:422615/71‑1 (MQ=255)
                                              gCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTg                         <  1:3286354/71‑1 (MQ=255)
                                              gCATCATCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTg                         <  1:2172188/71‑1 (MQ=255)
                                               catcatCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGa                        <  1:2528472/71‑1 (MQ=255)
                                               catcatCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGa                        <  1:1992104/71‑1 (MQ=255)
                                               catcatCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGa                        <  1:3239354/71‑1 (MQ=255)
                                               catcatCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGa                        <  1:2279645/71‑1 (MQ=255)
                                               catcatCGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGa                        <  1:947044/71‑1 (MQ=255)
                                                atcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGaa                       <  1:1652881/71‑1 (MQ=255)
                                                 tcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTaa                                                     <  1:909812/40‑1 (MQ=255)
                                                 tcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTa                                                      <  1:2494456/39‑1 (MQ=255)
                                                 tcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCt                                               <  1:1390788/46‑1 (MQ=255)
                                                 tcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGaa                       <  1:1200295/70‑1 (MQ=255)
                                                 tcatcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGaa                       <  1:4439608/70‑1 (MQ=255)
                                                  catcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgattga                                    <  1:2388035/56‑1 (MQ=255)
                                                   atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgat                                       <  1:4067205/52‑1 (MQ=255)
                                                   atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCtgat                                       <  1:4106111/52‑1 (MQ=255)
                                                   atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGata                                  >  1:2815054/1‑57 (MQ=255)
                                                   atcGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg                     >  1:4006463/1‑70 (MQ=255)
                                                   atcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTa                                                      <  1:3163705/37‑1 (MQ=255)
                                                   atcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTa                                                      <  1:2053007/37‑1 (MQ=255)
                                                   atcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCtgattga                                    >  1:1333370/1‑55 (MQ=255)
                                                   atcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATAtt                                >  1:4082219/1‑59 (MQ=255)
                                                   atcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg                     <  1:4407743/70‑1 (MQ=255)
                                                   atcGCACAGGGGGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAgg                    <  1:3357568/71‑1 (MQ=255)
                                                    tcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTaa                                                     <  1:481859/37‑1 (MQ=255)
                                                    tcGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGa                        <  1:736359/66‑1 (MQ=255)
                                                     cGCACAGGGTGGTGTGCCAAAACCTTCAATCAGTAAt                                                    >  1:4447979/1‑37 (MQ=255)
                                                          aGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTg             <  1:995930/71‑1 (MQ=255)
                                                            ggtggtGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGAt           >  1:4334674/1‑71 (MQ=255)
                                                                gtgtgtCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGaa                       <  1:3402387/55‑1 (MQ=255)
                                                                 tgtgtCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAgg                    <  1:3450088/57‑1 (MQ=255)
                                                                 tgtgtCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTAc      <  1:3082554/71‑1 (MQ=255)
                                                                   tgtCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAg                     <  1:2834470/54‑1 (MQ=255)
                                                                     tCAAAACCTTCAATCAGTAATCCGCTTTCCtgattga                                    <  1:837537/37‑1 (MQ=255)
                                                                     tCAAAACCTTCAATCAGTAATCCGCTTTCCtgatt                                      >  1:209960/1‑35 (MQ=255)
                                                                     tCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACtt    >  1:739633/1‑69 (MQ=255)
                                                                     tCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTg   >  1:333334/1‑70 (MQ=255)
                                                                     tCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGc  <  1:1843229/71‑1 (MQ=255)
                                                                     |                                                                      
TAACCCAGTTGTTCAATATGCGCGGCATACCACGGTTTGCCGTGTGGCATCATCGCACAGGGTGGTGTGTCAAAACCTTCAATCAGTAATCCGCTTTCCTGATTGATATTCAGGCTGAAAGGACCGCTGATCTTACTTGC  >  W3110S.gb/3128339‑3128478

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