Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 657,810 T→G 8.1% F33V (TTC→GTC)  cspE → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb657,8100TG 8.1% 253.9 / 18.0 99F33V (TTC→GTC) cspEDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (50/41);  new base G (6/2);  total (56/43)
Fisher's exact test for biased strand distribution p-value = 4.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 102 total aligned reads displayed.

GTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCAC  >  W3110S.gb/657740‑657877
                                                                      |                                                                   
gTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACt                                                                     >  1:4092656/1‑71 (MQ=255)
   gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                                  <  1:7325972/71‑1 (MQ=255)
   gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtct                                                                  <  1:1870735/71‑1 (MQ=255)
   gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct                                                                  <  1:6192711/71‑1 (MQ=25)
   gTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGAAGTGTACGAACACGtct                                                                  <  1:1673916/71‑1 (MQ=25)
      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)
             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)
               cAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTtctc                                                                 <  1:2329395/60‑1 (MQ=255)
                aaaGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGAc                                                     >  1:7511701/1‑71 (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)
                          gTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt                                           <  1:371197/71‑1 (MQ=255)
                           tttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                            >  1:4961038/1‑69 (MQ=255)
                            ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                                        >  1:6116692/1‑56 (MQ=255)
                            ttCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCa                                                        >  1:805729/1‑56 (MQ=255)
                             tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGtttt                                           >  1:1441973/1‑68 (MQ=255)
                             tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                         >  1:5446529/1‑70 (MQ=255)
                             tCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaa                                         >  1:5787280/1‑70 (MQ=255)
                              cATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                                       >  1:2539389/1‑71 (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)
                                  aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc                                    <  1:2718746/70‑1 (MQ=255)
                                  aCTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAActc                                    <  1:6520011/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:6065605/68‑1 (MQ=255)
                                           gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                             <  1:6056808/68‑1 (MQ=255)
                                           gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                             <  1:1332503/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:5720084/71‑1 (MQ=255)
                                                cagcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgc                     <  1:4612432/71‑1 (MQ=255)
                                                 agcaAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaaa                                       >  1:3387973/1‑52 (MQ=255)
                                                      aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                            >  1:2666809/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)
                                                          gtgtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                              <  1:2039999/52‑1 (MQ=255)
                                                          gtgtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTg                              <  1:7376272/52‑1 (MQ=255)
                                                            gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt                               <  1:6076128/49‑1 (MQ=255)
                                                            gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt                               <  1:2857452/49‑1 (MQ=255)
                                                            gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGtt           <  1:3086961/69‑1 (MQ=255)
                                                             ttCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCg         <  1:3209608/70‑1 (MQ=255)
                                                              tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                        >  1:5052724/1‑38 (MQ=255)
                                                              tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                        >  1:64293/1‑38 (MQ=255)
                                                              tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                        >  1:536051/1‑38 (MQ=255)
                                                              tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGa        >  1:6842191/1‑70 (MQ=255)
                                                               cGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTc          >  1:861075/1‑67 (MQ=255)
                                                               cGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGaaa      >  1:6533196/1‑71 (MQ=255)
                                                               cGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGaa       >  1:6160481/1‑70 (MQ=255)
                                                               cGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGaa       <  1:1274424/70‑1 (MQ=255)
                                                                gTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAg                           <  1:3768765/49‑1 (MQ=255)
                                                                gTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAg                           <  1:1651230/49‑1 (MQ=255)
                                                                gTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAAt     <  1:6464489/71‑1 (MQ=255)
                                                                gTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAAt     >  1:778411/1‑71 (MQ=255)
                                                                 tACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct                                     <  1:3874649/38‑1 (MQ=255)
                                                                 tACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAAt     >  1:527191/1‑70 (MQ=255)
                                                                  acacTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATc    >  1:4000257/1‑70 (MQ=255)
                                                                  acacTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCa   >  1:4915944/1‑71 (MQ=255)
                                                                  acacTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGAGTAGAGTTCg         >  1:3076572/1‑65 (MQ=255)
                                                                   cacTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAg                           >  1:2279828/1‑46 (MQ=255)
                                                                    acTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg                    <  1:1387147/52‑1 (MQ=255)
                                                                    acTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCAc  >  1:1451045/1‑70 (MQ=255)
                                                                      ttCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGc                                >  1:1338331/1‑38 (MQ=255)
                                                                      ttCTCTGCAATCCAGACTAATGGTTTTAAAAATCTTGc                                >  1:1937939/1‑38 (MQ=255)
                                                                      |                                                                   
GTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAATCAC  >  W3110S.gb/657740‑657877

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: