Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 657,805 T→A 7.8% V31E (GTA→GAA)  cspE → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb657,8050TA 7.8% 277.5 / 10.1 99V31E (GTA→GAA) cspEDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (50/41);  new base A (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 = 9.84e-01

Only 100 of 102 total aligned reads displayed.

TAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAAT  >  W3110S.gb/657737‑657874
                                                                    |                                                                     
tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                                      >  1:3956411/1‑70 (MQ=255)
tAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTa                                                                      >  1:2956402/1‑70 (MQ=255)
 aaGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTac                                                                     <  1:4326880/70‑1 (MQ=255)
   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)
         tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTacac                                                                   >  1:3305081/1‑64 (MQ=255)
         tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                                            >  1:2713369/1‑71 (MQ=255)
         tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                                            >  1:3708827/1‑71 (MQ=255)
         tAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTg                                                            >  1:2044733/1‑71 (MQ=255)
              aGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                       <  1:1736395/71‑1 (MQ=255)
               gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                       >  1:7063905/1‑70 (MQ=255)
               gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                       >  1:7035422/1‑70 (MQ=255)
               gTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATc                                                       >  1:5668175/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)
                     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)
                             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: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:6065605/68‑1 (MQ=255)
                                              gACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGa                          <  1:6056808/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)
                                                 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)
                                                         aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGttt                                         >  1:4568527/1‑42 (MQ=255)
                                                         aGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAAct                                  <  1:317109/49‑1 (MQ=255)
                                                         aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg                 >  1:7816748/1‑66 (MQ=25)
                                                         aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg                 >  1:487866/1‑66 (MQ=25)
                                                         aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg                 >  1:5940894/1‑66 (MQ=25)
                                                         aGAAGTGTACGAACACGTCTTTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAgcg                 >  1:497230/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:2857452/49‑1 (MQ=255)
                                                               gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCt                            <  1:6076128/49‑1 (MQ=255)
                                                               gtTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGtt        <  1:3086961/69‑1 (MQ=255)
                                                                ttCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCg      <  1:3209608/70‑1 (MQ=255)
                                                                 tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                     >  1:64293/1‑38 (MQ=255)
                                                                 tCGTACACTTCTCTGCAATCCAGACTAATGGTTTTaaa                                     >  1:5052724/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)
                                                                    |                                                                     
TAAGTGGTTTAATGAGTCCAAAGGATTCGGTTTCATTACTCCGGAAGACGGCAGCAAAGACGTGTTCGTACACTTCTCTGCAATCCAGACTAATGGTTTTAAAACTCTTGCTGAAGGTCAGCGCGTAGAGTTCGAAAT  >  W3110S.gb/657737‑657874

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: