Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,879,792 A→G 47.1% H78H (CAT→CAC 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,7920AG47.1% ‑0.5 / 42.5 51H78H (CAT→CACytfEpredicted regulator of cell morphogenesis and cell wall metabolism
Reads supporting (aligned to +/- strand):  ref base A (22/5);  new base G (12/12);  total (34/17)
Fisher's exact test for biased strand distribution p-value = 3.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGC  >  minE/2879730‑2879860
                                                              |                                                                    
gTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTc                                                                 <  1:1204648/68‑1 (MQ=255)
gTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTc                                                                 <  1:181041/68‑1 (MQ=255)
gTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTc                                                                 <  1:1607015/68‑1 (MQ=255)
gTCGCATGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTc                                                                 <  1:1725413/68‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:1645400/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:2549177/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:2673742/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:3159589/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:575671/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:616340/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:1617320/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:1497352/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:1438679/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:1329618/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:1299579/69‑1 (MQ=255)
   gCTTGCAGGATCAGCTCCGGCAGCTGCTCGCGATGGCGATCGTGATAGCGCACAATAATGTGGTCgatg                                                             <  1:1281608/69‑1 (MQ=255)
   gCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatg                                                             >  1:1324964/1‑69 (MQ=255)
   gCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatg                                                             >  1:3146866/1‑69 (MQ=255)
   gCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatg                                                             >  1:3035726/1‑69 (MQ=255)
   gCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatg                                                             >  1:2950904/1‑69 (MQ=255)
   gCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatg                                                             >  1:2886438/1‑69 (MQ=255)
   gCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatg                                                             >  1:2189597/1‑69 (MQ=255)
   gCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatg                                                             >  1:875324/1‑69 (MQ=255)
     ttGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCgatgat                                                           <  1:3149656/69‑1 (MQ=255)
             tCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGcc                                                    >  1:1527682/1‑68 (MQ=255)
             tCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCa                                                   >  1:1199181/1‑69 (MQ=255)
             tCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCa                                                   >  1:826492/1‑69 (MQ=255)
             tCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCAATGATTTCTGCCa                                                   >  1:305662/1‑69 (MQ=255)
                           tGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTAc                                       >  1:1721792/1‑67 (MQ=255)
                           tGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTAc                                       >  1:444233/1‑67 (MQ=255)
                           tGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTAc                                       >  1:3073388/1‑67 (MQ=255)
                           tGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTAc                                       >  1:2548437/1‑67 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:1786053/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:979180/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:1410127/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:733574/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:3342438/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:3308865/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:1782076/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:303441/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:291917/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:2108688/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:2678012/1‑69 (MQ=255)
                               cgcgATGGCGATCGTGATAGCGCACAATAATGTGGTCGATGATTTCCGCCAGCGGGGCGCTACGCCAGt                                 >  1:2381560/1‑69 (MQ=255)
                                         aTCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGc                                     >  1:2763736/1‑55 (MQ=255)
                                         aTCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAAt                       >  1:199337/1‑69 (MQ=255)
                                         aTCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAAt                       >  1:3081606/1‑69 (MQ=255)
                                         aTCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAAt                       >  1:2032793/1‑69 (MQ=255)
                                         aTCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAAt                       >  1:3395280/1‑69 (MQ=255)
                                           cGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCg                     >  1:200412/1‑69 (MQ=255)
                                                         atgatATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGc       >  1:739419/1‑69 (MQ=255)
                                                              atGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTg   <  1:3300432/68‑1 (MQ=255)
                                                              atGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTg   <  1:1664307/68‑1 (MQ=255)
                                                              atGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGc  <  1:1134217/69‑1 (MQ=255)
                                                              |                                                                    
GTCGCTTGCAGAATCAGCTCCGGCAGTTGCTCGCGGTGACGATCGTGGTAGCGCACGATGATATGGTCGATGATTTCTGCCAGCGGGGCGCTACGCCAGTCTTTCTCAATCGGTTGTTCAGCGAGCTTTGC  >  minE/2879730‑2879860

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