Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,552,318 T→C 63.2% 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,3180TC63.2% 13.6 / 15.7 19F67S (TTT→TCT) murPfused predicted enzyme IIBC components of PTS
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base C (3/9);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

AAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAG  >  W3110S.gb/2552278‑2552385
                                        |                                                                   
aaGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCaa                                        >  1:758279/1‑70 (MQ=255)
aaGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCaa                                        >  1:2077446/1‑70 (MQ=255)
aaGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCaa                                        >  1:2304790/1‑70 (MQ=255)
   ggCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    <  1:3441807/71‑1 (MQ=255)
   ggCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    <  1:2596676/71‑1 (MQ=255)
   ggCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    <  1:3244711/71‑1 (MQ=255)
   ggCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    <  1:3367010/71‑1 (MQ=255)
    gCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCg                                   <  1:1355731/71‑1 (MQ=255)
    gCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCg                                   <  1:3511124/71‑1 (MQ=255)
          ttttGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGc                                    <  1:367117/64‑1 (MQ=255)
            ttGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagc                             >  1:566864/1‑69 (MQ=255)
            ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagcg                            <  1:852047/70‑1 (MQ=255)
            ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagcg                            <  1:785017/70‑1 (MQ=255)
            ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagcg                            <  1:2218527/70‑1 (MQ=255)
            ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagcg                            <  1:2083819/70‑1 (MQ=255)
            ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagcg                            <  1:1710693/70‑1 (MQ=255)
            ttGACCAGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATgagc                             >  1:1524013/1‑69 (MQ=255)
                  aGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAgctgct                     >  1:2249488/1‑71 (MQ=255)
                  aGCGACCAGGTGCAGGTCGTTTCTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAgctgct                     >  1:175821/1‑71 (MQ=255)
                                     tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAACGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:1681776/71‑1 (MQ=255)
                                     tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:3433717/71‑1 (MQ=255)
                                     tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:689403/71‑1 (MQ=255)
                                     tttctGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAg  <  1:75709/71‑1 (MQ=255)
                                        |                                                                   
AAGGGCGTCATTTTGACCAGCGACCAGGTGCAGGTCGTTTTTGGGCCTGGTAAAGCGCATCGTGCCGCAAAAGCGATGAGCGAGCTGCTGGGAGAAGCACCGGTACAG  >  W3110S.gb/2552278‑2552385

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: