Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 2,552,318 T→C F67S (TTT→TCT)  murP → fused predicted enzyme IIBC components of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,552,3180TC97.4% 135.0 / ‑4.8 39F67S (TTT→TCT) murPfused predicted enzyme IIBC components of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (20/18);  minor base A (1/0);  total (21/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAG  >  W3110S.gb/2552252‑2552385
                                                                  |                                                                   
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:976181/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:734299/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:1636548/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:1759355/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:1771355/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:1811453/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:3340116/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:2209478/70‑1 (MQ=255)
ccAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTgg                                                                  <  1:2984386/70‑1 (MQ=255)
   tcatCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCt                                                              <  1:1529622/69‑1 (MQ=255)
   tcatCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCt                                                              <  1:3819069/69‑1 (MQ=255)
   tcatCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCt                                                              <  1:102441/69‑1 (MQ=255)
    catCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCt                                                              <  1:502819/69‑1 (MQ=255)
     atCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCt                                                              <  1:300084/69‑1 (MQ=255)
                                 tCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCg                                   <  1:1312826/68‑1 (MQ=255)
                                 tCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCg                                   <  1:2424215/68‑1 (MQ=255)
                                 tCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCg                                   <  1:2363747/68‑1 (MQ=255)
                                   aTTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAg                                     >  1:2514096/1‑64 (MQ=255)
                                   aTTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAg                                     >  1:736735/1‑64 (MQ=255)
                                   aTTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAg                                     >  1:142628/1‑64 (MQ=255)
                                   aTTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAg                                     >  1:3988719/1‑64 (MQ=255)
                                   aTTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAg                                     >  1:3568384/1‑64 (MQ=255)
                                   aTTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAg                                     >  1:213393/1‑64 (MQ=255)
                                   aTTTTGACCAGCGACCAGGTGCAGGTCGTTTATGGGCCTGGTAAAGCGCATCGTGCCGCaaaa                                      >  1:3668490/1‑63 (MQ=255)
                                      ttGACCAGCGCCCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCa                                                   >  1:481646/1‑47 (MQ=255)
                                      ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCaaaa                                      >  1:2439251/1‑60 (MQ=255)
                                      ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCaaaa                                      >  1:2587079/1‑60 (MQ=255)
                                      ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCaaaa                                      >  1:612277/1‑60 (MQ=255)
                                      ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCaaa                                       >  1:2535733/1‑59 (MQ=255)
                                      ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagc                             >  1:2113868/1‑69 (MQ=255)
                                      ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagc                             >  1:2149705/1‑69 (MQ=255)
                                      ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagc                             >  1:670527/1‑69 (MQ=255)
                                       tGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCg                                   <  1:3630831/62‑1 (MQ=255)
                                                             cGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAg                                     >  1:634722/1‑38 (MQ=255)
                                                             cGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    >  1:957519/1‑39 (MQ=255)
                                                             cGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    >  1:2455933/1‑39 (MQ=255)
                                                             cGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    >  1:3292290/1‑39 (MQ=255)
                                                             cGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTa     >  1:3699367/1‑70 (MQ=255)
                                                             cGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTa     >  1:256380/1‑70 (MQ=255)
                                                               tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:3921406/71‑1 (MQ=255)
                                                               tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:732332/71‑1 (MQ=255)
                                                               tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:3040059/71‑1 (MQ=255)
                                                               tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:2857640/71‑1 (MQ=255)
                                                               tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:917205/71‑1 (MQ=255)
                                                               tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:2790784/71‑1 (MQ=255)
                                                               tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:963093/71‑1 (MQ=255)
                                                                  |                                                                   
CCAATATCAAAACGCTGGAAGGCGTGAAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAG  >  W3110S.gb/2552252‑2552385

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