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%	103.2 / NA	29	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 (15/14);  total (15/14)

GGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816534‑3816625                                              |                                               ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                            >  1:2348751/1‑50 (MQ=255)      gctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        >  1:1700209/1‑49 (MQ=255)       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                       <  1:842770/49‑1 (MQ=255)       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      <  1:538854/50‑1 (MQ=255)       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      >  1:101810/1‑50 (MQ=255)           aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:1842191/1‑50 (MQ=255)               atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              <  1:1754365/50‑1 (MQ=255)                tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             >  1:1693734/1‑50 (MQ=255)                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                           >  1:14275/1‑50 (MQ=255)                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                           <  1:1830476/50‑1 (MQ=255)                     aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                         <  1:181231/49‑1 (MQ=255)                            cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                   <  1:1622250/48‑1 (MQ=255)                              gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc               >  1:1893353/1‑50 (MQ=255)                                ttGAAACAAATGTGACTGCGCATTATAAGGACTTCCGGCGGTTATTCCTa             >  1:1240143/1‑50 (MQ=255)                                 tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            >  1:622353/1‑50 (MQ=255)                                  gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTa             <  1:422582/48‑1 (MQ=255)                                     aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg        <  1:3627/50‑1 (MQ=255)                                      cAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa       >  1:254807/1‑50 (MQ=255)                                       aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      <  1:1720601/50‑1 (MQ=255)                                        aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      >  1:458188/1‑49 (MQ=255)                                        aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACa     <  1:2246698/50‑1 (MQ=255)                                         aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      >  1:1390160/1‑48 (MQ=255)                                         aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACaa    <  1:2205320/50‑1 (MQ=255)                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   <  1:1148084/50‑1 (MQ=255)                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   <  1:572682/50‑1 (MQ=255)                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   <  1:38968/50‑1 (MQ=255)                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   >  1:16356/1‑50 (MQ=255)                                           gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:57864/1‑50 (MQ=255)                                           gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:1332073/1‑50 (MQ=255)                                              |                                               GGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816534‑3816625

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