Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 522,599 A→C 45.6% P168P (CCT→CCG ybhK ← predicted transferase with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE522,5990AC45.6% 0.0 / 91.3 79P168P (CCT→CCGybhKpredicted transferase with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (19/24);  new base C (10/26);  total (29/50)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

AGGCCCAATAATGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCATG  >  minE/522533‑522667
                                                                  |                                                                    
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:2119934/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:115002/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:1199902/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:2275358/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:361002/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:3405586/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:3308832/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:1414820/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:1610618/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAgg                                                                    >  1:1607956/1‑69 (MQ=255)
aGGCCCAATAATGATGAGATCCGCTACACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCgg                                                                    >  1:163816/1‑69 (MQ=255)
       ataatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAAc                                                             <  1:2785996/69‑1 (MQ=255)
       ataatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAAc                                                             <  1:2224659/69‑1 (MQ=255)
       ataatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAAc                                                             <  1:2300774/69‑1 (MQ=255)
       ataatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAAc                                                             <  1:23348/69‑1 (MQ=255)
       ataatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAAc                                                             <  1:1597833/69‑1 (MQ=255)
       ataatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAAc                                                             <  1:2343927/69‑1 (MQ=255)
       ataatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAAc                                                             <  1:348429/69‑1 (MQ=255)
         aatGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACa                                                            >  1:479419/1‑68 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGc                                                                   >  1:1693886/1‑60 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACaa                                                           >  1:3117180/1‑68 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:2401996/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:2412350/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:1073183/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:3430968/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:370642/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:642318/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:1014387/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:1226887/1‑69 (MQ=255)
          atgatgAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:152525/1‑69 (MQ=255)
          atgatgAGATCCGCTTCACTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAAt                                                          >  1:3243171/1‑69 (MQ=255)
            gatgaGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATaa                                                        >  1:3040219/1‑69 (MQ=255)
            gatgaGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATaa                                                        >  1:1843811/1‑69 (MQ=255)
                             tGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAg                                       <  1:3289420/69‑1 (MQ=255)
                             tGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAg                                       <  1:2139079/69‑1 (MQ=255)
                              gATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAg                                       <  1:1953880/68‑1 (MQ=255)
                              gATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGt                                      <  1:2289435/69‑1 (MQ=255)
                              gATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGt                                      <  1:2186368/69‑1 (MQ=255)
                                  gCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAAc                                  >  1:983560/1‑69 (MQ=255)
                                  gCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAAc                                  >  1:1366319/1‑69 (MQ=255)
                                  gCGTGAACCGCCTCACGAGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAAc                                  >  1:453105/1‑69 (MQ=255)
                                             cTCACGCGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATAtt                       <  1:3359802/69‑1 (MQ=255)
                                              tCACGCGTTGCGGGTACATTAGGCGTTAACAATTACTCTTGAATCGGCGTAGTTAACTGGTCGATAtt                       <  1:1308791/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:3070207/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:3108549/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:3077256/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:3096695/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:3333105/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:1264845/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:1852004/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:395060/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:44154/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:1134679/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:808681/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:918797/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:2991561/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:2282553/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:1755787/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:2893745/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:1549913/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:2890173/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:2877514/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:1459213/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGGCGATATTGAcc                   <  1:1945781/68‑1 (MQ=255)
                                                  gcgTTGCGGGTACATTCCGCGTTAACAATAACTCTTGAATCGGCGAAGTTAACTGGTCGATATTGAcc                   <  1:1594955/68‑1 (MQ=255)
                                                   cgTTGCGGGTACATTCGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGAcc                   <  1:1322781/67‑1 (MQ=255)
                                                          ggTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTa           <  1:3173346/68‑1 (MQ=255)
                                                          ggTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTa           <  1:2241615/68‑1 (MQ=255)
                                                          ggTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTa           <  1:1284365/68‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGTGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:1768625/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:2812358/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:1268362/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:2640034/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:213541/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:2269459/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:162959/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:2307624/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:694079/69‑1 (MQ=255)
                                                              cATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACtt      <  1:858684/69‑1 (MQ=255)
                                                                  aGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCa    >  1:298203/1‑67 (MQ=255)
                                                                  aGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCAt   >  1:2906692/1‑68 (MQ=255)
                                                                  aGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCAt   >  1:2946268/1‑68 (MQ=255)
                                                                  aGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCATg  >  1:3121948/1‑69 (MQ=255)
                                                                  aGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCATg  >  1:897275/1‑69 (MQ=255)
                                                                  aGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCATg  >  1:2209359/1‑69 (MQ=255)
                                                                  |                                                                    
AGGCCCAATAATGATGAGATCCGCTTCATTGATAGCGTGAACCGCCTCACGCGTTGCGGGTACATTAGGCGTTAACAATAACTCTTGAATCGGCGTAGTTAACTGGTCGATATTGACCTCGCCGTAAACTTCATG  >  minE/522533‑522667

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