Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,119,294 C→T 23.2% Q416Q (CAG→CAA wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,119,2940CT23.2% 77.6 / 19.6 43Q416Q (CAG→CAAwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base C (18/15);  new base T (6/4);  total (24/19)
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.86e-01

GAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGC  >  W3110S.gb/2119235‑2119359
                                                           |                                                                 
gAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTg                                                        <  1:1602571/71‑1 (MQ=255)
gAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTg                                                        <  1:2695483/71‑1 (MQ=255)
 aaTTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTg                                                        >  1:1716309/1‑70 (MQ=255)
           aaaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCt                                             <  1:3788577/71‑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:1941588/71‑1 (MQ=255)
                    cGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCAc                                    <  1:3431443/71‑1 (MQ=255)
                         ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                    <  1:1633846/66‑3 (MQ=255)
                         ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                    <  1:413721/66‑3 (MQ=255)
                         ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                    <  1:3550499/66‑3 (MQ=255)
                         ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                    <  1:2618492/66‑3 (MQ=255)
                          tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCAc                                    <  1:2876005/65‑1 (MQ=255)
                          tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                  >  1:1283918/1‑67 (MQ=255)
                          tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                  >  1:487594/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:3230896/1‑70 (MQ=255)
                           cATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGa                              >  1:613996/1‑70 (MQ=255)
                           cATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGAc                             >  1:2671621/1‑71 (MQ=255)
                                 aCCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCtt                        <  1:2592283/70‑1 (MQ=255)
                                 aCCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCtt                        <  1:844057/70‑1 (MQ=255)
                                   cTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCAc                     >  1:3129392/1‑71 (MQ=255)
                                   cTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCAc                     >  1:1795821/1‑71 (MQ=255)
                                     ccccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACt                    >  1:2829198/1‑70 (MQ=255)
                                     ccccGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACt                    >  1:2849266/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)
                                                  aaGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc       <  1:3535425/70‑1 (MQ=255)
                                                  aaGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc       <  1:1726601/70‑1 (MQ=255)
                                                  aaGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc       <  1:2769581/70‑1 (MQ=255)
                                                      acaGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAgcgc  >  1:3317742/1‑71 (MQ=255)
                                                      acaGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgcgc  <  1:328746/71‑1 (MQ=255)
                                                           |                                                                 
GAATTTCAGTAAAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGC  >  W3110S.gb/2119235‑2119359

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: