Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,579	G→A	intergenic (‑30/‑97)	rph ← / → yicC	truncated RNase PH/UPF0701 family protein YicC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,579	0	G	A	100.0%	85.4 / NA	26	intergenic (‑30/‑97)	rph/yicC	truncated RNase PH/UPF0701 family protein YicC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/14);  total (12/14)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816533‑3816625                                               |                                               cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                             <  1:2676470/50‑1 (MQ=255)  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                            <  1:2248410/50‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      <  1:1418715/50‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      >  1:2169868/1‑50 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      <  1:2515144/50‑1 (MQ=255)         tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      >  1:236102/1‑49 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                   <  1:1750200/50‑1 (MQ=255)            aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  <  1:1491584/50‑1 (MQ=255)             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACgg                                 <  1:1152244/50‑1 (MQ=255)             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:553508/1‑49 (MQ=255)                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              >  1:273025/1‑50 (MQ=255)                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             >  1:1637942/1‑50 (MQ=255)                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                            <  1:407779/49‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:1631277/49‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:887577/49‑1 (MQ=255)                         tctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                     >  1:1342564/1‑50 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                     <  1:2514817/46‑1 (MQ=255)                               gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc               >  1:864854/1‑50 (MQ=255)                                  tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            >  1:736992/1‑50 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            <  1:1920340/49‑1 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc           >  1:1795175/1‑50 (MQ=255)                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         <  1:894596/50‑1 (MQ=255)                                       cAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa       >  1:278917/1‑50 (MQ=255)                                          aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACaa    <  1:1787614/50‑1 (MQ=255)                                          aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAACCTGAc      >  1:1723182/1‑48 (MQ=255)                                            gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:834491/1‑50 (MQ=255)                                               |                                               CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816533‑3816625

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑