Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 657,797 C→A 8.7% D28E (GAC→GAA cspE → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb657,7970CA 8.7% 292.7 / 7.5 91D28E (GAC→GAAcspEDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (46/37);  new base A (6/2);  total (52/39)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

AACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCG  >  W3110S.gb/657732‑657861
                                                                 |                                                                
aaCGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACgtgt                                                              >  1:4278237/1‑70 (MQ=255)
  cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                           <  1:3431917/71‑1 (MQ=255)
  cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                           <  1:4821509/71‑1 (MQ=255)
  cGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                           <  1:4435674/71‑1 (MQ=255)
   gTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                           >  1:3369993/1‑70 (MQ=255)
     tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                         >  1:2956402/1‑70 (MQ=255)
     tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                         >  1:3956411/1‑70 (MQ=255)
      aaGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac                                                        <  1:4326880/70‑1 (MQ=255)
        gTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACt                                                     >  1:4092656/1‑71 (MQ=255)
           gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                  <  1:1870735/71‑1 (MQ=255)
           gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                  <  1:7325972/71‑1 (MQ=255)
           gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct                                                  <  1:1673916/71‑1 (MQ=25)
           gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct                                                  <  1:6192711/71‑1 (MQ=25)
              tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTacac                                                      >  1:3305081/1‑64 (MQ=255)
              tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                               >  1:3708827/1‑71 (MQ=255)
              tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                               >  1:2044733/1‑71 (MQ=255)
              tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                               >  1:2713369/1‑71 (MQ=255)
                   aGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          <  1:1736395/71‑1 (MQ=255)
                    gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          >  1:5668175/1‑70 (MQ=255)
                    gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          >  1:7035422/1‑70 (MQ=255)
                    gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          >  1:7063905/1‑70 (MQ=255)
                     tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACgtgt                                                              <  1:7669629/49‑1 (MQ=255)
                     tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACgtgt                                                              <  1:4037709/49‑1 (MQ=255)
                     tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          >  1:1791814/1‑69 (MQ=255)
                     tCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                        <  1:2587352/71‑1 (MQ=255)
                      ccAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          >  1:5165310/1‑68 (MQ=255)
                       cAAAGTATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                 <  1:5409624/60‑1 (MQ=255)
                       cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTc                                                            >  1:4542572/1‑49 (MQ=255)
                       cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                 <  1:2329395/60‑1 (MQ=255)
                        aaaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc                                     >  1:7511701/1‑71 (MQ=255)
                          aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                           >  1:6676690/1‑47 (MQ=255)
                          aGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                         >  1:6912680/1‑49 (MQ=255)
                           ggATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac                                                        <  1:5762591/49‑1 (MQ=255)
                            gATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAAt                                 <  1:5627647/71‑1 (MQ=255)
                             aTTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          <  1:6614943/61‑1 (MQ=255)
                               tCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGt                              >  1:3941755/1‑71 (MQ=255)
                                cGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGt                              <  1:7902918/70‑1 (MQ=255)
                                 ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                  >  1:2361193/1‑49 (MQ=255)
                                 ggTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAAt                                 <  1:3558143/66‑1 (MQ=255)
                                  gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCg                                                           <  1:2895236/39‑1 (MQ=255)
                                  gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt                           <  1:371197/71‑1 (MQ=255)
                                   tttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                            >  1:4961038/1‑69 (MQ=255)
                                    ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                        >  1:805729/1‑56 (MQ=255)
                                    ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                        >  1:6116692/1‑56 (MQ=255)
                                     tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt                           >  1:1441973/1‑68 (MQ=255)
                                     tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                         >  1:5787280/1‑70 (MQ=255)
                                     tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                         >  1:5446529/1‑70 (MQ=255)
                                      cATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                       >  1:2539389/1‑71 (MQ=255)
                                       aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTaca                                                       >  1:7744232/1‑38 (MQ=255)
                                       aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                       <  1:3504256/70‑1 (MQ=255)
                                       aTTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc                      >  1:1308675/1‑71 (MQ=255)
                                        ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc                                     >  1:1283595/1‑55 (MQ=255)
                                        ttACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct                     >  1:5507068/1‑71 (MQ=255)
                                          aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTaca                                                       >  1:4893129/1‑35 (MQ=255)
                                          aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc                    <  1:6520011/70‑1 (MQ=255)
                                          aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc                    <  1:2718746/70‑1 (MQ=255)
                                           cTCCTGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAc                      >  1:1057902/1‑67 (MQ=255)
                                            tCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                 <  1:1127338/39‑1 (MQ=255)
                                              cGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCa                                             >  1:2385637/1‑41 (MQ=255)
                                               ggAAGACGGCAGCAAAGACGTGTTCTTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                         <  1:364742/60‑1 (MQ=255)
                                               ggAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                         <  1:2888888/60‑1 (MQ=255)
                                               ggAAGACGGAAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg              <  1:7858808/71‑1 (MQ=255)
                                                gAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                          >  1:2234123/1‑42 (MQ=255)
                                                gAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa             <  1:3585061/71‑1 (MQ=255)
                                                 aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGa                                      >  1:2249522/1‑45 (MQ=255)
                                                 aaGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg              >  1:3603404/1‑69 (MQ=255)
                                                  aGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg              >  1:6733109/1‑68 (MQ=255)
                                                   gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa             <  1:1332503/68‑1 (MQ=255)
                                                   gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa             <  1:6056808/68‑1 (MQ=255)
                                                   gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa             <  1:6065605/68‑1 (MQ=255)
                                                   gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAgg          >  1:1862221/1‑71 (MQ=255)
                                                    aCGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGt         <  1:6885882/71‑1 (MQ=255)
                                                     cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc                                     >  1:47654/1‑42 (MQ=255)
                                                     cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAg           >  1:5650000/1‑68 (MQ=255)
                                                     cGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTc        <  1:5249713/71‑1 (MQ=255)
                                                      ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGt         >  1:4800274/1‑69 (MQ=255)
                                                      ggCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCa       <  1:7136193/71‑1 (MQ=255)
                                                        cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTa                          <  1:6929864/50‑1 (MQ=255)
                                                        cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc     <  1:4612432/71‑1 (MQ=255)
                                                        cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc     <  1:5720084/71‑1 (MQ=255)
                                                         agcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                       >  1:3387973/1‑52 (MQ=255)
                                                              aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                            >  1:2666809/1‑42 (MQ=255)
                                                              aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                            >  1:4568527/1‑42 (MQ=255)
                                                              aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                            >  1:5670939/1‑42 (MQ=255)
                                                              aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct                     <  1:317109/49‑1 (MQ=255)
                                                              aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg    >  1:497230/1‑66 (MQ=25)
                                                              aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg    >  1:487866/1‑66 (MQ=25)
                                                              aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg    >  1:7816748/1‑66 (MQ=25)
                                                              aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg    >  1:5940894/1‑66 (MQ=25)
                                                              aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGcgcg  >  1:1229900/1‑68 (MQ=25)
                                                              aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTCCTGAAGGTCAGcgcg  >  1:6715600/1‑68 (MQ=25)
                                                                 |                                                                
AACGTTAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCG  >  W3110S.gb/657732‑657861

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: