Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,956,756	C→T	L476L (CTG→TTG)	ilvA →	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,956,756	0	C	T	96.7%	89.4 / NA	30	L476L (CTG→TTG)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/17);  minor base G (0/1);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCT  >  NC_000913/3956511‑3956985                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       aaaTTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTggcg                                                                                                                                                                                                                                   <  1:221043/251‑1 (MQ=255) aaaTTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTggcg                                                                                                                                                                                                                                   <  2:120868/251‑1 (MQ=255)     gttGCTGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAACTACACGCGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACATCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGgcggcg                                                                                                                                                                                                                                <  2:164418/249‑1 (MQ=255)         cAGATGCTCAAGGACGGCGTCTACAGCGTGGTTGAGCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGAGGTCCATCGCATCCGTTGCAGGACCGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCATCCATGGCACCGACTTCGGTCGCGTATTGGGGGCGGTCGa                                                                                                                                                                                                                           <  2:135725/251‑1 (MQ=255)                   acgacgGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGAc                                                                                                                                                                                                                 <  1:36190/251‑1 (MQ=255)                                                  cgacgacgaAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATg                                                                                                                                                                                  <  1:56003/251‑1 (MQ=255)                                                          aaaTGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGc                                                                                                                                                                          <  2:120751/251‑1 (MQ=255)                                                              ggCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTAc                                                                                                                                                                       <  1:44300/250‑1 (MQ=255)                                                                       gcacGTGCGCTATATGGTCGCGGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCac                                                                                                                                                              <  1:15144/249‑1 (MQ=255)                                                                          cGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTAGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGTGACCATGCACCGGGTTTCGAACCCCCGTTCAATGAGCTGGGCAACGATTGCCacgacg                                                                                                                                                          >  2:169568/1‑251 (MQ=255)                                                                            ggcgcTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGaa                                                                                                                                                        <  2:41346/250‑1 (MQ=255)                                                                                         cggcggACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGGATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAAACCgg                                                                                                                                           >  2:241449/1‑251 (MQ=255)                                                                                               aCGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCa                                                                                                                                     <  2:75571/251‑1 (MQ=255)                                                                                                     aTCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCt                                                                                                                               >  2:280271/1‑251 (MQ=255)                                                                                                      tCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGCATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACCTTTCTTTGTTTCACTATCGCAGCCATGGCCCCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGACTGGGCGGGGCTCCGATTGCCATGACGAAACCAATAACCCGGCGTTCAGGttatt                                                                                                                              >  2:15328/1‑248 (MQ=255)                                                                                                                                      cGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAgcg                                                                                              <  1:47820/251‑1 (MQ=255)                                                                                                                                      cGAATTCCCGGAATCACAGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGAAGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGACTggtca                                                                                                >  2:72296/1‑245 (MQ=255)                                                                                                                                                         ggCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTAc                                                                           <  2:23624/251‑1 (MQ=255)                                                                                                                                                                              aaCACGCTGGGTACGTACTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTtgt                                                      >  1:46791/1‑251 (MQ=255)                                                                                                                                                                                           cgtacTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGGTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAAt                                        >  1:246851/1‑251 (MQ=255)                                                                                                                                                                                             tacTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGGTAGGGAAAAAAGCCTTGTAGAGCTTCGCTTATTAGGCCTACCCGCGCGCCAACGTCATTT‑TGTGTTCGGCaaaaact                                      >  2:44632/1‑248 (MQ=255)                                                                                                                                                                                              ccTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGACCCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAGAATCt                                      <  2:157903/249‑1 (MQ=255)                                                                                                                                                                                              acTGGAACATTTCTTTGTTTCACTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGGGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGGTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGAGGGTTAGGGAAAAAAGCCTGATAGAGCTATGCTTAACAGGACTACCCGCGCGACAACGGCATTTGTG‑GTTCGGCAAAATCtt                                     >  1:118968/1‑251 (MQ=255)                                                                                                                                                                                                                 tcACTATCGCAGCCCTGGCACCGACTACGGGCGCGTCTTGGCGGCGTTCGAACTTCGCGACCATGAACCGGATTTCGAAACCCGGCTGCATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCTGGTTCTTTTTGGCGGGTTAGGGAAAAAATCCTGATTGCGCTTCTCTTTTCACGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGTCAAATCTTCCCGAATGCCTCACTTAGc                  >  2:30722/2‑251 (MQ=255)                                                                                                                                                                                                                   aCTATCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAAATCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATATTCCAGAATGCCTCAATTAACgg                >  1:105456/1‑251 (MQ=255)                                                                                                                                                                                                                      aTCGCAGCCATGGCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGGCAACGTCATTTGTG‑GTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTc             <  1:6337/251‑1 (MQ=255)                                                                                                                                                                                                                       acGCAGCCAGGGCACCGACTACGGGCGCGTATGGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCa            <  1:44632/250‑1 (MQ=255)                                                                                                                                                                                                                        cGCAGCCATGGCACCGACTACGGGCGCGTAGTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGTCTCAt           <  1:241449/251‑1 (MQ=255)                                                                                                                                                                                                                               aTGNCACCGACTACGGGTGCGTATTGGCGGCGTTCGAACTTTGCGACCATGATCCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACTATTGCCACGACGAAACCAATAACCCGGCGTGCAGGTTCTTTTTGGCGGGTTATGGAAAAATTCTTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCAGTTGTG‑GTTCGGCATAATCTTCCATAATCCCTCAATTACCGGCTCATGTAGCCg    <  2:255111/251‑1 (MQ=255)                                                                                                                                                                                                                                 ggCACCGACTACGGGCGCGTATTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCCTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCt  >  1:134680/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       AAATTTTGCAGATGCTCAACGACGGCGGCTACAGCGTGGTTGATCTCTCCGACGACGAAATGGCGAAGCTACACGTGCGCTATATGGTCGGCGGACGTCCATCGCATCCGTTGCAGGAACGCCTCTACAGCTTCGAATTCCCGGAATCACCGGGCGCGCTGCTGCGCTTCCTCAACACGCTGGGTACGTACTGGAACATTTCTTTGTTCCACTATCGCAGCCATGGCACCGACTACGGGCGCGTACTGGCGGCGTTCGAACTTGGCGACCATGAACCGGATTTCGAAACCCGGCTGAATGAGCTGGGCTACGATTGCCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTG‑GTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCT  >  NC_000913/3956511‑3956985

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: ‑