Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,879	T→C	Y239Y (TAT→TAC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,879	0	T	C	92.3%	62.9 / ‑3.5	26	Y239Y (TAT→TAC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (13/11);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCA  >  NC_000913/12651‑13111                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          ggAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGATAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTTGTGGTACTTTCGATATTTCTATTAACGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCGGGTTGAAGAATTCaagaaa                                                                                                                                                                                                                      <  1:291429/250‑1 (MQ=255)           cGTATCATCAACGAAACGACCGAAGCTGAGCTGGCTTACGGTCTGGACAAAGGAAATGGCAACCGTACTATCGCGGTTTATGACCTGGGAGGTGGAACTTTCGGTAATTCTATTAGCGAAATCGAAGAAGTTGACGGCGCAAAAACCTTCGAAGTTCTGGCAACAAACGGGGATAAACACCTGGGTGGTGCAGACTTCAACCACCGTCTGGTCAACGACCGTGTTGAAGAATTAAAGAAAGATCAGGGCa                                                                                                                                                                                                            >  2:331251/1‑250 (MQ=255)                                gCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCt                                                                                                                                                                                      <  1:344375/251‑1 (MQ=255)                                  aGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTgg                                                                                                                                                                                    >  1:328475/1‑251 (MQ=255)                                    ctgcaatGGCGTACGGTCTGGACACAGGTACTGGCAAACGTACTATCGAGGTTTATGATCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGCCGACGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGGGAACCCCACCTGGGTGGGGCAGCCTTCGACAGACGTCTGATCAACTACCTGGTTGAAGAATCAAAGAAAGATCAGGGCATAGACCTGCGAAACgtaccgcaggc                                                                                                                                                                                   >  2:328871/7‑240 (MQ=255)                                       cgcTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACCACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATg                                                                                                                                                                               <  1:206106/251‑1 (MQ=255)                                                      tGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAaaga                                                                                                                                                                 >  1:6803/1‑250 (MQ=255)                                                       ggACAAAGGTACTGTCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAaagaa                                                                                                                                                                <  2:6803/250‑1 (MQ=255)                                                              ggtaCTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGAGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGAGGGAGAAAAAACCTTCGAAGTTCTGGCAACCAACGGAGATACCCACCTGGGAGGAGAAGACTTCGTCAGCCGTCTGATCAACCACCTAGTTGAAGAATTCAAGAAAGATCAGGGCATAGACCTGCGCAACGCTCCGCTGGCAAAGCAGCGTCTGAaagacacggc                                                                                                                                                           >  2:58674/4‑242 (MQ=255)                                                                        aCCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAACGCGaa                                                                                                                                               >  1:265289/1‑250 (MQ=255)                                                                             aCTAGCGCGGTTTATGACCTGGGGGGTGGTACTTTCGATATGTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAAGTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCg                                                                                                                                          <  2:166570/250‑1 (MQ=255)                                                                                aTCGCGGGTGATGACCTGGGTGGTGGGATTTTCGATATTTCTATTAGCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTACTTGCAACCAACGGTGATACCCCCCTGGTTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCTCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACt                                                                                                                                      <  1:214768/251‑1 (MQ=255)                                                                                         tATGACCTGGGCGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGCAGTTCTGGCAACCAACGGTGATACCCACCTGGGTAGTGAAGCCTTAGAAAGCCGTATGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCATGGCCTGGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGAGCAAATCGAACCGTCTTCCGc                                                                                                                             >  2:375824/1‑251 (MQ=255)                                                                                                               tCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑t                                                                                                       <  1:258960/250‑1 (MQ=255)                                                                                                                        cTATTATCGAACTCGAAGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCGACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATATGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGAAATGCAGCGCCTGAAAGACGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGT‑AACCATGCCAAACATc                                                                                             >  2:200924/1‑251 (MQ=255)                                                                                                                                  aaTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATAACCACCTGGGTGGTGAAGACTTCGACCACCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGatg                                                                                   >  1:241709/1‑249 (MQ=255)                                                                                                                                   aTCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGatgc                                                                                  >  1:68465/1‑248 (MQ=255)                                                                                                                                          aatttGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGAGCATGAAAGAAGCGGCAGAAAAAGAGAACATCGAACTGTCTTCCGCTCAGCAGACCGACGTTACCC‑TGCCGTAcctaactgaaactgtgaccgtt                                                                           >  2:297256/4‑230 (MQ=255)                                                                                                                                                aCGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGGGAAGACTTCGAAAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGaaa                                                                     <  1:231965/251‑1 (MQ=255)                                                                                                                                                         aaaCCTGCGAAGTTCCGGCAACCAACGTGGATACCCACCTGGGGGTTGAAGACTGCGACTGCCTGCTGAACAACTGTCTGGTTGAAAAATACAAGAAAGATCAGTGCATGGACCGGAGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGACCTGCCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAAc                                                            <  2:34224/251‑1 (MQ=255)                                                                                                                                                                  aaGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTg                                                   >  1:277867/1‑251 (MQ=255)                                                                                                                                                                       cgggCAACCATCGGTGATACCCACCGGGGTGGTGAAGACTTCGACAGCCGTCTGTTCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTc                                               <  2:314493/248‑1 (MQ=255)                                                                                                                                                                           caaccaacGGTGAGACCCGCCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGGCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCg                                          <  2:59135/251‑1 (MQ=255)                                                                                                                                                                                                tGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGTGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACAGACGTTAACC‑TGCCGTACATAACTGCAGATGCGACCGGTCACAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTgc                      >  2:296050/1‑249 (MQ=255)                                                                                                                                                                                                         aaGACTTCGACAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGt             >  1:143963/1‑250 (MQ=255)                                                                                                                                                                                                                   cAGCCGTCTGATCAACTACCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCGTACATCACTGCAGATGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCa  <  1:329133/251‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                          GGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATCTGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGAAGCGGCAGAAAAAGCGAAAATCGAACTGTCTTCCGCTCAGCAGACCGACGTTAACC‑TGCCATACATCACTGCAGACGCGACCGGTCCGAAACACATGAACATCAAAGTGACTCGTGCGAAACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCA  >  NC_000913/12651‑13111

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