Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,537,470:1	+A	9.1%	coding (1569/2916 nt)	mrpAB ←	K(+)/H(+) antiporter subunit A/B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,537,470	1	.	A	9.1%	161.2 / 10.8	44	coding (1569/2916 nt)	mrpAB	K(+)/H(+) antiporter subunit A/B
Reads supporting (aligned to +/- strand):  ref base . (21/19);  new base A (2/2);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

TGCAGTCGGCGGGTGGTGAGCACGCGCTCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGCCAACAGGGGGCCGACC  >  NC_002947/2537327‑2537617                                                                                                                                                 |                                                                                                                                                    tGCAGTCGGCGGGTGGTGAGCACGCGCTCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCgc‑a                                                                                                                                                    <  4:217793/145‑1 (MQ=255) tGCAGTCGGCGGGTGGTGAGCACGCGCTCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCgc‑a                                                                                                                                                    >  3:217793/1‑145 (MQ=255)   cAGTCGGCGGGTGGTGAGCACGCGCTCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCgc‑ag                                                                                                                                                   <  8:17357/144‑1 (MQ=255)   cAGTCGGCGGGTGGTGAGCACGCGCTCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCgc‑ag                                                                                                                                                   >  7:17357/1‑144 (MQ=255)   cAGTCGGCGGGTGGTGAGCACGCGCTCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAAc                                                                                                                                               <  3:86907/148‑1 (MQ=255)                            tCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCAt                                                                                                                     >  4:243949/1‑149 (MQ=255)                                    aCGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGAc                                                                                                               >  2:193644/1‑147 (MQ=255)                                          cagcagCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCa                                                                                                             <  4:202367/143‑1 (MQ=255)                                          cagcagCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCa                                                                                                             >  3:202367/1‑143 (MQ=255)                                            gcagcaTCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGAGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCAt                                                                                                      >  5:267529/1‑148 (MQ=255)                                             cagcaTCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGGAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATgc                                                                                                    >  1:70722/1‑148 (MQ=255)                                               gcaTCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGat                                                                                                   >  7:237079/1‑148 (MQ=255)                                                  tCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGatcatc                                                                                               >  7:67763/1‑149 (MQ=255)                                                         aaCCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGt                                                                                        >  4:215798/1‑149 (MQ=255)                                                               gCCGTGTTGGCAAAGTCGCTTGCCATTGTAGCGCTCGATCACGGGCGGGGAGGGGAAGCGGGCGATGCGGAGTTGCTTGCGC‑AGCAACAGGTACAGGTCGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTGCAg                                                                                   <  6:267529/148‑1 (MQ=255)                                                                ccGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGcc                                                                                 <  3:243949/149‑1 (MQ=255)                                                                 cGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCa                                                                                >  5:182619/1‑149 (MQ=255)                                                                  gtgtTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCAt                                                                               <  7:138059/149‑1 (MQ=255)                                                                     ttCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGc                                                                             <  8:67763/148‑1 (MQ=255)                                                                        gAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCAt                                                                                                                     <  2:174259/105‑1 (MQ=255)                                                                        gAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCAt                                                                                                                     >  1:174259/1‑105 (MQ=255)                                                                                    gCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTc                                                                                                        <  7:296285/106‑1 (MQ=255)                                                                                    gCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTc                                                                                                        >  8:296285/1‑106 (MQ=255)                                                                                           aaGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCgt                                                      >  6:223696/1‑148 (MQ=255)                                                                                              cgcTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAg                                                   <  8:237079/149‑1 (MQ=255)                                                                                               gcTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGG‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAgg                                                  <  2:70722/149‑1 (MQ=255)                                                                                               gcTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAgg                                                  >  6:274086/1‑149 (MQ=255)                                                                                                 tCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAggg                                                 >  5:68502/1‑148 (MQ=255)                                                                                                     tCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCTACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCAAGATGGCCAGGCTGTATTCCGGCAGGGTACCg                                            <  1:193644/149‑1 (MQ=255)                                                                                                           gcgggTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCACATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAacc                                      >  5:270992/1‑149 (MQ=255)                                                                                                             gggTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAaccac                                    <  7:257579/149‑1 (MQ=255)                                                                                                               gTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTAc                                              >  1:130496/1‑137 (MQ=255)                                                                                                               gTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTAc                                              <  2:130496/137‑1 (MQ=255)                                                                                                                  ggggAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCACATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCa                               <  6:270992/149‑1 (MQ=255)                                                                                                                        gCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAggg                            >  8:194490/1‑146 (MQ=255)                                                                                                                        gCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAggg                            <  7:194490/146‑1 (MQ=255)                                                                                                                             ccAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCg                    <  6:68502/149‑1 (MQ=255)                                                                                                                              cAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAg                                                                                   <  8:228861/85‑1 (MQ=255)                                                                                                                              cAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAg                                                                                   >  7:228861/1‑85 (MQ=255)                                                                                                                               tgtacctgTTGC‑TGCGCAAGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGc                  >  3:206736/8‑149 (MQ=255)                                                                                                                               aGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCAAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGc                  >  6:3608/1‑149 (MQ=255)                                                                                                                                   cctgTTGC‑TGCGCAAGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAgg                                                  >  1:158851/4‑113 (MQ=255)                                                                                                                                   cctgTTGC‑TGCGCAAGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAgg                                                  <  2:158851/110‑1 (MQ=255)                                                                                                                                     aGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGCCAACAg            >  4:147411/1‑149 (MQ=255)                                                                                                                                       ttgcttgcGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCACCAACCACCGGCAGGGCCGCAGCGGCCAACAggg          <  3:215798/149‑1 (MQ=255)                                                                                                                                          cttgcGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGCCAACAggggg        <  3:147411/148‑1 (MQ=255)                                                                                                                                           ttgcGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGCCAACAGGGGGcc      >  2:249637/1‑149 (MQ=255)                                                                                                                                           ctgcGCAAGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGCCAACAGGGGGc       <  4:206736/148‑1 (MQ=255)                                                                                                                                               gc‑agcaACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGCCAACAGGGGGccgacc  >  5:145781/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TGCAGTCGGCGGGTGGTGAGCACGCGCTCGATGCGACGCGCCAGCAGCATCAGGTGAACCTGGCCGTGTTCGAACAGTCGCTTGCCATTGAAGCGCTCGATCACGGGCGGGTAGGGGAAGCGGCCCAGGCGCAGTTGCTTGCGC‑AGCAACAGGTACAGGACGATGCCGCCGCTCATGGCGACCAGGCTCATGATCATCGGCGCGTTCCAGCCATGCCAGATGGCCAGGCTGTATTCCGGCAGGGTACCGCCAACCACCGGCAGGGCCGCAGCGGCCAACAGGGGGCCGACC  >  NC_002947/2537327‑2537617

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: ‑