Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,331,394 A→G 25.0% N72N (AAT→AAC 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,3940AG25.0% 34.0 / 10.4 24N72N (AAT→AACobgEGTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base A (4/14);  new base G (3/3);  total (7/17)
Fisher's exact test for biased strand distribution p-value = 3.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GCGTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGT  >  W3110S.gb/3331327‑3331461
                                                                   |                                                                   
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:748071/1‑71 (MQ=255)
                                               gTACAGTCTCGGCTTGCACCGTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAg                   >  1:185073/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)
                                                                   aTTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCtcgt  >  1:2180140/1‑68 (MQ=255)
                                                                   |                                                                   
GCGTACCTACCGGCACTTTAATCGTCACGTCTTTACCGCGCTTACCGGTACAGTCGCGGCTTGCGCCATTCTGACCGCGCTCTGCACGGAAAGATTTTTCAAAACGATAATCGATAAGCGTGTTCAGGTTCTCGT  >  W3110S.gb/3331327‑3331461

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: