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

CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGG  >  NC_000913/3560394‑3560510                                                               |                                                        cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                    >  1:1063517/1‑68 (MQ=255) cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                    >  1:1349135/1‑68 (MQ=255) cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATCAGGGCGCAGGa                                                    >  1:2116168/1‑68 (MQ=255)      gATGGTGCCGCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                               <  1:2250645/68‑1 (MQ=255)      gATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAgg                                               <  1:2356310/68‑1 (MQ=255)             ccACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCt                                        >  1:808782/1‑68 (MQ=255)                   cGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCa                                  >  1:3105282/1‑68 (MQ=255)                    gtgtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAg                                 <  1:1164169/68‑1 (MQ=255)                     tgtTACCCGAGCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGc                                <  1:2532227/68‑1 (MQ=38)                        tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                             <  1:1954176/68‑1 (MQ=255)                        tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                             <  1:1841078/68‑1 (MQ=255)                         aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCAt                            >  1:2398273/1‑68 (MQ=255)                           ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                          >  1:1410133/1‑68 (MQ=255)                           ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                          >  1:8504/1‑68 (MQ=255)                               aCTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGt                      >  1:1621545/1‑68 (MQ=255)                                     aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCCGCTCATTGAGGTggcggc                >  1:1333391/1‑68 (MQ=38)                                     aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                >  1:2916958/1‑68 (MQ=255)                                            gccgcACCGCCATGATGGCGCAGGATCAGGGTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAgg         >  1:1810651/1‑67 (MQ=255)                                                   cGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCcgg  >  1:2391099/1‑68 (MQ=39)                                                               |                                                        CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGG  >  NC_000913/3560394‑3560510

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