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%	142.7 / NA	36	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 (24/12);  total (24/12)

ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816532‑3816625                                                |                                               aCGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                              >  1:746624/1‑50 (MQ=255)  cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                             >  1:590378/1‑50 (MQ=255)     ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                           >  1:700311/1‑49 (MQ=255)      ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                         >  1:1352952/1‑50 (MQ=255)          tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      >  1:2676222/1‑49 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                     <  1:13152/48‑1 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                     <  1:1020220/48‑1 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                   <  1:901713/50‑1 (MQ=255)            gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                   <  1:380139/50‑1 (MQ=255)              cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACgg                                 <  1:943142/50‑1 (MQ=255)              cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:421552/1‑49 (MQ=255)              cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:38240/1‑49 (MQ=255)               gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                >  1:2649067/1‑50 (MQ=255)                cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               >  1:248773/1‑50 (MQ=255)                 atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              >  1:1925928/1‑50 (MQ=255)                  tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             <  1:2195550/50‑1 (MQ=255)                    gtGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                           <  1:1897204/49‑1 (MQ=255)                     tGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          >  1:1092433/1‑50 (MQ=255)                     tGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          >  1:917765/1‑50 (MQ=255)                       aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                        >  1:120378/1‑50 (MQ=255)                              cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGtt                    <  1:2699962/47‑1 (MQ=255)                              cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAtt                 <  1:2809285/50‑1 (MQ=255)                                gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc               >  1:1860759/1‑50 (MQ=255)                                 cTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCt              <  1:1207753/50‑1 (MQ=255)                                   tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            >  1:87507/1‑50 (MQ=255)                                   tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            >  1:1507286/1‑50 (MQ=255)                                     aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc          >  1:1284389/1‑50 (MQ=255)                                       aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:793423/1‑49 (MQ=255)                                         aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa       >  1:2166031/1‑49 (MQ=255)                                          aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      >  1:1993750/1‑49 (MQ=255)                                           aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      >  1:1846336/1‑48 (MQ=255)                                            tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   <  1:1379940/50‑1 (MQ=255)                                            tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   <  1:1530294/50‑1 (MQ=255)                                            tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   >  1:1660568/1‑50 (MQ=255)                                             gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:1509944/1‑50 (MQ=255)                                             gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:2496330/1‑50 (MQ=255)                                                |                                               ACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816532‑3816625

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: ‑