Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,866	G→A	A568A (GCG→GCA)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,866	0	G	A	81.8%	38.8 / 2.6	23	A568A (GCG→GCA)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base G (1/3);  new base A (4/14);  total (5/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-02

CCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATA  >  NC_000913/13622‑14089                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  ccGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGTCGACCATCTGCTGCACCGCACCCGTAAGCAGGTTGAAGAAACAGGCGCAAAACTGCCGGATGACGACAAAACGGCTATCGAGTCTGCACTGACt                                                                                                                                                                                                                             >  2:284215/1‑251 (MQ=255)                 cAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCAAGGCTGCACTGACTGc                                                                                                                                                                                                                           <  1:112775/237‑1 (MQ=255)                 cAGCGGAAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGCTGAAATCCAGAAAATGGTACGCAACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGACGCAGGCGAAAAACTGACGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGc                                                                                                                                                                                                                           >  2:112775/1‑237 (MQ=255)                                      aaCCATCAAGGCTTTTTCTGGTCTGAACGAAGATTAAATCCAGAAAATGGTACGCGACGCCTAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCGGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACGGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCg                                                                                                                                                                                        <  2:242680/250‑1 (MQ=255)                                                              gAACGAAAATGAAACCCAGAAAATGGTACTCGTCGCAAAAGCTAACGCCGAAGCTAACCGTAAGGTTGAATAGCTGGTCCATACTCGCAACCAGGGCGACCAGCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCCCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGTCAAAGCCGCTATCGAATCGCAAAttcatgagc                                                                                                                                                                <  2:202242/251‑9 (MQ=255)                                                                 cGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTgg                                                                                                                                                             <  1:167294/251‑1 (MQ=255)                                                                                             gaagCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAAt                                                                                                                                 <  1:5858/248‑1 (MQ=255)                                                                                                         aaCGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGGCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATC‑AAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAgcagca                                                                                                                     <  1:154448/250‑1 (MQ=255)                                                                                                                     gACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCcagc                                                                                                          >  1:45812/1‑250 (MQ=255)                                                                                                                                   agagCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGTACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTg                                                                                           <  1:199161/251‑1 (MQ=255)                                                                                                                                                                       gctgcACAGCACCCCTAAGAAGGTTGAATAAGCATGCGACAAACTCCCGGCTGACGACAAAAATGCTATCGAGTCTGCACTGACTGCACTGGAAACGACTCTGAAAGGTGAAGACAATGCCGATATCGAAGCGAAAATGCAGGAGCTGGCACAGGATTCCCAGAAGCTGATGGAAATCGCCCATCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGaa                                                         <  2:141439/247‑2 (MQ=255)                                                                                                                                                                               gCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTTAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑TAACGCGAAAGATGTCGATGGGGTCGAc                                             <  2:161681/251‑1 (MQ=255)                                                                                                                                                                                   ccGTAAGCAGGTTGAAGAAGCTGGCGACAAACTGCCGGCGGACGACAAAACTTCGATCGAGTCTGCACTGACTGCACTGGAAACGGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGATCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACGGCCTGTGCTGATGCTTCTGTA‑AA‑CAACGCGAAAGAGGACGATGTTGTCgtcgatg                                         <  2:244576/251‑7 (MQ=255)                                                                                                                                                                                        agccgggtGAAGAGGCAGGCGACAAACTGCCGGCTGACGACAAAACTTCTATCGAGTCTGCACTGACTGCACTGGAAACGGCTCTGACAGGTCAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCTCAGGTTTCCCAGAATCTGATGTAAATCTCCCCGCAGCAACATGCCCAGCAGAAGGCTGCCGTTGCTGATGTTTCTGCA‑AA‑CAACGCGAAATATGACGATGTTTTCgtagttgaattt                                    <  2:47017/244‑12 (MQ=255)                                                                                                                                                                                          cAGGTTGAAGAAGCAGGCGACAAACGGCCGGCGGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTga                                  <  1:49605/251‑1 (MQ=255)                                                                                                                                                                                           aGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAGCGCGAAAGATGACGATGTTGTCGACGCTGAATTTgaa                                 <  1:295654/251‑1 (MQ=255)                                                                                                                                                                                               gaaagaagCAGGCGACAAACTGCCAGAGTACGAACAAACAGTTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATAGAAGCGACAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGTCGCTGAATTTGaagaag                             <  2:32605/249‑1 (MQ=255)                                                                                                                                                                                                       cAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGAc                     <  1:284215/251‑1 (MQ=255)                                                                                                                                                                                                              cAAACTGCCGGCTTACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTATGAAAGGTGAAGGCACAGCCGCTATCGAAGCGAAAATGCTGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTGCT‑CATAATCAA‑GCGAAAGAGGACGATGTTGTCGCCGCAGAAATGGAAGAAGTCAAAGaaaaaaca               >  2:2747/1‑248 (MQ=255)                                                                                                                                                                                                                aaCTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAaataa            <  1:362303/251‑1 (MQ=255)                                                                                                                                                                                                                aaCTGCAGGCTGACGACAAAACTGCTATCGAGTCTGCACTGACTGCAATGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAaata             >  1:333756/1‑250 (MQ=255)                                                                                                                                                                                                                       ggcggacgacAAAACTGCTATCTAGTCTGCACTTACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAGCTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGccca     <  2:4071/247‑2 (MQ=255)                                                                                                                                                                                                                           gacgacAAAACTGCTATCGAGTCTGCACATACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGGCGCTAACGAAGCGAAAATGCAGGAGGTGGCACAGGTGTCACAGAAACTGAGGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATTACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGTCCtaga  <  2:283729/250‑3 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                  CCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGACCATCTGCTGCACAGCACCCGTAAGCAGGTTGAAGAAGCAGGCGACAAACTGCCGGCTGACGACAAAACTGCTATCGAGTCTGCGCTGACTGCACTGGAAACTGCTCTGAAAGGTGAAGACAAAGCCGCTATCGAAGCGAAAATGCAGGAACTGGCACAGGTTTCCCAGAAACTGATGGAAATCGCCCAGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCA‑AA‑CAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATA  >  NC_000913/13622‑14089

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