Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	13,539	A→G	L459L (CTA→CTG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	13,539	0	A	G	91.2%	97.6 / NA	35	L459L (CTA→CTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (17/14);  minor base T (0/3);  total (17/18)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAC  >  NC_000913/13306‑13782                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    cTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGc                                                                                                                                                                                                                                    >  1:212316/1‑251 (MQ=255)   gACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCa                                                                                                                                                                                                                                   <  1:347458/250‑1 (MQ=255)     cTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcg                                                                                                                                                                                                                                >  1:304742/1‑251 (MQ=255)              tgataGACGTAATGCTGCTTGACATTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGTTGTGATGACGACACTGATCGCGAACAACACCACTATCCCGACCAAGCACAGCCAGGAGGTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGCTGAACGTAAACGTGCGGCTGATACCAAATCTCTGGCTAAGTTCAACCTGGCGGGGATCAACCCgactccgcccgtcc                                                                                                                                                                                                                         >  2:5858/5‑236 (MQ=255)               gaaggACGTACAGCTGCTGGACGTTACCCCGCTGTCTCTGGGTACCGAAACCATTGGCGGTGTGATGACTACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCATGTGGTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCAGCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGc                                                                                                                                                                                                                      <  2:143963/247‑1 (MQ=255)                aaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGcc                                                                                                                                                                                                                     >  1:244576/1‑251 (MQ=255)                aaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGGAAACGTGCGGCTGATAACAAAGCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGcc                                                                                                                                                                                                                     <  1:366959/251‑1 (MQ=255)                   gACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGc                                                                                                                                                                                                                   >  1:202242/1‑250 (MQ=255)                      gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGc                                                                                                                                                                                                                   >  1:359744/1‑247 (MQ=255)                                    cGTTACCCCGCTGTCTCTGGGTTGCGAAACCAGGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGAGGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCg                                                                                                                                                                                                 <  2:359744/251‑1 (MQ=255)                                               tGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCt                                                                                                                                                                                      >  1:242680/1‑251 (MQ=255)                                                    ctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGAc                                                                                                                                                                                   <  1:356191/249‑1 (MQ=255)                                                        gTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGCACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGTTATCGATGCTGACGGTATc                                                                                                                                                                             >  1:106447/1‑251 (MQ=255)                                                                aCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCCGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGt                                                                                                                                                                     >  1:381055/1‑251 (MQ=255)                                                                         ggTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGa                                                                                                                                                             >  1:225120/1‑250 (MQ=255)                                                                                    gacgCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCACGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGCTaaaaaa                                                                                                                                                  >  2:61014/1‑249 (MQ=255)                                                                                                      ccccacTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCTCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCagaaga                                                                                                                               <  1:193250/249‑1 (MQ=255)                                                                                                                  gACCAAGCACAGCCCGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGTTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCCCCATCAAgg                                                                                                                   <  1:288287/251‑1 (MQ=255)                                                                                                                    ccAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGct                                                                                                                 <  1:320708/251‑1 (MQ=255)                                                                                                                     cAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGctt                                                                                                                >  1:32605/1‑251 (MQ=255)                                                                                                                            aGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCTGCTGATAACAAATCTCTGGGTCAGTTCAACCTTGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCGTCGATATCGATGCTGACGGTATCCTGCACGTTTACGCTAAAGATAAAAACAGCGGTAAAGAGTAGAAGATCACCAGCAAGGCTTCTTCTgt                                                                                                         <  1:136233/251‑2 (MQ=255)                                                                                                                               cAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCt                                                                                                      >  1:313312/1‑251 (MQ=255)                                                                                                                               cAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCt                                                                                                      <  1:83514/251‑1 (MQ=255)                                                                                                                                             cgctgaagccccacAGTCGGCGGTAACCATCCATGTGCTGCAGGGCTAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTCGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTCCCTTCGATATCGACGCTGACTGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGTTAAAGAGCAGAAGATCACCATCAAGGTTTCTTCTGTTCTGAACGACGATTaaa                                                                                        <  2:112643/238‑1 (MQ=255)                                                                                                                                                  aaGACAACCAGTCTGCTGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGTTTAGTTCAACCTGGATGGTATCAACCCGGCTCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATTGATGCTGACGGTATCCTGCACGTTTCCGCGAAATATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTGCTTCTGGTCTGAACGAAGATGAAATCCa                                                                                    <  1:342954/250‑1 (MQ=255)                                                                                                                                                            gTCTGCGGTAACCATCCAAGTGCTGCAGGGTGAACGGAAACGTGCGGCTGAGAACAAATCTCTGGGTCAGTTCAACCTGTATGTGATCAACCCGGCGCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTGCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATTACATCCAGAAAATGGTAc                                                                         <  1:151571/251‑1 (MQ=255)                                                                                                                                                                 cGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATggatggttac                                                                    >  1:174807/1‑243 (MQ=255)                                                                                                                                                                   gTAACCATCCATGTGCTGAAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATggatggttaccg                                                                  <  2:174807/251‑11 (MQ=255)                                                                                                                                                                     aaCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCagaaga                                                                                                                               <  2:18764/188‑1 (MQ=255)                                                                                                                                                                     aaCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCagaaga                                                                                                                               >  1:18764/1‑188 (MQ=255)                                                                                                                                                                                                          gCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACAGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCCCGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAAAGGTACGCGACGCAGACGCTAACGCAGAAGCTGACCGTAAGTTTGAAGAGCt                           >  2:223601/1‑251 (MQ=255)                                                                                                                                                                                                              aTAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGt                        >  1:283729/1‑250 (MQ=255)                                                                                                                                                                                                               tAACAAAGCTCTGGTGCAGTTCAACCCTGATTGTATCAACCCGGCACCGCGCCGCCTGATGCAGATCGAAGTTACCTTCGATATCGATGCTTACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTACTGGTCTGAACGAAGATGAGATCCAGAAAATGGTACGGGTCGCAGAAGCTAACGCCGATGCTGACCGTAAGTTTGAAGAGCTggga                       <  2:313312/250‑3 (MQ=255)                                                                                                                                                                                                                  cAAATCACTGGGGCAGTTCAACCATGATGGCATCAACCCGGCTCCGCGCGGCATGCTGCAGATCGAAGGTACCGTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGAGAAAGATAAAAACAGCGGTAAAGAGCAGAATATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGAAGACGCTAACGCTGAAGCTGAACGTAAGTTTGAAGAGCTGTTACAGa                   <  2:212316/251‑1 (MQ=255)                                                                                                                                                                                                                                   ccAACCTGGATGGTATCAACCCGGCCCCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTAGCCTGCACGTTGCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCAGCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAc  <  2:381055/250‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    CTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGCATGCCGCAGATCGAAGTTACCTTCGATATCGATGCTGACGGTATCCTGCACGTTTCCGCGAAAGATAAAAACAGCGGTAAAGAGCAGAAGATCACCATCAAGGCTTCTTCTGGTCTGAACGAAGATGAAATCCAGAAAATGGTACGCGACGCAGAAGCTAACGCCGAAGCTGACCGTAAGTTTGAAGAGCTGGTACAGACTCGCAACCAGGGCGAC  >  NC_000913/13306‑13782

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