Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,507	T→C	I115I (ATT→ATC)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,507	0	T	C	85.3%	72.1 / 5.2	34	I115I (ATT→ATC)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (19/10);  total (21/13)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.48e-01

ACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACC  >  NC_000913/12268‑12724                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          acgctttcttccatgcccTATACCCAGGAGGGTGAACCGCTGGTTGGCCAGCGGGATAAACG‑CACGGCAGGGA‑AGAACCCCAAAAACACCCTG‑TTTGCGATTAAACGCCTGATGGGCCGCCGCTCCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTCAAAAATTTTTGCTGCGGATACCGGCGCCGCATGGGGCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGt                                                                                                                                                                                                                <  1:150574/235‑1 (MQ=255)  cGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCATCGCATGGGTCGAATTTAAAGGCCAGAAAATGGAACCGCCGCAGATTTCTGCTGAAggg                                                                                                                                                                                                               >  2:270393/1‑249 (MQ=255)    ccTTCTATCATTGCCTAGACCCAGGATGGTGAAACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGAAAAACCCCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCAGGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCt                                                                                                                                                                                                             <  1:388563/251‑1 (MQ=255)                        ccAGGATGGTGACACTCTGGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGAGCTGAAAAAAATGAAGAAACCCg                                                                                                                                                                                         >  2:443157/1‑251 (MQ=255)                              gggaaaaaCTCTGGTTGGTCGCCCGGCTAAACGTCG‑GACAGGAA‑CGAACCCGCAAAACCCCCTG‑TTTGCGATTAAACCCCGTATTGGCCGCCCCTTCAGGGACAAAAAAGTACAGCG‑TTATTTTTCCATCATGCCGTTAAAAATGATTGCTGTGGAAAACGGCGACGAAGGGGTCAAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAg                                                                                                                                                                                   <  1:227508/246‑1 (MQ=255)                                         tggTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGt                                                                                                                                                                        >  1:129539/3‑251 (MQ=255)                                          gtttGGCCAGCCGGCTAAACGTCA‑GGCAG‑GAGCCACCCCTAAAAACCCCCTG‑TTTGCGATTAACCGCCTGATGGGCCGCCGCTTCCAGAACGAAAAAGTACAGCG‑GGATGTTTCCATCATGCCGGTAAAATTTATTGCTGCTGATACCGGCGCCGCTTGGGCCAAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTg                                                                                                                                                                       <  1:326692/249‑1 (MQ=255)                                                     cGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACt                                                                                                                                                            >  1:88220/1‑251 (MQ=255)                                                         aaaaCGTGA‑GGCAGGGA‑CGACCCCGAAAAACACCCTG‑TTTGCGATAAAACGCCTGATTGGCCGCCTCTTCCAGGACGAAAAAGTACAGCG‑TGATTTTGCCATCATGCCTTTCAAAATTATTGCTGCTGATAACGGCTACGCATGGGTCGAAGTAAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAg                                                                                                                                                        <  1:627679/250‑1 (MQ=255)                                                          aaaCGTCA‑GGCAGTGA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGc                                                                                                                                                       >  1:678034/1‑251 (MQ=255)                                                              gTCA‑GGCAGTGA‑CGAACCCG‑AAAACCCTTTGGTTGGCGATAAACCGCCTGTTTGGCCGCCGCTTCAAGGACGAAGAAGTACAGGG‑TTATGTTTCAAACATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGAATGGTTCAAAGTAAAAGGCCAAAAAATGGCTCCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGtt                                                                                                                                                   <  1:270393/251‑1 (MQ=255)                                                                         gA‑CGAACCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTAc                                                                                                                                         >  2:314256/1‑251 (MQ=255)                                                                               aCCCGCAAAACACCCTG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATAACCGTACCGGc                                                                                                                                     >  2:280995/1‑250 (MQ=255)                                                                                         caccctG‑TTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACg                                                                                                                          >  1:130886/5‑251 (MQ=255)                                                                                                   tGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAAGATGCTCAGc                                                                                                                 >  2:180208/1‑251 (MQ=255)                                                                                                              gCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGAtt                                                                                                                                                                                >  1:536516/1‑177 (MQ=255)                                                                                                                 tGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCCGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaa                                                                                                    >  1:472739/1‑250 (MQ=255)                                                                                                                  ttggggccgccgtTTCCAGGTTGAATAATTACGGCG‑TGATGTTTACATTATGCGGTTCAAATTTATTGCTGTTGATAACGGAGTCGCATGGGTTGAAGTTAAAGGCAAAAAAATGTCTCCGCCGCAGATCCCTGCTGAAGTGCCGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaaa                                                                                                   <  2:472739/243‑1 (MQ=255)                                                                                                                      ggccgccgcTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGc                                                                                              >  1:626074/4‑251 (MQ=255)                                                                                                                       gcccccgcTTCAAGGACGAAGAAGT‑CAGGGGTGATGTTTCCACCATTCCTTTCAAAATAATAGCTGTGGAAAACGGCGACGCAGGGGTCGAAGTAAAAGGCAAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCa                                                                                             <  1:485011/247‑1 (MQ=255)                                                                                                                                   aGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCCCCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCt                                                                                 >  1:35120/1‑251 (MQ=255)                                                                                                                                                 cAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGCGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTa                                                                                                                                          >  2:335786/1‑180 (MQ=255)                                                                                                                                                 cAGCG‑TGATGTTTCCATCATGCCGTTAAAAATTATTGCTGCTAAAAACGGCGACGCATGGGTCAAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGCGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTa                                                                                                                                          <  1:335786/180‑1 (MQ=255)                                                                                                                                                         aTGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTatc                                                            >  1:322370/1‑251 (MQ=255)                                                                                                                                                                  tcatGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGCCCAGAAAATGGCACCGACGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGTGTGAACCGGTAACAGAAGCTGTTATCACCGTACCAGCATACCTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCCGGCCGTATCGCTGGTCTGGAAGTAACACGTATCCTCAACGa                                                    >  1:120422/1‑250 (MQ=255)                                                                                                                                                                  acagcccGTTAAAAATTATTGCTGCTAAAACCGGCGACGAGTGGGTGGAAGTTAAAGGCCAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGaa                                                   <  1:443157/246‑1 (MQ=255)                                                                                                                                                                          ataaaaaGTATTGCGGCTGATAACGGAGACGCATGGGTCGAAGTAAAAGGAAAGAAAAGGACACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAaccgacc                                             <  1:25964/246‑1 (MQ=255)                                                                                                                                                                               aaTTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTg                                      >  1:475883/1‑251 (MQ=255)                                                                                                                                                                                attattGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTgc                                     >  1:366354/1‑251 (MQ=255)                                                                                                                                                                                       ctgctgATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACATGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATAAACGAACCGACCGCAGCTGCGCTGGCt                              >  1:344280/1‑251 (MQ=255)                                                                                                                                                                                        tgctgATAACGGCGACGCATGGGGCGAAGTTAAAGGCCAGAAAATGGCTCCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCCAAGCCGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCCGCTGCGCGGGCtt                             >  2:650643/1‑251 (MQ=255)                                                                                                                                                                                                    gagagctgtggtagAAGTAAAAGGCAAAAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACa                 <  2:35120/238‑1 (MQ=255)                                                                                                                                                                                                      acgaaTGGGTCGAATTAAAAGGCAAAAAAAGGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACaaa               <  1:180208/247‑1 (MQ=255)                                                                                                                                                                                                                   aGTTAAAGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCaccc  >  2:7442/1‑248 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          ACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAGTTGGTCAGCCGGCTAAACGTCA‑GGCAGTGA‑CGAACCCGCAAAACACTCTG‑TTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCG‑TGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACC  >  NC_000913/12268‑12724

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