Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,973,590 T→G 54.4% intergenic (+6/‑15) radA → / → nadR predicted repair protein/bifunctional DNA‑binding transcriptional repressor and NMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,973,5900TG54.4% 14.3 / 43.4 46intergenic (+6/‑15)radA/nadRpredicted repair protein/bifunctional DNA‑binding transcriptional repressor and NMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (5/16);  new base G (21/4);  total (26/20)
Fisher's exact test for biased strand distribution p-value = 6.71e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGC  >  minE/2973536‑2973652
                                                      |                                                              
tGGCGTTAAAAAACTCTCCGACGCGCTTATCGTGTTCGACGACTTATAATGAGATATACGGAGGgaga                                                   <  1:1020587/68‑1 (MQ=255)
tGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGgaga                                                   <  1:2453238/68‑1 (MQ=255)
tGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGgaga                                                   <  1:3183153/68‑1 (MQ=255)
tGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGgaga                                                   <  1:2208495/68‑1 (MQ=255)
tGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGgaga                                                   <  1:563167/68‑1 (MQ=255)
tGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGgaga                                                   <  1:1859884/68‑1 (MQ=255)
       aaaaaaCTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACg                                                           <  1:2717479/53‑1 (MQ=255)
       aaaaaaCTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACg                                                           <  1:2927027/53‑1 (MQ=255)
       aaaaaaCTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATgtcgt                                           <  1:2936668/69‑1 (MQ=255)
             ctctCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTg                                     <  1:102522/69‑1 (MQ=255)
             ctctCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTg                                     <  1:431481/69‑1 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGTAGATATGTCGTCATTTGATTACCt                              >  1:2793820/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGAtt                                  >  1:1933413/1‑64 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGAtt                                  >  1:1323105/1‑64 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:490610/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:270583/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:565602/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:2504698/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:1123330/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:1209676/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:2450061/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:2334491/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:2313116/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:1238523/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:3194645/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:2074686/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:1565204/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:1908769/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:186220/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:1619861/1‑68 (MQ=255)
                     cgcgCTTAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:535226/1‑68 (MQ=255)
                     cgcgCTGAGCGTGTTCGACGACTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCt                              >  1:437638/1‑68 (MQ=255)
                       cgcTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGaa                           <  1:568197/69‑1 (MQ=255)
                       cgcTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGaa                           <  1:2290776/69‑1 (MQ=255)
                       cgcTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGaa                           <  1:2601494/69‑1 (MQ=255)
                       cgcTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGaa                           <  1:1086602/69‑1 (MQ=255)
                        gcTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGaa                           <  1:2164921/68‑1 (MQ=255)
                                  ttCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTAcc                               >  1:1377022/1‑54 (MQ=255)
                                  ttCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                >  1:2458703/1‑69 (MQ=255)
                                  ttCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                >  1:1062437/1‑69 (MQ=255)
                                          cTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAg        <  1:1002808/69‑1 (MQ=255)
                                          cTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAg        <  1:1689937/69‑1 (MQ=255)
                                          cTTATAATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAg        <  1:610860/69‑1 (MQ=255)
                                              taATGAGAGATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAg        <  1:1967294/65‑1 (MQ=255)
                                                aTGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGc  >  1:1929327/1‑69 (MQ=255)
                                                aTGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGc  >  1:1767860/1‑69 (MQ=255)
                                                      |                                                              
TGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGC  >  minE/2973536‑2973652

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: