Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,686	C→T	R173R (CGC→CGT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,686	0	C	T	90.3%	70.6 / ‑1.3	31	R173R (CGC→CGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (14/14);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

GTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATT  >  NC_000913/14451‑14932                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       gTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAg                                                                                                                                                                                                                                                     <  1:199722/239‑1 (MQ=255)                       tGGCGGCGGACGTGGGCGTCAACGTGTGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCGGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCacac                                                                                                                                                                                                                    <  1:54957/250‑1 (MQ=255)                             gggtcttgggagtcAACGTGCGGCGTGCGGTGGTGATGTATGCGATAACATGGACCTTAACCTCGAAGAAGCGGTACGTGGCGTGCCCAAAGAGATCCGCATTCCGACGCTGGAAGAGTGTGACGTTTGCCACGGTAGCGTTGCAAAACCAGGTACACGGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCAAACTGTCAgg                                                                                                                                                                                                             <  2:81250/240‑1 (MQ=255)                             cggACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCAGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAgg                                                                                                                                                                                                             >  1:27705/1‑251 (MQ=255)                                            aCGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTg                                                                                                                                                                                               >  1:252658/1‑250 (MQ=255)                                                     gcgcgGTGCTGATTTACGCTATAACATGGAGCTCTCCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGTCTCTGGAAGAGTGTGACGGTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGagc                                                                                                                                                                                     <  2:347539/251‑3 (MQ=255)                                                        cgGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGAACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGt                                                                                                                                                                                  >  1:176003/1‑251 (MQ=255)                                                                    aCGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtc                                                                                                                                                                      <  1:152248/251‑1 (MQ=255)                                                                       cTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtcatg                                                                                                                                                                   >  1:50088/1‑251 (MQ=255)                                                                            aCATGGAGTTCATCCGCGCAGAAGCTGTACGTGGCGTGACCAAAGAGATCTGCATTCCGACGCTGGTAGAGTTTGACGTTTGCCACGGTAGCGTTTCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGAGTCTCCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCATCAAATgtcatggtcat                                                                                                                                                              <  2:29117/251‑1 (MQ=255)                                                                             cATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAATGAGATCCGCATTCCGACTCTGGTAGAGCGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGCTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATgtcatggtcatg                                                                                                                                                             >  1:261615/1‑251 (MQ=255)                                                                                       aCCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTg                                                                                                                                                    >  1:82556/1‑250 (MQ=255)                                                                                             gaagaaGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAgcgc                                                                                                                                               >  1:200098/1‑249 (MQ=255)                                                                                                 aagCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCaaa                                                                                                                                         <  1:324289/251‑1 (MQ=255)                                                                                                                      aatagctccccattccctctctGGAAGTGTGTGCCGTTTGCCCCGTGCGCGGTGCAAAACCAGGTACCCCGCCGCCTACCTGTCCGACCTGTCATTGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCTCATCAAAACGCTGTCCGTTAAAATccct                                                                                                                    <  2:68999/233‑2 (MQ=255)                                                                                                                             ccGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAgggg                                                                                                             <  1:303326/251‑1 (MQ=255)                                                                                                                             ccGCATACCGACGCTGGAAGAGTGGTACGTTTGCCACGGGAGCGGTGCAAAACCTGGTCCAGAGCCGCAGACCGTTCCGACCTGTCATGGTTCTGGTCAGGGTCAGATGCGTCAGGGATTCTGTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGTTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAggg                                                                                                              <  2:27098/250‑1 (MQ=255)                                                                                                                                     cgattctGGAAGAGGTTGACTTTGGCCGTGGTTGCGGTGCAAATACAGGGACACAGCCGCATACTTGTCCGTCCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGAGCCGTGCAACAAATGTCATGGTCATGGGCGTGTTGAGCGCAGCAAAACGCTGTCCGGTAAAATCCCGGCATGGGTTGACACt                                                                                                     <  2:50088/247‑1 (MQ=255)                                                                                                                                                      gACGTTTGCCACGGTAGCGGTGCAAAGCCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACGGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCt                                                                                    <  2:108414/251‑1 (MQ=255)                                                                                                                                                                  ggTAGCGGTGCAAAACAAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGACAGGTGCAGATGCATCAGGGATTATTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAAAAAAACGATGTCCGTTAAAATCCCGGCAGGGGTGGAACCTGGAGCCCGCATCCGGCTTGCAGGCGaggg                                                                        >  2:315207/1‑248 (MQ=255)                                                                                                                                                                   gTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCCTGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGacggg                                                                       >  1:336343/1‑247 (MQ=255)                                                                                                                                                                       cGGTGCAAAACCAGGTACACAGCCGCCGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTATAATACCGGCAGGGGTGGCCACTGGAGACCGCATCCGTCTTGTGGGCGAAGGTGagg                                                                   >  2:86521/1‑249 (MQ=255)                                                                                                                                                                           gCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCggg                                                               <  1:112218/251‑1 (MQ=255)                                                                                                                                                                                                          cGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTAc                                >  1:258449/1‑251 (MQ=255)                                                                                                                                                                                                                gtcatggttctgtgtagttgtcgATGCGTCAGGGATTCTTGGCCGGGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCGGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAACCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttca                           <  1:23961/228‑1 (MQ=255)                                                                                                                                                                                                                     ggTTCTGGTCAGGTGCAGCTGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGAGCACCGGCAGGCGATCTGTACgttcaggttc                      >  1:95927/1‑250 (MQ=255)                                                                                                                                                                                                                     ggTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACgttcaggttc                      <  1:141314/250‑1 (MQ=255)                                                                                                                                                                                                                      gTTCTGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGCCCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCTCCGGCAGGCGATCTGTACgttcaggttcag                    <  1:68541/251‑1 (MQ=255)                                                                                                                                                                                                                          tGGTCAGGTGCAGATGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAAAGCTGTCCGTTAAACTCCCGGCAGGGGTGGACCCTGGCGACCGCATCCGTCTTGCGGGCGAAGGCGAAGCGGGCGAGCGTGGCGCACCGGCAGGCGATCTGTACGTTcaggttcaggtt                 >  2:222517/1‑250 (MQ=255)                                                                                                                                                                                                                                  tGCAGATGCGTCAGGGATTCTTTGCCGAGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCCGGTTCAGGTTAAACAGCAc        >  2:162003/1‑251 (MQ=255)                                                                                                                                                                                                                                        tGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTAGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGCAGCGGGCGAGCATGGCGCACCGGCAGGCGGTCTGTAAGTTCCGGTTACGGTTAAACAGAACCCgctt  >  2:236094/1‑248 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                       GTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATT  >  NC_000913/14451‑14932

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