Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,564,141	T→G	A183A (GCA→GCC)	yjiL ←	putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,564,141	0	T	G	90.3%	69.8 / ‑0.8	31	A183A (GCA→GCC)	yjiL	putative ATPase, activator of (R)‑hydroxyglutaryl‑CoA dehdratase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (15/13);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

AAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCA  >  NC_000913/4563907‑4564369                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     aaaataaaaaCTCTTTCATCGGCGTGTCCTCACTCGTTGAACAATTACAGCCGCGCCAATTGCGCCAGCAAATTGCGCCTCAGGATGGGTATTTACCGGCATTCGCAGGGGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCTCCAGTAAACAGAAACGGCGCTTCACAGGAGAGACGGGCAATGAAATTGGCACTCCTCCGCGCCATCTCGTTAATCCCTCCTGAGAGAATAGCttatg                                                                                                                                                                                                                      >  2:12893/1‑248 (MQ=255)            cctTTCACCGGCGTGTCCTCACTCTTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGccc                                                                                                                                                                                                           <  1:27627/250‑1 (MQ=255)                     tgCGTGTCCTCACTCGTTGAACAATTACCGCCGCGACAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTAGCAGGTGAGATTACAGCATCAGGGAAAACTTCTGGCAATGACTAACGCCACCAGTAACCAGAATCGGCGCTTCACCGGAGAGGCGAGACATGAAATTGGCAGTCCTCCGCGCCACAGCGTACATCACTCCGGCGAGAATCACTTCTGGCGCGACGCCCGCTGAgcgc                                                                                                                                                                                                  >  2:336949/2‑251 (MQ=255)                     ggCGTGGCCTCACCCGTTAACCAATTACCTCCGCTCCAATTGCGCCAGCAAATTGCGCATTAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTTCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGAACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAgcgc                                                                                                                                                                                                  <  1:134588/251‑1 (MQ=255)                      gCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAgcgc                                                                                                                                                                                                  >  1:226640/1‑250 (MQ=255)                              tCACTCGTTGACCAATTACCGCCTCGACAATTGCGCCAGCAAATTGCCCATCAGGCTGGGTCTTTACCGGCATTCGCAGGTGCGAATCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGACTCGGCGCTTCACATGAGAGCCGAGCAATGAAATTGGCACTCCTCCTCGCCATCGCGTTATTCACCCCGGCGAGAAACACATCTGACGCGACGCCCGCTGAGCGCCGGCTGATc                                                                                                                                                                                         >  2:302833/1‑251 (MQ=255)                                  cggTTGACAAATTCCCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCtt                                                                                                                                                                                     <  2:142545/249‑1 (MQ=255)                                              tACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCa                                                                                                                                                                          >  1:43868/1‑250 (MQ=255)                                                                       aaTTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGc                                                                                                                                                <  1:148812/251‑1 (MQ=255)                                                                           gcgcATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATGCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCTGCTAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTgc                                                                                                                                            <  1:50577/251‑1 (MQ=255)                                                                           gagcATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGAATCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAAGTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTgc                                                                                                                                            <  2:26625/249‑1 (MQ=255)                                                                                     aTGGGTAGTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGGGCGGCGTGACAt                                                                                                                                  <  1:61524/251‑1 (MQ=255)                                                                                               aCCGGCATTCGCAGGGGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAAt                                                                                                                        >  1:320729/1‑251 (MQ=255)                                                                                                             gTGAGATTCCAGCTTCCGGGCAACCTTCTGGCAATGACTAACGCCACCAGAAAACAGCAGCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGAGCCATCGCGTTAATCACTCAGGCGAGAATCGCTTATGGCACGACGCCCGCCGAGCGCAGGCTGACCTCTTCAGATTCAGCAAACAATGTGCACATAGTCGAGATGGCGTGCGGCGTGACATTCTAGGCAATGCGGTCTAGTTgac                                                                                                          >  2:203435/1‑249 (MQ=255)                                                                                                             gTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCt                                                                                                          <  1:177932/251‑1 (MQ=255)                                                                                                                       aGCATCCGGGCAAACTTATGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGt                                                                                                                                       >  1:293799/1‑212 (MQ=255)                                                                                                                              gggCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGc                                                                                         <  1:331275/251‑1 (MQ=255)                                                                                                                                         tGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCAGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGCGAt                                                                              >  1:248666/1‑251 (MQ=255)                                                                                                                                                    aaCGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAgg                                                                    >  1:318701/1‑250 (MQ=255)                                                                                                                                                    aaCGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAg                                                                     <  1:224503/249‑1 (MQ=255)                                                                                                                                                    aaCGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGa                                                                   <  1:6558/251‑1 (MQ=255)                                                                                                                                                      cGCCTCCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGaaa                                                                 <  1:12893/251‑1 (MQ=255)                                                                                                                                                         caccaGTAAACAGAATCGGCGCTTCACAGGAGAGACGCGCAATGAAATTGGCACTCCTACGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGCGAGCAGGCTGATCGCTTCTGATTAAGCCACCACTGTGCACATACTCGTGATGGCGTGCGGCGTGACCTTTTCGGTAATGATGTCGAGTTGCTCGACGCTGGTGAACAGCGTGCGAGAGATCCCCTCCAGGAAACGc                                                              >  2:102074/1‑251 (MQ=255)                                                                                                                                                                           tcgcTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCCTCGCGTTAATCACTCCGGCGAGACTCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCCCATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCTCGAGATCACCTCCCGGAAACGCCCGGATCCCGCCgcgctt                                            >  2:199918/2‑249 (MQ=255)                                                                                                                                                                                               gCAAAGAAATTGCCATTCCTCCGCCCTCTCTCGTTAATCCCTCCGGCGATATTCGTTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGTTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCTTTACTTTTTCGGTAATGCTGTCGAGTTGCTCGGCGCTGGTGCTAAGCGTGCTCGAGATCGCCTCCAGGAACCGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATc                        <  2:226640/251‑1 (MQ=255)                                                                                                                                                                                                cAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCg                       >  1:54367/1‑251 (MQ=255)                                                                                                                                                                                                        ttGGCACTCCTCAGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACCTTTTCGGTAATACAGTAGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCAGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCCCAGGtt               >  2:240198/1‑251 (MQ=255)                                                                                                                                                                                                           gCACTCCTCCGCGCCATCGCGTTAATCACCCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGGGCCCGCCGCGCATTGGTCATTCATCAGGAAATCGCACAGGTTAcc            <  1:102074/251‑1 (MQ=255)                                                                                                                                                                                                             aCTCCTCCGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGt          >  1:314561/1‑251 (MQ=255)                                                                                                                                                                                                                   ccGCGCCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGtcatc      <  1:82800/249‑1 (MQ=255)                                                                                                                                                                                                                       gcCATCGCGTTAATCACTCCGGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCCGGAAATCGCACAGGTTACCGtcttcatca  >  1:16691/1‑249 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     AAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAGATTCCAGCATCCGGGCAAACTTCTGGCAATGACTAACGCCACCAGTAAACAGAATCGGCGCTTCACAGGAGAGACGAGCAATGAAATTGGCACTCCTCCGCGCCATCGCGTTAATCACTCCTGCGAGAATCGCTTCTGGCGCGACGCCCGCTGAGCGCAGGCTGATCGCTTCTGATTCAGCAAACACTGTGCACATACTCGTGATGGCGTGCGGCGTGACATTTTCGGTAATGCTGTCGAGTTGCTCGACGCTGGTGCCAAGCGTGCGCGAGATCACCTCCAGGAAACGCCCGGTGCCCGCCGCGCATTTGTCATTCATCAGGAAATCGCACAGGTTACCGTCATCATCA  >  NC_000913/4563907‑4564369

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