Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,698	C→T	F177F (TTC→TTT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,698	0	C	T	90.6%	85.0 / ‑1.6	34	F177F (TTC→TTT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (16/13);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

TGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTG  >  NC_000913/14473‑14942                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     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)                                                                                                                                                                                                              cagatgctgtaGGGATTCTTGGCCGTGAAGCAGACTTGTCCACACTTGCAGGGCCGCGGGACGCTGAGCAAAGATACGTGCAAAAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGTAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATTGCGCACCGGCAGGCGATCGGTCTGTGCAGGTTCAGGTTAAACAGCatac                <  2:46451/241‑4 (MQ=255)                                                                                                                                                                                                                  tGCGTCAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTAGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGCAGCGGGCGAGCATGGCGCACCGGCAGGCGGTCTGTAAGTTCCGGTTACGGTTAAACAGAACCCgctt            >  2:236094/1‑248 (MQ=255)                                                                                                                                                                                                                      tcAGGGATTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCg        >  1:226625/2‑251 (MQ=255)                                                                                                                                                                                                                            aTTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTg  <  1:49571/251‑1 (MQ=255)                                                                                                                                                                                                                            aTTCTTTGCCGTGCAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTg  >  1:166403/1‑251 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                     TGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGGTCAGGTGCAGATGCGCCAGGGATTCTTCGCTGTACAGCAGACCTGTCCACACTGTCAGGGCCGCGGTACGCTGATCAAAGATCCGTGCAACAAATGTCATGGTCATGGTCGTGTTGAGCGCAGCAAAACGCTGTCCGTTAAAATCCCGGCAGGGGTGGACACTGGAGACCGCATCCGTCTTGCGGGCGAAGGTGAAGCGGGCGAGCATGGCGCACCGGCAGGCGATCTGTACGTTCAGGTTCAGGTTAAACAGCACCCGATTTTCGAGCGTG  >  NC_000913/14473‑14942

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