Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	47,776	G→A	*177* (TAG→TAA)  S3N (AGC→AAC)	kefF → kefC →	potassium‑efflux system ancillary protein for KefC, glutathione‑regulated, quinone oxidoreductase, FMN‑dependent potassium:proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	47,776	0	G	A	98.0%	152.6 / NA	49	*177* (TAG→TAA)  S3N (AGC→AAC)	kefF kefC	potassium‑efflux system ancillary protein for KefC, glutathione‑regulated, quinone oxidoreductase, FMN‑dependent potassium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/23);  minor base T (0/1);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGAACTGGCTGCCACCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGGCCCGTGGGGG  >  NC_000913/47653‑47903                                                                                                                            |                                                                                                                                gCTGAACTGGCTGCCACCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCa                                                                                                                  <  1:39814/139‑1 (MQ=255)   tGAACTGGCTGCCACCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCa                                                                                                                              <  2:340995/125‑1 (MQ=255)   tGAACTGGCTGCCACCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCa                                                                                                                              >  1:340995/1‑125 (MQ=255)     aaCTGGCTGCCACCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGgcg                                                                                                              <  1:989128/139‑1 (MQ=255)      aCTGGCTGCCACCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGgcgc                                                                                                             <  1:389002/139‑1 (MQ=255)                aCCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTAtct                                                                                                   <  1:331805/139‑1 (MQ=255)                            gcactgcacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAgc                                                                                       <  2:995413/139‑1 (MQ=255)                                tgcacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAgcgc                                                                                     <  2:777544/137‑1 (MQ=255)                                tgcacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAgcgc                                                                                     >  1:777544/1‑137 (MQ=255)                                tgcacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTg                                                                                   >  2:959014/1‑139 (MQ=255)                                  cacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGgcg                                                                                                              >  1:289788/1‑110 (MQ=255)                                  cacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGgcg                                                                                                              <  2:289788/110‑1 (MQ=255)                                  cacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCCGCGCGGAt                                                                                 >  2:83308/1‑139 (MQ=255)                                  cacCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGAt                                                                                 >  1:278559/1‑139 (MQ=255)                                                acgacgAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCgg                                                                   >  2:512424/1‑139 (MQ=255)                                                      aaaCCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTc                                                             <  2:404199/139‑1 (MQ=255)                                                       aaCCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCt                                                            <  2:278559/139‑1 (MQ=255)                                                            tCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTc                                                       >  1:135844/1‑139 (MQ=255)                                                              gAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTg                                                     >  1:530111/1‑139 (MQ=255)                                                                    cAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCg                                               >  1:718197/1‑139 (MQ=255)                                                                          cgTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACtt                                         <  1:959014/139‑1 (MQ=255)                                                                                       aaCGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCAt                                                                        >  1:633487/1‑95 (MQ=255)                                                                                       aaCGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCAt                                                                        <  2:633487/95‑1 (MQ=255)                                                                                       aaCGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCg                            >  1:1058879/1‑139 (MQ=255)                                                                                         cGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGcc                          <  2:135844/139‑1 (MQ=255)                                                                                         cGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGcc                          >  2:222333/1‑139 (MQ=255)                                                                                          gTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGc                                                                                          <  2:572459/74‑1 (MQ=255)                                                                                          gTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGc                                                                                          >  1:572459/1‑74 (MQ=255)                                                                                          gTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCAt                                                                        <  1:154583/92‑1 (MQ=255)                                                                                          gTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCAt                                                                        >  2:154583/1‑92 (MQ=255)                                                                                          gTCTGCTGGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCg                         >  1:542927/1‑139 (MQ=255)                                                                                            ctgctgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGc                       <  2:718197/139‑1 (MQ=255)                                                                                              gctgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAgcg                                                                                      <  2:822540/74‑1 (MQ=255)                                                                                              gctgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAgcg                                                                                      >  1:822540/1‑74 (MQ=255)                                                                                              gctgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGt                                                                              <  2:800580/82‑1 (MQ=255)                                                                                              gctgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGt                                                                              >  1:800580/1‑82 (MQ=255)                                                                                              gctgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTg                     <  2:530111/139‑1 (MQ=255)                                                                                              gctaGAATGGCAGGGGAGCCATCATGGATATCCATACGCTGATTCAAGCGCTGATTTATGTCGGTTCGGCAGCGCTGATTGTACCCTTTGCGGTACGTCTTGGTGTGGGATCGGTAGTTGGCTACCTGATCGCCGGCTg                     <  2:646585/139‑1 (MQ=255)                                                                                                tgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGcc                          <  2:555640/132‑1 (MQ=255)                                                                                                tgGAATGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGcc                          >  1:555640/1‑132 (MQ=255)                                                                                                   aaTGGCAGGAGGCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCatta                >  1:1258731/1‑139 (MQ=255)                                                                                                            aggCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGt                                                                  <  2:1269366/80‑1 (MQ=255)                                                                                                            aggCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGt                                                                  >  1:1269366/1‑80 (MQ=255)                                                                                                            aggCCCATCATGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGGCCCGt       >  1:1174975/1‑139 (MQ=255)                                                                                                                 catcatGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGcc                          >  1:90713/1‑115 (MQ=255)                                                                                                                 catcatGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGcc                          <  2:90713/115‑1 (MQ=255)                                                                                                                 catcatGGATAACCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGGCCCGTggggg  >  2:1019957/1‑139 (MQ=255)                                                                                                                  atcatGGATAACCATACGCTGATTCAGGCGCTGAtt                                                                                                        >  1:988959/1‑36 (MQ=37)                                                                                                                  atcatGGATAACCATACGCTGATTCAGGCGCTGAtt                                                                                                        <  2:988959/36‑1 (MQ=37)                                                                                                                            |                                                                                                                                GCTGAACTGGCTGCCACCGTTTGCCATGCACTGCACCTTTATTTGTGACGACGAAACCCTCGAAGGGCAGGCGCGTCACTATAAGCAACGTCTGCTGGAATGGCAGGAGGCCCATCATGGATAGCCATACGCTGATTCAGGCGCTGATTTATCTCGGTTCGGCAGCGCTGATTGTACCCATTGCGGTACGTCTTGGTCTGGGATCGGTACTTGGCTACCTGATCGCCGGCTGCATTATTGGCCCGTGGGGG  >  NC_000913/47653‑47903

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑