Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 113,035 T→C 17.4% D62G (GAT→GGT)  coaE ← dephospho‑CoA kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb113,0350TC17.4% 49.9 / 5.7 23D62G (GAT→GGT) coaEdephospho‑CoA kinase
Reads supporting (aligned to +/- strand):  ref base T (14/5);  new base C (2/2);  total (16/7)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GTTAAGCCAGTTTTTCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAGCGCAGGTGCACCTGGTTCAAC  >  W3110S.gb/112965‑113102
                                                                      |                                                                   
gTTAAGCCAGTTTTTCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCAt                                                                     >  1:1024525/1‑71 (MQ=255)
   aaGCCAGTTTTTCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCACcag                                                                  >  1:1319933/1‑71 (MQ=255)
   aaGCCAGTTTTTCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCACcag                                                                  >  1:213250/1‑71 (MQ=255)
     gCCAGTTTTTCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCATcag                                                                  >  1:1541587/1‑69 (MQ=255)
        aGTTTTTCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCagcagc                                                              >  1:807076/1‑70 (MQ=255)
          tttttCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCACCAGCAGCAAt                                                           <  1:670536/71‑1 (MQ=255)
          tttttCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCACCAGCAGCAAt                                                           <  1:1736118/71‑1 (MQ=255)
                    ccGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAgcg                                                 >  1:588241/1‑71 (MQ=255)
                           ggCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGt                                          >  1:7539/1‑71 (MQ=255)
                                      gCTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAAt                                >  1:15329/1‑70 (MQ=255)
                                       cTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAg                                           <  1:394216/58‑1 (MQ=255)
                                       cTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATgg                              <  1:1276050/71‑1 (MQ=255)
                                        tCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCAAAAGTGAt                                       >  1:1744413/1‑61 (MQ=255)
                                             cAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCAt                           >  1:908183/1‑68 (MQ=255)
                                                ggCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAgcg                     >  1:1376974/1‑71 (MQ=255)
                                                 gCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAgcgc                    <  1:700912/71‑1 (MQ=255)
                                                  cccGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAgcg                     >  1:1412521/1‑69 (MQ=255)
                                                  cccGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAgcg                     >  1:1212348/1‑69 (MQ=255)
                                                               tGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAGAGCAGGTGCACCTGGtt      <  1:1051743/71‑1 (MQ=255)
                                                                gTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAg                                    >  1:116719/1‑40 (MQ=255)
                                                                 ttCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAGCGCAGGTGCACCTGGTTc     >  1:1580062/1‑70 (MQ=255)
                                                                   ccATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAGCGCAGGTGCACCTGGTACAac  >  1:743591/1‑71 (MQ=255)
                                                                     aTCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATg                          <  1:1306197/45‑1 (MQ=255)
                                                                     aTCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAGCGCAGGTGCACCTGGTTCAac  >  1:23876/1‑69 (MQ=255)
                                                                      |                                                                   
GTTAAGCCAGTTTTTCTCTTCCGGGTTGGCGAAGATCCGCTCGCGCAAGGCCCGGCGCTGCAATGTTCCATCAGCAGCAATCATGTTAGCGCCAAAGTGATCAGCAATGGCATGTAGCGCAGGTGCACCTGGTTCAAC  >  W3110S.gb/112965‑113102

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