Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,119,312 T→C 14.7% T410T (ACA→ACG wcaK ← predicted pyruvyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,119,3120TC14.7% 143.2 / 16.1 61T410T (ACA→ACGwcaKpredicted pyruvyl transferase
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base C (5/4);  total (33/28)
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.94e-01

AAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCC  >  W3110S.gb/2119246‑2119381
                                                                  |                                                                     
aaaaGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCt                                                                   <  1:3788577/71‑1 (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)
              ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                                          <  1:3550499/66‑3 (MQ=255)
              ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                                          <  1:413721/66‑3 (MQ=255)
              ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                                          <  1:2618492/66‑3 (MQ=255)
              ttCATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTcgc                                                          <  1:1633846/66‑3 (MQ=255)
               tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCAc                                                          <  1:2876005/65‑1 (MQ=255)
               tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                        >  1:1283918/1‑67 (MQ=255)
               tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                        >  1:2121304/1‑67 (MQ=255)
               tCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGa                                                        >  1:487594/1‑67 (MQ=255)
                cATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGAc                                                       >  1:1978661/1‑67 (MQ=255)
                cATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGAc                                                       >  1:10062/1‑67 (MQ=255)
                cATTTCACCTCCCCGATGCGTTCAAGCACAGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGAc                                                       >  1:92340/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: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)
                                       aaGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                             <  1:3535425/70‑1 (MQ=255)
                                       aaGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                             <  1:2769581/70‑1 (MQ=255)
                                       aaGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                             <  1:1726601/70‑1 (MQ=255)
                                           acaGATTGCACCATCTGCATTCCCGTCTGACGTTCGCGACTGACAGCTTCGTTAAGTCGTGCGTTCAgcgc                        >  1:3317742/1‑71 (MQ=255)
                                           acaGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAgcgc                        <  1:328746/71‑1 (MQ=255)
                                                   cACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCa                                            >  1:2568473/1‑43 (MQ=255)
                                                   cACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCt                >  1:3432483/1‑71 (MQ=255)
                                                   cACCATATGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTcg                                    >  1:492601/1‑51 (MQ=255)
                                                     ccATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                >  1:2539554/1‑53 (MQ=255)
                                                     ccATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCa                                >  1:402197/1‑53 (MQ=255)
                                                         cTGCATTCCTGTCTGACGCTCACGACTGACGGCTTc                                             >  1:16700/1‑36 (MQ=255)
                                                         cTGCATTCCTGTCTGACGCTCACGACTGACGGCTTc                                             >  1:2217601/1‑36 (MQ=255)
                                                         cTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAg                           >  1:2241975/1‑54 (MQ=255)
                                                          tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                             <  1:2792994/51‑1 (MQ=255)
                                                          tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTc                             <  1:1015431/51‑1 (MQ=255)
                                                          tGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCt           >  1:1087228/1‑69 (MQ=255)
                                                            cATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTa          <  1:3241606/68‑1 (MQ=255)
                                                            cATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTa          <  1:1914233/68‑1 (MQ=255)
                                                            cATTCCTGTCTGACGCTCACGACTGACGGCCTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAgg       <  1:2458027/71‑1 (MQ=255)
                                                             aTTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTa          <  1:2126944/67‑1 (MQ=255)
                                                             aTTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt      <  1:3285587/71‑1 (MQ=255)
                                                             aTTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt      <  1:828416/71‑1 (MQ=255)
                                                             aTTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGt      <  1:1626519/71‑1 (MQ=255)
                                                              ttCCTGTCTGACGCTCACGACTGACGGCTTCACTAAg                                       <  1:688330/37‑1 (MQ=255)
                                                              ttCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTa     <  1:3541453/71‑1 (MQ=255)
                                                                ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATc   >  1:235272/1‑71 (MQ=255)
                                                                ccTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATc   >  1:2120416/1‑71 (MQ=255)
                                                                 cTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATcc  >  1:1651465/1‑71 (MQ=255)
                                                                  |                                                                     
AAAAGAAGCCGACCTTCATTTCACCTCCCCGATGCGCTCAAGCACAGACTGCACCATCTGCATTCCTGTCTGACGCTCACGACTGACGGCTTCACTAAGTCGCGCATTCAGCGCCGGAAGCTGGCCTAAGGTATCC  >  W3110S.gb/2119246‑2119381

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: