Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,608,013 G→A 25.0% 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,0130GA25.0% 57.3 / 15.4 32G67E (GGA→GAA) lsrFpredicted aldolase
Reads supporting (aligned to +/- strand):  ref base G (15/9);  new base A (5/3);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAATGTGTACGC  >  W3110S.gb/1607943‑1608078
                                                                      |                                                                 
tttAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCgg                                                                   <  1:1954599/71‑1 (MQ=255)
       cGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGAc                                                                 >  1:197347/1‑66 (MQ=255)
       cGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGaa                                                            >  1:2640511/1‑71 (MQ=255)
       cGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGaa                                                            >  1:1917399/1‑71 (MQ=255)
                gTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACg                                                          <  1:2256466/64‑1 (MQ=255)
                gTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGa                                                    >  1:1058645/1‑70 (MQ=255)
                gTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGa                                                    >  1:3869582/1‑70 (MQ=255)
                gTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGa                                                    >  1:2150126/1‑70 (MQ=255)
                 tAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGata                                                  >  1:2125234/1‑71 (MQ=255)
                  aaaaCCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGata                                                  >  1:2142860/1‑70 (MQ=255)
                  aaaaCCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGata                                                  >  1:3848850/1‑70 (MQ=255)
                   aaaCCTTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGAACTTGAACGCATTGAtata                                                <  1:1464079/71‑1 (MQ=255)
                           aTGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGaa                                                            >  1:3113936/1‑51 (MQ=255)
                           aTGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGc                                        <  1:601695/71‑1 (MQ=255)
                             gCTGGCTTTTGACCATGGTTATTTTCAGGGACCGCCTACCGAACTTGAACGCATTGATATAAATATCGccc                                      >  1:682205/1‑71 (MQ=255)
                              cTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGcccc                                     >  1:687462/1‑71 (MQ=255)
                              cTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGcccc                                     >  1:2297350/1‑71 (MQ=255)
                              cTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGcccc                                     >  1:2386238/1‑71 (MQ=255)
                                  cTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTg                                 >  1:3699934/1‑71 (MQ=255)
                                       gACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTg                             <  1:762681/70‑1 (MQ=255)
                                        aCCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGaa                           <  1:3170769/71‑1 (MQ=255)
                                        aCCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGaa                           <  1:1596726/71‑1 (MQ=255)
                                         ccATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGaa                           <  1:3091174/70‑1 (MQ=255)
                                           aTGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGAtata                                                <  1:3919406/47‑1 (MQ=255)
                                           aTGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGAtata                                                <  1:1679997/47‑1 (MQ=255)
                                              gTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGcc                     >  1:3486876/1‑71 (MQ=255)
                                              gTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGcc                     >  1:4038446/1‑71 (MQ=255)
                                              gTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGcc                     >  1:2078286/1‑71 (MQ=255)
                                              gTTATTTTCAGGGACCGACTACCGAACTTGAAGGCATTGATATAAATATCGCCCCGCTg                                 <  1:3943934/59‑1 (MQ=255)
                                                 aTTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAat                                             <  1:2515195/44‑1 (MQ=255)
                                                          gACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAAtg         >  1:714790/1‑71 (MQ=255)
                                                               aCTACCGAACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACa                         <  1:793564/50‑1 (MQ=255)
                                                                 tACCGAACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAATGTGTAcgc  >  1:1283769/1‑71 (MQ=255)
                                                                      |                                                                 
TTTAATCCGAAAACGGGTAAAACCGTGATGCTGGCTTTTGACCATGGTTATTTTCAGGGACCGACTACCGGACTTGAACGCATTGATATAAATATCGCCCCGCTGTTTGAACATGCCGATGTATTAATGTGTACGC  >  W3110S.gb/1607943‑1608078

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: