Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,695,607 Δ1 bp 61.5% intergenic (+100/+66) 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,6070G.61.5% 28.0 / 66.1 39intergenic (+100/+66)nudE/mrcAADP‑ribose diphosphatase/fused penicillin‑binding protein 1a murein transglycosylase and murein transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (6/8);  new base . (9/15);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

GGGTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTA  >  minE/2695551‑2695671
                                                        |                                                                
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                       <  1:1098527/68‑1 (MQ=255)
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                       <  1:533532/68‑1 (MQ=255)
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                       <  1:1315646/68‑1 (MQ=255)
gggTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGc                                                       <  1:208468/68‑1 (MQ=255)
  gTTTCTGCAATTTATTGTGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATa                                                    <  1:2572830/69‑1 (MQ=255)
  gTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATa                                                    >  1:2350435/1‑69 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                   >  1:887537/1‑67 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                   >  1:2825651/1‑67 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                   >  1:2390426/1‑67 (MQ=255)
     tCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATaa                                                   >  1:1350806/1‑67 (MQ=255)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:606870/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:2081369/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:2174877/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:778255/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:813354/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:2999090/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:2405166/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:1046574/2‑66 (MQ=2)
     ccTGCAATTTATTGAATTTTCAATTTTTTGTTGGCCGGATAACGCGTTCACCCGCATCTGGCATaa                                                   >  1:3129678/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:2415982/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:3377377/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:407320/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:85103/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:536930/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:229493/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:1830146/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:1558139/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:1305918/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:1156988/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:1111050/40‑1 (MQ=21)
                                            tAACGCGTTCACCCGCATCTGGCATAAACAAGTGCACttt                                      <  1:1054931/40‑1 (MQ=21)
                                            tAACGCGTTCAACCGCATCTGGCATAAACAAGTGCACttt                                      <  1:2775946/40‑1 (MQ=19)
                                                     cACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGAAAACTGAAAAGGCGCCGAAGCGCCTTTTTa  >  1:2753902/1‑68 (MQ=255)
                                                        |                                                                
GGGTTTCTGCAATTTATTGAGTTTTCACGATTTTGTTAGCCGGATAACGCGTTCACGCCGTATCCGGCATAAACAAGTGCACTTTGTCAGCAAACTGAAAAGGCGCCGAAGCGCCTTTTTA  >  minE/2695551‑2695671

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: