Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,119,282 C→T 21.9% E420E (GAG→GAA wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,119,2820CT21.9% 75.8 / 16.8 41E420E (GAG→GAAwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base C (18/14);  new base T (5/4);  total (23/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 = 8.84e-01

CGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCG  >  W3110S.gb/2119213‑2119349
                                                                     |                                                                   
cGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCt                                                                    >  1:3121411/1‑71 (MQ=255)
        tCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGca                                                              >  1:324214/1‑69 (MQ=255)
            acgacAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAg                                                           <  1:1066905/68‑1 (MQ=255)
                      gAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTg                                              <  1:1602571/71‑1 (MQ=255)
                      gAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTg                                              <  1:2695483/71‑1 (MQ=255)
                       aaTTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTg                                              >  1:1716309/1‑70 (MQ=255)
                         tttCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAg                                                           <  1:2040386/55‑1 (MQ=255)
                                 aaaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCt                                   <  1:3788577/71‑1 (MQ=255)
                                  aaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCaa                                                                 <  1:750480/40‑1 (MQ=255)
                                    agaagCCGACCTTCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTccc                                       >  1:1952566/1‑63 (MQ=255)
                                        gCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTc                            <  1:731362/71‑1 (MQ=255)
                                          cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCAc                          <  1:925139/71‑1 (MQ=255)
                                          cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCAc                          <  1:3431443/71‑1 (MQ=255)
                                          cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCAc                          <  1:1941588/71‑1 (MQ=255)
                                            aCCTTCATTTCACCTCCCCGATGCGCTCAAGCACAg                                                           <  1:413277/36‑1 (MQ=255)
                                               ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                          <  1:2618492/66‑3 (MQ=255)
                                               ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                          <  1:413721/66‑3 (MQ=255)
                                               ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                          <  1:3550499/66‑3 (MQ=255)
                                               ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                          <  1:1633846/66‑3 (MQ=255)
                                                tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCAc                          <  1:2876005/65‑1 (MQ=255)
                                                tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                        >  1:487594/1‑67 (MQ=255)
                                                tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                        >  1:1283918/1‑67 (MQ=255)
                                                tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                        >  1:2121304/1‑67 (MQ=255)
                                                 cATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGAc                       >  1:10062/1‑67 (MQ=255)
                                                 cATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGAc                       >  1:92340/1‑67 (MQ=255)
                                                 cATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGAc                       >  1:1978661/1‑67 (MQ=255)
                                                 cATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGATCGCGAc                       >  1:1403891/1‑67 (MQ=255)
                                                 cATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGa                    >  1:613996/1‑70 (MQ=255)
                                                 cATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGa                    >  1:3230896/1‑70 (MQ=255)
                                                 cATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                   >  1:2671621/1‑71 (MQ=255)
                                                       aCCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCtt              <  1:844057/70‑1 (MQ=255)
                                                       aCCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCtt              <  1:2592283/70‑1 (MQ=255)
                                                         cTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCAc           >  1:1795821/1‑71 (MQ=255)
                                                         cTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCAc           >  1:3129392/1‑71 (MQ=255)
                                                           ccccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACt          >  1:2849266/1‑70 (MQ=255)
                                                           ccccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACt          >  1:2829198/1‑70 (MQ=255)
                                                            cccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGATTCACt          >  1:3253535/1‑69 (MQ=255)
                                                            cccGATGCGCTCAAGCACAGACTGAACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACt          >  1:845578/1‑69 (MQ=255)
                                                              cGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGa                                 >  1:2474880/1‑44 (MQ=255)
                                                              cGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGAc                       >  1:2637940/1‑54 (MQ=255)
                                                                 tGCGCTCAAGCACAGACTGCACCATCTGCATTCCt                                       >  1:705444/1‑35 (MQ=255)
                                                                  gcgcTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgcg  >  1:2339735/1‑71 (MQ=255)
                                                                   cgcTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACg                  >  1:1504743/1‑54 (MQ=255)
                                                                   cgcTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcgcg  <  1:3489712/70‑1 (MQ=255)
                                                                     |                                                                   
CGAAGGTTTCTGACGACAGCGGGAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCG  >  W3110S.gb/2119213‑2119349

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: