Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,705,414	T→A	21.0%	intergenic (‑152/+25)	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_002947	4,705,414	0	T	A	21.0%	73.6 / 16.8	43	intergenic (‑152/+25)	kdpF/PP_4162	potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (21/12);  new base A (7/2);  total (28/15)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.23e-01

CCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGC  >  NC_002947/4705273‑4705557                                                                                                                                               |                                                                                                                                                ccccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCAAGCCcgc                                                                                                                                           >  1:97362/1‑148 (MQ=255)     gtgtTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAg                                                                                                                                      >  6:384020/1‑149 (MQ=255)            ggCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGt                                                                                                                                >  1:298276/1‑148 (MQ=255)                  ggCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                                                                                         >  1:7462/1‑149 (MQ=255)                   gCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCt                                                                                                                        >  7:301932/1‑149 (MQ=255)                             ggAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCACGCGAAAGGGCCGTTACAGGCTTATtcttgtt                                                                                                              >  4:252629/1‑148 (MQ=255)                                             agagGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGcgccgc                                                                                              >  8:48445/1‑149 (MQ=255)                                                     aaaGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGAGGGCTTGCCCGCGAAATGGCCGATACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACgg                                                                                       >  2:325031/1‑148 (MQ=255)                                                         gggCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTATTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc                                                                                    >  2:319905/1‑147 (MQ=255)                                                              tctgtggggggggggTGTGCGTCATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGAGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGTTGCTGCGCCGCGTTACGGCtccagtcc                                                                              <  8:301932/144‑1 (MQ=255)                                                              agaTCGGGGGGTGGGGGTGGGT‑ATAAAGAAAGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCATCGCGGGCTTGCCCGCGAAAGGGCAGTTAAAGGCTTATTGATGTGGCGGCTCCTGCTCCGCCTTACGGCtcccctcac                                                                             >  4:389497/1‑143 (MQ=255)                                                                aTCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtccag                                                                            >  3:131567/1‑148 (MQ=255)                                                                        gtgggggtgggT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCt                                                                   >  4:105888/1‑149 (MQ=255)                                                                                 gtt‑attAAGAATGCGTAATAATCTGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTGGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCGAGCAGCAGGCTGTAGCCTAc                                                          <  1:421824/147‑1 (MQ=255)                                                                                        aaGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGTGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCATTCCAGCAGCAGGCTGTAGCCTACCGCCa                                                     <  6:180072/148‑1 (MQ=255)                                                                                        aaGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAg                                                    >  3:152693/1‑149 (MQ=255)                                                                                                   aaaaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTGCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGcc                                         <  4:152693/149‑1 (MQ=255)                                                                                                   aaaaTCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGcc                                         <  7:48445/149‑1 (MQ=255)                                                                                                      aTCAGCGTGGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGa                                      <  3:252629/149‑1 (MQ=255)                                                                                                       tCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCCCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                     <  1:319905/149‑1 (MQ=255)                                                                                                        cAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGAt                                     <  2:427020/148‑1 (MQ=255)                                                                                                         aGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATg                                    >  3:307458/1‑148 (MQ=255)                                                                                                         aGCGTTGTGGCCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATg                                    >  3:262272/1‑148 (MQ=255)                                                                                                                tgGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGc                             >  1:31523/1‑148 (MQ=255)                                                                                                                        cTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGc                    >  5:65776/1‑149 (MQ=38)                                                                                                                          gTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGACTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc                   >  1:367073/1‑148 (MQ=255)                                                                                                                           tAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaa                                  >  6:127406/1‑132 (MQ=255)                                                                                                                           tAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGaa                                  <  5:127406/132‑1 (MQ=255)                                                                                                                             aCGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTTGGCCCGATGAAAAGGCCAATGAAGCCAAAg               >  4:391942/1‑149 (MQ=255)                                                                                                                             aCGGCCCTTTCGCGGGCAAGCCCGCGAAGGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAg               >  5:279978/1‑149 (MQ=39)                                                                                                                               ggCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAgg              >  3:390148/1‑148 (MQ=255)                                                                                                                                gccctttCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCCGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAACGGCCAATGAAGCCAAAGgc             >  2:419357/6‑148 (MQ=38)                                                                                                                                 ccTCTTCGCGGGCTTGCCCGCTAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaa           <  4:390181/149‑1 (MQ=255)                                                                                                                                  ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGcc                                                      <  3:342994/105‑1 (MQ=255)                                                                                                                                  ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGcc                                                      >  4:342994/1‑105 (MQ=255)                                                                                                                                  ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCAAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaag          >  6:359339/1‑149 (MQ=255)                                                                                                                                   ctttCGCGGGCAAGCCCGCGAAGGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATg                        >  1:228952/3‑134 (MQ=255)                                                                                                                                   ctttCGCGGGCAAGCCCGCGAAGGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATg                        <  2:228952/132‑1 (MQ=255)                                                                                                                                    tttCGCGGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaagc         >  4:297691/2‑148 (MQ=255)                                                                                                                                        gcgGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc                                                                                    >  4:324425/1‑69 (MQ=255)                                                                                                                                        gcgGGCAAGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtc                                                                                    <  3:324425/69‑1 (MQ=255)                                                                                                                                          gggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccgc  <  3:105888/149‑1 (MQ=255)                                                                                                                                          gggCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCAGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAgccgc  <  2:298276/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                CCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGT‑ATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCAAGCCGC  >  NC_002947/4705273‑4705557

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