Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,213,568 C→T 52.9% intergenic (+24/+55) rpoD → / ← ygjF RNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,213,5680CT52.9% ‑5.1 / 21.3 17intergenic (+24/+55)rpoD/ygjFRNA polymerase, sigma 70 (sigma D) factor/G/U mismatch‑specific DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (4/5);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

CCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAC  >  W3110S.gb/3213508‑3213637
                                                            |                                                                     
ccGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCg                                                             >  1:184531/1‑71 (MQ=255)
     cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCgg                                                            <  1:1060672/67‑1 (MQ=255)
     cTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCAc                                                         >  1:482220/1‑70 (MQ=255)
               gcGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCctc                                                >  1:1507895/1‑69 (MQ=255)
                    gCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGcaca                                          <  1:842535/70‑1 (MQ=255)
                            gACGATTAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCAcc                                 >  1:1640759/1‑71 (MQ=255)
                               gATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCttt                              <  1:651420/71‑1 (MQ=255)
                                 ttAATCGGTAGGCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACaa                                         <  1:1335181/58‑1 (MQ=255)
                                         tAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCTGTGGCGt                    <  1:241469/71‑1 (MQ=255)
                                            gCCGGATCAGGCGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttt                 <  1:1113823/71‑1 (MQ=255)
                                                   cAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGc          >  1:1309481/1‑71 (MQ=255)
                                                      gCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGtttt                 >  1:1634285/1‑61 (MQ=255)
                                                       cGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAAc                                       >  1:663253/1‑38 (MQ=255)
                                                       cGTTATGCCGCACCCGGCACTAGGCCCTCTGCACAAAc                                       >  1:109173/1‑38 (MQ=255)
                                                       cGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCAc      <  1:1294459/71‑1 (MQ=255)
                                                         ttATGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTc                            <  1:551632/47‑1 (MQ=255)
                                                           aCGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGccc  >  1:1476049/1‑69 (MQ=255)
                                                            |                                                                     
CCGTTCTGAAGTGCTGCGTAGCTTCCTGGACGATTAATCGGTAGGCCGGATCAGGCGTTACGCCGCACCCGGCACTAGGCCCTCTGCACAAACGCCACCTTTTCGGTGGCGTTTTTTATCGCCCACGCAC  >  W3110S.gb/3213508‑3213637

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: