Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,690,469	A→G	R210R (CGT→CGC)	glrK ←	sensor protein kinase regulating glmY sRNA in two‑component system with response regulator GlrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,690,469	0	A	G	97.7%	139.4 / ‑6.3	43	R210R (CGT→CGC)	glrK	sensor protein kinase regulating glmY sRNA in two‑component system with response regulator GlrR
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (16/26);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CATGAGATAAATGTCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGACGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGCCCACGTTCGGCGATTTCA  >  NC_000913/2690321‑2690601                                                                                                                                                     |                                                                                                                                     cATGAGATAAATGTCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCg                                                                                                                                    >  2:621745/1‑151 (MQ=255)    gagaTAAATGTCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGtt                                                                                                                                 >  2:94761/1‑151 (MQ=255)     agaTAAATGTCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTg                                                                                                                                <  1:620822/151‑1 (MQ=255)       aTAAATGTCTTAAAAAGTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGAt                                                                                                                              <  1:567454/151‑1 (MQ=255)         aaaTGTCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATc                                                                                                                             >  1:220210/1‑150 (MQ=255)           aTGTCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAAGACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATg                                                                                                                          <  1:408517/151‑1 (MQ=255)              tCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCAt                                                                                                                           >  1:542942/1‑147 (MQ=255)               cTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCt                                                                                                                      >  1:433353/1‑151 (MQ=255)                  aaaaaTTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCg                                                                                                                    >  1:28409/1‑150 (MQ=255)                   aaaaTTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGat                                                                                                                  <  2:470146/151‑1 (MQ=255)                          aTGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTAGATAttttt                                                                                                            >  1:688217/1‑150 (MQ=255)                                ttGGGATTCCAGCCATGACAGGCGCGCACTTAACCAAAGAATACGTTGCCCAACCGAGCGAAACTCGCTCGGTCCACTGAAAGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGACGGTTGATCCTGCGCTCGAAATTTTTCACCGGc                                                                                                     >  1:29502/1‑151 (MQ=255)                                    gATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGa                                                                                                 <  2:224480/151‑1 (MQ=255)                                        ccAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATaa                                                                                              <  1:286129/150‑1 (MQ=255)                                            ccATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAatcatc                                                                                         <  2:28167/151‑1 (MQ=255)                                                        cTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGc                                                                             >  2:436159/1‑151 (MQ=255)                                                          cACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAg                                                                           <  2:49773/151‑1 (MQ=255)                                                              tAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTa                                                                         <  1:367876/149‑1 (MQ=255)                                                              tAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTAc                                                                        <  2:152267/150‑1 (MQ=255)                                                               aaCCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTAcc                                                                       >  2:15680/1‑150 (MQ=255)                                                                 ccAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCAtt                                                                    <  1:78210/151‑1 (MQ=255)                                                                           cGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCa                                                          <  2:282704/151‑1 (MQ=255)                                                                               gCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCa                                                       <  1:600430/150‑1 (MQ=255)                                                                                   aaCCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAAt                                                  <  2:685076/151‑1 (MQ=255)                                                                                    aCCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAAt                                                  >  1:313281/1‑150 (MQ=255)                                                                                    aCCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAAt                                                  >  2:134259/1‑150 (MQ=255)                                                                                             aaCTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGc                                         <  1:69571/150‑1 (MQ=255)                                                                                                 cgctcgGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGAtt                                     <  2:513438/150‑1 (MQ=255)                                                                                                 cgctcgGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTg                                    >  2:661841/1‑151 (MQ=255)                                                                                                      ggTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATCTTTTTCAACGGCCCGATACTCATCCGCGTGAAAAGCTGTACCATTACCAGACTCACCAGAAATAACACCAGCGATTGccaa                                >  1:486568/1‑150 (MQ=255)                                                                                                       gTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGccaacc                              <  2:670723/151‑1 (MQ=255)                                                                                                       gTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTAACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGccaacc                              >  2:654352/1‑151 (MQ=255)                                                                                                            cTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGccaacca                             <  1:4083/147‑1 (MQ=255)                                                                                                              gAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAata                        <  2:114028/150‑1 (MQ=255)                                                                                                               aaCGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAatat                       <  1:581971/150‑1 (MQ=255)                                                                                                                aCGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTg                     <  1:248062/151‑1 (MQ=255)                                                                                                                 cGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGc                    >  1:20232/1‑151 (MQ=255)                                                                                                                   agaCGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGccc                  >  2:542688/1‑151 (MQ=255)                                                                                                                     aCGCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGCCCAc                <  2:49972/151‑1 (MQ=255)                                                                                                                       gCTATTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGCCCACGt              <  2:507578/151‑1 (MQ=255)                                                                                                                          aTTGCCCAGAGAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGCCCACGTTc            <  1:484069/150‑1 (MQ=255)                                                                                                                                 agagAACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGCCCACGTTCGGCGAttt    <  2:445830/151‑1 (MQ=255)                                                                                                                                    gaACGCCCTTCCCCCAGGCGGTTGATCATGCGCTCGATATTTGTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTGCCAGACTCTCCAGAAATAGCTCCAGCGATTGCCAACCAATATATTGCCCACGTTCGGCGATTTCa  <  1:669123/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                     CATGAGATAAATGTCTTAAAAATTGATGGCGTTGGGATTCCAGCCATGACAGGCGCTCACTTAACCAAAGAATACGTTGCCCAACCGAGCGTAACTCGCTCGGTCCACTGAACGAGACGCTATTGCCCAGAGAACGCCCTTCCCCCAGACGGTTGATCATGCGCTCGATATTTTTCACCGGCCCGATAATCATCCGCGTGAAAAGCAGTACCATTACCAGACTCACCAGAAATAGCACCAGCGATTGCCAACCAAAATATTGCCCACGTTCGGCGATTTCA  >  NC_000913/2690321‑2690601

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