Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,637,001	C→T	G102G (GGC→GGT)	creC →	sensory histidine kinase CreC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,637,001	0	C	T	96.8%	83.1 / NA	31	G102G (GGC→GGT)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/12);  total (19/12)
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.81e-01

CGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGC  >  NC_000913/4636758‑4637235                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           cGTACTGGCTATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCAGGTCTATATGACCGATTCGCAGGGTAAAGTAt                                                                                                                                                                                                                                     <  2:44931/251‑1 (MQ=255)                      agaagTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAg                                                                                                                                                                                                                <  1:24358/251‑1 (MQ=255)                          gTTAAACCGGGAGTGCGAAGAGCAACGGAGGGGACGTTAAGCGATACCGCAACGTTGCTGGCGGAGCTGGCGTGTCCCGATTTGCTCTCTGGGTACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGt                                                                                                                                                                                                            <  2:77076/251‑1 (MQ=255)                                    gCGTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCCTTAACAAAGGGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGat                                                                                                                                                                                                  >  1:149215/1‑251 (MQ=255)                                      gTGCGAAGAGCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGatt                                                                                                                                                                                                 >  1:207078/1‑250 (MQ=255)                                               gCAACGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCCTGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGCTTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTg                                                                                                                                                                                       >  2:28596/1‑251 (MQ=255)                                                 aaCGGAGGGGACGTTAATCGATACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTgg                                                                                                                                                                                      >  1:73579/1‑250 (MQ=255)                                                              ttaaTCGATACCGCAACGATGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTAGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGAACGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCt                                                                                                                                                                        >  1:52499/4‑251 (MQ=255)                                                                    gatacCGCAACGTTGCATGCGGAGCTGGCGCGTCCCGATGTGCTCTCTGTGGATCCAACGCATGGGCAACTTGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACTAAGTGCGCAATGAATAGCATGTCTATATGACCGAGTCGCAGGGTAAATTTTTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGtt                                                                                                                                                                  <  1:236748/248‑1 (MQ=255)                                                                      tacCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTc                                                                                                                                                                             <  2:127685/237‑1 (MQ=255)                                                                      tacCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTc                                                                                                                                                                             >  1:127685/2‑238 (MQ=255)                                                                                        ggAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTgcgc                                                                                                                                              >  2:129016/1‑251 (MQ=255)                                                                                              ggCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCgcag                                                                                                                                          <  1:210142/249‑1 (MQ=255)                                                                                                cgcgTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACg                                                                                                                                      >  2:71069/1‑251 (MQ=255)                                                                                                     cccGATTTGCTCTCTGGGGACCCAAAGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGATTCGCGCCACTATCGGTGGCATTAAAAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGATGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGc                                                                                                                                  >  1:130703/1‑250 (MQ=255)                                                                                                                      ggACCCAACGCATGGGCACCTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAAGGAATATCATGTCTATATAACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCTGTTGGGCAGGATTAGTCGCGCTGGAATGACGTCTGGCTATCGTTGCGTGGTCAGTATGGTGCGCTCAGCACGTTTCAAAATCCTACCAATCCCg                                                                                                                <  2:114376/251‑1 (MQ=255)                                                                                                                               gCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTcgg                                                                                                       >  2:49835/1‑249 (MQ=255)                                                                                                                                cATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATGGTTCGATGCGGCAAATAAAGCCGTTGGGCAGGATGAGTCGCGCTGGAATGGCGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTcggt                                                                                                      <  1:115606/251‑4 (MQ=255)                                                                                                                                       aaCTGGCGCAGGAGTTTACTGCGCTACAACATCGCCCGTTTAGAGCCAAGCGCGGTGGCATTAACACCGTGCGCACTGAATATCCTGTCTATATGCCCGATTCGCAGGGTAATGCATTGTTCGATTCGGCACATAAAGCCGTTGGGCAGGATTATTGGCGCTGGAATGAGGTATGGATAACGTTGCGTGGTCAGTACGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCACGAAAGTTcgacgatggat                                                                                               >  2:195733/1‑241 (MQ=255)                                                                                                                                              gcAGGCGTTTTATCAGCTACAACATCGCCCGTTTAGCGCCAATATCGGTGGCATTAACAGAGTACGCAATGAATATCATGTCTATATGACCGATTAGCAGGGCAAAGTCTTGTTCGATTAGGCAAATAACGACGATGGGCAGGATTATTCGCGCTGGTATGACGTCTGGCTACCGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGATCGGTGAGATATGTTGCCg                                                                                        >  2:145956/1‑251 (MQ=255)                                                                                                                                                  gCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATAGCGGTGGCATTAACAAAGTGCGCAATGAATATCAGGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGTGCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGcgcc                                                                                    <  1:29484/251‑4 (MQ=255)                                                                                                                                                         aTCAGCTGCAACATCGCCCGTTTCGCGCAATTATCGGTGGCATTAACAAAGAGCGCAATGTATATCATGTATATCTGACCGATTCGCAGGGTAACGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGAGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTACCGCGACGATTATg                                                                             >  2:25510/1‑251 (MQ=255)                                                                                                                                                              cTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGACTATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAAAGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACgg                                                                        >  2:214686/1‑251 (MQ=255)                                                                                                                                                                              ttCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTgg                                                         <  1:71069/250‑1 (MQ=255)                                                                                                                                                                                           cGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGt                                             >  2:254356/1‑249 (MQ=255)                                                                                                                                                                                                       cAAAGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAAcg                               <  1:49835/251‑1 (MQ=255)                                                                                                                                                                                                        aaaGTGCGCAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGAGATGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGCTGTATGTTGCCGCGACGATTATGGACGTCTCGCGGCTTCTTGGCGTTTTGAGCGTAGGCAAACAGAAcgc                              >  2:140776/1‑251 (MQ=255)                                                                                                                                                                                                                cAATGAATATCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTAGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGAGCAGCACGTTGAAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATgg                      >  1:120553/1‑251 (MQ=255)                                                                                                                                                                                                                      atatCATGTCTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCgg                >  1:98519/1‑251 (MQ=255)                                                                                                                                                                                                                               cTATATGACCGATTCGCAGGGTAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGGTCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGc       >  1:2193/1‑251 (MQ=255)                                                                                                                                                                                                                                     gACCGATTCGCAGGGTAAAGTATGGTGCGATTCGGCAAATAAAGCCGTTGGGCAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTGCGGTGATGTATGTTGCCGCGCCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTagc  <  2:207078/250‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                           CGTACTGGCCATTTTTGTCAAAGAAGTTAAACCGGGCGTGCGAAGAGCAACGGAGGGGACGTTGATCGACACCGCAACGTTGCTGGCGGAGCTGGCGCGTCCCGATTTGCTCTCTGGGGACCCAACGCATGGGCAACTGGCGCAGGCGTTTAATCAGCTACAACATCGCCCGTTTCGCGCCAATATCGGTGGCATTAACAAAGTGCGCAATGAATATCATGTCTATATGACCGATGCGCAGGGCAAAGTATTGTTCGATTCGGCAAATAAAGCCGTTGGACAGGATTATTCGCGCTGGAATGACGTCTGGCTAACGTTGCGTGGTCAGTATGGTGCGCGCAGCACGTTGCAAAATCCTGCCGATCCCGAAAGTTCTGTGATGTATGTTGCCGCACCGATTATGGACGGCTCGCGGCTTATTGGCGTTTTGAGCGTAGGCAAACCGAACGCGGCGATGGCTCCGGTCATTAAGCGTAGC  >  NC_000913/4636758‑4637235

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