Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	883,396	C→T	intergenic (‑8/‑277)	ybjG ← / → mdfA	undecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	883,396	0	C	T	92.3%	75.5 / ‑0.4	26	intergenic (‑8/‑277)	ybjG/mdfA	undecaprenyl pyrophosphate phosphatase/multidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (11/13);  total (12/14)
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.85e-01

GTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAA  >  NC_000913/883163‑883639                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    gTCCCATCGGCCAGAACACAAACAGTCTGACCGCCAGCGCGATAGCGATTGTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGATAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCATGACGCAAAGTGATTaa                                                                                                                                                                                                                                     <  2:266415/250‑1 (MQ=255)   cccATCGTCCACGACACAAACAGGCTGACCGCCAGCGCGATAGCGCTTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAACGTACCACGGCCAGCAAAGGCACCACGGAAATCAAATCTTTAGCAATACAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGGGATTCAAATTTTACAGACTAACTTCCTGACGCAAAGGGATTaaacg                                                                                                                                                                                                                                  >  2:194139/1‑249 (MQ=255)                                                      aTCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGt                                                                                                                                                                                <  2:323870/250‑1 (MQ=255)                                                        cACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTaa                                                                                                                                                                              <  1:45983/250‑1 (MQ=255)                                                             gTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGc                                                                                                                                                                        >  2:99912/1‑251 (MQ=255)                                                               tGCCGTTGGGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAGAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCAt                                                                                                                                                                      >  2:320646/1‑251 (MQ=255)                                                                                 ccccACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAg                                                                                                                                                    <  2:377173/251‑1 (MQ=255)                                                                                      caACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGcc                                                                                                                                               <  1:431811/251‑1 (MQ=255)                                                                                        aCCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTg                                                                                                                                             <  1:99912/251‑1 (MQ=255)                                                                                                        cAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGAtt                                                                                                                             <  2:410008/251‑1 (MQ=255)                                                                                                              cGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACg                                                                                                                       <  1:470074/251‑1 (MQ=255)                                                                                                                        gTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGtt                                                                                                             <  2:377996/251‑1 (MQ=255)                                                                                                                                             aaaaTCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCg                                                                                        >  2:249882/1‑251 (MQ=255)                                                                                                                                                               atcCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGtaataa                                                                      <  2:270244/251‑1 (MQ=255)                                                                                                                                                                 ccACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATg                                                                    >  2:316031/1‑251 (MQ=255)                                                                                                                                                                    cGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTa                                                                  >  1:466408/1‑250 (MQ=255)                                                                                                                                                                       aGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGaa                                                              <  1:203912/251‑1 (MQ=255)                                                                                                                                                                        gCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACGTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAAt                                                             <  2:466408/251‑1 (MQ=255)                                                                                                                                                                         ccGAGTCTGGCGTCGCGCTAATATGAGCGAATAGAGAGAGATTCAAATTTTCCAGCATCACTTCCTGACGCAAAGCGATTAAAAGGGCGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAAtg                                                            >  2:121234/1‑251 (MQ=255)                                                                                                                                                                                           taataaGAGAGAATAGAGAGAGATTCAACTTTTCCAGCATAACTTACTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACAGTGGttta                                           >  2:455771/1‑249 (MQ=255)                                                                                                                                                                                                    agaATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACaaa                                  >  2:83416/1‑250 (MQ=255)                                                                                                                                                                                                                    cAAATTTTCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTa                  <  2:403304/250‑1 (MQ=255)                                                                                                                                                                                                                           tCCAGCATAACTTCCTGACGCAAAGTGATTAAAAGGGCGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCt          >  1:28194/1‑251 (MQ=255)                                                                                                                                                                                                                                   aaCTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAATTCGTTGAGAATAAATTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTAGTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGGTAAATAAAACTTATAGAAAGCTTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCCCTTTTATTaa  <  2:55344/251‑1 (MQ=255)                                                                                                                                                                                                                                   aaCTTCCTGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGAAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATATAACTTAAGGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTACACGTGGTTTAAGCTACAAAATTCGCTTTCTCGTTAGCTGCACTTttatta   >  2:215762/1‑250 (MQ=255)                                                                                                                                                                                                                                   aaCTTCCTGACGCAAAGGGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTAAGCTGCGCTTttatta   >  1:208308/1‑250 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    GTCCCATCGTCCAGGACACAAACAGGCTGACCGCCAGCGCGATAGCGATTTTTATCACCAGTTGCCGTTGTGCTGTAAGCCCCCACAACCAAAGTACCACGGCCAGCAACGGCACCACGGTAATCAAATCTTTAGCAATAAAAATCGCCAACGAGATCATCCACGGAGCCGAGTCTGGCGTCGCGTTAATAAGAGAGAATAGAGAGAGATTCAAATTTTCCAGCATAACTTCCCGACGCAAAGTGATTAAAAGGGAGCCAATACAGGCAAGTCGTTGAGAATAAAGTGCAGGTTAAACTGGGTAAAGCGGCATCGTCTTATTTCCCTCAAGCGGCCTGTTTACGGTGGGTGATTGTAACGGGCATAGGTTAAATAAAACTTAAAGAAAGCGTAGCTATACTCGTAATAATGTAAGAATGTGCTTAACCGTGGTTTCAGCTACAAAATTCGCTTTCTCGTTAGCTGCGCTTTTATTAA  >  NC_000913/883163‑883639

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