Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,255,904	+GT	100%	intergenic (+74/+168)	RCAP_rcc01184 → / ← RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,255,904	1	.→G	100.0%	105.3	45	intergenic (+74/+168)	RCAP_rcc01184/RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer
Reads supporting (aligned to +/- strand):  new base (25/19):  ref base (0/0):  total (26/19)
*	NC_014034	1,255,904	2	.→T	100.0%	109.7	45	intergenic (+74/+168)	RCAP_rcc01184/RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer
Reads supporting (aligned to +/- strand):  new base (26/19):  ref base (0/0):  total (26/19)

GGACCGGAACCGGGGGCCTTCGG..CT..G.TCTCGGTCCGCCCGAAGAGGGGGGGCGCTGTGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCAAATATCCCGGGGGGGGGGAACTTCCC  >  NC_014034/1255843‑1256003                                                                    ||                                                                                                    ggaccggaACCGGGGGCCTTCGGGCCTCCGCTTTCGGTCAGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCAACCCGGGCTATTTCCGcc                                                                     >  2:694836/1‑101 (MQ=2)                         cccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                             >  1:109536/10‑101 (MQ=17)                         cccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                             <  2:2429170/92‑1 (MQ=17)                         cccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGcgcg                                             >  1:3169166/10‑101 (MQ=17)                          ccccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCgg                                            >  1:1974942/9‑101 (MQ=17)                           cccgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGc                                           >  1:973919/8‑101 (MQ=17)                             cgctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCCTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGGtt                                         >  1:3498132/6‑101 (MQ=18)                              gctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCttt                                        >  2:606584/5‑101 (MQ=18)                              cctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCttt                                        <  2:1532806/97‑1 (MQ=18)                               ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc                                       <  1:1192617/98‑1 (MQ=18)                               ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc                                       >  1:1481668/4‑101 (MQ=18)                               ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc                                       <  2:2843341/98‑1 (MQ=18)                                tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                      >  2:83680/3‑101 (MQ=18)                                tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                      <  1:3884891/99‑1 (MQ=18)                                tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa                                      >  1:3484318/3‑101 (MQ=18)                                 ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                     <  2:795150/100‑1 (MQ=18)                                 ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                     >  2:619885/2‑101 (MQ=18)                                 ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc                                     >  2:1653994/2‑101 (MQ=18)                                 ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGCCGCGGCTTTCAc                                     >  1:487296/2‑101 (MQ=18)                                  tcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAcc                                    <  1:3355168/101‑1 (MQ=18)                                  tcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAcc                                    >  2:248591/1‑101 (MQ=18)                                   cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt                                   >  1:1151652/1‑101 (MQ=18)                                   cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt                                   <  1:1559726/101‑1 (MQ=18)                                   cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt                                   >  1:4140493/1‑101 (MQ=18)                                   cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt                                   <  2:3362930/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt                                  <  2:4073225/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt                                  <  1:2823861/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt                                  <  1:2092326/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGAAAATGAACCCGAGGCGCGGCTTTCACCtt                                  >  2:713376/1‑101 (MQ=17)                                     gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg                                 >  1:2408705/1‑101 (MQ=18)                                     gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg                                 >  2:3639583/1‑101 (MQ=18)                                     gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg                                 <  2:410544/101‑1 (MQ=18)                                      tCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTgg                                >  2:2761157/1‑101 (MQ=18)                                      tCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTgg                                >  2:3226681/1‑101 (MQ=18)                                         gCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGcc                             >  2:2895051/1‑101 (MQ=18)                                         gCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGcc                             <  1:2211676/101‑1 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa                            >  1:1365399/1‑101 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa                            >  2:4091236/1‑101 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa                            <  1:1817064/101‑1 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCCTGGGCCa                            >  2:3925807/1‑101 (MQ=18)                                           ccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaa                           <  1:2210259/101‑1 (MQ=18)                                           ccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaa                           <  2:3542341/101‑1 (MQ=18)                                            cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGGAATGAACCCGGGGCGCGGCTTTCCCCTTGGGCCaaa                          >  1:995546/1‑101 (MQ=17)                                            cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGNTTTCACCTTGGGCCaaa                          >  2:674688/1‑101 (MQ=17)                                            cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCTaaa                          <  2:3583557/101‑1 (MQ=18)                                            cGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaaa                          <  1:2372339/101‑1 (MQ=18)                                                                  gc..cccGGCGGGCCCCCCCCCGGGATATTTCCCCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCCCCTTGGGCCAAATATCCCGGGGGGGGGGAACTTccc  <  1:1654109/100‑1 (MQ=21)                                                                    ||                                                                                                    GGACCGGAACCGGGGGCCTTCGG..CT..G.TCTCGGTCCGCCCGAAGAGGGGGGGCGCTGTGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCAAATATCCCGGGGGGGGGGAACTTCCC  >  NC_014034/1255843‑1256003 Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG