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%	101.2 / NA	27	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 (19/8);  total (19/8)

GGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAA  >  NC_000913/3816534‑3816623                                              |                                             ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                           <  1:1322093/49‑1 (MQ=255) ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                           <  1:1477326/49‑1 (MQ=255) ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                           >  1:1846682/1‑49 (MQ=255)   ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                         >  1:727809/1‑49 (MQ=255)   ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                         >  1:554247/1‑49 (MQ=255)      gctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                      >  1:1414343/1‑49 (MQ=255)       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                     <  1:2029040/49‑1 (MQ=255)       ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                      >  1:1734579/1‑48 (MQ=255)         ggACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                    >  1:744567/1‑48 (MQ=255)          gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                  >  1:2289751/1‑49 (MQ=255)           aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                 >  1:2256360/1‑49 (MQ=255)            cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                >  1:638752/1‑49 (MQ=255)             gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                >  1:287311/1‑48 (MQ=255)              cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                              >  1:460438/1‑49 (MQ=255)                tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                            <  1:1469648/49‑1 (MQ=255)                 atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                           >  1:984518/1‑49 (MQ=255)                  ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                           <  1:2329827/48‑1 (MQ=255)                     aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                       <  1:2097790/49‑1 (MQ=255)                     aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                       >  1:458397/1‑49 (MQ=255)                      aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                      >  1:1657790/1‑49 (MQ=255)                       aTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                     >  1:1121369/1‑49 (MQ=255)                       aTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                     >  1:1819926/1‑49 (MQ=255)                              gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc              >  1:2078539/1‑49 (MQ=255)                                       aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa     >  1:2041711/1‑49 (MQ=255)                                        aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc    <  1:137433/49‑1 (MQ=255)                                         aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc    >  1:1063708/1‑48 (MQ=255)                                          tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACaa  <  1:2175442/49‑1 (MQ=255)                                              |                                             GGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAA  >  NC_000913/3816534‑3816623

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