Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,481	A→C	E389D (GAA→GAC)	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,628,481	0	A	C	100.0%	74.6 / NA	25	E389D (GAA→GAC)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (13/12);  total (13/12)

ATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCA  >  NC_000913/4628236‑4628722                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  attgacactgATTTTGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTGCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGAtt                                                                                                                                                                                                                                               <  1:98672/241‑1 (MQ=255)             tttGTCACTACCCAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGttt                                                                                                                                                                                                                                  >  2:70569/1‑251 (MQ=255)                 tcactacccAGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCctg                                                                                                                                                                                                                              >  2:91059/9‑251 (MQ=255)                          aGGCGTTCTGCAAAAAGTACGAAGGGCGTGAGCATCCGTTCGTACAGGCGTTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTg                                                                                                                                                                                                                     >  2:88501/1‑251 (MQ=255)                                                       gagcATCCGTTCGTACAGGCGTTGATTGCTGAATACCGTTGCGAGCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGGTGAAGAGGACGATTATGACAGCCGTGTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCa                                                                                                                                                                                        <  2:138130/248‑1 (MQ=255)                                                                          cgttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACACAACCCCGTGGGTGGCGGATGGTTTAAGCAGCCTAGGCCGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATAGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGCTGGGGGAGCAGAGATAACCGCGCTGaaa                                                                                                                                                                      >  1:121419/4‑250 (MQ=255)                                                                          cgttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAAc                                                                                                                                                                     <  1:192368/248‑1 (MQ=255)                                                                           gttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTAAGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGAGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAAcg                                                                                                                                                                    >  1:120180/3‑250 (MQ=255)                                                                            ttgattgatGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGAGGATCGCCAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAAcgg                                                                                                                                                                   >  1:176433/2‑249 (MQ=255)                                                                                     gaATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGc                                                                                                                                                          <  1:70569/251‑1 (MQ=255)                                                                                            gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGtat                                                                                                                                                   <  1:91059/251‑1 (MQ=255)                                                                                            gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGtat                                                                                                                                                   >  1:264393/1‑251 (MQ=255)                                                                                                      ggggATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGGGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATGATGACAGCCGTTTCCTGCGCTGCTTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAAta                                                                                                                                             <  1:252470/244‑2 (MQ=255)                                                                                                                  ggAGAACAACACGCCGTGGGGGGCGGATTGTTTACGGAGCCTCGGCAGGTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCTTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGATATAACCGCGATGAAACGGCTCAAAGGCGAGTGATAAAATAAGTTTTTT‑GATGAGtt                                                                                                                            <  2:249719/251‑1 (MQ=255)                                                                                                                            acGCCGAGGGTGGAGGATGGTTAACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGATGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGAGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCGGAGATAACCGCGATGAAACGGCCCAAAGGCGAGGAATAAAATAAGTTTTTT‑GATGAGTTGAAAACctct                                                                                                                  >  2:191165/1‑251 (MQ=255)                                                                                                                                 gtgggtggCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGCGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTgg                                                                                                              >  2:21223/1‑250 (MQ=255)                                                                                                                                         cGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTg                                                                                                     >  1:101854/1‑251 (MQ=255)                                                                                                                                                                 gTTCGGTTGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGATCTGTTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGATTTGAAAACCTCTCTGGAAGAGGCTGTCGAGAGTAATGAAGGGTTGaaaa                                                                             <  2:159911/251‑1 (MQ=255)                                                                                                                                                                     ggTGGATCGCAAAGAGTTCCAGAACTTGCTTGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAATAAAcc                                                                         <  2:120180/251‑1 (MQ=255)                                                                                                                                                                     ggTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAAcc                                                                         <  1:88501/251‑1 (MQ=255)                                                                                                                                                                                                                    agaACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATAAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCtgttgt                            >  1:74131/1‑249 (MQ=255)                                                                                                                                                                                                                        cAATATCGGATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGCGCTGGTGCGGGAGATGCTGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAACGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGCGGCTGTCGGGATTAAACAAGGTTTAAAAAACCCTGCACGGGTGACCCGCTACGAACTTGCGGATACTAAGGCTGTTGTAGTCCg                      >  2:32893/1‑251 (MQ=255)                                                                                                                                                                                                                               gaattcgggCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGaaaagg               <  1:241303/243‑4 (MQ=255)                                                                                                                                                                                                                                        cGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCAcgc      >  1:30078/1‑251 (MQ=255)                                                                                                                                                                                                                                             tGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCGCGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCAcgccgta  <  1:21223/250‑3 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCA  >  NC_000913/4628236‑4628722

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