Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,608,013 G→A 42.9% G67E (GGA→GAA)  lsrF → predicted aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,608,0130GA42.9% 6.8 / 22.4 21G67E (GGA→GAA) lsrFpredicted aldolase
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base A (3/6);  total (7/14)
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.31e-01

TCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAATGT  >  W3110S.gb/1607948‑1608072
                                                                 |                                                           
tCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTg                                                        <  1:521621/71‑1 (MQ=255)
    aaaaCGGGTAAAACCTTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACg                                                    <  1:1249768/71‑1 (MQ=255)
              aaaCCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGAtata                                          <  1:195109/71‑1 (MQ=255)
              aaaCCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGAtata                                          <  1:2617463/71‑1 (MQ=255)
               aaCCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGAtata                                          <  1:195112/70‑1 (MQ=255)
                 ccGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAat                                       <  1:746364/71‑1 (MQ=255)
                         cTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGcccc                               >  1:2018107/1‑71 (MQ=255)
                           ggCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGATATAAATATCGCCCCGc                             >  1:1966555/1‑71 (MQ=255)
                             cTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGaa                                                      >  1:126723/1‑44 (MQ=255)
                                     cATGGTTATTTTCAGGGACCGACTACCGGACTTGAACg                                                    <  1:860513/38‑1 (MQ=255)
                                     cATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACa                   <  1:138332/71‑1 (MQ=255)
                                         gTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGcc               >  1:2115999/1‑71 (MQ=255)
                                          ttATTTTCAGGGACCGACTACCGAACTTGAACGCATTGAtata                                          <  1:1559632/43‑1 (MQ=255)
                                          ttATTTTCAGGGACCGACTACCGAACTTGAACGCATTGAtata                                          <  1:925889/43‑1 (MQ=255)
                                           tATTTTCAGGGACCGACTACCGAACTTGAACGCAATGATATAAATATCGCCCCGCTGTTTGAAc                    >  1:530398/1‑64 (MQ=255)
                                                  aGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTAtt       <  1:2254895/70‑1 (MQ=255)
                                                  aGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTAtt       <  1:1022517/70‑1 (MQ=255)
                                                  aGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTAtt       <  1:156525/70‑1 (MQ=255)
                                                      aCCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAAtgt  >  1:2654449/1‑71 (MQ=255)
                                                          aCTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAAc                    <  1:721982/49‑1 (MQ=255)
                                                              ccGAACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACa                   >  1:2603746/1‑46 (MQ=255)
                                                                 gACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACa                   >  1:2170448/1‑43 (MQ=255)
                                                                 gACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACa                   >  1:1841700/1‑43 (MQ=255)
                                                                 |                                                           
TCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAATGT  >  W3110S.gb/1607948‑1608072

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