Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	11,366	C→T	intergenic (‑10/+16)	yaaW ← / ← yaaI	putative enzyme‑specific chaperone YaaW/DUF2541 domain‑containing protein YaaI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	11,366	0	C	T	96.9%	87.2 / ‑4.2	32	intergenic (‑10/+16)	yaaW/yaaI	putative enzyme‑specific chaperone YaaW/DUF2541 domain‑containing protein YaaI
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (14/17);  total (14/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 = 7.84e-01

ATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCAT  >  NC_000913/11124‑11605                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 aTACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTatat                                                                                                                                                                                                                                          <  1:198579/251‑1 (MQ=255) aTACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAAACCTTCAGGTAATTAACATTCATTTTTTGTGGTTTCtatct                                                                                                                                                                                                                                          >  1:371293/1‑249 (MQ=255)              caAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGGGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCCGGTAATTAAAATTCATTTTTTGTGGTTTCTATTTTCTGGCGTTAgtc                                                                                                                                                                                                                             >  2:479964/1‑251 (MQ=255)               aaaaTGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAgtcg                                                                                                                                                                                                                            >  1:165626/1‑251 (MQ=255)                               ccagCAATCAGCTGCCAGTTGCGGCGATGGTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCa                                                                                                                                                                                                            <  1:15305/248‑1 (MQ=255)                                                             tGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTAt                                                                                                                                                                             >  1:275931/1‑251 (MQ=255)                                                              gCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATg                                                                                                                                                                            <  2:291873/251‑1 (MQ=255)                                                                   ggATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCg                                                                                                                                                                       <  2:167415/251‑1 (MQ=255)                                                                         ccTTCCATCGATTTAAACATTTCGTTGCGTATCAGTCCGCTGGAGATGCGAGTGTTGCCGTTTTCATTATGGGTGAGCAATCGGGCGAAATTGGCCAACTGTTCCTCTCTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGTCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTGCACCGTATTACCTGAAAAgg                                                                                                                                                                 <  2:382089/251‑1 (MQ=255)                                                                                               cGTTGCGTATCAGTACGCTGTAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCATATCTTAATCATTCAGGTAATTAACATGCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCGAAGAGGTGATTTTTGAGACGCAGCGtt                                                                                                                                            <  1:429831/250‑1 (MQ=255)                                                                                                gTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTa                                                                                                                                          <  1:41473/251‑1 (MQ=255)                                                                                                  tGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGAcc                                                                                                                                                                         >  1:176676/1‑218 (MQ=255)                                                                                                  tGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGAcc                                                                                                                                                                         <  2:176676/218‑1 (MQ=255)                                                                                                   gcgtaTCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAAGCGGGCGACATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAAGTTTCAGCGTG‑GCCATATCCGATGATTTCACCCGATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTAtt                                                                                                                                        <  1:427005/246‑1 (MQ=255)                                                                                                     gtaTCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTc                                                                                                                                     <  1:455425/249‑1 (MQ=255)                                                                                                         cAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTa                                                                                                                                 <  1:55003/251‑1 (MQ=255)                                                                                                               gCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTg                                                                                                                           <  1:353144/251‑1 (MQ=255)                                                                                                                                 gCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCcagtca                                                                                                         <  1:479964/251‑1 (MQ=255)                                                                                                                                       ttCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAg                                                                                                                                                                                                                               >  2:343754/1‑128 (MQ=255)                                                                                                                                       ttCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAg                                                                                                                                                                                                                               <  1:343754/128‑1 (MQ=255)                                                                                                                                              tGGGTGAGCAATCGGGCGAAATTTGCCAACTGATCCTCACTACAATGCTGCAGAAAATCCAGATCTGAGTCCTTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCctt                                                                                            >  2:186508/1‑251 (MQ=255)                                                                                                                                                 gTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGAcc                                                                                                                                                                         >  1:135059/1‑171 (MQ=255)                                                                                                                                                 gTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGAcc                                                                                                                                                                         <  2:135059/171‑1 (MQ=255)                                                                                                                                                                          aCTGTTCCTCACTACAATGCTGAAGAAACTCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGAGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAgact                                                                >  1:349491/1‑248 (MQ=255)                                                                                                                                                                                     cTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTCTATTATGGCGTTAGTCGTCGCCGTTAATTTTCCGCGTG‑GCCATATCCGATGAGTTCACCGTCTGACCGGAAAAGGTGATTTTTGAGACGAAGCGTTTATTGTCGTTATCCCTGTTACTGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTCTTTCTGAAGGATTATTCAGACTCTGACTGgcgc                                                     >  2:482718/1‑251 (MQ=255)                                                                                                                                                                                                aagaaAAGCCAGATCTGAATCATGCAGGTAATTAACATTCATTTTTTGTGGTTGCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGGG‑GCCATATCCGATGAGTTCACCGTATGAACGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCt                                          <  2:161444/251‑1 (MQ=255)                                                                                                                                                                                                 agaaAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCtt                                         >  1:174731/1‑251 (MQ=255)                                                                                                                                                                                                     aaTCCAGATCTGAATCATTCAGGTACTTAAATTTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCTTCGCCGATAAGTTTCAGCGTG‑GCCATAGCCGATGAGTTCACCGTATGACCGTAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTTATCCAGGCAGTGGTTTGCCCTTCTTTTAGTTCGGAAGGAATATGCAGACTCTGACTGGCGCTACGGGCGGCTTTGaa                                     <  1:398725/251‑1 (MQ=255)                                                                                                                                                                                                     aaTCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGTCGGCTTTGa                                      >  1:169598/1‑250 (MQ=255)                                                                                                                                                                                                                 aaTCATTCAGGGAATTAACATTCATTTTTTGTGGTTTCTATATTCTTGCGTTAGTCGTCGCCGATAATTTTCAGC‑TGTGCCATATCCGATTAGTTCACCGTATGACCGGAAAATGTGATTTTTGAGACGCAGCGTTTATTTTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGTGCGGCTTTGAAATAAACCGATGc                         <  1:454853/251‑1 (MQ=255)                                                                                                                                                                                                                                 aaCATTCATTTTTTGTGGTTTATATATTCTGGCGTTCGTCGTAGCCGATACTTTTCAGCGTG‑GCCATATCCGATGAGTTAACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCATTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAGCCGATGCACTGCTTAACTgcaa          >  1:389602/1‑247 (MQ=255)                                                                                                                                                                                                                                        attTTTTGTGGTTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCGGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCACTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGACTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCAt  >  1:402309/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 ATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGTTAGTCGTCGCCGATAATTTTCAGCGTG‑GCCATATCCGATGAGTTCACCGTATGACCCGAAAAGGTGATTTTTGAGACGCAGCGTTTATTGTCGTTATCGCTGTTAATGTTGATCCAGTCAGTGGTTTGCCCTTCTTTTATTTCTGAAGGAATATTCAGGCTCTGACTGGCGCTACGGGCGGCTTTGAAATAAACCGATGCACCGCTTAACTGTAAATCGCCAT  >  NC_000913/11124‑11605

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