Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,636,796	C→A	Q205H (CAG→CAT)	yhiN ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,636,796	0	C	A	95.7%	165.3 / 0.8	46	Q205H (CAG→CAT)	yhiN	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (23/21);  total (24/22)
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.31e-01

TGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCG  >  NC_000913/3636651‑3636940                                                                                                                                                     |                                                                                                                                                 tGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTa                                                                                                                                              >  1:585784/1‑150 (MQ=255)   aCAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACCTGTAAct                                                                                                                                           >  2:705019/1‑151 (MQ=255)      aaTTCCCCCGGGTGCCAGGTGGTTGAAATCCGCAACCC‑CCCCCGGTCCAGACAAGCCGGGGTGGGGAAAGAGTAAGTTCTCACGAAAACCGGTGCCTTTTTCAGCGGTAAGCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGt                                                                                                                                               <  1:140229/145‑1 (MQ=255)       aTTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTc                                                                                                                                       >  1:113903/1‑151 (MQ=255)           ccccGGTTTCCAAGTGGTTAAAATCCTCAA‑CCCCCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCg                                                                                                                                      <  2:355809/149‑1 (MQ=255)            cccGGTTGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGc                                                                                                                                   >  1:684853/1‑150 (MQ=255)              cGGTTGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAAc                                                                                                                                >  2:761449/1‑151 (MQ=255)                  tGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGtt                                                                                                                            >  2:144931/1‑151 (MQ=255)                          gCTTGAAATTCGCAACCC‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCCCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGt                                                                                                                             <  2:38201/143‑1 (MQ=255)                                       aa‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCAcc                                                                                                        >  2:386901/1‑150 (MQ=255)                                         acccccgccgGTCCAGCCAACCGGCGGTGGGTAAAGAGTAGTTCCCCCCGAAAAAGGGTCCCGTTTTCAGGGGTAATCCGGGAG‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCa                                                                                                       <  1:629571/147‑1 (MQ=255)                                             ccgccgGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGAcc                                                                                                   >  2:468526/1‑151 (MQ=255)                                                           aggCCGCGGTGGGTGAAGATTAAGTTCCCACGAAAAACGGTCCCGTTTTCGGCGGAAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGgca                                                                                     <  2:765104/149‑1 (MQ=255)                                                           aggCCGCGGTGGGTAAAGAGAAAGTCCCCACGAAAAACGGTGCCTTTTTAGGCGGTAATCCGGGAG‑GGCCCCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGgca                                                                                     <  1:116073/149‑1 (MQ=255)                                                           aaGCCGGGGGGGGTAAAGATTAAGTTCTCCCGAAAACCGGTCCCTTTTCCGGGGTTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGt                                                                                          <  1:681298/146‑1 (MQ=255)                                                              ccGGGGGGGGTAAAGAGAAGGTTCCCACGAAAAAGGGTGCCTTTTTCAGCGGT‑ATCCCGGAAGGGCCCCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCg                                                                                        <  1:386852/145‑1 (MQ=255)                                                                gcgGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACg                                                                                >  1:11681/1‑151 (MQ=255)                                                                 cgGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGgca                                                                                     >  1:226489/1‑145 (MQ=255)                                                                  ggtgggtgAAGAGTAATTTCCCACGAAAAACGGTGCCTTTTTCGGGGGTAACCCCGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGgca                                                                                     <  2:226522/144‑1 (MQ=255)                                                                    tgggtgAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTg                                                                             >  1:234382/1‑150 (MQ=255)                                                                      ggtgAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAgg                                                                          >  2:297089/1‑151 (MQ=255)                                                                              gTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTg                                                                  >  2:17148/1‑151 (MQ=255)                                                                               tAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTg                                                                  >  1:460305/1‑150 (MQ=255)                                                                                 aGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTc                                                               >  2:366485/1‑151 (MQ=255)                                                                                   ttCTCCCGAAAAACGGGCCCGTTTTAGGGGGTAATCCGGGAG‑GGCCCCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAAtt                                                                   <  1:571786/145‑1 (MQ=255)                                                                                       cACGAAAAACGGTGCTTTTTTCAGCGGTAACCACGGAG‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCgg                                                             <  2:168252/147‑1 (MQ=255)                                                                                           aaaaccggcccgtttttCAGCGTTAACCCCGAAA‑GGCCCCCCCCCGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCCCCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCa                                                           <  1:43697/132‑1 (MQ=255)                                                                                             aaccGGTCCCTTTTTCAGGGTTAATCCCGAAG‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTtata                                                   <  2:653630/148‑1 (MQ=255)                                                                                                cGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAAcc                                                >  1:98835/1‑151 (MQ=255)                                                                                                 gggccctttttCAGGGGAAACCACGGAG‑GCCCCCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATaa                                                  <  1:762095/141‑1 (MQ=255)                                                                                                 gatGCCGTTTTCAGCGGTAATCATGGAA‑GGCCCCGCCACGCCCGCAAGCACATGTAACTCTTTGAGCAACGGCTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAAcc                                                >  1:394267/3‑150 (MQ=255)                                                                                                     ccGTTTTCAGGGGTAACCCCGGAG‑GGCCCCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCaa                                              <  2:354350/148‑1 (MQ=255)                                                                                                     acGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCGCATGTAACTCTTCGAGCAACGGTTTATGCCGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCCATCTTATACCCAAAc                                            >  1:522464/2‑150 (MQ=255)                                                                                                       gTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATCGTTCGGCAATCTTATAACCAAACGGc                                         >  2:30913/1‑151 (MQ=255)                                                                                                          ttCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGAc                                      >  1:287761/1‑151 (MQ=255)                                                                                                                       aCGGAA‑GGCACCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCg                         >  2:129515/1‑151 (MQ=255)                                                                                                                        cGGAG‑GGCCCCCCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCa                               <  1:144898/144‑1 (MQ=255)                                                                                                                              ggCCCCGCCACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCAtt                    <  1:241801/150‑1 (MQ=255)                                                                                                                                     ccACGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGccc            <  2:132210/151‑1 (MQ=255)                                                                                                                                       aCGCCCGCCAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCccac          >  2:468790/1‑151 (MQ=255)                                                                                                                                          cccggcAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGATGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCACAGCCCCGGCATTGACAGCccaccg        <  1:431211/150‑2 (MQ=255)                                                                                                                                          cccgccAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCccacca        <  2:616613/150‑1 (MQ=255)                                                                                                                                          cccgccAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCccacc         >  1:228115/1‑149 (MQ=255)                                                                                                                                             gccAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTcg    >  2:295353/1‑151 (MQ=255)                                                                                                                                             gccAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTcg    <  1:68598/151‑1 (MQ=255)                                                                                                                                               cAGCACATGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTcgcg  <  2:91996/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 TGACAAATTCCCCCGGTTGCCAGTAGCTTGAAATCTGCAA‑CA‑CCGCCGGTCCAGACAAGCCGCGGTGGGTGAAGAGTAAGTTCTCACGGAAAACGGTGCCGTTTTCAGCGGTAATCACGGAA‑GGCACCGCCACGCCCGCCAGCACCTGTAACTCTTCGAGCAACGGTTTATGCAGAGTGAATGGCACCAGACCCGCGCGGGTCGGCAGCACGTTGAGGCCAAATTGTTCGGCAATCTTATAACCAAACGGCGACGCGCCCAGCCCCGGCATTGACAGCCCACCAGTCGCG  >  NC_000913/3636651‑3636940

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