Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,749	A→G	A18A (GCA→GCG)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,636,749	0	A	G	97.0%	87.6 / ‑4.2	33	A18A (GCA→GCG)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (13/19);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

TTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAT  >  NC_000913/4636502‑4636977                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      ttACTCAAGTCGCCAGGCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGcc                                                                                                                                                                                                                                    >  1:30696/1‑251 (MQ=255)   attcaaatcgccaggctgcgTCTGGCCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGTCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTg                                                                                                                                                                                                                                  <  1:263395/234‑1 (MQ=255)             ccaggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGGTCGTACTg                                                                                                                                                                                                                        >  1:173635/4‑251 (MQ=255)              caggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACCCCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGGCTGGTTCGGACTgg                                                                                                                                                                                                                       >  2:80363/3‑251 (MQ=255)              caggCCGCGTCTGGTCCCGCCAGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTgg                                                                                                                                                                                                                       <  1:231802/249‑1 (MQ=255)                                   aGCAACTAATGGATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTa                                                                                                                                                                                                 <  1:80363/251‑1 (MQ=255)                                     cAACTGATGTATAGCATATGGGAAGACGCGCAGTACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGGCAAGCTGCGCGCCAGCAGCCCCGATCTTTCCCCGAATAATACTCAACGCTGCATGGGATATATTCTGAGGGGCCTGTAATGGTTATCCGCAGGCGGTTGTTGCTGGGCTAGTTGTTACTGGTGGCGGGGGCAGCCGGGTTCTTACTGGCCATTTTTGTCAAAG‑AAGTTaaa                                                                                                                                                                                               <  1:15879/251‑1 (MQ=255)                                             ggATAGCGTGTGGGAAGATGCGCAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCgg                                                                                                                                                                                       >  1:137866/1‑250 (MQ=255)                                                       tGGGAAGATGCGCAGGACAACTACGATCGCACCGATGATACCCCCATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGAATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTTAAAGTAAGT‑AAACCGGGCGTGCGAGGAGc                                                                                                                                                                              >  2:253636/1‑250 (MQ=255)                                                                 cgcAGGACACCTACGATCGCACCGTTGATACCCACATTAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAggggg                                                                                                                                                                   >  1:250591/1‑250 (MQ=255)                                                                         ccctccgaccgccccgtgtaTACCCACATTAAAACGTGGCGTGCCAAGTTGCGCGCCATCAACCCCGATCTGTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGTGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGTGGACGTTAATc                                                                                                                                                           <  2:250591/232‑1 (MQ=255)                                                                                          tgATACCCACATTAAAACG‑TTCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTg                                                                                                                                         >  1:74948/2‑251 (MQ=255)                                                                                                      tAAAACGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGCTGTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTgg                                                                                                                               <  2:67123/250‑1 (MQ=255)                                                                                                         aaCGTTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGcgcg                                                                                                                           >  1:245171/1‑251 (MQ=255)                                                                                                                cGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTTGTACTGGCTATTTTTGTCAAAG‑AAGTTTAAACGGGCGTGCGAAGGGCAACGGAGGGGGCGGTAAATGATACCGCCACGTTGCTGGGGGGGCTGGggcgggccgct                                                                                                                    >  2:364210/1‑244 (MQ=255)                                                                                                                           ccgcgcCATCAACCCCGATCTTTCACCGTTTAATACTCATCGCGGCATGGGATACAGCCTGAGGGGCCTGTAATGCGTAGCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTGCGTACTGGCTATTGTTGTCAAAG‑AAGTTAAACCGGTCATGCTAAGAGCAACGGATGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGctctct                                                                                                           <  1:216340/248‑1 (MQ=255)                                                                                                                                    aaaacccGATCTGTCACCTAGTAATTCTATTCGCGGCATGGGATATAGCCTGAGGGGCCTTGAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTGTTACTGGTTGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGGCAAAG‑AAGGTAAACCGGGCGGGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAAc                                                                                                <  2:302516/247‑1 (MQ=255)                                                                                                                                                  cACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCTTGTAATGCGTATCGGCATGCGGTTGCTGCTTGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTg                                                                                   <  1:379396/250‑1 (MQ=255)                                                                                                                                                             gttgaTCGCGGTTTGGGATATAGTCTGGGGTGGCTGTAATGCGTATCGGCAGTCGGTTGCTGCCGGGATAGTTTTGACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGGCCAAG‑AAGTTAAACCGGGCGGGCGAAGAGCAACTGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCGGGCGTGTCCCGAGTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGAGCTGGCGttt                                                                       <  2:207666/247‑1 (MQ=255)                                                                                                                                                                   cgcgGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTTCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATc                                                                   >  1:373166/1‑249 (MQ=255)                                                                                                                                                                      ggCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTa                                                              <  1:364210/251‑1 (MQ=255)                                                                                                                                                                       gCATGGGATATATCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGTCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTac                                                             <  1:194777/251‑1 (MQ=255)                                                                                                                                                                       gCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAAATGGCGCAGGCGTTTAATCAGcaac                                                             >  1:109948/1‑248 (MQ=255)                                                                                                                                                                                   gCCTTAGGTGCCTGTTATGCGTGTCGGCATGCGGTTGTTGCTGGGCTATTTTTGACTGGTGGCGGTGGCGGCCTGGGTTGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGATCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGTGCTGGCGCGTCCCGATTTGCTCTCTGGGTACCCACCGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGt                                                 <  2:9772/251‑1 (MQ=255)                                                                                                                                                                                         ggggCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCTGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcg                                           <  1:209500/251‑1 (MQ=255)                                                                                                                                                                                         ggggCCTGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCCACGCAAGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcc                                           >  1:16459/1‑250 (MQ=255)                                                                                                                                                                                         ggggCCTGTAAGGCGTACCGGCATGCGGTTGCTGCTTGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCGGGCGGAGCTGGCGCGGCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTcgcg                                           <  1:52470/251‑1 (MQ=255)                                                                                                                                                                                               tGTAATGCGTATCGGCATGCGGTTGCTGCTGGGCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACt                                                                                    <  1:223895/204‑1 (MQ=255)                                                                                                                                                                                                  aaTGCGTATCGGCATGCGGTTGCTGCTGGGCTAGTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCg                                  <  1:323345/251‑1 (MQ=255)                                                                                                                                                                                                                 catgttgctGCTTGGCTATTATTTACGGGAGGCGGTGGCGGCCTGGTTCGTACTGGCTATTGTTGTCAAAG‑AAGTTAAACCGGTCGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCTCAACGTTGCTGGCGGAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTGCGCGCCAATATCGGTGGCATTAACaaa                    <  2:282371/247‑1 (MQ=255)                                                                                                                                                                                                                   ggTTGCTGCTGGGCTCTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGGACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATAGGTGGCATTAACaaa                    >  1:282371/1‑248 (MQ=255)                                                                                                                                                                                                                               gCTATTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAatat     <  1:118215/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTTTTTTACTGGTGGCGGTGGCGGCCTGGTTCGTACTGGCTATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAt  <  1:114126/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      TTACTCAAGTCACCGGGCCGCGTCTGGTCCCGCCAGCAACTGATGGATAGCGTATGGGAAGATGCGCAGGACACCTACGATCGCACCGTCGATACCCACATTAAAACGCTGCGTGCCAAGCTGCGCGCCATCAACCCCGATCTTTCACCGATTAATACTCATCGCGGCATGGGATATAGCCTGAGGGGCCTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAG‑AAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAT  >  NC_000913/4636502‑4636977

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