Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,381,341 T→C 9.3% S248S (AGT→AGC degQ → serine endoprotease, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,381,3410TC 9.3% 115.7 / 4.4 43S248S (AGT→AGCdegQserine endoprotease, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base C (2/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

AACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGCATCA  >  W3110S.gb/3381282‑3381411
                                                           |                                                                      
aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgcgc                                                              >  1:3122528/1‑70 (MQ=255)
aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgcgc                                                              >  1:464645/1‑70 (MQ=255)
aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgcgc                                                              >  1:2418505/1‑70 (MQ=255)
aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGgcgc                                                              >  1:3646384/1‑70 (MQ=255)
aaCACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGcgcg                                                             <  1:3422909/71‑1 (MQ=255)
  cacTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGaa                                                           <  1:1675760/71‑1 (MQ=255)
   acTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGACATCCCCAGTAATATGGCGCGAac                                                          <  1:2264188/71‑1 (MQ=255)
    cTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:1441645/58‑1 (MQ=255)
    cTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:1561023/58‑1 (MQ=255)
    cTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:3440094/58‑1 (MQ=255)
    cTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:3466969/58‑1 (MQ=255)
      gCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGCAata                                                                    >  1:517748/1‑58 (MQ=255)
      gCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGCAata                                                                    >  1:418151/1‑58 (MQ=255)
        aaTCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:153497/54‑1 (MQ=255)
        aaTCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:1725132/54‑1 (MQ=255)
        aaTCATTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:1415377/54‑1 (MQ=255)
          tCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTaa                                                                      <  1:2544808/52‑1 (MQ=255)
            cTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAata                                                                    >  1:3634196/1‑52 (MQ=255)
                                    gggATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAAt                         >  1:2347191/1‑71 (MQ=255)
                                            attTGCCATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAAcgcg                 <  1:3518317/71‑1 (MQ=255)
                                             tttGCCATCCCCAGTAATATGGCGGGAACACTGGCgca                                                 <  1:3641210/38‑1 (MQ=255)
                                             tttGCCATCCCCAGTAATATGGCGCGAACACTGGCgca                                                 <  1:29948/38‑1 (MQ=255)
                                               tGCCATCCCCAGTAATATGGCGCGAACACTGGCgcagc                                               >  1:2970124/1‑38 (MQ=255)
                                               tGCCATCCCCAGTAATATGGCGCGAACACTGGCgcagc                                               >  1:1704145/1‑38 (MQ=255)
                                                  cATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGt           <  1:1224439/71‑1 (MQ=255)
                                                  cATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGt           <  1:491170/71‑1 (MQ=255)
                                                  cATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGt           <  1:2137983/71‑1 (MQ=255)
                                                  cATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGt           <  1:1747428/71‑1 (MQ=255)
                                                  cATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGt           <  1:1560213/71‑1 (MQ=255)
                                                  cATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGt           <  1:78829/71‑1 (MQ=255)
                                                    tCCCCAGCAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCa                       <  1:2309288/57‑1 (MQ=255)
                                                    tCCCCAGCAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCa                       <  1:677008/57‑1 (MQ=255)
                                                     ccccAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTa         >  1:751419/1‑70 (MQ=255)
                                                      cccAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAgg       >  1:624110/1‑71 (MQ=255)
                                                      cccAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAgg       >  1:1079200/1‑71 (MQ=255)
                                                        cAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCa                       <  1:1689782/53‑1 (MQ=255)
                                                         aGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCaaa                     >  1:972690/1‑54 (MQ=255)
                                                         aGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCaaa                     >  1:3627156/1‑54 (MQ=255)
                                                         aGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCaaa                     >  1:2772507/1‑54 (MQ=255)
                                                         aGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCaaa                     >  1:2039847/1‑54 (MQ=255)
                                                          gTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCaaa                     >  1:2864394/1‑53 (MQ=255)
                                                          gTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCaaa                     >  1:1494627/1‑53 (MQ=255)
                                                          gTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCaaa                     >  1:1350653/1‑53 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACtt                                  >  1:3179436/1‑39 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGc      <  1:3622454/67‑1 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGc      <  1:3008081/67‑1 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGc      <  1:1678063/67‑1 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGc      <  1:141194/67‑1 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGCATc   >  1:3415299/1‑70 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGCATc   >  1:912729/1‑70 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGCATCa  <  1:2865610/71‑1 (MQ=255)
                                                           tAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGGGTTTGTTAGGCATCa  >  1:1580064/1‑70 (MQ=255)
                                                           |                                                                      
AACACTGCAATCCTTGCGCCTGGCGGCGGGAGCGTCGGGATTGGATTTGCCATCCCCAGTAATATGGCGCGAACACTGGCGCAGCAGCTTATCGACTTTGGTGAAATCAAACGCGGTTTGTTAGGCATCA  >  W3110S.gb/3381282‑3381411

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: