Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,740 T→G 36.9% G118G (GGT→GGG araC → DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb70,7400TG36.9% ‑2.7 / 22.4 19G118G (GGT→GGGaraCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (6/6);  new base G (5/2);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

TACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCC  >  W3110S.gb/70675‑70806
                                                                 |                                                                  
tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctt                                                               >  1:1585433/1‑71 (MQ=255)
tACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttct                                                                <  1:18467/70‑1 (MQ=255)
 aCTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttcttt                                                              >  1:1648538/1‑71 (MQ=255)
  cTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGgttctttc                                                             >  1:297577/1‑71 (MQ=255)
     tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctt                                                               >  1:806029/1‑66 (MQ=255)
     tCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctt                                                               >  1:1884384/1‑66 (MQ=255)
       gTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCg                                                        <  1:1988477/71‑1 (MQ=255)
       gTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCg                                                        <  1:685264/71‑1 (MQ=255)
           gcgcgcCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATg                                                    <  1:2204909/71‑1 (MQ=255)
                  tACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttcttt                                                              >  1:1431912/1‑54 (MQ=255)
                    cTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCAc                                            >  1:1167841/1‑70 (MQ=255)
                         aTGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                      <  1:1304526/71‑1 (MQ=255)
                              tGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCa                                             >  1:1941449/1‑59 (MQ=255)
                                         gCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGt                      <  1:1976573/71‑1 (MQ=255)
                                           cGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                    <  1:470186/71‑1 (MQ=255)
                                              cAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTg                       <  1:520584/65‑1 (MQ=255)
                                               aaTATTTGCCAATACGGGTTTCTTTCGCCCGGATGaa                                                  >  1:39535/1‑37 (MQ=255)
                                                  atTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAg                              >  1:1765306/1‑54 (MQ=255)
                                                  atTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaaa             >  1:727489/1‑71 (MQ=255)
                                                   tttGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCaa              >  1:2040916/1‑69 (MQ=255)
                                                      gCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCa          <  1:981078/70‑1 (MQ=255)
                                                            aCGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGAc                          >  1:1978711/1‑48 (MQ=255)
                                                            aCGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                     >  1:1769901/1‑53 (MQ=255)
                                                             cGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGcc  >  1:2210342/1‑71 (MQ=255)
                                                             cGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGcc  >  1:1604972/1‑71 (MQ=255)
                                                               ggTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCa                               >  1:261110/1‑40 (MQ=255)
                                                                 |                                                                  
TACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCC  >  W3110S.gb/70675‑70806

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: