Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,824,650 A→G 37.5% Y107Y (TAT→TAC yjdC ← predicted transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,824,6500AG37.5% 8.2 / 11.4 24Y107Y (TAT→TACyjdCpredicted transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (6/9);  new base G (4/5);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAATCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGAT  >  minE/2824586‑2824713
                                                                |                                                               
tGAAATCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATaa                                                              <  1:3235615/68‑1 (MQ=255)
tGAAATCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATaa                                                              <  1:2862634/68‑1 (MQ=255)
    aTCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACg                                                          <  1:1220428/68‑1 (MQ=255)
    aTCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACg                                                          <  1:3430676/68‑1 (MQ=255)
    aTCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACg                                                          <  1:323730/68‑1 (MQ=255)
    aTCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACg                                                          <  1:2726351/68‑1 (MQ=255)
       gTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTAc                                                       >  1:554145/1‑68 (MQ=255)
                  ttttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGGTAAAACGTACAGGCAGCGATa                                            >  1:1058724/1‑68 (MQ=255)
                  ttttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGGTAAAACGTACAGGCAGCGATa                                            >  1:958995/1‑68 (MQ=255)
                  ttttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGGTAAAACGTACAGGCAGCGATa                                            >  1:101911/1‑68 (MQ=255)
                   tttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGGTAAAACGTACAGGCAGCGATa                                            >  1:1991208/1‑67 (MQ=255)
                   tttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATa                                            >  1:2288615/1‑67 (MQ=255)
                   tttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATa                                            >  1:2241157/1‑67 (MQ=255)
                   tttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATa                                            >  1:1847578/1‑67 (MQ=255)
                   tttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATa                                            >  1:525908/1‑67 (MQ=255)
                   tttttGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATa                                            >  1:715317/1‑67 (MQ=255)
                                     ttGATGAATAGGGTGGCCAGGATCGGGATAAAACg                                                          <  1:687682/35‑1 (MQ=255)
                                                   ggCCATGATCGGGGTAAAACGTACAGGCAGCGATAAACAGACAGCCCGGATAg                          <  1:485893/53‑1 (MQ=255)
                                                   ggCCAGGATCGGGGTAAAACGTACAGGCAGCGATAAACAGACAGCCCGGATAg                          <  1:2574662/53‑1 (MQ=255)
                                                   ggCCAGGATCGGGGTAAAACGTACAGGCAGCGATAAACAGACAGCCCGGATAg                          <  1:312814/53‑1 (MQ=255)
                                                   ggCCAGGATCGGGGTAAAACGTACAGGCAGCGATAAACAGACAGCCCGGATAg                          <  1:118020/53‑1 (MQ=255)
                                                    gCCAGGATCGGGGTAAAACGTACAGGCAGCGATAAACAGACAGCCCGGATAg                          <  1:552399/52‑1 (MQ=255)
                                                    gCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCa         <  1:756419/69‑1 (MQ=255)
                                                           tCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGAt  <  1:3070091/69‑1 (MQ=255)
                                                                |                                                               
TGAAATCGTAGGCCGCGCTTTTTTGCTGATCGGCCAGTTGATGAATAGGGTGGCCAGGATCGGGATAAAACGTACAGGCAGCGATAAATAGACAGCCCGGATAGCGGTTGTTTTTAACGCACTCCGAT  >  minE/2824586‑2824713

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: