Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,771,979	G→T	intergenic (‑235/‑302)	yibI ← / → mtlA	DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,771,979	0	G	T	96.9%	110.1 / ‑2.6	32	intergenic (‑235/‑302)	yibI/mtlA	DUF3302 domain‑containing protein YibI/mannitol‑specific PTS enzyme IICBA component
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (19/12);  total (20/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 = 9.30e-01

ACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATT  >  NC_000913/3771741‑3772222                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     aCATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAAcc                                                                                                                                                                                                                                          <  1:31972/251‑1 (MQ=255)   atatCTGCCGCGGATTCACGTTTATCAGTGTTGGTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCAt                                                                                                                                                                                                                                        <  2:354904/251‑1 (MQ=255)     atCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATACCACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACACATATTCAACAAACCATc                                                                                                                                                                                                                                       >  1:418825/1‑250 (MQ=255)      tCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATca                                                                                                                                                                                                                                      >  1:337375/1‑250 (MQ=255)                        ttagcAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACact                                                                                                                                                                                                                    <  1:319304/246‑1 (MQ=255)                              gttttgttTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTcaca                                                                                                                                                                                                             <  1:49455/248‑1 (MQ=255)                                                                     gATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCttgtt                                                                                                                                                                        <  2:354336/249‑1 (MQ=255)                                                                        atTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGGAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGGGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTGGTTG‑Cttt                                                                                                                                                                  <  2:223030/251‑1 (MQ=255)                                                                               cAATAATCCGGGTGCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAGCCAGCAAAACAAAAATGTGACTCTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTacac                                                                                                                                                           <  1:418652/251‑1 (MQ=255)                                                                                   aatCCGGGTTCGCGATTCTCGCCATAACACCAAAGCATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGACGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGc                                                                                                                                                       >  2:184960/1‑251 (MQ=255)                                                                                       cGGGTTCGCGCATCGCGCCATACCACAAAAGAATAATTTTTAGCGGTGATGAGTTGCTTCGTTACAAAACGATTGTATGACGAAGGCATAACATGCTGTAGCTCACATCAGGTGAACGACGTAAGAACATATCTTGTGATTCAGATCACAAATATTCACAAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGATG‑CTTTTTTACACAAGAGttt                                                                                                                                                   >  2:178477/1‑251 (MQ=255)                                                                                                 aTTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAac                                                                                                                                         <  2:139836/251‑1 (MQ=255)                                                                                                    ctcGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAacgtc                                                                                                                                      >  1:446497/1‑251 (MQ=255)                                                                                                                  aaaGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACctc                                                                                                                        >  2:345067/1‑251 (MQ=255)                                                                                                                       ataatTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTAc                                                                                                                   >  1:241614/1‑251 (MQ=255)                                                                                                                             tttAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCt                                                                                                              >  1:139703/1‑250 (MQ=255)                                                                                                                                agagGTGTTGCGTTGATTAGTTCGATACCGATTATATGACGAAGGCGTAACATGCTGAAGATCACCTCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCAAAAATATTCACCAAACAATCAACACAAAAATGTGACAATACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGACATGTTG‑CTTTTTTCCACAAGCGTTTTGTGATGAACGTCACGTCACTTACCTCTCTACCCCCTatat                                                                                                          >  2:401986/1‑251 (MQ=255)                                                                                                                                        tgaGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCAGTGAGACAATATCTTGTGATTCAGCTCACAAATATTCAACAAACAATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGAATTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTCTATTTAtgtg                                                                                                   >  1:250885/1‑250 (MQ=255)                                                                                                                                                    ttACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATcacac                                                                                      >  1:418529/1‑251 (MQ=255)                                                                                                                                                     tACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATcacac                                                                                      <  1:320903/250‑1 (MQ=255)                                                                                                                                                      aCATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACaa                                                                                    <  2:337375/251‑1 (MQ=255)                                                                                                                                                                           ggCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGCCAGtt                                                               >  2:279822/1‑251 (MQ=255)                                                                                                                                                                           ggCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGtt                                                               >  1:77846/1‑251 (MQ=255)                                                                                                                                                                              aTAACCTGCTGTAGATAACATCAGGTGAACGCCGTCAGAAAATATCTTGTGATTCAGCTCAAAAATATTCAACAAACCATCAAAACACAAATGTGACACTACTCACATTTAACTGCCATTTTTAGCGCAAATCGCAGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGACACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGCCTggcca                                                                  >  2:36267/1‑242 (MQ=255)                                                                                                                                                                              aTAACATGCTGTAGATCACATCATGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAAGTTGACACTACTCACATTTAAATGCCATGTTTAGCGAAAATCGCCGTCTT‑TTGTCTTTTTTACACAAGCGTTTTGTGATGAACGGCACGTCAATTACCTCGCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACa                                                                  <  1:36267/245‑1 (MQ=255)                                                                                                                                                                                    tGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTc                                                      >  2:28788/1‑251 (MQ=255)                                                                                                                                                                                            aTCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGcc                                                >  1:320431/1‑249 (MQ=255)                                                                                                                                                                                                        gAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTAt                                  >  2:281976/1‑251 (MQ=255)                                                                                                                                                                                                           cGCCGTAAGAAAATATCTTGTGATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCt                                >  1:267027/1‑250 (MQ=255)                                                                                                                                                                                                                             tgtgATTCAGATCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACctc              <  1:50315/250‑1 (MQ=255)                                                                                                                                                                                                                                 aTTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTgccg         >  1:30371/1‑251 (MQ=255)                                                                                                                                                                                                                                        tCACAAATATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACAtt  >  2:11166/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     ACATATCTGCCTCGGATTCACGTTTATCAGTGTTGTTTTTGGGCTGGCAGCCAGAAGGGAGTCAGGCTGATATTTTGACAATAATCCGGGTTCGCGATTCTCGCCATAACACCAAAGAATAATTTTTAGAGGTGATGAGTTGCTTAGTTACATAACGATTGTATGACGAAGGCATAACATGCTGTAGATCACATCAGGTGAACGCCGTAAGAAAATATCTTGTGATTCAGATCACAAAGATTCAACAAACCATCAAAACAAAAATGTGACACTACTCACATTTAAATGCCATTTTTAGCGAAAATCGCCGCCTTGTTG‑CTTTTTTACACAAGCGTTTTGTGATGAACGTCACGTCAATTACCTCTCTACCCCCTATATTTATGTGATTGATATCACACAAAAGGCCGTCGACTGGACAGTTAACCGATTCAGTGCCAGATTTCGCAGTATCTACAAGGTCCGGCTACCTCTGCCGCCACATT  >  NC_000913/3771741‑3772222

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