Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,095,057	G→A	V175M (GTG→ATG)	hisH →	imidazole glycerol phosphate synthase subunit HisH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,095,057	0	G	A	93.3%	97.2 / ‑5.8	31	V175M (GTG→ATG)	hisH	imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (11/17);  total (13/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

CGACTTGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCGTGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCT  >  NC_000913/2094813‑2095298                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  cGACTTGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATCTGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAAACGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGt                                                                                                                                                                                                                                             >  1:56887/1‑251 (MQ=255) cGACTTGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGt                                                                                                                                                                                                                                             >  2:165194/1‑251 (MQ=255)       gCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCAcc                                                                                                                                                                                                                                       >  1:315323/1‑251 (MQ=255)                                                 tGGTCTGCCCCTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAAc                                                                                                                                                                                              <  1:392982/250‑1 (MQ=255)                                                  gacgtgccaacgccaaaaaggggCGGGAACCGGGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTGCTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACtt                                                                                                                                                                                            <  2:341025/231‑1 (MQ=255)                                                                          tGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCg                                                                                                                                                                    >  1:124256/1‑251 (MQ=255)                                                                                   gTTTACCAGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACTGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCACAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCGTGCAGTTCCACCCGGAGCGCTCTAGTGCCGCTGGTGGTCAGTTGCTGAAACACTTCCTGGAGATGTGCTTCTTATTCCGGCATTCGAt                                                                                                                                                           >  2:387712/1‑251 (MQ=255)                                                                                        cccGCAGGCGGGGAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTTAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAAt                                                                                                                                                      <  1:59897/251‑1 (MQ=255)                                                                                             agtcaggcaACCGCCTGTTTCAGGGGATTGAAGACGGAGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCTCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTGAATCGACg                                                                                                                                                 <  1:264546/248‑1 (MQ=255)                                                                                                           cTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGt                                                                                                                                   >  1:288414/1‑251 (MQ=255)                                                                                                                       aTTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTGCGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGACGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGTTTCCGCTGGTGCTAAGGTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGCCGGCACTGTGGTGCGTCTTCATCAGGGc                                                                                                                       <  2:271699/251‑1 (MQ=255)                                                                                                                       aTTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGCACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGCTTTAATAGACGGCACTGTGGTGCGTCTCCATCAGGGc                                                                                                                       >  1:236938/1‑251 (MQ=255)                                                                                                                          gAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGAt                                                                                                                    <  2:182133/251‑1 (MQ=255)                                                                                                                           aaGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGAtt                                                                                                                   <  1:35156/251‑1 (MQ=255)                                                                                                                              aCGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCAGCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACg                                                                                                                <  2:337563/251‑1 (MQ=255)                                                                                                                                     gTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTGCCACCCGGAGCGTTCTGGTGCCGCGGGTGGTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACTGCACTGTGGTGCTTCTCCATCAGGGCGATTACGGCAAACa                                                                                                         <  1:141756/251‑1 (MQ=255)                                                                                                                                                               cAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGc                                                                               >  1:157183/1‑251 (MQ=255)                                                                                                                                                               cAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGCTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGc                                                                               >  1:46787/1‑251 (MQ=255)                                                                                                                                                                 aTGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCcg                                                                             <  1:154152/251‑1 (MQ=255)                                                                                                                                                                      ggtcatcccgctgACCATAGCCAAGTGCAATTACGGAGAACAGATCACCGCGGCGGTAAAAAACGACAACTTCTATGGCCTGCAGTTACACCCGGCGCGTTCTGGTACCGCTGGTGCTAAGCTGCTGAAAAACTTCCTGGGGCTGAGATGATTATTCCGGCCTTAGAGTTAATCGACGGCACTGTGGTGCGTCTCCATGAGGGCGATTAAGGCCAACAGCGCGATTACGGTAACGACACGCTTGCGCGAtg                                                                        >  2:243767/13‑250 (MQ=255)                                                                                                                                                                              cGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTcc                                                                                                                                                                                          <  2:16683/129‑1 (MQ=255)                                                                                                                                                                              cGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTcc                                                                                                                                                                                          >  1:16683/1‑129 (MQ=255)                                                                                                                                                                                                    ttACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGt                                          >  1:204960/1‑251 (MQ=255)                                                                                                                                                                                                     tACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGt                                          <  2:83213/250‑1 (MQ=255)                                                                                                                                                                                                                    aCCGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGAt                          <  2:88067/251‑1 (MQ=255)                                                                                                                                                                                                                     ccGCGGCGGTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATc                         <  1:93306/251‑1 (MQ=255)                                                                                                                                                                                                                          gcggTACAAAAAGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGAcc                    <  2:315323/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaGATAACTTCTACGGCATGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAg          <  2:179814/251‑1 (MQ=255)                                                                                                                                                                                                                                       aTAACTTCTACGGCGTGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATc       <  1:387712/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaCTTCTACGGCATGCAGTTCCACCCGGAGCTTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGAT‑ACGGTAACGTCCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGATGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCgg    <  1:312184/251‑1 (MQ=255)                                                                                                                                                                                                                                            gtCTACGGCATGTAGTTCCACCCGGGGTGTTCTGGTGCCGCTGGCGCTAAGTTGCTGAATAACTTCCTGGAGATGTGATGATTATTCCGGAATTAGATTTAATAGAAGGCACTGTGGTGCGTCTCCATCAGGGCGATTATGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGGGCCGAAGGGATGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCt  >  2:43186/2‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  CGACTTGCTGGGCATCATCGACGAAGACGTGCCGAAAATGACCGACTTTGGTCTGCCACTGCCACATATGGGCTGGAACCGCGTTTACCCGCAGGCAGGCAACCGCCTGTTTCAGGGGATTGAAGACGGCGCGTACTTTTACTTTGTTCACAGCTACGCAATGCCGGTCAATCCGTGGACCATCGCCCAGTGTAATTACGGCGAACCGTTCACCGCGGCGGTACAAAAAGATAACTTCTACGGCGTGCAGTTCCACCCGGAGCGTTCTGGTGCCGCTGGTGCTAAGTTGCTGAAAAACTTCCTGGAGATGTGATGATTATTCCGGCATTAGATTTAATCGACGGCACTGTGGTGCGTCTCCATCAGGGCGATTACGGCAAACAGCGCGATTACGGTAACGACCCGCTGCCGCGATTGCAGGATTACGCCGCGCAGGGTGCCGAAGTGCTGCACCTGGTGGATCTGACCGGGGCAAAAGATCCGGCT  >  NC_000913/2094813‑2095298

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