Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,710,918 A→T 11.3% V168E (GTG→GAG)  yfiC ← predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,710,9180AT11.3% 144.7 / 10.2 57V168E (GTG→GAG) yfiCpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base T (4/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 4.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGAT  >  W3110S.gb/2710851‑2710974
                                                                   |                                                        
cAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                                                         <  1:2686299/69‑1 (MQ=255)
cAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                                                         <  1:1440635/69‑1 (MQ=255)
cAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                                                         <  1:799450/69‑1 (MQ=255)
cAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                                                         <  1:1482297/69‑1 (MQ=255)
cAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCAc                                                         <  1:1659408/69‑1 (MQ=255)
      tAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGc                                                    >  1:2083363/1‑68 (MQ=255)
          tGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa                                               >  1:1855775/1‑69 (MQ=255)
          tGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa                                               >  1:1807105/1‑69 (MQ=255)
                    tGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCtt                                     >  1:3676744/1‑69 (MQ=255)
                    tGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCtt                                     >  1:2503947/1‑69 (MQ=255)
                    tGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCtt                                     >  1:2956527/1‑69 (MQ=255)
                    tGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCtt                                     >  1:457603/1‑69 (MQ=255)
                    tGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCtt                                     >  1:1223684/1‑69 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGTAGCACTACGCAGAAAAATCCct                                        >  1:2789230/1‑64 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTa                                                       >  1:431770/1‑49 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa                                               <  1:496613/57‑1 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGaaa                                               <  1:2534696/57‑1 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCctc                                       <  1:3475901/65‑1 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCctc                                       <  1:195662/65‑1 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCctc                                       <  1:2358701/65‑1 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCctc                                       <  1:2469760/65‑1 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTc                                    <  1:1486978/68‑1 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTc                                    >  1:2185724/1‑68 (MQ=255)
                      cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTc                                    >  1:3578750/1‑68 (MQ=255)
                        tAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGt                                 <  1:3214179/69‑1 (MQ=255)
                        tAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGt                                 <  1:1171125/69‑1 (MQ=255)
                        tAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGt                                 <  1:795749/69‑1 (MQ=255)
                        tAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGt                                 <  1:1471544/69‑1 (MQ=255)
                           tGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTaa                               >  1:557559/1‑68 (MQ=255)
                             ccAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGaaaa                                              >  1:437251/1‑51 (MQ=255)
                             ccAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATg                             >  1:2199360/1‑68 (MQ=255)
                             ccAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATg                             >  1:3190279/1‑68 (MQ=255)
                             ccAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATg                             >  1:2097993/1‑68 (MQ=255)
                                      tAAAACCATTACCAATCTGCTCTGGCAGGACTACGCAGAAAAATCCCTCTTCGGTAATGCACTccgccg                   <  1:743351/69‑1 (MQ=255)
                                        aaaCCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATGCACTccgc                     <  1:2952759/65‑1 (MQ=255)
                                        aaaCCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATGCACTccgc                     <  1:3617603/65‑1 (MQ=255)
                                        aaaCCATTACCAATCTGCTCTGGCAGCCCTACGCAGAAAAATCCCTCTTCGGTAATGCACTccgc                     <  1:885363/65‑1 (MQ=255)
                                        aaaCCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTccgccg                   >  1:3476887/1‑67 (MQ=255)
                                        aaaCCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTACGGTAATGCACTccgccg                   >  1:940910/1‑67 (MQ=255)
                                          aCCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCg                                   >  1:2768520/1‑49 (MQ=255)
                                          aCCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTa                                >  1:87952/1‑52 (MQ=255)
                                                aCCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCa          >  1:291378/1‑68 (MQ=255)
                                                aCCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCa          >  1:2460764/1‑68 (MQ=255)
                                                  cAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGc        <  1:1031257/68‑1 (MQ=255)
                                                    aTCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:3714839/69‑1 (MQ=255)
                                                    aTCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:295816/69‑1 (MQ=255)
                                                    aTCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:1727274/69‑1 (MQ=255)
                                                    aTCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:1254557/69‑1 (MQ=255)
                                                       tggtctGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGAt  >  1:3253591/4‑69 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAgaac      <  1:3403521/65‑4 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCCAc     <  1:1593860/66‑1 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:3620790/66‑1 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:348532/66‑1 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:759512/66‑1 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc     <  1:970315/66‑1 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGAt  >  1:2137217/1‑69 (MQ=255)
                                                       tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGAt  >  1:2132114/1‑69 (MQ=255)
                                                                   |                                                        
CAAACGTAAATGCCAGCCCATGCTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGAT  >  W3110S.gb/2710851‑2710974

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: