Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,255,899	5 bp→CCCCC	100%	intergenic (+69/+169)	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,899	0	T→C	100.0%	108.0	45	intergenic (+69/+173)	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)
*	NC_014034	1,255,900	0	G→C	100.0%	107.0	45	intergenic (+70/+172)	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)
*	NC_014034	1,255,901	0	G→C	100.0%	107.9	45	intergenic (+71/+171)	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)
*	NC_014034	1,255,902	0	T→C	100.0%	109.0	45	intergenic (+72/+170)	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)
*	NC_014034	1,255,903	0	A→C	100.0%	109.7	45	intergenic (+73/+169)	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..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCAAA  >  NC_014034/1255843‑1255979                                                              |||||                                                                               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               <  2:2843341/98‑1 (MQ=18)                               ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc               <  1:1192617/98‑1 (MQ=18)                               ctttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTc               >  1:1481668/4‑101 (MQ=18)                                tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa              >  2:83680/3‑101 (MQ=18)                                tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa              >  1:3484318/3‑101 (MQ=18)                                tttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCa              <  1:3884891/99‑1 (MQ=18)                                 ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc             <  2:795150/100‑1 (MQ=18)                                 ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc             >  2:1653994/2‑101 (MQ=18)                                 ttcGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCAc             >  2:619885/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:4140493/1‑101 (MQ=18)                                   cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt           >  1:1151652/1‑101 (MQ=18)                                   cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt           <  1:1559726/101‑1 (MQ=18)                                   cGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCt           <  2:3362930/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt          <  1:2092326/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt          <  1:2823861/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCtt          <  2:4073225/101‑1 (MQ=18)                                    ggTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGAAAATGAACCCGAGGCGCGGCTTTCACCtt          >  2:713376/1‑101 (MQ=17)                                     gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg         >  1:2408705/1‑101 (MQ=18)                                     gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg         <  2:410544/101‑1 (MQ=18)                                     gTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTg         >  2:3639583/1‑101 (MQ=18)                                      tCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTgg        >  2:3226681/1‑101 (MQ=18)                                      tCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTgg        >  2:2761157/1‑101 (MQ=18)                                         gCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGcc     >  2:2895051/1‑101 (MQ=18)                                         gCCCGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGcc     <  1:2211676/101‑1 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa    >  2:4091236/1‑101 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa    >  1:1365399/1‑101 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCa    <  1:1817064/101‑1 (MQ=18)                                          cccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCCTGGGCCa    >  2:3925807/1‑101 (MQ=18)                                           ccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaa   <  2:3542341/101‑1 (MQ=18)                                           ccGAAGAGGGGGGGCGCTGCCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCaa   <  1:2210259/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)                                                              |||||                                                                               GGACCGGAACCGGGGGCCTTCGG..CT..G.TCTCGGTCCGCCCGAAGAGGGGGGGCGCTGTGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAAGGGGAAATGAACCCGGGGCGCGGCTTTCACCTTGGGCCAAA  >  NC_014034/1255843‑1255979 Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG