Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	14,424	G→A	G86D (GGC→GAC)	dnaJ →	chaperone protein DnaJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	14,424	0	G	A	93.8%	35.9 / ‑4.1	17	G86D (GGC→GAC)	dnaJ	chaperone protein DnaJ
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (4/11);  total (6/11)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

CGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGG  >  NC_000913/14214‑14670                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       cGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTggaatcggaagagc                                                                                                                                                                                                                 >  1:11239/1‑238 (MQ=255)                                aCGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAAGTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGGGCGGCATACTATCAGGATGTTCATGCTGCGTTGGAGCAAGGTGTCATGTGCGGCGTCGGTTTTGGCGACTGCGCATACTTCTGCGATAGTTTTGGTTACGTTTTCTGCGATAGTTTTGGCGGCGGACGTGGTCGTCAAc                                                                                                                                                                                 <  1:454721/251‑1 (MQ=255)                                      ggCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTgc                                                                                                                                                                             <  1:175594/249‑1 (MQ=255)                                       gCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTgcggc                                                                                                                                                                          <  1:353372/251‑1 (MQ=255)                                                                                 gCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCGGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGTTCATGCTGCGTTTGAGCAAGGTGGCATGGGCTGCGGCTGTTTTGGCGACTGCGCAGACTTCAGCGATATGTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACg                                                                                                                                                                                <  1:346558/203‑1 (MQ=255)                                                                                       gCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTggn                                                                                                                                                                                                    <  1:95518/177‑2 (MQ=255)                                                                                                           ggAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAag                                                                                                      <  1:289625/251‑1 (MQ=255)                                                                                                                         ttCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGACGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCTGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGc                                                                                                                <  1:492398/227‑1 (MQ=255)                                                                                                                           cTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGCCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGc                                                                                                                >  1:485123/1‑225 (MQ=255)                                                                                                                               ccGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCCTTCCGACTCTg                                                                                  >  1:64564/1‑251 (MQ=255)                                                                                                                                 gACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTgg                                                                                 >  1:24933/1‑250 (MQ=255)                                                                                                                                     cGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTa                                                                                                                                                          >  1:420793/1‑173 (MQ=255)                                                                                                                                                   gCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCGGCGATATTTTTGGTGACGGTTTAGGCGTTATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGGGTGACg                                                                     >  1:114392/1‑244 (MQ=255)                                                                                                                                                                        gctgcGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTGTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATAGTTTTGGTGGCGGACGTGGTCGGCAACGTGCGGCGCGCGGTGCTGATTGACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGGGCGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGt                                                          <  1:283073/234‑1 (MQ=255)                                                                                                                                                                                           gCATGGGCGGCGGCGGTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATAtttttn                                                                                                                                                                                                      <  1:93407/75‑2 (MQ=255)                                                                                                                                                                                                      ggcggTTTTGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAgtg                                                                          <  1:413918/188‑1 (MQ=255)                                                                                                                                                                                                          gtgttGGCGACGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTGTGGCGACGGACGTTGCCGTTACCGTGCCGCGCGCGGTGCAAATTTACGCTATAGCATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTg                                                                                                             <  1:448146/146‑1 (MQ=255)                                                                                                                                                                                                               ggcgacggcGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCCCCCTCGAAGAAGCTGTACGAGGCGTGCCCATAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCACAGACTTGTCCGACCTGTCATGGTTCTgg  >  1:485798/1‑251 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       CGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAGGGTGACAAAGAGGCCGAGGCGAAATTTAAAGAGATCAAGGAAGCTTATGAAGTTCTGACCGACTCGCAAAAACGTGCGGCATACGATCAGTATGGTCATGCTGCGTTTGAGCAAGGTGGCATGGGCGGCGGCGGTTTTGGCGGCGGCGCAGACTTCAGCGATATTTTTGGTGACGTTTTCGGCGATATTTTTGGCGGCGGACGTGGTCGTCAACGTGCGGCGCGCGGTGCTGATTTACGCTATAACATGGAGCTCACCCTCGAAGAAGCTGTACGTGGCGTGACCAAAGAGATCCGCATTCCGACTCTGGAAGAGTGTGACGTTTGCCACGGTAGCGGTGCAAAACCAGGTACACAGCCGCAGACTTGTCCGACCTGTCATGGTTCTGG  >  NC_000913/14214‑14670

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