Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,191	C→T	Q293* (CAA→TAA)	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,191	0	C	T	100.0%	104.0 / NA	32	Q293* (CAA→TAA)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/17);  total (15/17)

AAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGA  >  NC_000913/4628062‑4628314                                                                                                                                  |                                                                                                                            aaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAAt                                                                                                                          >  1:537395/1‑133 (MQ=255)            tatCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGAtttt                                                                                                               <  2:70048/133‑1 (MQ=255)             atCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGATGGGCCACGCTTAATACATTGATTTTg                                                                                                              <  2:357540/133‑1 (MQ=255)                ttCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCgg                                                                                                           <  1:474024/133‑1 (MQ=255)                   aaCACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGa                                                                                                        >  1:338957/1‑133 (MQ=255)                     caccacCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGaaa                                                                                                      <  2:584266/133‑1 (MQ=255)                       ccaccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAata                                                                                                    <  1:600355/133‑1 (MQ=255)                        caccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAatat                                                                                                   <  1:8546/133‑1 (MQ=255)                          ccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGc                                                                                                 >  2:380232/1‑133 (MQ=255)                               gCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATa                                                                                            >  2:55261/1‑133 (MQ=255)                                cGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGGTGGGGTCGCATTGCAGTATTTTGACCACGATAAAAACGCGCTGGGCCACGGTTAATACATTTGTTTTTCGGTGAAATATGCCAATaa                                                                                           >  2:157161/1‑133 (MQ=255)                                  tGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAg                                                                                         >  2:20411/1‑133 (MQ=255)                                  tGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAg                                                                                         >  1:100847/1‑133 (MQ=255)                                  tGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGATGGAATATGGCAATAAAg                                                                                         >  1:171146/1‑133 (MQ=255)                                   gggAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTATGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGACATATGCCAATAAAGt                                                                                        >  1:226457/1‑133 (MQ=255)                                                           tttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGt                                                                <  1:159635/133‑1 (MQ=255)                                                            ttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTc                                                               >  1:516353/1‑133 (MQ=255)                                                             tCACACCTCGGCGGTGATGAGATCTCATTGCATTATTATGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTca                                                              <  2:226457/133‑1 (MQ=255)                                                                   cTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGGGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTc                                                        >  1:623908/1‑133 (MQ=255)                                                                             atgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTTCGGTGAAATATGCCAATAAAGGCGCATTTATCGATGCCGACTTTGTCATCACTCCGGCGTTcgg                                              >  2:296816/1‑131 (MQ=255)                                                                                 gaTCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCaaa                                          <  2:338957/133‑1 (MQ=255)                                                                                    cGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAg                                       >  2:134948/1‑133 (MQ=255)                                                                                        ttGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACg                                   <  1:55261/133‑1 (MQ=255)                                                                                        ttGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACg                                   <  2:100847/133‑1 (MQ=255)                                                                                                ttCTGACTTCGATAATATTGCGCTGGGCCAGGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCgg                           <  1:221136/133‑1 (MQ=255)                                                                                                ttCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCgg                           <  1:157161/133‑1 (MQ=255)                                                                                                 tCTGACTACGATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTGAGGCGTTGTGCAAATAGTACGGAGGGCggg                          >  1:226441/1‑133 (MQ=255)                                                                                                          gATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGt                 >  1:7371/1‑133 (MQ=255)                                                                                                          gATAAAATCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGt                 >  2:455146/1‑133 (MQ=255)                                                                                                               aaTCGCGCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTg            <  2:638155/133‑1 (MQ=255)                                                                                                                  cgcgCTGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAg         <  2:516353/133‑1 (MQ=255)                                                                                                                       tGGGCCACGCTTAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCt    <  1:380232/133‑1 (MQ=255)                                                                                                                         ggCCACGCTTAATACGTTGATTTTGCGGTGAACTATGCCAAGAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtga  <  2:226441/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                            AAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGA  >  NC_000913/4628062‑4628314

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