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%	140.6 / NA	35	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 (17/18);  total (17/18)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACA  >  NC_000913/3816533‑3816622                                               |                                            cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                           <  1:1105618/49‑1 (MQ=255)  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                          <  1:1834882/49‑1 (MQ=255)   gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                         <  1:2450159/49‑1 (MQ=255)     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcg                                        >  1:1469162/1‑48 (MQ=255)      tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                      >  1:1603920/1‑49 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                    <  1:768340/49‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                     <  1:438401/48‑1 (MQ=255)         tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                   >  1:210915/1‑49 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                 >  1:1943058/1‑49 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtata                                 >  1:1938603/1‑49 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                  <  1:864299/48‑1 (MQ=255)            aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                <  1:692880/49‑1 (MQ=255)             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                               <  1:1005993/49‑1 (MQ=255)              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACgg                              <  1:350848/49‑1 (MQ=255)              gCATATTGAAATC‑CCGGCTTGAAACAAATGTAACTGCGCATTATACGGa                             <  1:50091/49‑1 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                             >  1:1473033/1‑49 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                             >  1:1851891/1‑49 (MQ=255)                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                            >  1:861148/1‑49 (MQ=255)                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                           >  1:2780574/1‑49 (MQ=255)                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                      >  1:2661231/1‑49 (MQ=255)                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                       <  1:1137849/48‑1 (MQ=255)                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                       <  1:1484555/48‑1 (MQ=255)                         tctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcgg                   >  1:1140253/1‑49 (MQ=255)                          ctcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                  >  1:241760/1‑49 (MQ=255)                              ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAtt              >  1:1300840/1‑49 (MQ=255)                                cTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc            >  1:1939096/1‑49 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt         <  1:535546/48‑1 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt         <  1:266612/48‑1 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc        >  1:839962/1‑49 (MQ=255)                                      aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt      <  1:2751385/49‑1 (MQ=255)                                       cAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg     <  1:1362401/49‑1 (MQ=255)                                       cAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg     <  1:1757549/49‑1 (MQ=255)                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGa    >  1:52521/1‑49 (MQ=255)                                          aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACa  <  1:1584078/49‑1 (MQ=255)                                           tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc   >  1:1808342/1‑47 (MQ=255)                                               |                                            CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACA  >  NC_000913/3816533‑3816622

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