Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,647	T→C	T160T (ACT→ACC)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,647	0	T	C	87.1%	67.2 / 0.8	32	T160T (ACT→ACC)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base C (13/14);  total (15/17)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

CGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGG  >  NC_000913/14410‑14871                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 cggcggTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATGTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGGGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGAGCTGTc                                                                                                                                                                                                                     <  2:175680/251‑1 (MQ=255)           ggcggcggcGCAGACTGCAGCGATATTTTGGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTg                                                                                                                                                                                                            <  2:341772/250‑1 (MQ=255)                 ggcgCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTCCGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGCCTCTGGCAGAGTGTGCCGTTTGCCACGGTAGCGGTGCCAAACCAGGTCCACAGCCGCAAACCTGTCCGCCCTGTCATGGTTCTgcgacc                                                                                                                                                                                                       >  2:196928/1‑244 (MQ=255)                 ggcgCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTg                                                                                                                                                                                                            <  1:196928/244‑1 (MQ=255)                   cgcAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCCTGGTTCTGGTCAGGTGc                                                                                                                                                                                                   >  1:71578/1‑251 (MQ=255)                     cAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGc                                                                                                                                                                                                   >  1:68999/1‑249 (MQ=255)                        aCTTCAGAGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGAGCGGCGCGCGGTGCTTATTTACGCTCTAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTTACCAAAGAGATCCGAATTCCGAACCGGGAAGAGGGTGACCAATGACACGGTAGCGGTGCAAAACCAGGTACAAATCCGCAGGCATGTACGACCCGTCATGGTTCTggaatgtgtcagatg                                                                                                                                                                                              >  2:141862/1‑238 (MQ=255)                             aGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATgcgtc                                                                                                                                                                                          >  1:173501/1‑248 (MQ=255)                                          gTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACCTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGTCAg                                                                                                                                                                                        <  1:199722/239‑1 (MQ=255)                                                  ttCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGt                                                                                                                                                                                                     <  1:97488/218‑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)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                 CGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGG  >  NC_000913/14410‑14871

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: ‑