Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,705,448	T→C	17.0%	intergenic (‑140/+37)	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,448	0	T	C	17.0%	81.6 / 5.7	39	intergenic (‑140/+37)	kdpF/PP_4162	potassium‑transporting ATPase subunit F/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (15/17);  new base C (2/5);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCA  >  NC_002947/4705308‑4705597                                                                                                                                             |                                                                                                                                                      ggtgggaaTCTCCCCGTGTGACGGCTTGGGCGGTTTTGGGGAAGTGGGCGGAGGGAGGGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCggg                                                                                                                                              <  2:153992/143‑1 (MQ=255)         tctcCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGcccc                                                                                                                                      >  2:626577/1‑149 (MQ=255)             cggatcttAGGGCTTGGGCCGCTTTGGGGAAATGGTCGGGGGGAGAGGGGAAAGGGGGCGAGATGGGGGGGTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCATTCGCGGGCAAGCCCGCGaa                                                                                                                                  <  2:420232/144‑1 (MQ=255)             cccGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaa                                                                                                                                  >  1:600947/1‑150 (MQ=255)              ccGTGTTCCGGCTGGGGCCGCTTTGGGGAAATTGTCCGAGGGAGGGGCTAAAAGGGGCGAGATGGGGGGGTGGGGCTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                 <  1:99281/150‑1 (MQ=255)                 tgtTACGGCTTGGGCCGCTTTGGGGAATTTGCCCGAGGGAGAGGCGTAAAGGGGCGAGAGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGaagagg                                                                                                                              <  1:80496/150‑5 (MQ=37)                   ttAGGGCGTGGGCGGTGTTGGGGAAATTGCCCGAAGGAGGGGTGTAAGGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCTTAAAAATTAGCTTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGcc                                                                                                                            <  2:275727/150‑1 (MQ=255)                     cggccTTGGGCCGTTTGGGGGAAATTGTCGGAGGAAGAGGCGAAAAGGGGCGAGAGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGt                                                                                                                          <  1:442719/146‑1 (MQ=255)                            gggCGGCTGGGGGGACGTTGGCAGGGGGGGCGGGGTAGAGGGGGGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGc                                                                                                                   <  1:1118083/150‑1 (MQ=255)                             ggcggcTTGGGGGAAATTGGCCGGGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGca                                                                                                                  <  2:600947/150‑2 (MQ=25)                             ggCCGGTTGGGGGAAATTGTCGGGGGGGGAGGGGAAAAGGGGCGAGAGGGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCt                                                                                                                  <  2:1336674/150‑1 (MQ=255)                              gCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGGTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCtt                                                                                                                 >  1:989490/1‑150 (MQ=255)                                      ggagAAACTGTGCGCAGGAGTGGGGTATAGGGGTGTGTTCGGGGTGGGGAGGTCGGTATAAAGAATGCGTAAACATCCGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATtgttgt                                                                                                         <  1:103304/147‑1 (MQ=255)                                         gAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCgg                                                                                                                                               >  1:544167/1‑109 (MQ=255)                                         gAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCgg                                                                                                                                               <  2:544167/109‑1 (MQ=255)                                         gAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTAttgagatcgg                                                                                                      >  1:1045948/1‑143 (MQ=255)                                                tCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctgct                                                                                               >  1:188440/1‑150 (MQ=255)                                                tCCGAGGGAGAGGAGTAAAGGGGCGACGGCGGGGGGTGGGGGTGGGAATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgctgct                                                                                               <  1:551396/150‑1 (MQ=255)                                                   gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAg                                                                                                                     <  2:963133/125‑1 (MQ=25)                                                   gAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAg                                                                                                                     >  1:963133/1‑125 (MQ=25)                                                     gggAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGATGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCCTATTGTTGTGGCGGCTGCTGcgcc                                                                                          >  2:675948/1‑150 (MQ=255)                                                                              ggggggTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagcagca                                                                 <  2:188440/150‑1 (MQ=255)                                                                              ggggggTGGGGGGGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCcagcagca                                                                 <  2:737451/150‑1 (MQ=255)                                                                                ggggggggggggTGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAgg                                                               <  1:1064434/150‑1 (MQ=255)                                                                                       gggtgggtATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGc                                                        <  2:1063553/150‑1 (MQ=255)                                                                                                taAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAgg                                                               >  1:590573/1‑134 (MQ=255)                                                                                                taAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGCCGGCTGCTGCTCCGCGTTACGGCTCCAGTCCAGCAGCAgg                                                               <  2:590573/134‑1 (MQ=255)                                                                                                      aTGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgct                                                                                                  <  2:870097/93‑1 (MQ=255)                                                                                                      aTGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGgct                                                                                                  >  1:870097/1‑93 (MQ=255)                                                                                                                              cTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGGTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGGCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGa                 >  2:762257/1‑150 (MQ=255)                                                                                                                                  cTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGcc             >  1:1217710/1‑150 (MQ=255)                                                                                                                                      aaCGGCCCTTTCGCGGGCAAGCCCGCGAAGAGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTCACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAg         <  2:813802/150‑1 (MQ=38)                                                                                                                                         ggCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATAGTTGTGGCGGCTGCCGCACCGGATTACGGCTCCAATCCAGCAGCACGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATCAAGCCAAAGgcc      >  2:520187/1‑149 (MQ=255)                                                                                                                                          gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTtgtg                                                                                                        <  1:1139958/51‑1 (MQ=255)                                                                                                                                          gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTtgtg                                                                                                        >  2:1139958/1‑51 (MQ=255)                                                                                                                                          gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTTGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGGCGGCCCGATGAAAAGGCCCATGAAGCCAAAGgcaa     >  2:566893/1‑150 (MQ=255)                                                                                                                                          gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                >  1:390813/1‑139 (MQ=255)                                                                                                                                          gCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGaa                <  2:390813/139‑1 (MQ=255)                                                                                                                                           ccTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGgcaag    >  2:591651/1‑150 (MQ=255)                                                                                                                                            ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                       <  1:1207376/82‑1 (MQ=255)                                                                                                                                            ctctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCtccagtcca                                                                       >  2:1207376/1‑82 (MQ=255)                                                                                                                                             tctTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCCGCTGCTACGCCGCGTTACGGCTCCAGTCAAGCAGCAGGCGGTAGCCTACCGCCAGAAAGGTCGGGCCGCTGAAAAGGCCAATGAAGCCAAAGgcaagca  >  2:431954/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                                      TGTGTGGATCTCCCCGTGTTACGGCTTGGGCCGCTTTGGGGAAATTGTCCGAGGGAGAGGCGTAAAGGGGCGAGATCGGGGGGTGGGGGTGGGTATAAAGAATGCGTAAAAATCAGCGTTGTGGTCTGCCTGTAACGGCCTCTTCGCGGGCTTGCCCGCGAAAGGGCCGTTACAGGCTTATTGTTGTGGCGGCTGCTGCGCCGCGTTACGGCTCCAGTCCAGCAGCAGGCTGTAGCCTACCGCCAGCAAGGTCGGCCCGATGAAAAGGCCAATGAAGCCAAAGGCAAGCA  >  NC_002947/4705308‑4705597

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