Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,401	C→T	G78G (GGC→GGT)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,401	0	C	T	93.3%	89.7 / ‑6.6	30	G78G (GGC→GGT)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (20/8);  minor base G (1/1);  total (21/9)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACT  >  NC_000913/14180‑14647                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       gATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGgcgc                                                                                                                                                                                                                           >  1:114814/1‑251 (MQ=255)    tattaCGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGgcgc                                                                                                                                                                                                                           >  1:62718/1‑248 (MQ=255)                              aaCAGCGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGTCGGAGCAGACTTCAGCGATATTTTTGGTGACGttt                                                                                                                                                                                              >  2:102667/1‑251 (MQ=255)                                   cGGAAGAGCATGAAATCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGCTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGc                                                                                                                                                                                         >  1:127553/1‑251 (MQ=255)                                        gagcatgcaaTCAAAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGatat                                                                                                                                                                                    >  1:72871/9‑251 (MQ=255)                                                  tcaaaaaGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTtgcggc                                                                                                                                                                          >  1:4871/5‑251 (MQ=255)                                                  acaaaaaGGCCTACAACCGCCTGGCCATGAAATACCACCCGCACCGTAGCCAGGGTGACAAAGAGGCCGAGGAGAAATTGAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGATCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCGGCGATATTTGTGGTGACGTTTTCGGCGATATTTTTggcggc                                                                                                                                                                          <  2:24527/247‑1 (MQ=255)                                                   caaaaaGTCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGATCAAGGTGGTATGGGCGGCTGCTGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTggcggc                                                                                                                                                                          <  1:164873/247‑1 (MQ=255)                                                         ggCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTg                                                                                                                                                                   >  1:54913/1‑251 (MQ=255)                                                                  aCGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAAc                                                                                                                                                          <  1:250546/251‑1 (MQ=255)                                                                   cGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGCCGTGGTCGTCAACg                                                                                                                                                         >  1:60997/1‑251 (MQ=255)                                                                      cTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTgc                                                                                                                                                      >  1:150469/1‑251 (MQ=255)                                                                        ggCCATGAAATACCACCCGGGCCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTgcgg                                                                                                                                                    >  1:58043/1‑251 (MQ=255)                                                                                aaTACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCTAGGCGAAATGTAAAGATATCAAGGAAGCTTATGAGGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGt                                                                                                                                            <  1:131573/251‑1 (MQ=255)                                                                                   accaccCGGACCGTAACCCGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCCGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGGCAACGTGCGGCGCGCGGTGCt                                                                                                                                         >  2:239893/1‑251 (MQ=255)                                                                                          ggACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGCGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTAc                                                                                                                                  >  1:235832/1‑251 (MQ=255)                                                                                            aCCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTGCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGc                                                                                                                                <  2:114814/251‑1 (MQ=255)                                                                                                                     cGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCga                                                                                                         >  1:224028/1‑249 (MQ=255)                                                                                                                               tttAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGAAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTCTAACATGGAGCTCACCCTCGAAGCAGCTGTACg                                                                                             >  1:123256/1‑251 (MQ=255)                                                                                                                                   aaGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGc                                                                                         >  2:250099/1‑251 (MQ=255)                                                                                                                                            aGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCACGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTa                                                                                                                                   >  1:256648/1‑200 (MQ=255)                                                                                                                                            aGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGGATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTGGGTGACGTTTGCGGCGATATTTTGGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTa                                                                                                                                   <  2:256648/200‑1 (MQ=255)                                                                                                                                                 gCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTTTAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAag                                                                               >  1:267458/1‑247 (MQ=255)                                                                                                                                                       gacgTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGCTTTACGCTATAACATGGAGCTCCCCCTCGAAGAAGCTGTACGTGGCGTTACCAAAGAGATCCGcct                                                                     >  2:215603/4‑249 (MQ=255)                                                                                                                                                       gAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCGCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCAt                                                                     <  2:60997/251‑1 (MQ=255)                                                                                                                                                                                    ggCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGcc                                        <  1:63748/251‑1 (MQ=255)                                                                                                                                                                                    ggCATAAGATCAGAATGGTCATGCTGCGTTTGAGCAAGGAGGGATGGGCGGCGGCGGTTCTGGCGGCGGCGCAGACTTCAGCGATAATTTTGGGGACGTTTTAGGCGATATTTTTGGCGGCGGACGTGGAAGTCAACGTGCGGGGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGGGGCGTGACCAAAGAGAACCGCAATACGACTATGGAAGAGGGTGACGTTTGc                                         >  2:178294/1‑250 (MQ=255)                                                                                                                                                                                               aGTATGGTCATGCTGCGTTTGAGCCAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGAACCGCATTCCGACTCTGGAAGAGTGTGCCGTTTGCCACGGTAGCGGt                             >  2:200153/1‑251 (MQ=255)                                                                                                                                                                                                                      cAAGGTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGGCTTCAGCGATATTTTTGGTGACGTTTTAGGCGCTATTTTTGGCGGCGGACGTGGTAGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGGGCTCACCCTCGAAGTAGCTGTACGTGGCGTGACCAAAGAGAACCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGGAGCGGAGCAAAACCAGGTACACAGACGc       >  1:207122/1‑250 (MQ=255)                                                                                                                                                                                                                         ggtggtATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGAc   <  2:4871/251‑1 (MQ=255)                                                                                                                                                                                                                          gTGGTATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGAcc  <  2:224028/251‑2 (MQ=255)                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       GATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACT  >  NC_000913/14180‑14647

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