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%	66.0 / NA	20	pseudogene (151/758 nt)	glpR	DNA‑binding transcriptional repressor GlpR
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (14/6);  total (14/6)

CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCT  >  NC_000913/3560398‑3560505                                                           |                                                   cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAt                                              <  1:330954/65‑1 (MQ=255) cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGAt                                              >  1:2010544/1‑65 (MQ=255)    tCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAg                                           <  1:1675957/65‑1 (MQ=255)     cGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                          <  1:3094724/65‑1 (MQ=255)      gTGCCGCGGCGTGTTAACCGACCTGGAAGGCGGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                         <  1:1356570/65‑1 (MQ=255)          cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                     <  1:3245352/65‑1 (MQ=255)             ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGGTCAGGTTTTGCTc                                  >  1:2958735/1‑65 (MQ=255)             ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                  >  1:999316/1‑65 (MQ=255)             ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTc                                  >  1:98968/1‑65 (MQ=255)                  gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCa                             >  1:1910338/1‑65 (MQ=255)                       ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTcc                        >  1:1614353/1‑64 (MQ=39)                       ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                        >  1:1936073/1‑65 (MQ=255)                       ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                        >  1:2081152/1‑65 (MQ=255)                       ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                        >  1:3003697/1‑65 (MQ=255)                       ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                        >  1:1902204/1‑65 (MQ=255)                                aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcg               >  1:2833795/1‑65 (MQ=255)                                aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCACTGAGGAcg               >  1:1898404/1‑65 (MQ=38)                                 aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTcgc              >  1:459098/1‑65 (MQ=255)                                   gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg            >  1:1987096/1‑65 (MQ=255)                                   gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg            <  1:3068757/65‑1 (MQ=255)                                   gCAGCGCCGCACCGCCATGACGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGCGGGCgcgg            >  1:1961736/1‑65 (MQ=37)                                     aGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg          >  1:1424468/1‑65 (MQ=255)                                      gcgcCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa         >  1:973920/1‑65 (MQ=255)                                        gccgcACCGCCATGATGGCGCCGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt       >  1:1107994/1‑65 (MQ=255)                                        gccgcACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt       >  1:272421/1‑65 (MQ=255)                                             aCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt  <  1:2921693/65‑1 (MQ=255)                                                           |                                                   CGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCT  >  NC_000913/3560398‑3560505

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