Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,705,413 T→A 20.7% intergenic (‑151/+26) kdpF ← / ← PP_4162 potassium‑transporting ATPase subunit F/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,705,4130TA20.7% 56.8 / 8.7 29intergenic (‑151/+26)kdpF/PP_4162potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (17/6);  new base A (4/2);  total (21/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGC  >  NC_002947/4705268‑4705557
                                                                                                                                                 |                                                                                                                                                 
gATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGc                                                                                                                                                >  3:99606/1‑149 (MQ=255)
     ccccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTAGCCGGCTTGCCcgc                                                                                                                                            >  6:41686/1‑148 (MQ=255)
                                 ggAAATTGTCCGAGCGAGAGGCGAAAAGGGGCGAGATCGGCGGGTGCGCGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCGTAATGTTGTG                                                                                                               >  8:184784‑M1/1‑149 (MQ=255)
                                  gAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGG                                                                                                              >  7:67858‑M1/1‑149 (MQ=255)
                                  gAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTCCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAggg                                                                                                                                     >  1:68323/1‑126 (MQ=255)
                                                                       ggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagcag                                                                         >  7:80645/1‑148 (MQ=255)
                                                                                             gAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGtaggctacagcctgctgctggactggagccgtaacgcggcgcagcagccgccacaa                                                  <  8:67858‑M1/149‑57 (MQ=255)
                                                                                                     aaaaaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGc                                          >  7:198323/1‑149 (MQ=255)
                                                                                                       aaaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGccc                                        >  3:144377/1‑149 (MQ=255)
                                                                                                          tCAGCGTTGTGGTCTGCCTGTAACCGCCTCTTCGCGGGCTAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTTCGCCGCGTTACGGCTCCAGTCCCGCAGCATGCTGTAGCCTCCCGCCAGCCAGGTCGGCCCGAt                                     >  5:7277/1‑149 (MQ=255)
                                                                                                                   tgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGcc                            <  8:80645/149‑1 (MQ=255)
                                                                                                                      tCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGt                                                                 <  7:177979/109‑1 (MQ=255)
                                                                                                                      tCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGt                                                                 >  8:177979/1‑109 (MQ=255)
                                                                                                                      tCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAAt                         >  1:15338/1‑149 (MQ=255)
                                                                                                                             gTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGAGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCAACCGCCAGCAAGGTCGGCCCGATGAAAAGGACAATGAAGcc                   >  6:176798/1‑148 (MQ=37)
                                                                                                                                aCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                >  7:120640/1‑132 (MQ=255)
                                                                                                                                aCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaaaa                                <  8:120640/132‑1 (MQ=255)
                                                                                                                                  ggccctttCGCGGGCAAGCCCGCGAAGGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAgg              >  7:100394/7‑148 (MQ=39)
                                                                                                                                  ggCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgc             <  2:15338/149‑1 (MQ=255)
                                                                                                                                   gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgc             >  5:107319/1‑148 (MQ=255)
                                                                                                                                    ccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgc             >  1:206082/5‑148 (MQ=38)
                                                                                                                                      tctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaag          >  2:102416/1‑148 (MQ=255)
                                                                                                                                        ttCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGtgggctacagcctgctgctggactggagccgtaacgcggcgcagcagccgccacaacaa                                               >  1:244722‑M1/1‑50 (MQ=255)
                                                                                                                                        ttCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGtgggctacagcctgctgctggactggagccgtaacgcggcgcagcagccgccacaacaa                                               <  2:244722‑M1/109‑60 (MQ=255)
                                                                                                                                         tCGCGGGCAAGCCCGCGAAGGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagcaa       <  8:100394/148‑1 (MQ=255)
                                                                                                                                          cgcgGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagca                                                                          <  1:238669/80‑1 (MQ=255)
                                                                                                                                          cgcgGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagca                                                                          >  2:238669/1‑80 (MQ=255)
                                                                                                                                            cgGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccg   >  3:57229/1‑149 (MQ=255)
                                                                                                                                             cggCTTGCCCGCGAAAGGTCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGTCTCCAGAACAGCAGCATGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATTAAAACGCCAATGAATCCAAATGCAAGCAAgccgc  >  2:101314/2‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGC  >  NC_002947/4705268‑4705557

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