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%	109.1 / NA	28	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 (14/14);  total (14/14)

ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAC  >  NC_000913/3816532‑3816621                                                |                                           aCGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                          >  1:2090626/1‑50 (MQ=255)  cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                         <  1:1273190/50‑1 (MQ=255)      ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                     >  1:1631816/1‑50 (MQ=255)         ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                    <  1:2068839/48‑1 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                <  1:2175196/49‑1 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                 <  1:2187646/48‑1 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                 <  1:2182704/48‑1 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                               >  1:1819906/1‑50 (MQ=255)             aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                              >  1:1934007/1‑50 (MQ=255)              cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACgg                             <  1:2141631/50‑1 (MQ=255)                   atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                        >  1:2095044/1‑50 (MQ=255)                      gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                     >  1:1001788/1‑50 (MQ=255)                      gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                     >  1:2024831/1‑50 (MQ=255)                         aTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcgg                  >  1:1679029/1‑50 (MQ=255)                           ctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                 >  1:1047872/1‑49 (MQ=255)                           ctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                 >  1:1983735/1‑49 (MQ=255)                               ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc            <  1:1816321/50‑1 (MQ=255)                               ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc            <  1:1142854/50‑1 (MQ=255)                               ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc            >  1:807647/1‑50 (MQ=255)                                 cTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCt          <  1:822415/50‑1 (MQ=255)                                   tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt        >  1:785289/1‑50 (MQ=255)                                    gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTa         <  1:2175986/48‑1 (MQ=255)                                    gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt        <  1:1303551/49‑1 (MQ=255)                                     aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc      <  1:1108501/50‑1 (MQ=255)                                      aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAttcctttcct     >  1:1467319/1‑50 (MQ=255)                                       aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt     >  1:1625348/1‑49 (MQ=255)                                         aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc  <  1:1798140/50‑1 (MQ=255)                                         aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc  <  1:896681/50‑1 (MQ=255)                                                |                                           ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAC  >  NC_000913/3816532‑3816621

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