Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 17,076 C→A 21.4% A349A (GCC→GCA nhaA → sodium‑proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE17,0760CA21.4% 145.2 / 28.0 79A349A (GCC→GCAnhaAsodium‑proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base A (4/13);  total (31/48)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

TATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTC  >  minE/17013‑17141
                                                               |                                                                 
tATGGTGGTGGTGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCttt                                                                >  1:56810/1‑67 (MQ=255)
tATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCttt                                                                >  1:2102275/1‑67 (MQ=255)
      ggttgGTATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCg                                                         <  1:2440428/68‑1 (MQ=255)
      ggttgGTATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCg                                                         <  1:2108057/68‑1 (MQ=255)
      ggttgGTATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCg                                                         <  1:677473/68‑1 (MQ=255)
      ggttgGTATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCg                                                         <  1:651897/68‑1 (MQ=255)
      ggttgGTATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCg                                                         <  1:560874/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:1037460/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:207653/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:1850918/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:1780430/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:389363/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:521896/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:1449988/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:1392922/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:1282032/68‑1 (MQ=255)
                       gTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:747742/68‑1 (MQ=255)
                        tATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGAtt                                        <  1:3311563/67‑1 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:1650526/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:249347/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:2254709/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:3201694/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:3235694/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:2157640/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:710750/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:3337231/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:1907821/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:3406761/1‑69 (MQ=255)
                          tCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACt                                    >  1:665066/1‑69 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTATATCCAGAACTGATTAACTgg                                  >  1:2349659/1‑69 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTgg                                  >  1:1758035/1‑69 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:758659/68‑1 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:3418841/68‑1 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:2762865/68‑1 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:2638026/68‑1 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:938155/68‑1 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   >  1:2320536/1‑68 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:1155744/68‑1 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:1464057/68‑1 (MQ=255)
                            ggTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:2204968/68‑1 (MQ=255)
                             gTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTg                                   <  1:3201854/67‑1 (MQ=255)
                               tttACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCGTAg                                                      <  1:1035529/46‑1 (MQ=255)
                               tttACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCGTAg                                                      <  1:936350/46‑1 (MQ=255)
                               tttACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCGTAg                                                      <  1:811007/46‑1 (MQ=255)
                               tttACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCGTAg                                                      <  1:2701732/46‑1 (MQ=255)
                               tttACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCGTAg                                                      <  1:1583024/46‑1 (MQ=255)
                               tttACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGCGTAg                                                      <  1:3173079/46‑1 (MQ=255)
                               tttACTATGTCTATCTTTATTGCCAGCCTGGCATTTGGTAGAGTAg                                                      <  1:3099310/46‑1 (MQ=39)
                                      tGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTa                                       >  1:2353757/1‑54 (MQ=255)
                                      tGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTa                                       >  1:1222041/1‑54 (MQ=255)
                                            tCTTTATTGCCAGCCTGGCCTTTGGTAGCGTATATCCAGAACTGATTAACTGGGCGAAACTCGGTATcc                  >  1:911802/1‑69 (MQ=255)
                                            tCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATcc                  >  1:848803/1‑69 (MQ=255)
                                            tCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATcc                  >  1:266258/1‑69 (MQ=255)
                                            tCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATcc                  >  1:681836/1‑69 (MQ=255)
                                            tCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATcc                  >  1:412377/1‑69 (MQ=255)
                                            tCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATcc                  >  1:1218646/1‑69 (MQ=255)
                                             cTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCt                 >  1:1972356/1‑69 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTc                         <  1:3201980/60‑1 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTc                         <  1:791394/60‑1 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTc                         <  1:780661/60‑1 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTc                         <  1:2727489/60‑1 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTc                         <  1:1733523/60‑1 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTc                         <  1:2171134/60‑1 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCt                 >  1:179990/1‑68 (MQ=255)
                                              tttATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTg                >  1:2386022/1‑69 (MQ=255)
                                              tttATTGCCAGCCTCGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTc                         <  1:913905/60‑1 (MQ=255)
                                                tATTGCCAGCCTGGCATTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCAAAATTAGGTATCCTTGt              >  1:1435507/1‑69 (MQ=255)
                                                tATTGCCAGCCTGGCATTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCAAAATTAGGTATCCTTGt              >  1:1108233/1‑69 (MQ=255)
                                                tATTGCCAGCCTGGCATTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCAAAATTAGGTATCCTTGt              >  1:1091382/1‑69 (MQ=255)
                                                tATTGCCAGCCTGGCATTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCAAAATTAGGTATCCTTGt              >  1:2513817/1‑69 (MQ=255)
                                                       aGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTc        <  1:2149433/68‑1 (MQ=255)
                                                          ttGGCATTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTc        <  1:2331157/64‑1 (MQ=255)
                                                            ggCCTTTGGTAGCGTAGATCCAGAACTTATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTAtctc  <  1:1447197/69‑1 (MQ=255)
                                                            ggCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTAtctc  <  1:1892985/69‑1 (MQ=255)
                                                            ggCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTAtctc  <  1:3131275/69‑1 (MQ=255)
                                                            ggCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTAtctc  <  1:855350/69‑1 (MQ=255)
                                                            ggCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTAtctc  <  1:304654/69‑1 (MQ=255)
                                                            ggCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTAtctc  <  1:2258609/69‑1 (MQ=255)
                                                            ggCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTAtctc  <  1:648438/69‑1 (MQ=255)
                                                               |                                                                 
TATGGTGGTGGGGATCCTGTGCGGTATCGGTTTTACTATGTCTATCTTTATTGCCAGCCTGGCCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTC  >  minE/17013‑17141

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: