Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,331,391 G→A 27.2% G73G (GGC→GGT obgE ← GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,331,3910GA27.2% 30.4 / 8.9 22G73G (GGC→GGTobgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base G (2/14);  new base A (3/3);  total (5/17)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTC  >  W3110S.gb/3331324‑3331459
                                                                   |                                                                    
cacGCGTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGcc                                                                    <  1:1614150/70‑1 (MQ=255)
   gcgTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTc                                                                <  1:24210/71‑1 (MQ=255)
    cgTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCt                                                               <  1:3536940/71‑1 (MQ=255)
                cACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTg                                                   <  1:1206169/71‑1 (MQ=255)
                  cTTTAATCGTGACGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCa                                                 <  1:898852/71‑1 (MQ=255)
                  cTTTAATCGTGACGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCa                                                 <  1:448761/71‑1 (MQ=255)
                  cTTTAATCGTGACGTCTTTACCGCGCTTACCGGTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCa                                                 <  1:1728068/71‑1 (MQ=255)
                  cTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTg                                                              <  1:1087841/58‑1 (MQ=255)
                             aCGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGAttt                                      <  1:277639/71‑1 (MQ=255)
                                    tACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCaa                                 <  1:620691/69‑1 (MQ=255)
                                      ccGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAAc                              <  1:2356464/70‑1 (MQ=255)
                                       cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                            <  1:2751402/71‑1 (MQ=255)
                                       cgcgCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGAGCTCTGCACGGAAAGAt                                        <  1:1395431/59‑1 (MQ=255)
                                       cgcgCTTACCGGTACAGTCGCGGCTCGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGa                            <  1:2997427/71‑1 (MQ=255)
                                                  gTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATaa                  >  1:1695475/1‑70 (MQ=255)
                                                  gTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                 >  1:185073/1‑71 (MQ=255)
                                                  gTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                 >  1:748071/1‑71 (MQ=255)
                                                        tCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGa                     <  1:2477209/61‑1 (MQ=255)
                                                          gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGtt           >  1:2826448/1‑69 (MQ=255)
                                                          gcgGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGtt           >  1:1768686/1‑69 (MQ=255)
                                                               ttGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATc                       <  1:197596/52‑1 (MQ=255)
                                                               ttGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATc                       <  1:2661749/52‑1 (MQ=255)
                                                                 gcgcCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATa                          >  1:468427/1‑47 (MQ=255)
                                                                 gcgcCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtctc  >  1:3119331/1‑71 (MQ=255)
                                                                  cgcCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTtctc  <  1:133859/70‑1 (MQ=255)
                                                                   |                                                                    
CACGCGTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTC  >  W3110S.gb/3331324‑3331459

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: