Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,637,857	T→C	C388R (TGC→CGC)	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,857	0	T	C	93.8%	74.1 / ‑4.3	32	C388R (TGC→CGC)	creC	sensory histidine kinase CreC
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (10/20);  total (10/22)
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.92e-01

GCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGC  >  NC_000913/4637611‑4638089                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         gCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCAt                                                                                                                                                                                                                                      <  1:209287/251‑1 (MQ=255)       gTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCt                                                                                                                                                                                                                                <  1:197893/251‑1 (MQ=255)                      cGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGa                                                                                                                                                                                                                 <  1:41988/251‑1 (MQ=255)                        ccAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTgggtc                                                                                                                                                                                                               >  1:104248/1‑248 (MQ=255)                                                           tGACTGTTGTTGATGTGGCGGCATTATTTCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGt                                                                                                                                                                            >  1:202727/1‑251 (MQ=255)                                                                              ggcATTATTTCGCCGCGGCAGCGAATCGATCACCGTGCAGTTGGCAGAAAAAAATATCCCTTTACATGTTATGCCCTCTGAGGTTAATTTTGATGCTGAACCGGCTTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCTCCCTTAAGGTGCTGTAGACCGGTAGTGGTATTCCTGACTAcg                                                                                                                                                         <  2:63590/251‑1 (MQ=255)                                                                                    atttcgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAac                                                                                                                                                                                                         <  1:365738/193‑1 (MQ=255)                                                                                    atttcgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCAGAAGTGGATCAGGAac                                                                                                                                                                                                         >  2:365738/5‑197 (MQ=255)                                                                                        cgccgcGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACg                                                                                                                                               <  1:248712/251‑1 (MQ=255)                                                                                          ctgcgTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCTACCGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCTCTGGGGAATTTACTGGATAACGCCAGCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTTACTACGCGCTGTCACg                                                                                                                                               <  2:37153/247‑1 (MQ=255)                                                                                                  gCGAAGCGCGCACCGTGCAGTTTGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGGTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGATGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGaa                                                                                                                                     <  2:104248/251‑1 (MQ=255)                                                                                                    gAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATGCCGGTAGAGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACg                                                                                                                                   >  1:263088/1‑251 (MQ=255)                                                                                                    gAAGCGCGCACCGTGCAGTTGGCAGAAAAAAACATCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACg                                                                                                                                   >  1:371013/1‑251 (MQ=255)                                                                                                        cgcgcACCGTGCAGTTGGCAGAAAAAAAAATCAATCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCttt                                                                                                                               <  1:59037/251‑1 (MQ=255)                                                                                                         gcgcTCCGGGCAGTTGGCTGAAAAAAGCATCACGCGACATGTTATGCCCGCTGAGGTTAATGTGGCTGCTGACCCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCtttt                                                                                                                              <  2:114885/251‑1 (MQ=255)                                                                                                                             aaaaaaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCTGTCGCATAACTCTAAGCGCCGAAGTGGATCATGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAAt                                                                                                          <  1:240021/251‑1 (MQ=255)                                                                                                                                aaaaaaTCACTCTACATGTTATGCCCACTGAGGTTAATTTTGCTGCTGAACCGTCGTTACTTGAGCAGGCTCTGTGGAATTTACTTGATAACGCCATCGATTTTACCCCCGAGAGCGTTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAATGTGCTGGATACCGGTAGTGGTATTCCTGACTACTCGCTGTCACATATTTTTGAACGCTTTTACTCTTTGCCGCATGCAAATggg                                                                                                       <  2:173717/251‑1 (MQ=255)                                                                                                                                 aaacaTCACTCTACATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCACATGGGc                                                                                                      >  1:143313/5‑251 (MQ=255)                                                                                                                                          tctacATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAagc                                                                                              >  1:4952/5‑250 (MQ=255)                                                                                                                                             acATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCg                                                                                          >  1:151729/2‑251 (MQ=255)                                                                                                                                             acATGTTATGCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCg                                                                                          <  1:159609/250‑1 (MQ=255)                                                                                                                                              cATGTTATTCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGGTACTGGAGCAGGCGCTGGGGAATTTACTGGAGAACGCCATCGATTTTACCCCCGAGAGCGGTCGCAGAACGCTAAGCGCCGAAGTGTATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTGTGAACGCTTTTACGCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCg                                                                                          <  2:236118/250‑1 (MQ=255)                                                                                                                                                    atgCCCACTGAGGTTAATGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTgg                                                                                    <  1:93393/248‑1 (MQ=255)                                                                                                                                                                         gctgctGACCCGGCGTTATTGTCGCAGGCGCTGGGGAATTTACTGGCTAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCGAAGCGCCGAAGTGGATCCGGAACACGTCGCCCTTAAGGTGTTGGATTCCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAATCTGCCGTCTGGGGTTGGCGTTCGGCAGTGAGGt                                                              <  2:276643/251‑1 (MQ=255)                                                                                                                                                                           gtcggAACCGTCGTTACTGGAGCAGGCGCTTGGGAATTTTCTGTATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATAACGCTAAGTGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCTCTGTCACGTATTTTTGAATGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTc                                                             <  2:186163/246‑1 (MQ=255)                                                                                                                                                                            gctgAACCGGCGTTACTGGAGCAGGCGCTGGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGc                                                           <  1:195276/251‑1 (MQ=255)                                                                                                                                                                             cagAACCGGCGTTACTGGAGCAGGCGCGGGGGAATTTACTGTAGAACTCCTTCCATTTTCCCCCCTAGAGCGGGCGCAGAACGCTAAGCGACGAAGTGGACGAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCGGTCACGTATTTTTGAACGCTTGTACTCTTTGCCGCGTGCATGTGGTCAAAAAAGCATCGTTCTGGGGGTGGCGTTCGTCAGTGTGGTCgac                                                          <  2:249851/249‑3 (MQ=255)                                                                                                                                                                               gAACCGGCGTTACGGGAGCGGGCGCTTGGGAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCGGTCACGGATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATTGGCAAAAAAGCAGCTGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCg                                                        <  2:179516/251‑1 (MQ=255)                                                                                                                                                                                                            gaaTTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTgcgc                            >  1:346976/2‑250 (MQ=255)                                                                                                                                                                                                                 tACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTg                      >  1:10534/1‑251 (MQ=255)                                                                                                                                                                                                                           cGCCATCGATTTTACCCCCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGTCGCCCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTATTTTTGAACGCTTTTACTCTTTGCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTg            <  1:301558/251‑1 (MQ=255)                                                                                                                                                                                                                                      atACCACCGAGAGCGGTCGCATAACGCTAAGCGCCGAAGTGGATCAGGAACACGACGCCCTTAAGGTGCGGGATACCGGTAGTGGTATTCCTGACTACGCGCTGTCACGTAATTTTGAACGCTTGTACTCTTTTCCGCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGTGTTCGTCAGTGAGGTCGGCCGTTTGTTGAACGGCGAAGTCAAGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGc  <  2:188819/249‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         GCATTGGTAGAAACGTTACTACGCCAGGCAAGACTGGAGAATCGTCAGGAAGTCGTTCTGACTGCTGTTGATGTGGCGGCATTATTCCGCCGCGTCAGCGAAGCGCGCACCGTGCAGTTGGCAGAAAAAAAAATCACTCTGCATGTTACGCCCACCGAGGTTAACGTTGCTGCTGAACCGGCGTTACTGGAGCAGGCGCTGGGAAATTTACTGGATAACGCCATCGATTTTACCCCCGAGAGCGGTTGCATCACGCTAAGCGCCGAAGTGGATCAGGAACACGTCACGCTTAAGGTGCTGGATACCGGTAGTGGTATTCCTGACTACGCGCTTTCACGTATTTTTGAACGCTTTTACTCTTTGCCTCGTGCAAATGGGCAAAAAAGCAGCGGTCTGGGGTTGGCGTTCGTCAGTGAGGTCGCCCGTTTGTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGC  >  NC_000913/4637611‑4638089

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