Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 70,244 T→C 30.4% intergenic (‑196/‑143) araB ← / → araC L‑ribulokinase/DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb70,2440TC30.4% 59.8 / 35.5 46intergenic (‑196/‑143)araB/araCL‑ribulokinase/DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (11/21);  new base C (10/4);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 2.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAG  >  W3110S.gb/70184‑70299
                                                            |                                                       
tgtgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGt                                               >  1:1381078/1‑71 (MQ=255)
  tgGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTa                                             >  1:1500581/1‑71 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:983994/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:1134196/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:1159838/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:1213849/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:59844/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:500423/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:1400295/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:359158/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:313618/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:198188/71‑1 (MQ=255)
     aCTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAg                                          <  1:498101/71‑1 (MQ=255)
              ccGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCttt                                  <  1:1065690/70‑1 (MQ=255)
              ccGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCttt                                  <  1:1695801/70‑1 (MQ=255)
              ccGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCttt                                  <  1:1027763/70‑1 (MQ=255)
                          acacTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAg                                          <  1:1534755/50‑1 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:1699471/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:926910/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:165926/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:27286/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:187167/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:1818580/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:1741888/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg                >  1:89794/1‑71 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCg                >  1:477549/1‑71 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGTCACGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTa                   >  1:1305801/1‑68 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCtt                                   >  1:975174/1‑52 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCtt                                   >  1:1661567/1‑52 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCtt                                   >  1:628161/1‑52 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCtt                                   >  1:1268494/1‑52 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCtt                                   >  1:376256/1‑52 (MQ=255)
                               tttGTTACGCGTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCtt                                   >  1:420022/1‑52 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:1662336/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:905312/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:68604/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:521535/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:431537/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:175408/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:225722/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:202425/70‑1 (MQ=255)
                                     aCGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTgg           <  1:193407/70‑1 (MQ=255)
                                        cgTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGt         >  1:1769251/1‑69 (MQ=255)
                                        cgTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGt         >  1:1890985/1‑69 (MQ=255)
                                        agTTTTTGTCATGGCTTTGGCCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGt         >  1:1610962/2‑69 (MQ=255)
                                              tGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCgag  <  1:1725818/70‑1 (MQ=255)
                                                            |                                                       
TGTGGACTTTTCTGCCGTGATTATAGACACTTTTGTTACGCGTTTTTGTCATGGCTTTGGTCCCGCTTTGTTACAGAATGCTTTTAATAAGCGGGGTTACCGGTTGGGTTAGCGAG  >  W3110S.gb/70184‑70299

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: