Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,636,925	G→C	R77P (CGC→CCC)	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,925	0	G	C	97.0%	97.8 / NA	33	R77P (CGC→CCC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (17/15);  minor base A (1/0);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTC  >  NC_000913/4636691‑4637168                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     cTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCcccc                                                                                                                                                                                                                                                   <  1:282213/238‑1 (MQ=255) cTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGCGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGCGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAat                                                                                                                                                                                                                                      >  1:328546/1‑251 (MQ=255)   gTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGAACCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGGTTCGCGCCAatat                                                                                                                                                                                                                                    >  1:156527/1‑251 (MQ=255)     aaTGCGTATCGGCATGCGGTTGTTGCTGGGATATTTTTTACCGGTGTCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGGCAAAGAAGTTAAACCGGGCGTGCGAATAGCAACGGAGGGGACGTTGATCGACACCGCCACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAAGATCg                                                                                                                                                                                                                                  <  1:125320/251‑1 (MQ=255)        gCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTg                                                                                                                                                                                                                               <  1:77048/251‑1 (MQ=255)           tATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGGGGCGGTGGAAGCATGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGATACCGCAACGTTGCTGGAGGAGCTGGCGCGGCCCGATTTTCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTCAAAAATCCCCCGTTTCGCGCCAATATCGGTGGcc                                                                                                                                                                                                                            >  1:203861/1‑250 (MQ=255)            aTCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCa                                                                                                                                                                                                                            <  1:261346/250‑1 (MQ=255)              cGGCATGCGGTTGTTGCTGGGGTATTTTTTACTTGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCCACGGAGTGGACGTTGATCGACACCGCAACGGTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCTCCAATATCGGTGGCagt                                                                                                                                                                                                                          <  2:21409/250‑3 (MQ=255)              cGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAAACGGGCGTGCGAAGAGCAACGGAGGGGAAGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGGCCCGCTTTGTTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAAAATCACCCGTTTCGCGCCAAAATCGGCGGCAtt                                                                                                                                                                                                                          >  2:72026/1‑250 (MQ=255)                 cATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCCGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGGGGCATTAAc                                                                                                                                                                                                                       >  1:228447/1‑250 (MQ=255)                  aTGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACa                                                                                                                                                                                                                      >  1:313988/1‑250 (MQ=255)                                            actgggtgcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGACGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCtata                                                                                                                                                                                           <  2:179043/244‑1 (MQ=255)                                                   gcggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTTTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGa                                                                                                                                                                                    <  2:242140/251‑1 (MQ=255)                                                    cggtggcAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGAt                                                                                                                                                                                   <  1:60562/251‑1 (MQ=255)                                                                     gTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGCGCGCAATGAATATCATGTCTATATGACCGGTGCGCAGGGCAAAGTAtt                                                                                                                                                                  >  2:309615/1‑251 (MQ=255)                                                                      tACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGAAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGAATTg                                                                                                                                                                 >  1:157338/1‑251 (MQ=255)                                                                          ggCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCTATGCGCAGGGCAAAGTATTGtt                                                                                                                                                               <  1:123116/249‑1 (MQ=255)                                                                              agtgttggcAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGAGGTTGACCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGACCCGATTTGCTCTCTGGGGACCCAACGGATGGGCAACTGGCGCAGGGGTTTAATCAGCTACAGCATCCCCCGTTTCGCTCCAATATCGGTTGCATTAACAAAGTGCGCAATGAAGATCATGTCTATATGACCGATGCGCTGGGCAAAGTATTGTTCGATTc                                                                                                                                                         <  2:312748/243‑1 (MQ=255)                                                                                                aaaCCGAGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCACGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAGGTCTATATGACCGCTGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGGTgg                                                                                                                                       >  2:308477/1‑251 (MQ=255)                                                                                                             gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGCTTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGCTTATTcg                                                                                                                          >  2:116011/1‑251 (MQ=255)                                                                                                             gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGtt                                                                                                                                                               >  2:48853/1‑214 (MQ=255)                                                                                                             gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGtt                                                                                                                                                               <  1:48853/214‑1 (MQ=255)                                                                                                             gAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTcg                                                                                                                          >  1:253036/1‑251 (MQ=255)                                                                                                                          ggggACGCTGATCGCCACCGCAACGTTGATGGCGGCGCTGCCGCGACCCGATTTGCACTCTGAGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTCCTCAGCTACAACCTCCCCCGTTTCGCCACAATCTCGGTGGAATTAACAAAGTACGCCATGAATATCATGTCTATATCACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAACGCCCTTGGACAGGATTCTTCTCGCTGGAATGACg                                                                                                             >  2:315874/1‑251 (MQ=255)                                                                                                                                                 cGTTGCGGGGGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGAGTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAg                                                                                     <  2:313988/251‑1 (MQ=255)                                                                                                                                                   ttGCTGGCGGCGCTGGCGCGTTTAGCTTTGCTCTATGGGGACCCAACGCATGGGCAACAAGCGCAGGCGTTTAATCAGCTACAACATCCCCCTGTTCGCGCCCATCTCGGTGGCATTAAAAAAGTGCGCAATGAATATCATGTCTATATACCCGATTCCCAGGGCAACGGATTGTTCGATTCGGCAAATACAACAGTTGGACAGGATTATTCGCGCGTGTATGACGTTCTG‑CTAACGTTACGCGGTctgt                                                                                    >  2:219619/1‑247 (MQ=255)                                                                                                                                                     gCTGGCGGAGCGGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAGTAt                                                                                  <  1:309615/250‑1 (MQ=255)                                                                                                                                                              gctggaccgTCCCGATTTTCTCTCTGGGGACCCAACGGATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATACCCCGTTTCGCGACAAGATCGGTGGCAATAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGCCGT‑CTGGCTAACGGTTCGTGGGCAGtttagtgctcgc                                                                         >  2:317874/8‑239 (MQ=255)                                                                                                                                                                                                 cGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTAAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCt                                     >  2:258319/1‑251 (MQ=255)                                                                                                                                                                                                   cATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCtgt                                   >  2:126624/1‑251 (MQ=255)                                                                                                                                                                                                      gggCAACTGGCGCAGGCGTTTAATCAGCTACAACATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGACAGTTCTGGGat                                >  2:244993/1‑251 (MQ=255)                                                                                                                                                                                                                           aaTCAGCTACAACTTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATTAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGGATGTTGCCGCACCGATTa            <  2:228447/250‑1 (MQ=255)                                                                                                                                                                                                                               aGCTACATCATCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATAGCATGTCTATATGACCGATGCGCAGGGCAAAGTATGGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CGGGCTAACGTTGCGTGGTCAGTAGGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCTAAAGTTCTGTGATGTATGTTGCCGCTCCGATTATGGAc       <  2:157338/251‑1 (MQ=255)                                                                                                                                                                                                                                    caacaTCCCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTTCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTc  <  1:27181/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     CTGTAATGCGTATCGGCATGCGGTTGTTGCTGGGCTATTTTTTACTGGTGGCGGTGGCAGCCTGGTTCGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGT‑CTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTC  >  NC_000913/4636691‑4637168

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