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%	264.3 / NA	66	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 (32/34);  total (32/34)

CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816533‑3816625                                               |                                               cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                              <  1:776697/49‑1 (MQ=255) cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                              <  1:3849465/49‑1 (MQ=255) cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACt                                              <  1:4040958/49‑1 (MQ=255) cGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTg                                             <  1:830087/50‑1 (MQ=255)  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                            >  1:826395/1‑50 (MQ=255)  ggCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTgc                                            >  1:2964785/1‑50 (MQ=255)     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                         >  1:2769139/1‑50 (MQ=255)     ctgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCa                                         >  1:332964/1‑50 (MQ=255)      tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        >  1:74859/1‑50 (MQ=255)      tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        >  1:2093257/1‑50 (MQ=255)      tgctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        >  1:4432161/1‑50 (MQ=255)       gctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        >  1:2467601/1‑49 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                       <  1:2856988/49‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAtt                                       <  1:3734155/49‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        <  1:2425641/48‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCAt                                        <  1:3219527/48‑1 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      >  1:4595606/1‑50 (MQ=255)        ctgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      <  1:1796789/50‑1 (MQ=255)         tgGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATta                                      >  1:1507801/1‑49 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATtat                                     <  1:4357274/48‑1 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                   <  1:2865535/50‑1 (MQ=255)           gACGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATAc                                   <  1:1243571/50‑1 (MQ=255)            aCGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:4589428/1‑50 (MQ=255)             cGCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACg                                  >  1:53227/1‑49 (MQ=255)              gCATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGa                                <  1:4817864/50‑1 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               <  1:4987306/50‑1 (MQ=255)               cATATTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGAc                               >  1:2614134/1‑50 (MQ=255)                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              <  1:425455/50‑1 (MQ=255)                atatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACt                              >  1:2135512/1‑50 (MQ=255)                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             >  1:494330/1‑50 (MQ=255)                 tatTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACtt                             >  1:4462010/1‑50 (MQ=255)                  atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                            >  1:2592517/1‑50 (MQ=255)                  atTGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTc                            >  1:5836178/1‑50 (MQ=255)                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                           >  1:4706266/1‑50 (MQ=255)                   ttGAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTcc                           >  1:2114645/1‑50 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                         >  1:306715/1‑50 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:1473547/49‑1 (MQ=255)                     gAAATCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcg                          <  1:4450985/49‑1 (MQ=255)                      aaaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcgg                         <  1:3149508/49‑1 (MQ=255)                       aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                       >  1:2645972/1‑50 (MQ=255)                       aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                       >  1:5805084/1‑50 (MQ=255)                       aaTCTCCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCcggcg                       <  1:2586451/50‑1 (MQ=255)                           tcCGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                   >  1:1052376/1‑50 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGtt                    <  1:5935010/47‑1 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGt                     <  1:3906516/46‑1 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTa                   <  1:2950182/48‑1 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTAtt                 <  1:5681229/50‑1 (MQ=255)                             cGGCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGACGGt                     <  1:5540793/46‑1 (MQ=255)                              ggCTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTc                <  1:4225759/50‑1 (MQ=255)                                cTTGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTcc               >  1:837999/1‑49 (MQ=255)                                  tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            >  1:2581677/1‑50 (MQ=255)                                  tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            >  1:3668635/1‑50 (MQ=255)                                  tGAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            >  1:1592663/1‑50 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            <  1:4630594/49‑1 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTAt            <  1:2620418/49‑1 (MQ=255)                                   gAAACAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATc           >  1:2806088/1‑50 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc          <  1:5391708/50‑1 (MQ=255)                                    aaaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATcc          <  1:1849386/50‑1 (MQ=255)                                     aaCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:5326735/1‑50 (MQ=255)                                      aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCt         >  1:5839125/1‑49 (MQ=255)                                      aCAAATGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTg        <  1:2437230/50‑1 (MQ=255)                                        aaaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGAc      <  1:4908527/50‑1 (MQ=255)                                         aaTGTGACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACa     <  1:1290016/50‑1 (MQ=255)                                           tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   <  1:5078069/50‑1 (MQ=255)                                           tgtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAg   <  1:3255793/50‑1 (MQ=255)                                            gtgACTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAgg  >  1:5131637/1‑50 (MQ=255)                                               |                                               CGGCCTGCTGGACGCATATTGAAATCTCCGGCTTGAAACAAATGTGGCTGCGCATTATACGGACTTCCGGCGGTTATTCCTATCCTGACAAGG  >  NC_000913/3816533‑3816625

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