Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 449,567 C→T 48.5% A257A (GCG→GCA kdpB ← potassium translocating ATPase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE449,5670CT48.5% 1.0 / 66.7 77A257A (GCG→GCAkdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base C (16/23);  new base T (21/16);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGC  >  minE/449501‑449629
                                                                  |                                                              
cATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAgcgc                                                              <  1:2258462/69‑1 (MQ=255)
cATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAgcgc                                                              <  1:2796453/69‑1 (MQ=255)
cATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAgcgc                                                              <  1:1346284/69‑1 (MQ=255)
cATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAgcgc                                                              <  1:1133990/69‑1 (MQ=255)
  tCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAgcgc                                                              <  1:1637855/67‑1 (MQ=255)
   cccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccac                                                           <  1:2725633/69‑1 (MQ=255)
   cccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccac                                                           <  1:1674581/69‑1 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:2006427/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:2695091/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:3403152/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:568612/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:1623512/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:1617626/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:1544792/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGccacc                                                          >  1:1525188/1‑69 (MQ=255)
    ccGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACAAGCAGCGccacc                                                          >  1:901211/1‑69 (MQ=255)
      ggCGACGCCAATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAg                                                        >  1:691611/1‑69 (MQ=255)
      ggCGACGCCAATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAg                                                        >  1:338682/1‑69 (MQ=255)
      ggCGACGCCAATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAg                                                        >  1:2999544/1‑69 (MQ=255)
      ggCGACGCCAATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAg                                                        >  1:2858898/1‑69 (MQ=255)
      ggCGACGCCAATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGACACCAg                                                        >  1:70600/1‑69 (MQ=255)
            gCCGATCGCTGACAACTGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGt                                                  <  1:2544364/69‑1 (MQ=255)
            gCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGt                                                  <  1:1707024/69‑1 (MQ=255)
            gCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGt                                                  <  1:589841/69‑1 (MQ=255)
            gCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGt                                                  <  1:2882378/69‑1 (MQ=255)
            gCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGt                                                  <  1:2634947/69‑1 (MQ=255)
            gCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGt                                                  <  1:2230346/69‑1 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACGGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:642869/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGtt                                                 >  1:1482350/3‑68 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGtt                                                 >  1:2417207/3‑68 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:3381938/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:3412433/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:365642/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:41147/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:2498712/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:1222982/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:1745705/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:895915/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:1875607/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:1716504/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:2138529/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:1557914/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGTGCCACCAGTACCGTTa                                                >  1:1950079/3‑69 (MQ=255)
              caaTCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGAGCCACCAGTACCGTTa                                                >  1:2684785/3‑69 (MQ=255)
                  cGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                            <  1:2455865/69‑1 (MQ=255)
                  cGCTGACAACAGGCCGCCAATAGTGGTGGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCt                                            <  1:2610469/69‑1 (MQ=255)
                         aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                       <  1:2256316/67‑1 (MQ=255)
                         aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                       <  1:1876275/67‑1 (MQ=255)
                         aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                       <  1:3370588/67‑1 (MQ=255)
                         aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                       <  1:157497/67‑1 (MQ=255)
                         aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                       <  1:3122401/67‑1 (MQ=255)
                         aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                       <  1:1982292/67‑1 (MQ=255)
                         aacaGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTg                                       <  1:537413/67‑1 (MQ=255)
                                       aGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCAcg                        >  1:1955312/1‑68 (MQ=255)
                                           gTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGaa                   >  1:3100562/1‑69 (MQ=255)
                                                   cagacagaCCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             >  1:1957928/1‑67 (MQ=255)
                                                   cagacagaCCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             >  1:1535839/1‑67 (MQ=255)
                                                   cagacagaCCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             >  1:2335016/1‑67 (MQ=255)
                                                   cagacagaCCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             >  1:2677750/1‑67 (MQ=255)
                                                   cagacagaCCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             >  1:1319743/1‑67 (MQ=255)
                                                    agacagacCTGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:849678/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACTCGGAAAACGGc             <  1:88517/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGACAACGGc             <  1:204/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:43047/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:378994/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:1403140/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:3094908/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:1890880/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:3281792/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:3144886/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:3058913/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:2974709/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:2781632/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:2039374/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:2169983/66‑1 (MQ=255)
                                                    agacagacCAGCAGTGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGc             <  1:2171889/66‑1 (MQ=255)
                                                            cagcagCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:2561665/69‑1 (MQ=255)
                                                            cagcagCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:430457/69‑1 (MQ=255)
                                                            cagcagCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:413852/69‑1 (MQ=255)
                                                            cagcagCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:2554969/69‑1 (MQ=255)
                                                            cagcagCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:2232206/69‑1 (MQ=255)
                                                            cagcagCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:1481602/69‑1 (MQ=255)
                                                             agcagcGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:3221848/68‑1 (MQ=255)
                                                             agcagcGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGc  <  1:915555/68‑1 (MQ=255)
                                                                  |                                                              
CATCCCGGCGACGCCGATCGCTGACAACAGGCCGCCAATAGTGGTTGGGATCAGACAGACCAGCAGCGCCACCAGTACCGTTACGCTGACTGCATTACCGCCCCACGCGGAAAACGGCCACAGCGTGGC  >  minE/449501‑449629

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: