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%	129.2 / NA	33	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 (12/21);  total (12/21)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAG  >  NC_000913/3816533‑3816624                                               |                                              cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                             <  1:1914069/49‑1 (MQ=255) cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                             <  1:424625/49‑1 (MQ=255)  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                            >  1:1123180/1‑49 (MQ=255)   gCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                           <  1:1330912/49‑1 (MQ=255)     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgcgc                                         <  1:1355805/49‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                        <  1:1446840/47‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                      <  1:780109/49‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                       <  1:134175/48‑1 (MQ=255)         tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                     >  1:1008529/1‑49 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                    <  1:590448/48‑1 (MQ=255)            aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                  <  1:409373/49‑1 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                               <  1:1183333/49‑1 (MQ=255)                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                              <  1:303642/49‑1 (MQ=255)                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                              >  1:27764/1‑49 (MQ=255)                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                            <  1:614042/48‑1 (MQ=255)                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                            <  1:1955724/48‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                          <  1:1581444/48‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                          <  1:1517035/48‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                         >  1:1078639/1‑49 (MQ=255)                           tcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                    >  1:57744/1‑48 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGtt                   <  1:1972450/47‑1 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGtt                   <  1:1165109/47‑1 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                    <  1:704071/46‑1 (MQ=255)                                  tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTa            >  1:991039/1‑49 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCt             <  1:866750/47‑1 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt           >  1:212187/1‑49 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt           <  1:706594/48‑1 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc          >  1:994434/1‑49 (MQ=255)                                         aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc     <  1:655531/49‑1 (MQ=255)                                          aTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc     >  1:2077341/1‑48 (MQ=255)                                           tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc     >  1:507095/1‑47 (MQ=255)                                           tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc     >  1:175680/1‑47 (MQ=255)                                            gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg  >  1:763017/1‑49 (MQ=255)                                               |                                              CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAG  >  NC_000913/3816533‑3816624

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