Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,254,557	Δ1 bp	coding (630/942 nt)	rihB ←	pyrimidine‑specific ribonucleoside hydrolase RihB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,254,554	0	G	.	100.0%	78.0 / NA	19	coding (633/942 nt)	rihB	pyrimidine‑specific ribonucleoside hydrolase RihB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (12/7);  total (12/7)

ACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGA  >  NC_000913/2254474‑2254621                                                                                 |                                                                    aCCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCcggcg                                                                     <  1:354451/81‑1 (MQ=255)   cGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCggg                                                                   >  2:354491/1‑81 (MQ=255)     ggCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCgggcc                                                                 <  2:354452/81‑3 (MQ=255)        cgccAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGgccgcc                                                             >  2:354488/1‑81 (MQ=255)         gccAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCt                                                            <  1:354453/81‑1 (MQ=255)         gccAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCt                                                            >  2:354486/1‑81 (MQ=255)              cAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCt                                                       <  1:354454/81‑1 (MQ=255)               aaGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCtt                                                      <  2:354455/81‑1 (MQ=255)                 ggCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTc                                                    >  2:354476/1‑81 (MQ=255)                              gAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAAt                                        >  1:354503/1‑80 (MQ=255)                              gAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAAt                                        >  2:354485/1‑80 (MQ=255)                                 cTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACg                                    >  1:354490/1‑81 (MQ=255)                                     gTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCg                                >  1:354487/1‑81 (MQ=255)                                       tttGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCggg                              >  1:354499/1‑81 (MQ=255)                                           agagTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGc                          <  2:354456/81‑1 (MQ=255)                                               tGAAGTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAAc                      >  2:354481/1‑81 (MQ=255)                                                  aGTTCATGATGTCGCTGAACAGCCCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGt                   >  2:354506/1‑81 (MQ=255)                                                   gTTCATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTc                  >  1:354484/1‑81 (MQ=255)                                                      cATGATGTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCtg                <  1:354457/80‑1 (MQ=255)                                                          atgTCGCTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCtggttgg           <  1:354458/81‑1 (MQ=255)                                                                cTGAACAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCtggttggt          >  1:354495/1‑76 (MQ=255)                                                                   aaCAGCTCTCCGGC‑GGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGa  <  2:354459/81‑1 (MQ=255)                                                                                 |                                                                    ACCGGGCCGCCAGCAAGGCCGTAGTTTTCGAACTGCGTTTTGAGAGTGAAGTTCATGATGTCGCTGAACAGCTCTCCGGCGGGGCCGCCTGCCCTTTCCATCCGAGCAATCACGTCCGGGGTGCAAACGGTCTGGTTGGTGAGATCGA  >  NC_000913/2254474‑2254621

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