Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,500,114	G→A	20.4%	S112S (AGC→AGT)	ydcK ←	putative enzyme YdcK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,500,114	0	G	A	20.4%	100.1 / 18.6	54	S112S (AGC→AGT)	ydcK	putative enzyme YdcK
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base A (5/6);  total (24/30)
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.90e-01

CCCTTGCGCGATCATAAATTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGATACGCGTATTACCGGAAATCACCGCTCCGCCGAACGCAA  >  NC_000913/1499972‑1500250                                                                                                                                               |                                                                                                                                         cccTTGCGCGATCATAAATTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACg                                                                                                                                          <  1:665987/143‑1 (MQ=255)  ccTTGCGCGATCATAAATTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATAGGCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGc                                                                                                                                         >  1:1046194/1‑143 (MQ=255)       cgcgATCATAAATTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCa                                                                                                                                    >  1:765332/1‑143 (MQ=255)          gATCATAAATTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGAGACACTGTCACta                                                                                                                                 >  2:358648/1‑142 (MQ=255)               tAAATTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTtatata                                                                                                                            >  1:1067065/1‑143 (MQ=255)                aaaTTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTatatat                                                                                                                           <  1:1001675/143‑1 (MQ=255)                  cttgggAGCAGGAGTTGGTGTTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGc                                                                                                                         <  2:614879/139‑1 (MQ=255)                       gAGCAGGGGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCt                                                                                                                    <  1:242085/143‑1 (MQ=255)                         gCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTgg                                                                                                                  >  1:236973/1‑143 (MQ=255)                            ggAGTTCATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTg                                                                                                               <  1:216286/143‑1 (MQ=255)                                ttGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGAttt                                                                                                           >  2:635652/1‑143 (MQ=255)                                  aaggggcGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCa                                                                                                         <  1:502197/137‑1 (MQ=255)                                   aTGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGCCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCAc                                                                                                        <  1:608035/143‑1 (MQ=255)                                       tCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTCATATATGCACCCTGGCTGATGTCACtatt                                                                                                    >  2:395608/1‑143 (MQ=255)                                       gcGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATAAAATAAGTCTTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACtatt                                                                                                    <  1:754034/142‑1 (MQ=255)                                           gCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCAGCCGTGCTTGTTTCACTATTATCg                                                                                                >  2:718036/1‑143 (MQ=255)                                              ctAAGCCTTGTGCGGCCACGATCATAGAATGTTGACCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATc                                                                                             <  2:820886/142‑1 (MQ=255)                                                      tgtgCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGtt                                                                                     >  1:611526/1‑143 (MQ=255)                                                          cggcgACGATCATAGAATGTTGATCAATTAAGGCGAGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATcc                                                                                 <  1:423552/143‑1 (MQ=255)                                                           ggcgACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCg                                                                                >  1:1009448/1‑143 (MQ=255)                                                            tgtacgaTCATAGAATGTTGATCAATTGAGGCGTGACCAAAGATTCTACAGTGCCCGTAAACTAGTGATTCATGAATGGTGACGCTGTCACGTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGt                                                                               <  1:348276/140‑1 (MQ=255)                                                             cgacgaTCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGtt                                                                              >  2:766682/1‑143 (MQ=255)                                                             cgacgaTCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGtt                                                                              >  2:820217/1‑143 (MQ=255)                                                             cgacgaTCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGtt                                                                              <  1:453010/143‑1 (MQ=255)                                                              gacgaTCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTg                                                                             <  1:637124/143‑1 (MQ=255)                                                               acgaTCATAGTATGTTGATCAATTCAGGTGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGCCACTCCAATCTTCATCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGc                                                                            <  2:237545/143‑1 (MQ=255)                                                               acgaTCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGc                                                                            <  1:929743/143‑1 (MQ=255)                                                                  attatgtaaTGTTGATCAATTAAGGCGTCACGAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTa                                                                         <  2:794430/136‑1 (MQ=255)                                                                     aTAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAAc                                                                      <  2:257174/143‑1 (MQ=255)                                                                     aTAGAATGTTGAGCAAGTAAGGCGTGACCAAATATGCTACACTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATCTATGCACCCTGGCTGTTTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAAc                                                                      <  1:210363/143‑1 (MQ=255)                                                                         aaTGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACGTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACctct                                                                  <  2:742403/143‑1 (MQ=255)                                                                            gTTGATCGATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTccc                                                               <  2:278589/143‑1 (MQ=255)                                                                            gTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACATCTccc                                                               >  1:721735/1‑143 (MQ=255)                                                                              ggATCAATTAGGGCGTGACCAAATATTCTACACTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCca                                                             <  2:36737/142‑1 (MQ=255)                                                                               gATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCcac                                                            >  2:751449/1‑143 (MQ=255)                                                                                aTCGATTAAGGCCTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCcaca                                                           >  1:570176/1‑143 (MQ=255)                                                                                aTCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCcaca                                                           <  2:765332/143‑1 (MQ=255)                                                                                 tCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAACGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAg                                                          <  1:167406/143‑1 (MQ=255)                                                                                    aTTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGccc                                                       >  2:746488/1‑141 (MQ=255)                                                                                       aaGGCGTGACCAAATATTCGACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCt                                                    <  2:1029302/143‑1 (MQ=255)                                                                                          gCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCCCCGCCCGCTgtt                                                 >  2:366069/1‑141 (MQ=255)                                                                                          gCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATAGATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTCGCGTAAACCTCTCCCCACAGCACGCTAGt                                                 >  2:677326/1‑143 (MQ=255)                                                                                                 cAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCCCGCTGGTGCCGGTg                                          >  1:192425/1‑143 (MQ=255)                                                                                                  aaaTATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGa                                         <  2:236973/143‑1 (MQ=255)                                                                                                               tGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATGCCGGCGTTATACCTCGCGTAAGCCTCTCCCCACAGCACGCTGGTGCAGGTGATACGCGTATTAcg                            >  1:206939/1‑142 (MQ=255)                                                                                                               tGCCCGTAAACTAATGACACCAGCATGTTGACCCGGTCAATTAGAGGTGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGATACGCGTATTAcc                            <  1:983565/143‑1 (MQ=255)                                                                                                                gCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTCTCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGATACGCGTATTACCg                           <  2:166088/143‑1 (MQ=255)                                                                                                                 cccGTAAACTAATGAGTCATGAATGGTGACGCTGGCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGGGATAAGCGTATTACCgg                          >  2:799100/1‑143 (MQ=255)                                                                                                                       aaCTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCACTCCTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGATACGCGTATTACCGGAAATCa                    <  1:695984/143‑1 (MQ=255)                                                                                                                        aCTAATGTGTCATGAATGGTGACACTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGATACGCGTATTACCGGAAATCAc                   <  1:103207/143‑1 (MQ=255)                                                                                                                          tAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGTGTAAACCTCTCCCCACAGCACGATGGTGCCGGTGATACGCGTATTACCAGAAaccacca                 >  2:545005/1‑142 (MQ=255)                                                                                                                            aTGAGTCATGAATGGTGACGCTGCCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTCATCCGTTGCGTAAACCTCTCCCCAAACCACGCCGGTGCCACCGACACGCGAACTACCGGAAAACAccccc               >  2:1015894/1‑140 (MQ=255)                                                                                                                                    tGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCTCAGCACGCTGGTGCCGGTGATACGCGTATTACCGGAAATCACCGCTCCGCCGaa       <  1:635652/143‑1 (MQ=255)                                                                                                                                      taTGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGATACGCGTATTACCGGAAATCACCGCTCCGCCGAACg     >  2:576119/2‑143 (MQ=255)                                                                                                                                         ggTGACGCTGTGACTGATACATGCAGGCTGGCTGCTTGCGCTATTATCGATGCAGACGTTAGCCGTTGCGTATGCGTGTCCCCTAAGCACGCTGGTGCCGGTGATACGCGTATTACCGGAAATCACCGCTCCGCCGAACGCaa  <  1:304213/143‑1 (MQ=255)                                                                                                                                               |                                                                                                                                         CCCTTGCGCGATCATAAATTTGGAGCAGGAGTTGATGGTCGGGCGTTAAGCCTTGTGCGGCGACGATCATAGAATGTTGATCAATTAAGGCGTGACCAAATATTCTACATTGCCCGTAAACTAATGAGTCATGAATGGTGACGCTGTCACTTATATATGCACCCTGGCTGATTTCACTATTATCGATCCAGACGTTATCCGTTGCGTAAACCTCTCCCCACAGCACGCTGGTGCCGGTGATACGCGTATTACCGGAAATCACCGCTCCGCCGAACGCAA  >  NC_000913/1499972‑1500250

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