Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,560,455	+G	pseudogene (151/758 nt)	glpR ←	DNA‑binding transcriptional repressor GlpR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,560,455	1	.	G	100.0%	64.1 / NA	19	pseudogene (151/758 nt)	glpR	DNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/3);  total (16/3)

GGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACG  >  NC_000913/3560388‑3560516                                                                     |                                                              ggTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGc                                                                <  1:1421702/68‑1 (MQ=255)       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                          >  1:521337/1‑68 (MQ=255)       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                          >  1:2931243/1‑68 (MQ=255)       cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                          >  1:1526327/1‑68 (MQ=255)            gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                     >  1:2670822/1‑68 (MQ=255)            gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                                     >  1:3911286/1‑68 (MQ=255)              tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCCGGtt                                                   >  1:2512800/1‑68 (MQ=255)                    cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                             >  1:492303/1‑68 (MQ=255)                       ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                          >  1:3179932/1‑68 (MQ=255)                          gtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAg                                       <  1:2741378/68‑1 (MQ=255)                              tAACCGAGCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                   <  1:3373028/68‑1 (MQ=39)                              tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGCTTTGCTCCGCCACCTCa                                   >  1:1584309/1‑68 (MQ=255)                                 ccGAACTGGAAGGCAGCGCCGCACCGCCATGCTGGCGCAGGATCAGGTTTTGCTCCGCCCGCTAATTg                                >  1:2623572/1‑68 (MQ=37)                                 ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                >  1:1015026/1‑68 (MQ=255)                                 ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                >  1:1367933/1‑68 (MQ=255)                                 ccGAACTGGAAGGCAGCGCCGCACCGCCAAGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                                >  1:248946/1‑68 (MQ=39)                                           aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                      >  1:4016547/1‑68 (MQ=255)                                             gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                    >  1:1040389/1‑68 (MQ=255)                                             gCAGCGCCGCACCGCCACGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                    >  1:2815824/1‑68 (MQ=39)                                                 cgccgcCCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAg                <  1:2692248/68‑1 (MQ=39)                                                  gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt               >  1:2402470/1‑68 (MQ=255)                                                   ccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTc              >  1:1162741/1‑68 (MQ=255)                                                             atgatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGa    >  1:3726390/1‑68 (MQ=255)                                                               gatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGaag  >  1:1738429/1‑66 (MQ=39)                                                               gatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg  >  1:3007716/1‑68 (MQ=255)                                                               gatgGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACg  >  1:4043415/1‑68 (MQ=255)                                                                     |                                                              GGTGGCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTGACG  >  NC_000913/3560388‑3560516

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