Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 666,009 C→T 36.3% N489N (AAC→AAT ycbY → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE666,0090CT36.3% 35.1 / 29.9 55N489N (AAC→AATycbYpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (23/12);  new base T (8/12);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 9.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TCGGTACTGGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCACTTGTG  >  minE/665944‑666074
                                                                 |                                                                 
tCGGTACTGGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCaa                                                                >  1:2351498/1‑69 (MQ=255)
        ggggATTGCGCCAAACCAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGaa                                                        >  1:3081614/1‑69 (MQ=255)
         gggATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGaaa                                                       >  1:1057641/1‑69 (MQ=255)
         gggATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGaaa                                                       >  1:1116678/1‑69 (MQ=255)
         gggATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGaaa                                                       >  1:676706/1‑69 (MQ=255)
         gggATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGaaa                                                       >  1:3248649/1‑69 (MQ=255)
         gggATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGaaa                                                       >  1:271542/1‑69 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCATAAACCAATACCAGAAACTGGGCgag                                              <  1:1388151/64‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:2394194/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:698880/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:2857786/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:2794798/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:275337/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:2661191/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:2624317/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:2254052/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:1750830/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:1583891/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:1276321/57‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATAc                                                             >  1:3008751/1‑49 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCgag                                              <  1:715736/64‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCgag                                              <  1:650899/64‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCgag                                              <  1:1496804/64‑1 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAgg                                          >  1:1400185/1‑68 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAgg                                          >  1:2828095/1‑68 (MQ=255)
                       caaacTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAgg                                          >  1:1457084/1‑68 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:125201/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:389034/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:1370176/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:3121879/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:2962307/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:2625293/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:1758300/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:258422/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:2573042/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:1908241/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:2124716/1‑69 (MQ=255)
                           cTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAg                                     >  1:2322741/1‑69 (MQ=255)
                            tGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAATCAATACCAGAAACt                                                     <  1:3079042/52‑1 (MQ=255)
                             ggTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGtt                                   <  1:165083/69‑1 (MQ=255)
                             ggTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGtt                                   <  1:3333229/69‑1 (MQ=255)
                                              aaCGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCg                     <  1:21376/66‑1 (MQ=255)
                                              aaCGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGaa                   >  1:148974/1‑68 (MQ=255)
                                              aaCGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGaa                   >  1:2458317/1‑68 (MQ=255)
                                              aaCGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGaa                   >  1:142890/1‑68 (MQ=255)
                                                  ccAGAAGGGCAAAAATCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTGGaa                            >  1:1538493/1‑55 (MQ=255)
                                                  ccAGAAGGGCAAAAATCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTGGAAGTTACCGAAtata               >  1:633504/1‑68 (MQ=255)
                                                  ccAGAAGGGCAAAAATCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTGGAAGTTACCGAATATaa              >  1:1832428/1‑69 (MQ=255)
                                                         ggCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCAc       <  1:1055511/69‑1 (MQ=255)
                                                           cAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCACtt     >  1:2922856/1‑69 (MQ=255)
                                                           cAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCACtt     >  1:1817717/1‑69 (MQ=255)
                                                              aaaCCAATACCAGAAACTGGGCGAGAAGGGCGAGtt                                   <  1:2298132/36‑1 (MQ=255)
                                                              aaaCCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCACTtgtg  <  1:1685942/69‑1 (MQ=255)
                                                              aaaCCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCACTtgtg  <  1:2518083/69‑1 (MQ=255)
                                                               aaCCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCACTtgtg  <  1:82880/68‑1 (MQ=255)
                                                                 |                                                                 
TCGGTACTGGGGATTGCGCCAAACAAACTGGTGCTGAAAACCCGTGAACGCCAGAAGGGCAAAAACCAATACCAGAAACTGGGCGAGAAGGGCGAGTTTCTTGAAGTTACCGAATATAACGCTCACTTGTG  >  minE/665944‑666074

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: