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%	202.5 / NA	50	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 (28/22);  total (28/22)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816533‑3816625                                               |                                               cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                             <  1:564905/50‑1 (MQ=255)  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                            <  1:1231898/50‑1 (MQ=255)  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                            <  1:3027980/50‑1 (MQ=255)   gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                           <  1:2636109/50‑1 (MQ=255)   gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                           >  1:3271009/1‑50 (MQ=255)   gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                           <  1:4101450/50‑1 (MQ=255)    ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                           >  1:1942195/1‑49 (MQ=255)    ccTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                           >  1:1092496/1‑49 (MQ=255)     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                         >  1:2167115/1‑50 (MQ=255)     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                         >  1:987438/1‑50 (MQ=255)     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                         >  1:2253683/1‑50 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                       <  1:577494/49‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        <  1:2935287/48‑1 (MQ=255)         tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      >  1:783288/1‑49 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                    <  1:2440546/49‑1 (MQ=255)            aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  <  1:1754463/50‑1 (MQ=255)             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:1960929/1‑49 (MQ=255)             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:1301700/1‑49 (MQ=255)              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                >  1:3971076/1‑50 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               >  1:645416/1‑50 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               >  1:912858/1‑50 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               <  1:134367/50‑1 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               <  1:3436644/50‑1 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               <  1:1453900/50‑1 (MQ=255)                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              >  1:547155/1‑50 (MQ=255)                  atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                            >  1:1269135/1‑50 (MQ=255)                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                            <  1:1997642/49‑1 (MQ=255)                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                            <  1:498596/49‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:3699986/49‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:855013/49‑1 (MQ=255)                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggc                        >  1:2170080/1‑50 (MQ=255)                        aTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcgg                      >  1:1375931/1‑50 (MQ=255)                        aTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcgg                      >  1:1750512/1‑50 (MQ=255)                           tcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                   >  1:3751626/1‑50 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                   <  1:1691404/48‑1 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAtt                 >  1:2138045/1‑50 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTGCCGGCGGTTa                   <  1:3702332/48‑1 (MQ=255)                               gCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc               <  1:3612799/50‑1 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            <  1:2630271/49‑1 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc           >  1:3175612/1‑50 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc          >  1:2866229/1‑50 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc           <  1:405791/49‑1 (MQ=255)                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         <  1:4339432/50‑1 (MQ=255)                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:1845698/1‑50 (MQ=255)                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:1635858/1‑50 (MQ=255)                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa       >  1:4328089/1‑49 (MQ=255)                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa       >  1:1886738/1‑49 (MQ=255)                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa       >  1:88285/1‑49 (MQ=255)                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa       >  1:890524/1‑49 (MQ=255)                                            gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:855113/1‑50 (MQ=255)                                               |                                               CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816533‑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: ‑