Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 666,267 A→G 36.1% E575E (GAA→GAG ycbY → predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE666,2670AG36.1% 17.7 / 18.9 36E575E (GAA→GAGycbYpredicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (2/21);  new base G (7/6);  total (9/27)
Fisher's exact test for biased strand distribution p-value = 4.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTG  >  minE/666202‑666335
                                                                 |                                                                    
gTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATg                                                                    >  1:849425/1‑68 (MQ=255)
          gATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCa                                                         <  1:1089118/69‑1 (MQ=255)
          gATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCa                                                         <  1:2739775/69‑1 (MQ=255)
          gATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCa                                                         <  1:1418267/69‑1 (MQ=255)
                  ggTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGc                                                   >  1:973914/1‑67 (MQ=255)
                  ggTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGt                                                 >  1:3048705/1‑69 (MQ=255)
                  ggTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGt                                                 >  1:1811463/1‑69 (MQ=255)
                  ggTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGt                                                 >  1:2154503/1‑69 (MQ=255)
                  ggTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGt                                                 >  1:2724159/1‑69 (MQ=255)
                  ggTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGt                                                 >  1:487370/1‑69 (MQ=255)
                  ggTGCACGCAGCACCACCACCGTGGATATGTCGCATACTTATCTGGAGTGGGCAGAACGCAACCTGCGt                                                 >  1:2658847/1‑69 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGTCTg                                              <  1:3182939/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGTCTg                                              <  1:463484/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGTCTg                                              <  1:3037507/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGTCTg                                              <  1:2031731/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGTCTg                                              <  1:220603/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAGTGGGCAGAACGCAACCTGCGTCTg                                              <  1:2869666/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:3410220/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:295477/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:2877865/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:1144464/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:2271543/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:1861517/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:1737240/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:934201/69‑1 (MQ=255)
                     gCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTg                                              <  1:1303202/69‑1 (MQ=255)
                                          gATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTgcg                         >  1:1246056/1‑69 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:3377355/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:446330/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:1185422/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:520759/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:800893/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:2494776/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:1729291/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGc                   <  1:1102748/69‑1 (MQ=255)
                                                tCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGGGCGCATCGc                   <  1:75150/69‑1 (MQ=255)
                                                                 aTGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTg  <  1:3131876/69‑1 (MQ=255)
                                                                 |                                                                    
GTGCATGCGGGATTAGGCGGTGCACGCAGCACCACCACCGTGGATATGTCGCGTACTTATCTGGAATGGGCAGAACGCAACCTGCGTCTGAATGGCCTGACCGGGCGTGCGCATCGCCTGATTCAGGCCGATTG  >  minE/666202‑666335

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