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

GCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGG  >  NC_000913/3560392‑3560510                                                                  |                                                        gCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cg                                                         >  1:2534754/1‑64 (MQ=255)   cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑cgca                                                       >  1:3632387/1‑62 (MQ=255)       cGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGAt                                                   >  1:2107913/1‑65 (MQ=255)            gTGCCACGGCGTGTTAACCGAACTGGACGGCAGCGCCGCGCCGCCATGATGG‑CGCAGGATCAGGt                                              <  1:3163162/65‑1 (MQ=255)                cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTg                                          >  1:3152175/1‑65 (MQ=255)                   ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGAG‑CGGGGTCAGGTTTTgggc                                       >  1:2262746/1‑62 (MQ=255)                         ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATTATGG‑CGCAGGATCAGGTTTTGCTCCGCCAg                                 <  1:381100/65‑1 (MQ=39)                          tAACCGGACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                <  1:95161/65‑1 (MQ=39)                          tAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGc                                <  1:3579479/65‑1 (MQ=255)                           agccGACCTGGGAGGCAGCGCCGCTCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCt                               <  1:3125882/63‑1 (MQ=37)                           aaCCGAACTGGAAGGCGGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCt                               <  1:3100218/65‑1 (MQ=39)                           aaCCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCt                               >  1:184385/1‑65 (MQ=255)                             ccGAACTGGAAGGCAGCGCCGCCCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCa                             <  1:1411383/65‑1 (MQ=39)                             ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCa                             <  1:1770897/65‑1 (MQ=255)                             ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCa                             >  1:378524/1‑65 (MQ=255)                             ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCa                             >  1:3698204/1‑65 (MQ=255)                                   tGGAAGGCAGCGCCGCGCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAgg                       <  1:1812584/65‑1 (MQ=39)                                         gCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                 >  1:3181668/1‑65 (MQ=255)                                         gCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                 >  1:2759636/1‑65 (MQ=255)                                         gCAGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                 >  1:2655363/1‑65 (MQ=255)                                           aGCGCCGCACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggcg               >  1:1945598/1‑65 (MQ=255)                                              gccgcACCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAAt            >  1:1918759/1‑65 (MQ=255)                                                   aCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCt       >  1:339396/1‑65 (MQ=255)                                                   aCCGCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCAATGAGGTCGCGGGGAATAGGCt       >  1:2949332/1‑65 (MQ=37)                                                      gCCATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCg    >  1:999799/1‑65 (MQ=255)                                                        cATGATGG‑CGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCggg  >  1:2078798/1‑65 (MQ=255)                                                                  |                                                        GCCTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGG‑C‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGG  >  NC_000913/3560392‑3560510

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