Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 218,450 C→T 26.4% V188V (GTG→GTA 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,4500CT26.4% 54.4 / 15.7 49V188V (GTG→GTAyafVpredicted C‑N hydrolase family amidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (22/14);  new base T (5/8);  total (27/22)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

ATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCGG  >  minE/218385‑218512
                                                                 |                                                              
aTCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCcaca                                                             >  1:2335597/1‑69 (MQ=255)
aTCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCcaca                                                             >  1:413118/1‑69 (MQ=255)
aTCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCcaca                                                             >  1:907866/1‑69 (MQ=255)
aTCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCcaca                                                             >  1:2651595/1‑69 (MQ=255)
 tCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCcaca                                                             <  1:804981/68‑1 (MQ=255)
                      aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                       >  1:242187/1‑69 (MQ=255)
                      aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                       >  1:935476/1‑69 (MQ=255)
                      aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                       >  1:935116/1‑69 (MQ=255)
                      aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATcgcgc                                       >  1:3127260/1‑69 (MQ=255)
                      aaTGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTATCAATcgcgc                                       >  1:681695/1‑69 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:2432832/1‑57 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:370127/1‑57 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:471631/1‑57 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:1643731/1‑57 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:92414/1‑57 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:1392255/1‑57 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:110981/1‑57 (MQ=255)
                          gCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATtctc                                               >  1:107064/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:2756266/69‑1 (MQ=255)
                                 ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                            <  1:2545246/69‑1 (MQ=255)
                                 ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                            <  1:2124375/69‑1 (MQ=255)
                                 ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                            <  1:3145941/69‑1 (MQ=255)
                                 ttGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                            <  1:1816085/69‑1 (MQ=255)
                                 gttccATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGc                            <  1:1832120/66‑1 (MQ=255)
                                   gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCat                          <  1:2387548/69‑2 (MQ=255)
                                   gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                          <  1:2586860/69‑1 (MQ=255)
                                   gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                          <  1:97370/69‑1 (MQ=255)
                                   gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                          <  1:2949390/69‑1 (MQ=255)
                                   gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                          <  1:2254609/69‑1 (MQ=255)
                                   gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                          <  1:1585186/69‑1 (MQ=255)
                                   gCCATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                          <  1:203660/69‑1 (MQ=255)
                                     cATCGCTGCCGACGCGATTGCATCCCGCTACATACGTCTGATTCTCAATCGCGCGGGCCGTCAGCaa                          <  1:3363340/67‑1 (MQ=255)
                                        cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                     >  1:270197/1‑69 (MQ=255)
                                        cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                     >  1:2332408/1‑69 (MQ=255)
                                        cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                     >  1:2089563/1‑69 (MQ=255)
                                        cGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcc                      >  1:172949/1‑68 (MQ=255)
                                         gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                     >  1:2757042/1‑68 (MQ=255)
                                         gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                     >  1:1184784/1‑68 (MQ=255)
                                         gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                     >  1:2223624/1‑68 (MQ=255)
                                         gCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcct                     >  1:478954/1‑68 (MQ=255)
                                           tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                  <  1:991211/69‑1 (MQ=255)
                                           tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                  <  1:119713/69‑1 (MQ=255)
                                           tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                  <  1:1533080/69‑1 (MQ=255)
                                           tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                  <  1:424266/69‑1 (MQ=255)
                                           tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                  <  1:2311301/69‑1 (MQ=255)
                                           tGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAAtgcctgcc                  <  1:2448645/69‑1 (MQ=255)
                                                       tGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAgc      >  1:124377/1‑69 (MQ=255)
                                                       tGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAgc      >  1:3072250/1‑69 (MQ=255)
                                                       agCATCCCGCCACATACGCCTGATTCTCAATCGCGCggg                                    >  1:616046/2‑39 (MQ=255)
                                                             ccGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCgg  >  1:1013467/1‑67 (MQ=255)
                                                                 |                                                              
ATCACCCGGCTGTCACCGCGATAATGGCAGCCGTTGCCATCGCTGCCGACGCGATTGCATCCCGCCACATACGCCTGATTCTCAATCGCGCGGGCCGTCAGCAATGCCTGCCAGTGCAGAGAGCGCGG  >  minE/218385‑218512

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: