Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	11,925	A→C	intergenic (‑139/‑238)	yaaI ← / → dnaK	DUF2541 domain‑containing protein YaaI/chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	11,925	0	A	C	91.2%	90.4 / ‑0.5	34	intergenic (‑139/‑238)	yaaI/dnaK	DUF2541 domain‑containing protein YaaI/chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (16/15);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.45e-01

TGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGAT  >  NC_000913/11688‑12164                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                tGCAATGACGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCGGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCGCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGGTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGAGTCATT‑CTTTATCTGAATaaaactg                                                                                                                                                                                                                                     <  2:238904/250‑4 (MQ=255)   cAATGACGCCGAGGAGATTATGGACGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGGTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGCTGAATGCCTTGGCTGCGATTCATT‑CTTTTTCTGAATAAAATTGCTg                                                                                                                                                                                                                                  >  2:376013/1‑251 (MQ=255)                        ggTCGGTTGCCCGCGCCGTGCAGCCAAGCATCAGGGTAATCGCCCGGCTGGGGGAAATCGGAAAAACGGATTTCATAAGGATTTTCTTAGTGGGAAGAGGTAGGGGGAAGAATTCCCAATAGTTTTCTGCTGATAAAGAGAAGATTCAGGGACCTAATCTTTTTTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTg                                                                                                                                                                                                             <  2:290778/251‑1 (MQ=255)                                         gTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAAAGGATTTCATAAGGATTCTCTTAGTGGGCAGAGGTAGAGGGATGCATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGTATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAAtttt                                                                                                                                                                                            >  1:282808/1‑250 (MQ=255)                                               cacagcaTCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACATTTTTTTGATGAATGCCTTGGCTGCGTTTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTAAGTCTTGTCCTGCCATATCGCGCAATTTCTgcgc                                                                                                                                                                                       >  2:318996/1‑250 (MQ=255)                                                         ggCTAATCGCCAGGCTGGCGTAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGTTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCCGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACaaa                                                                                                                                                                            <  1:174847/251‑1 (MQ=255)                                                          gCTAATAGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACaaaa                                                                                                                                                                           <  1:198659/251‑1 (MQ=255)                                                             aaTCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACaaaaaa                                                                                                                                                                         <  1:261247/250‑1 (MQ=255)                                                                cGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGACTAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAAtttt                                                                                                                                                                     >  1:292306/1‑251 (MQ=255)                                                                     ggCTGGCGGAAATCGTCAAACCGGATTTAATACGGATTCTCTTAGTGGGAAGAGGTAGAGGATTGAATACCCACTAGTTTACTGCTGCTAAAGAGAAGGTTCCGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGTTTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTATGCGAAACAGCACAAAAAATTTTTGCAt                                                                                                                                                                >  2:313007/1‑251 (MQ=255)                                                                                  agTAAAAACGGTGTTCATAGGGATTATCTTAGTGGGAAGAGGGAGAGGGATGAATACCCACTATATTTCTACAGATAAAGAGAAGGTTCAGGTACGTAATCTTTTTTTTTTATTACAATGATTTAATGAATGCCTTGGCGGCGATCCATT‑CTTTATCTGAATAAAATTGCTGTAAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTtgat                                                                                                                                                    <  2:80548/249‑1 (MQ=255)                                                                                        accGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAAt                                                                                                                                                                                                                             >  2:351719/3‑171 (MQ=255)                                                                                        aaCGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAAt                                                                                                                                                                                                                             <  1:351719/171‑1 (MQ=255)                                                                                            gATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGACTGCGATTCATT‑CTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAA‑TTTTTGCATCTCCCCCTTGATGACGTGGTTTAc                                                                                                                                        >  1:16881/1‑251 (MQ=255)                                                                                              tttCATAAGGATTCTCTTAGTGGGAAGAGGTAGAGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACg                                                                                                                                       >  1:198251/1‑251 (MQ=255)                                                                                                  aTAAGGATTTTCTTAGTGGGAAGAGGTAGAGGGAAGAATACCCACTAGTTTACTGCGGATAAAGAGAATTGTCAGTTACGTAATCTTTTCTTTTTATTACAATTTTTTGAGGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGCCAATTTTACTTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCGCCCCCTTGATGACGTGGGTTACGAccc                                                                                                                                   <  2:301696/251‑1 (MQ=255)                                                                                                                           gtatagggATTAATACCCACTAGTTTACTTCTGAAAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCagtgagt                                                                                                          <  1:405988/246‑1 (MQ=255)                                                                                                                              gagggaggaaacccccttagagtACTCCTGATAAATAGAAGGTTCTGTTACGTAACCTTTCCGTTTTATTACAATTTTTTGAGGAATGACTCGGCTGTGATTCCTT‑CTTTTTCTGAATAAAATTGTGGTCATTTTTCCGTATTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCTCCTTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGtgagtgag                                                                                                       <  2:141778/229‑1 (MQ=255)                                                                                                                                  gATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGGTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTg                                                                                                   >  1:70318/1‑251 (MQ=255)                                                                                                                                                           gataaagagaaggttagggtactgaaTCTTTTCTTTTTATTACAATTTTTTGATTAATGCTTCGGCTGTGCATCATTTCTTT‑TCTGAATAAAATTGCTGTCAATTTTACGCTTTGTCCGGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAAGTGTTGTATCTCCACCTTGATGCCGTGGGTTCCGACCCCATTTAGTAGTTAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTg                                                                           <  1:39028/226‑1 (MQ=255)                                                                                                                                                                  agaaggTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCAtt                                                                                                                               >  1:53886/1‑191 (MQ=255)                                                                                                                                                                  agaaggTTCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCAtt                                                                                                                               <  2:53886/191‑1 (MQ=255)                                                                                                                                                                    aaggttCAGGTACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACg                                                                 <  1:259878/247‑1 (MQ=255)                                                                                                                                                                                               ttattaCAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAAttttt                                                                                                                                                                    >  1:257176/1‑125 (MQ=255)                                                                                                                                                                                               ttattaCAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAAcc                                      >  1:204541/1‑251 (MQ=255)                                                                                                                                                                                               ttattaCAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTATTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAAttttt                                                                                                                                                                    <  2:257176/125‑1 (MQ=255)                                                                                                                                                                                                  ttaCAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACcac                                    >  1:395683/1‑250 (MQ=255)                                                                                                                                                                                                      aaTTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACatgat                               >  2:217678/1‑251 (MQ=255)                                                                                                                                                                                                         tttttGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAGAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAACATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCCCAACCACatgatgcc                            >  1:66905/1‑249 (MQ=255)                                                                                                                                                                                                            ttGATGAATGCCTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCgc                         >  1:332322/1‑250 (MQ=255)                                                                                                                                                                                                                   gtGCCTTGGCGGCGATTCATT‑CTTTATATGAATACATTTGCTGTCAATTTTACGTCTTGTCCTGCCATTTCGAGACATTTCTGCGCAAAAGAACAAAAAATTTTTGCATCTCCCCCTTGCTGTCGTGGGTTACGACCCAATTTAGTAGTCAACCGCAGTGAGTGAGTCTGAAAAAACATGAAATTGGGCAGTTGAACCCAGACGTTTCGACCATATTACAAACTCAAAACCACATTATGACCGAatatat                   >  2:49330/2‑250 (MQ=255)                                                                                                                                                                                                                      ccTTGGCTGCGATTCATT‑CTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTg               >  1:322606/1‑251 (MQ=255)                                                                                                                                                                                                                           gcgtcGATTCATT‑CTTTATCTGAATAAAATTGCGGTCAAGTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTGTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGAc          <  2:292306/247‑1 (MQ=255)                                                                                                                                                                                                                                    catt‑cTTTATCTGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGGTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGAt  <  1:293793/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                TGCAATGGCGCCGAGGAGTTTATGGTCGTTTGCCTGCGCCGTGCAGCACAGCATCAGGCTAATCGCCAGGCTGGCGGAAATCGTAAAAACGGATTTCATAAGGATTCTCTTAGTGGGAAGAGGTAGGGGGATGAATACCCACTAGTTTACTGCTGATAAAGAGAAGATTCAGGCACGTAATCTTTTCTTTTTATTACAATTTTTTGATGAATGCCTTGGCTGCGATTCATT‑CTTTATATGAATAAAATTGCTGTCAATTTTACGTCTTGTCCTGCCATATCGCGAAATTTCTGCGCAAAAGCACAAAAAATTTTTGCATCTCCCCCTTGATGACGTGGTTTACGACCCCATTTAGTAGTCAACCGCAGTGAGTGAGTCTGCAAAAAAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGAT  >  NC_000913/11688‑12164

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