Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	49,035	+G	coding (1267/1863 nt)	kefC →	K(+) : H(+) antiporter KefC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	49,035	1	.	G	95.5%	58.5 / ‑2.3	22	coding (1267/1863 nt)	kefC	K(+) : H(+) antiporter KefC
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (10/11);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

CTTTGTGGTATTTGGCGCGGCGCAGATGGCGAATGTGCTGGAGCCGGAGTGGGCGAAATCGCTGACCCTGGCGGTGGCGCTGTCGATGGCAGCAACGCCGATTCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GG‑TGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGATGGTGAAAGAACATTTCCCGCATTTGCAGATTATTGCCCGCGCCCGCGATGTCGACCACTACATTC  >  NC_000913/48785‑49281                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       ctttgtgccatctGGCGCGGCGCAGATGGCGAATGCGCTGGAGACGGAGTGGGCGAAATCGCTGAACCTGGCGGTGGCGCTGTCGATGGCAGCAACGCCGATTCTGCTGGTGATCCTCACTCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGACGCCGATGAGATCGCCGAAGAACCGCCGCGAGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGCGGTTTACTGCTCTCCAGCgg‑gg                                                                                                                                                                                                                                                         >  2:122759/13‑251 (MQ=255)       ggTATGTGGCGCGGCGCAGATGGCGAATGTGCTTGAGCCGGAGGGGGCGAAATCGCTGCCCCGGGCGGTTGCGCTGTCGATGGCAGCAACGCCGATTCTTCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGATATCGACGAAGAACATCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGaaa                                                                                                                                                                                                                                                   <  2:73919/251‑1 (MQ=255)                         gATGGCGAATGGGCTGGGGCCGGAGTGGGCGAAATCGCTGACCCTGGCGGTGGCGCTGTCGATGGCAGCAACGCCGATTCTGCTGGTGATCCGCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCAc                                                                                                                                                                                                                                 <  2:25742/251‑1 (MQ=255)                                        ggAGCCGGAGTGGGCGAAATCGCGGACCCTGGCGGTGGCGCTGTCGATGGCAGCAACGCCGATTCGGCTGGTGATCCTCAAGCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACGCGATCACTATCCGGACCATATc                                                                                                                                                                                                                  <  2:42608/251‑1 (MQ=255)                                                                      tgccgttgcgcTGTCGATGGCAGCATCGCCGATTCTGCTGGTGATTCTCAATCGCCTTGAGCAATCCTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGTGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGTACGTTTACTGCTCTCCAGCGGTGG‑TGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAgc                                                                                                                                                                                     <  1:13682/243‑2 (MQ=255)                                                                            ggCGCTGTCGATGGCAGCAACGCCGATTCTGCTGGTGATCCTCAATCGCCTTGAGAAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTAt                                                                                                                                                                              >  1:152107/1‑251 (MQ=255)                                                                                             caACGCCGATTCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCggat                                                                                                                                                             <  1:100228/251‑1 (MQ=255)                                                                                                                     tCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGCATCTGCCGGAgc                                                                                                                                     >  2:272744/1‑251 (MQ=255)                                                                                                                                           tACTGAGGCAGCGGGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGAtt                                                                                                               >  2:176314/1‑251 (MQ=255)                                                                                                                                                                cGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCa                                                                                           >  1:106484/1‑250 (MQ=255)                                                                                                                                                                    gagaTCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAaacc                                                                                       >  1:158735/1‑250 (MQ=255)                                                                                                                                                                         cgacgaAGAACAGCCGCGCGTGATTATCGCCGGCTTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGCCGATCCGCAAACCAACCTg                                                                                 >  2:179014/1‑251 (MQ=255)                                                                                                                                                                             gaagaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCgg                                                                                                                                                                                                                          <  2:108639/110‑1 (MQ=255)                                                                                                                                                                             gaagaaCAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCgg                                                                                                                                                                                                                          >  1:108639/1‑110 (MQ=255)                                                                                                                                                                                       cgcgcgTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCTTATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGAt                                                                   <  2:167263/251‑1 (MQ=255)                                                                                                                                                                                                      ccGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGATGGTGAAAGAACAttt                                                    <  1:176314/251‑1 (MQ=255)                                                                                                                                                                                                          aTTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGATGGTGAAAGAACATTTCCCg                                                <  2:276519/251‑1 (MQ=255)                                                                                                                                                                                                             cGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGATGGTGAAAGAACATTTCCCGCAt                                             >  1:269376/1‑251 (MQ=255)                                                                                                                                                                                                              ggTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGATGGTGAAAGAACATTTCCCGCAtt                                            <  1:209621/251‑1 (MQ=255)                                                                                                                                                                                                               gtggTTTTGGGCATATTACCGGACGTTTTCGGCTCTCCAGCGG‑GGGAGAAAATGGTGGTGCTCGATCACGTGCCGGACCATATCGAAACCTTGTGTAAATTTGGGATGATAGTGTTTTATGTCGATGCCACGCGGTTGGATTTTCTGGAATCTGCCGGAGCGGCTAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCTCATACCATCCTGCAACTGACAGAGATGGTGAAAGACCATTTCCCGCAttt                                           <  2:50693/248‑1 (MQ=255)                                                                                                                                                                                                                                  ctgacgataacagCTCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATACGTAACACATCGAAACATTGCCTAAATATGGGATGAAAGTGTTTTAAGGTGATGTCACGCGGATGGATTTAATGGAATCTGCCGGATAGGCGAAAGAGTAAGTGTTGATTAACGCCATCGACGATCCGCAAACCAACCTCCAACTGCCAGAGATGGTGAAAGTACATTTCCCGCATTTGCAGATTATTgctcccagc                        >  2:223431/13‑244 (MQ=255)                                                                                                                                                                                                                                         ttACTGCGCTCCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCTGGACCATATCGAAACCTTGCGTAAATTTGGTATGCAAGCGATTTATGGCGATGCCACGCGGATGGATTTACTGGAAACTGCCGGCGCGGCGAAAGAGGAAGTGCTGATTTACGCCATCGAAGATCCGCAAACCAACCGACAACTGACAGGGATGGGTAAAGAACATTTCCCGCACTTGCAAATTATTGCCCGCGCCCCCGATg                 >  2:217198/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ctcCAGCGG‑GGGTGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGATGGTGAAAGAACATTTCCCGCATTTGCAGATTATTGCCCGCGCCCGCGATGTCGAcc           <  2:106484/249‑1 (MQ=255)                                                                                                                                                                                                                                                       ggg‑gg‑TGAAAATGGTGGTACTCGATCACGATCCGGAACATATCGGAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGTCGCGGCGAAAGCGGAAGTGCTGATTAACGCAATCGACGATCCGCAAACCAACCTGCCACTGCCAGAGATGGTGAAAGAACATTTCCCGCATTTTCAGATTATTGCACGCGCCCGCGATGTAGTGCACTACATTc  >  2:217430/2‑251 (MQ=255)                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                       CTTTGTGGTATTTGGCGCGGCGCAGATGGCGAATGTGCTGGAGCCGGAGTGGGCGAAATCGCTGACCCTGGCGGTGGCGCTGTCGATGGCAGCAACGCCGATTCTGCTGGTGATCCTCAATCGCCTTGAGCAATCTTCTACTGAGGAAGCGCGTGAAGCCGATGAGATCGACGAAGAACAGCCGCGCGTGATTATCGCCGGATTCGGTCGTTTTGGGCAGATTACCGGACGTTTACTGCTCTCCAGCGG‑GG‑TGAAAATGGTGGTACTCGATCACGATCCGGACCATATCGAAACCTTGCGTAAATTTGGTATGAAAGTGTTTTATGGCGATGCCACGCGGATGGATTTACTGGAATCTGCCGGAGCGGCGAAAGCGGAAGTGCTGATTAACGCCATCGACGATCCGCAAACCAACCTGCAACTGACAGAGATGGTGAAAGAACATTTCCCGCATTTGCAGATTATTGCCCGCGCCCGCGATGTCGACCACTACATTC  >  NC_000913/48785‑49281

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