Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,582,272	C→T	G169E (GGG→GAG)	ypfG ←	DUF1176 domain‑containing protein YpfG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,582,272	0	C	T	89.7%	60.6 / ‑1.1	31	G169E (GGG→GAG)	ypfG	DUF1176 domain‑containing protein YpfG
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (14/12);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

ATCAATGGTGTTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAG  >  NC_000913/2582028‑2582513                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  aataaTGGCGTTATAGGCCCCTGCCTCACATCTAATCAGCATCAGCGCTTTACTAACAGTCAGCGCCGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCATCGCAGACCATTCACCCGCCAGTTTCCATAATCCAGCAAAGCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCCCCGCGACCTCTGTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCt                                                                                                                                                                                                                                             <  2:24159/248‑1 (MQ=255)           ttATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCATTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCCcgc                                                                                                                                                                                                                                   >  1:147148/1‑251 (MQ=255)               aGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGAGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGAGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGGGCGGCGGTTCGCACTCTTTCTTGATCCACGCGGtt                                                                                                                                                                                                                               >  1:179217/1‑251 (MQ=255)                 gCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCgggttt                                                                                                                                                                                                                             >  1:181784/1‑247 (MQ=255)                       gccccacaacaaaTCAACATCAGCGCTTTATCATCAGGCAGCGCAGTGCCATTCACCTCACGCCGCAATGTATCAAGCGGGCAGCGCAGCCCAGTTATCCGCCAGTGTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTACCCACCGCGACCTCTTTCAGCGCAGGCTCGTGCTGTACGCTGAGCGGCTGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTg                                                                                                                                                                                                                         <  2:135449/238‑1 (MQ=255)                               gCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCTGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCt                                                                                                                                                                                                               <  1:107808/251‑1 (MQ=255)                                   atcatcatcAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTg                                                                                                                                                                                                           <  1:348364/251‑1 (MQ=255)                                            aGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTTATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATc                                                                                                                                                                                                  >  1:22924/1‑251 (MQ=255)                                            aGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGAACCACGCGGTTTCAATGCCAACGCGCTTTTGCTGAGCATc                                                                                                                                                                                                  >  1:253505/1‑251 (MQ=255)                                                         cagtcagCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAaacaacaac                                                                                                                                                                                     <  1:381753/251‑1 (MQ=255)                                                             cagCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCCCTGCCAAAGCGCTTTTGCTGAACATCAATAAACAAAAACgctg                                                                                                                                                                                 >  1:322081/1‑251 (MQ=255)                                                             cagCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACgctg                                                                                                                                                                                 <  1:35385/251‑1 (MQ=255)                                                              agCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACgctgc                                                                                                                                                                                <  1:288442/251‑1 (MQ=255)                                                                                                 aaGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAACCACTCCGCGAAATGGTCTGATCACCAt                                                                                                                                             >  1:156744/1‑251 (MQ=255)                                                                                                                 ccATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGGTTCACTGCCAACCCGCTTTTGCTGAGCATCAATAAACAACAAAGCTGCTTTTAAGCCACTCAGAGAAATGGTTTTATCGCCATCGCGTAAGGTGATTGc                                                                                                                             >  1:19249/1‑251 (MQ=255)                                                                                                                                 ttCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGGGGGTTAACTGCCACCGCGCTTTTGCTGAGCATCa                                                                                                                                                                                                 >  2:60412/1‑167 (MQ=255)                                                                                                                                 ttCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCa                                                                                                                                                                                                 <  1:60412/167‑1 (MQ=255)                                                                                                                                    cATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGCGAAATGGTCTGATCGCCATCCCGTAAGGTGATTGCCTTCCCTTCCTGAatcctc                                                                                                          >  1:99464/1‑248 (MQ=255)                                                                                                                                       aaTCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTg                                                                                                        <  1:122844/250‑1 (MQ=255)                                                                                                                                            aGCAAATCTTTGCGTTCTTCGAGTGAGAGTGTCGTCTGCGTTGTGTTTACCGTCGCTCTCTCTGTTATCTCTTGCGCGGGTGGTACTGTGAGTGGCGGTGCGTCCTCTATCTTGAGCCACGCGGTTTCACTGTCAACGCTCTTTTGCTGAGCCTCTATATACAACAACGCTGCTTTTAATCCACTTATTGAAATGGTCTGTTCGCTATCTCGTACGGTGATTGCCTTCCCTTCCTGAATCATCTGCaaaaa                                                                                                  <  2:96854/251‑1 (MQ=255)                                                                                                                                                       gcggtcttcggttGAGAGTGGCGTCGGCGTGGGGTTAACCCTCGCGGCCTCTTTCAGCGCATGCGCGTGCGGTACGCTGAGCGGCGGTTCGTCCGTTGTCTTGATCCACGCGGTTGCAATGACAACGCGCTTTTTCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTGCCCTGCCTGAATCATCTGCAAAAACGCGGGGATgg                                                                                       <  2:322106/239‑1 (MQ=255)                                                                                                                                                           tcttcGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCATGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGGAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTTATGGTTg                                                                                    <  2:189791/250‑1 (MQ=255)                                                                                                                                                                                         ttAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCctct                                                                                                         >  1:300305/1‑199 (MQ=255)                                                                                                                                                                                         ttAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCt                                                                                                         <  2:300305/199‑1 (MQ=255)                                                                                                                                                                                                    gACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACt                                          >  1:304414/1‑251 (MQ=255)                                                                                                                                                                                                             cAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGACCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAAAGCTTCTTTTAAACCACTCAGAGAAATGGTCTGATAGCCATCGCGTAAGGGGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGGGATTGTTGCCGGATCATCGGTTACTAATAAACATGGTGAAATCCGCCACTTGTccccac                                 >  1:270854/1‑251 (MQ=255)                                                                                                                                                                                                                              gTACGCTGAGCGGCGGTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCg                <  1:179120/251‑1 (MQ=255)                                                                                                                                                                                                                                         ggcggTTCGTCCTCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAAACATGGTGAAATCAGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAAc       >  1:126166/1‑249 (MQ=255)                                                                                                                                                                                                                                          ccctatcttCCCTTTTCTAGATCCCCCCGGTTCAACTGCCAACGCGCGTTTGCTGAGCATCAATAAACAACAACGTTGCTTGTAAGCCACTCAGAGAAAGGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCAGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGGGAAATCCGCCACTGGTCACCACTTAACGCCAACGGCTCGCCATCTAAca      <  2:187671/241‑1 (MQ=255)                                                                                                                                                                                                                                           cggatCGTCCTCTTTCTGGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGGCTGATCGCCATCGCGTAAGGGGATTGACTTACCTTCCTGAGTCATCTGCAAAAACGCGGTGATGATTGCCGTATCATCGGTTCCTAAAAACAAAGGCGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCACGCCCACTAAcagcc   >  2:98380/5‑250 (MQ=255)                                                                                                                                                                                                                                            gtttCGTCCGCTTTTTTTTTCCCCGCGGTTTCCCTGCCAAAGCCCTTTTGCTGAGCATCAATATACAACAACGCTGCGTTTAAGCCACTCAGAGAAATGGTCTGATCGCCACCGCGTAAGGTGATTGCCTGCCCTTCCTGAATCATCTGCAAAAACGCGGTGCTGGTTGCCGTATCATCGTTTACTAATAACCATGTTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTTGCCATCTAAcatcaa  <  2:15608/249‑5 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  ATCAATGGTGTTATAGGCCCCTGCCTCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGTTGCGTTCTTCGAGTGAGAGTGGCGTCGGCGTTGGGTTAACCACCGCGACCTCTTTCAGCGCAGGCGCGGGCGGTACGCTGAGCGGCGGTTCGTCCCCTTTCTTGATCCACGCGGTTTCACTGCCAACGCGCTTTTGCTGAGCATCAATAAACAACAACGCTGCTTTTAAGCCACTCAGAGAAATGGTCTGATCGCCATCGCGTAAGGTGATTGCCTTCCCTTCCTGAATCATCTGCAAAAACGCGGTGATGGTTGCCGTATCATCGGTTACTAATAACCATGGTGAAATCCGCCACTTGTCACCACTTAACGCCAACGGCTCGCCATCTAACAGCAG  >  NC_000913/2582028‑2582513

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