Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,879,768 A→G 100% R86R (CGT→CGC ytfE ← predicted regulator of cell morphogenesis and cell wall metabolism

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,879,7680AG78.1% 104.8 / 24.3 56R86R (CGT→CGCytfEpredicted regulator of cell morphogenesis and cell wall metabolism
Reads supporting (aligned to +/- strand):  ref base A (12/0);  new base G (43/0);  total (56/0)
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
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GTCGGCGTGAACGCGCTCGACTTTAGTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGC  >  minE/2879705‑2879820
                                                               |                                                    
gTCGGCGTGAACGCGCTCGACTTTAGTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGCCGa                                                   >  1:145192/1‑67 (MQ=255)
gTCGGCGTGAACGCGCTCGACTTTAGTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGa                                                   >  1:76874/1‑67 (MQ=255)
gTCGGCGTGAACGCGCTCGACTTTAGTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGa                                                   >  1:172776/1‑67 (MQ=255)
gTCGGCGTGAACGCGCTCGACTTTAGTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGa                                                   >  1:781424/1‑67 (MQ=255)
gTCGGCGTGAACGCGCTCGACTTTAGTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGa                                                   >  1:1635513/1‑67 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATa                                           >  1:44865/1‑50 (MQ=19)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:244659/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1889407/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1959911/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1986797/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:202219/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:2039297/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:2148252/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:730004/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:2251580/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1820744/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:245453/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:312702/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:394086/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:436891/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:490217/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:572145/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:631092/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:698223/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:712033/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:725782/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1886634/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:188910/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1812674/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:961212/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1728656/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1320238/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1727953/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1630042/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:162659/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1403626/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:148234/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:931288/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1502084/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1564976/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1028106/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1584537/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:843082/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:159038/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:2273134/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:2130481/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:179644/1‑66 (MQ=255)
                         gTCGCTTGCAGGATCAGCTCCGGAAGCTGCTCGCTATGGCGATCGTGATAGCGCACAATAATGTgg                           >  1:1476116/1‑66 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:961803/1‑67 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:2162449/1‑67 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:1742928/1‑67 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:1702255/1‑67 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:1587692/1‑67 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:1502082/1‑67 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:1212517/1‑67 (MQ=255)
                                                 aGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGgcgc  >  1:1146448/1‑67 (MQ=255)
                                                               |                                                    
GTCGGCGTGAACGCGCTCGACTTTAGTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGC  >  minE/2879705‑2879820

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