Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,818	A→G	K168K (AAA→AAG)	nadR →	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,818	0	A	G	100.0%	93.6 / NA	30	K168K (AAA→AAG)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/13);  total (17/13)

GCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAAT  >  NC_000913/4627593‑4628047                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      acTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTaca                                                                                                                                                                                                              <  1:176881/250‑1 (MQ=255)       tatCATTATGGGTTTTGAAGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATAAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAAAGGCATCAAAAAGTTTATAGCTGAAAAGGGGGTTCAGCCGGAACTGATCTACACCTCgg                                                                                                                                                                                                        >  2:170786/1‑251 (MQ=255)                    tttGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATgc                                                                                                                                                                                           >  1:67170/1‑251 (MQ=255)                       gACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGcc                                                                                                                                                                                        <  2:212187/251‑1 (MQ=255)                                      cgcgcgTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGCTCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGCGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGGTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAAc                                                                                                                                                                          >  1:41643/1‑250 (MQ=255)                                        cgcgTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATc                                                                                                                                                                       >  1:102343/1‑251 (MQ=255)                                                               tGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGAGCCACAGTATATGGAACATCTGGGGATCGATACGGAGCTGGtg                                                                                                                                                >  2:228165/1‑250 (MQ=255)                                                                                 ggCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCttt                                                                                                                              <  2:197469/250‑1 (MQ=255)                                                                                  gCCGGATCGTATGCGCTGGTTATTACAAACTTGTAAATATCAGAAAAATATGCGCATGCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTa                                                                                                                             <  2:62854/251‑1 (MQ=255)                                                                                     ggATCGTCTGCGCTAGTTATTACAATCTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAATGAAGAGGGCATAGAGCCTTATCCACACGGCTGGGATGTGTGGAGCAACGGCATCAAAAATTTTATGTCTGAAAAGGGGATTCAGCCGGCCCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATTGAACATCTGGGTATCGATACGGTGCTGGTCGATCCGAAACGTACCATTATGa                                                                                                                          <  2:54259/251‑1 (MQ=255)                                                                                      gATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATg                                                                                                                           >  2:199446/1‑249 (MQ=255)                                                                                                 ctGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCt                                                                                                                                                    <  2:199479/212‑1 (MQ=255)                                                                                                 ctGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCt                                                                                                                                                    >  1:199479/2‑213 (MQ=255)                                                                                                          acAAACTTATAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGACGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATc                                                                                                      <  2:74942/249‑1 (MQ=255)                                                                                                              aCTTTTACGTATAAGAAAAATATGCGCATTCATGCTTTCAACGAAGAAGACATAGACCCGTATCCGCACGGCTGGGATGTGTGGAGCACCGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTAAGCCGGACCTGATCTACACCTCGGAAGAAGCCGAAGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGAACCTTTATGAGTATCAGCGGTTCGCAGCTCcg                                                                                                    >  2:131525/1‑248 (MQ=255)                                                                                                               cTTTTAAATATCAGAAACATATTCGCATTCATGCTTTCATCGAAGAGGACATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAAAGGCATCAAAAAGTTTATGGATGAAACGGGGATTCAGCCGGACCTGCTCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTATATGAGTATCCTCGGTGCGCATATCCGCGa                                                                                                 >  2:236532/1‑250 (MQ=255)                                                                                                                                aaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTg                                                                               <  1:79857/251‑1 (MQ=255)                                                                                                                                aaTATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTg                                                                               >  1:77642/1‑251 (MQ=255)                                                                                                                                 atatTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACCGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTgg                                                                              >  1:166949/1‑251 (MQ=255)                                                                                                                                       gCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAatat                                                                        >  2:69449/1‑251 (MQ=255)                                                                                                                                             aTGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCt                                                                  >  1:90035/1‑251 (MQ=255)                                                                                                                                             aTGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGAGGTGTGGAGCATCGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCGCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGTGATCGATACGGTGCTGGTCGATCCGAAACGTTCCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCt                                                                  <  1:117632/251‑1 (MQ=255)                                                                                                                                                             gggCATGGAGCCGTATCCGCACTGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGt                                                  <  1:152736/251‑1 (MQ=255)                                                                                                                                                                cATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGttct                                               >  2:148872/1‑249 (MQ=255)                                                                                                                                                                                cACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGc                               <  2:90035/251‑1 (MQ=255)                                                                                                                                                                                            gtgtgGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAAACCGTTCCGCTACTGGGAATATATTCCTACGGCAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTggtgg                   >  1:124563/1‑248 (MQ=255)                                                                                                                                                                                                   gCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCg            >  1:73447/1‑251 (MQ=255)                                                                                                                                                                                                     aaCGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAg          <  2:71434/251‑1 (MQ=255)                                                                                                                                                                                                             cAAAAAGTTTATGGCTGAAAAGGGGCTTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGAGTCGAGCGGTAAAt  >  2:189118/1‑251 (MQ=255)                                                                                                                                                                                                             cAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACCCCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGATACGGTGCTTGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACGGAAGTGAAGCCGTTCTTTGTACGTACCGTGGCGATCCTTGGTGGCGATTCGAGCTGTAAAt  <  1:51971/251‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      GCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCATGGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATATGGAACATCTGGGGATCGAGACGGTGCTGGTCGATCCGAAACGTACCTTTATGAGTATCAGCGGTGCGCAGATCCGCGAAAACCCGTTCCGCTACTGGGAATATATTCCTACCGAAGTGAAGCCGTTTTTTGTGCGTACCGTGGCGATCCTTGGCGGCGAGTCGAGCGGTAAAT  >  NC_000913/4627593‑4628047

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