Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,710,918 A→T 100% 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,9180AT80.9% 38.9 / 6.0 21V168E (GTG→GAG) yfiCpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base T (2/15);  total (6/15)
Fisher's exact test for biased strand distribution p-value = 2.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGATGGGTGATCAAGCG  >  W3110S.gb/2710873‑2710987
                                             |                                                                     
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:671476/71‑1 (MQ=255)
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:489208/71‑1 (MQ=255)
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:127636/71‑1 (MQ=255)
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:1475032/71‑1 (MQ=255)
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:276793/71‑1 (MQ=255)
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:1633557/71‑1 (MQ=255)
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:1025400/71‑1 (MQ=255)
cTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGt                                              <  1:1871813/71‑1 (MQ=255)
   aaTGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAAt                                           >  1:292739/1‑71 (MQ=255)
   aaTGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAAt                                           >  1:1800038/1‑71 (MQ=255)
    aTGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATg                                          <  1:351969/71‑1 (MQ=255)
    aTGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATg                                          <  1:336977/71‑1 (MQ=255)
    aTGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATg                                          <  1:1890281/71‑1 (MQ=255)
     tGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATg                                          <  1:1298118/70‑1 (MQ=255)
       ccAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACACAGAAAAATCCCTCTTCGGTAATGCAc                                       >  1:1762977/1‑71 (MQ=255)
                          aCCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGc                     >  1:1418438/1‑70 (MQ=255)
                          aCCAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGc                     >  1:1133660/1‑70 (MQ=255)
                            cAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc                  <  1:159183/71‑1 (MQ=255)
                            cAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc                  <  1:130671/71‑1 (MQ=255)
                            cAATCTGCTCTGGCAGCTCTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAAc                  <  1:403859/71‑1 (MQ=255)
                                 tGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGATgg             >  1:1138629/1‑71 (MQ=255)
                                            caCTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGATGGGTGATCAAGCg  <  1:1883934/71‑1 (MQ=255)
                                            caCTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGATGGGTGATCAAGCg  <  1:1484082/71‑1 (MQ=255)
                                             |                                                                     
CTTAATGCCAGCTCCGTAAAACCATTACCAATCTGCTCTGGCAGCACTACGCAGAAAAATCCCTCTTCGGTAATGCACTCCGCCGCACAGGTCAGCAACGATGGGTGATCAAGCG  >  W3110S.gb/2710873‑2710987

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: