Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,151,144 G→A 51.9% A295A (GCC→GCT ftsH ← protease, ATP‑dependent zinc‑metallo

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,151,1440GA51.9% ‑1.3 / 70.9 79A295A (GCC→GCTftsHprotease, ATP‑dependent zinc‑metallo
Reads supporting (aligned to +/- strand):  ref base G (5/33);  new base A (21/20);  total (26/53)
Fisher's exact test for biased strand distribution p-value = 3.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGAAACGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGT  >  minE/2151087‑2151209
                                                         |                                                                 
gTCGAAACGGCCAGGACGCTGCAGGGCGGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgat                                                        <  1:3242589/69‑1 (MQ=255)
gTCGAAACGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgat                                                        <  1:2646747/69‑1 (MQ=255)
gTCGAAACGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgat                                                        <  1:804876/69‑1 (MQ=255)
gTCGAAACGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgat                                                        <  1:2143470/69‑1 (MQ=255)
gTCGAAACGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgat                                                        <  1:1968894/69‑1 (MQ=255)
 tCGAAACGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgat                                                        <  1:2689507/68‑1 (MQ=255)
      aCGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgatgat                                                     <  1:272643/66‑1 (MQ=255)
      aCGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgatgat                                                     <  1:800242/66‑1 (MQ=255)
      aCGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACgatgat                                                     <  1:1744842/66‑1 (MQ=255)
           cAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:2614649/69‑1 (MQ=255)
           cAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:273054/69‑1 (MQ=255)
           cAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:2835496/69‑1 (MQ=255)
           cAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:300962/69‑1 (MQ=255)
           cAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:2112300/69‑1 (MQ=255)
           cAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:3322638/69‑1 (MQ=255)
           cAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:1323785/69‑1 (MQ=255)
            aGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGt                                             <  1:2200323/68‑1 (MQ=255)
              gACGCAGCAGGGCGGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATTATACCTTCGTTAc                                          <  1:136941/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:2423278/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:139755/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:2130268/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:1427180/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:2877698/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:1126349/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:3053490/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:1886916/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:1824979/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:1688755/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:1795909/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:3363403/69‑1 (MQ=255)
              gACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:2468513/69‑1 (MQ=255)
               aCGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:857692/68‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:98417/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:939091/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:1107435/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:2618203/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:492574/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:3027200/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:3036588/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:775725/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:312784/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:623572/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:3373111/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:3420295/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:241984/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:1814686/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:1514390/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:2033551/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:2081177/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:1325717/54‑1 (MQ=255)
                   agcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:113039/54‑1 (MQ=255)
                    gcagGGCTGGGTCGAGAACGTCCGGACGGTTAGTCGCAGCGATAACGATGATa                                                    <  1:2413949/53‑1 (MQ=255)
                      agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                   >  1:2841183/1‑68 (MQ=255)
                      agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                   >  1:219704/1‑68 (MQ=255)
                      agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                   >  1:88011/1‑68 (MQ=255)
                      agGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGaa                                   >  1:1945970/1‑68 (MQ=255)
                               tCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGccatccat                          >  1:3138939/1‑68 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:855419/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1164299/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1052954/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:787225/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1502162/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1747011/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1750831/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:9883/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1751417/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1769921/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:3320813/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1835606/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:3237479/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:3188291/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:1966813/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:2182857/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:2246261/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:2735019/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:2381855/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTACCTTCAAAACCATCCATCTCAAcc                   >  1:2717140/1‑69 (MQ=255)
                                     aCGTCCGGACGGTTAGTCGCAGCGATAACGATGATACCTTCGTTAACTTCAAAACCATCCATCTCAAcc                   >  1:723580/1‑69 (MQ=255)
                                              cGGTTAGTCGCGGCGATAACGATGATACCTTCGTTAc                                          <  1:2445456/37‑1 (MQ=255)
                                                      cgcgGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGt  >  1:1132630/1‑69 (MQ=255)
                                                      cgcgGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGt  >  1:2283953/1‑69 (MQ=255)
                                                         |                                                                 
GTCGAAACGGCCAGGACGCAGCAGGGCCGGGTCGAGAACGTCCGGACGGTTAGTCGCGGCGATAACGATGATACCTTCGTTACCTTCGAAGCCATCCATCTCAACCAGCATCTGGTTCAGAGT  >  minE/2151087‑2151209

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: