Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,695,615 C→T 55.6% intergenic (+108/+58) nudE → / ← mrcA ADP‑ribose diphosphatase/fused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,695,6150CT55.6% 3.2 / 33.8 45intergenic (+108/+58)nudE/mrcAADP‑ribose diphosphatase/fused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (11/9);  new base T (9/16);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-01

GGGTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAACAA  >  minE/2695551‑2695681
                                                                |                                                                  
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                                 <  1:1315646/68‑1 (MQ=255)
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                                 <  1:208468/68‑1 (MQ=255)
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                                 <  1:533532/68‑1 (MQ=255)
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                                 <  1:1098527/68‑1 (MQ=255)
  gTTTCTGCAATTTATTGTGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATa                                                              <  1:2572830/69‑1 (MQ=255)
  gTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATa                                                              >  1:2350435/1‑69 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                             >  1:1350806/1‑67 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                             >  1:2825651/1‑67 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                             >  1:2390426/1‑67 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                             >  1:887537/1‑67 (MQ=255)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:778255/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:606870/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:813354/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:3129678/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:2405166/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:2999090/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:1046574/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:2174877/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                             >  1:2081369/2‑66 (MQ=2)
                         cACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCaa                                        <  1:2259720/68‑1 (MQ=255)
                         cACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCaa                                        <  1:1071653/68‑1 (MQ=255)
                            gATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCaa                                        <  1:2509534/65‑1 (MQ=255)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:1054931/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:85103/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:1111050/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:1156988/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:536930/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:407320/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:3377377/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:1305918/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:2956307/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:2815865/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:2778066/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:1558139/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:2415982/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:1830146/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                                <  1:229493/40‑1 (MQ=21)
                                            tAACGCGTTCAACCGCATCTGGCATAAACAAGTGCACttt                                                <  1:2775946/40‑1 (MQ=19)
                                                     cACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGAAAACTGAAAAGGCGCCGAAGCGCCTTTTTa            >  1:2753902/1‑68 (MQ=255)
                                                            tATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAAc    <  1:1399803/69‑1 (MQ=255)
                                                              tCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAACaa  >  1:1347920/1‑69 (MQ=255)
                                                              tCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAACaa  >  1:474258/1‑69 (MQ=255)
                                                              tCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAACaa  >  1:1235807/1‑69 (MQ=255)
                                                              tCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAACaa  >  1:2705330/1‑69 (MQ=255)
                                                              tCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAACaa  >  1:1914497/1‑69 (MQ=255)
                                                                |                                                                  
GGGTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTAATCAGAACAA  >  minE/2695551‑2695681

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: