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

CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG  >  NC_000913/3560394‑3560513                                                               |                                                           cTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGa                                                       >  1:1832095/1‑68 (MQ=255)       aTCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGt                                                 <  1:3132699/68‑1 (MQ=255)           tGCCGCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTg                                             <  1:1096608/68‑1 (MQ=255)           tGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCTCCGCCATGATGGCGCAGGATCAGGTTTTg                                             <  1:4180257/68‑1 (MQ=255)              cACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTCGGCTc                                          >  1:4069713/1‑68 (MQ=255)               aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                         >  1:2623226/1‑68 (MQ=255)               aCGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTcc                                         <  1:1158715/68‑1 (MQ=255)                 ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                       >  1:920713/1‑68 (MQ=255)                 ggCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGc                                       >  1:3945343/1‑68 (MQ=255)                    gtgtTAACCAAACTTGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAg                                    <  1:1673961/68‑1 (MQ=38)                      gtTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCt                                  <  1:169809/68‑1 (MQ=255)                       ttAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTc                                 <  1:3755202/68‑1 (MQ=255)                        tAACCAAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCa                                <  1:2153062/68‑1 (MQ=39)                           ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                             >  1:1060812/1‑68 (MQ=255)                           ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                             >  1:1338003/1‑68 (MQ=255)                           ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                             >  1:3858488/1‑68 (MQ=255)                           ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTg                             >  1:1494832/1‑68 (MQ=255)                           ccGAACTGGAAGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGGTTTGCTCCGCCAGCTCATTg                             >  1:1476930/1‑68 (MQ=255)                                    aaGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcgg                    >  1:3737504/1‑68 (MQ=255)                                     aGGCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCgcggc                   >  1:3309768/1‑68 (MQ=255)                                       gCAGCGCCGCACCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGa                 >  1:3313589/1‑68 (MQ=255)                                            gccgccCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGt            <  1:2681502/68‑1 (MQ=39)                                                 aCCGCCATGATGGCGCAGGATCAGGTTTTGCTCCCCCAGCTCATTGAGGTCGCGGCGAATAGTATGCg       >  1:1714199/1‑68 (MQ=255)                                                    gCCATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGc    >  1:1290360/1‑68 (MQ=255)                                                      cATGATGGCGCAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTg  >  1:3461839/1‑68 (MQ=255)                                                               |                                                           CTTGCGATCGTGCCACGGCGTGTTAACCGAACTGGAAGGCAGCGCCGCACCGCCATGATGGC‑CAGGATCAGGTTTTGCTCCGCCAGCTCATTGAGGTCGCGGCGAATAGTCTGCGGGCTG  >  NC_000913/3560394‑3560513

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