Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,458 C→T 25.3% A186T (GCG→ACG)  yafV ← predicted C‑N hydrolase family amidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE218,4580CT25.3% 59.8 / 15.0 51A186T (GCG→ACG) yafVpredicted C‑N hydrolase family amidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (24/14);  new base T (5/8);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

AATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCGGAGCAGGCCAGTT  >  minE/218407‑218524
                                                   |                                                                  
aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                                   >  1:935476/1‑69 (MQ=255)
aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                                   >  1:935116/1‑69 (MQ=255)
aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                                   >  1:3127260/1‑69 (MQ=255)
aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                                   >  1:242187/1‑69 (MQ=255)
aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTATCAATcgcgc                                                   >  1:681695/1‑69 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:107064/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:92414/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:471631/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:370127/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:2432832/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:1643731/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:110981/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:1392255/1‑57 (MQ=255)
    gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATACCGCCACATACGCCTGATtctc                                                           >  1:3130115/1‑57 (MQ=255)
    gCAGCCGTTGCCATCG‑‑GCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                                           >  1:3333230/1‑55 (MQ=255)
          gTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAg                                         <  1:1814653/69‑1 (MQ=255)
          gTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAg                                         <  1:352880/69‑1 (MQ=255)
           ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                                        <  1:1816085/69‑1 (MQ=255)
           ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                                        <  1:2756266/69‑1 (MQ=255)
           ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                                        <  1:2124375/69‑1 (MQ=255)
           ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                                        <  1:2545246/69‑1 (MQ=255)
           ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                                        <  1:3145941/69‑1 (MQ=255)
           gttccATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                                        <  1:1832120/66‑1 (MQ=255)
             gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCat                                      <  1:2387548/69‑2 (MQ=255)
             gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                                      <  1:2254609/69‑1 (MQ=255)
             gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                                      <  1:97370/69‑1 (MQ=255)
             gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                                      <  1:2586860/69‑1 (MQ=255)
             gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                                      <  1:1585186/69‑1 (MQ=255)
             gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                                      <  1:2949390/69‑1 (MQ=255)
             gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                                      <  1:203660/69‑1 (MQ=255)
               cATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                                      <  1:3363340/67‑1 (MQ=255)
                  cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                                 >  1:2332408/1‑69 (MQ=255)
                  cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                                 >  1:270197/1‑69 (MQ=255)
                  cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                                 >  1:2089563/1‑69 (MQ=255)
                  cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcc                                  >  1:172949/1‑68 (MQ=255)
                   gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                                 >  1:478954/1‑68 (MQ=255)
                   gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                                 >  1:1184784/1‑68 (MQ=255)
                   gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                                 >  1:2757042/1‑68 (MQ=255)
                   gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                                 >  1:2223624/1‑68 (MQ=255)
                     tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                              <  1:424266/69‑1 (MQ=255)
                     tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                              <  1:2448645/69‑1 (MQ=255)
                     tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                              <  1:119713/69‑1 (MQ=255)
                     tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                              <  1:1533080/69‑1 (MQ=255)
                     tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                              <  1:2311301/69‑1 (MQ=255)
                     tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                              <  1:991211/69‑1 (MQ=255)
                                 tGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAgc                  >  1:3072250/1‑69 (MQ=255)
                                 tGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAgc                  >  1:124377/1‑69 (MQ=255)
                                 agCATCCCGCCACATACGCCTGATTCTCAATCGCGCggg                                                >  1:616046/2‑39 (MQ=255)
                                       ccGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCgg              >  1:1013467/1‑67 (MQ=255)
                                                 cGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCGGAGCAGGCCAGtt  >  1:1296386/1‑69 (MQ=255)
                                                 cGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCGGAGCAGGCCAGt   >  1:1272310/1‑68 (MQ=255)
                                                   |                                                                  
AATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCGGAGCAGGCCAGTT  >  minE/218407‑218524

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: