Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,255,868	+C	100%	intergenic (+38/+204)	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,868	1	.→C	100.0%	144.7	57	intergenic (+38/+204)	RCAP_rcc01184/RCAP_rcc01185	CsbD family protein/methyl‑accepting chemotaxis sensory transducer
Reads supporting (aligned to +/- strand):  new base (23/33):  ref base (0/1):  total (23/34)

AAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GG.....CTG.TC.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..TGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAA  >  NC_014034/1255795‑1255938                                                                                 |                                                                            aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggg                                                      >  2:895453/1‑101 (MQ=255) aagaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggg                                                      >  1:3086497/1‑101 (MQ=255)  agaagGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggg                                                     >  1:2999180/1‑101 (MQ=255)   gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggggg                                                    <  1:662126/101‑1 (MQ=255)   gaaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGGGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggggg                                                    >  2:4591084/1‑101 (MQ=255)    aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   <  2:2126450/101‑1 (MQ=255)    aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   >  2:1777990/1‑101 (MQ=255)    aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   <  1:651559/101‑1 (MQ=255)    aaggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAggggggg                                                   >  1:2622892/1‑101 (MQ=255)     aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAgggggggg                                                  >  2:3268860/1‑100 (MQ=255)     aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgc                                                  <  2:4528672/101‑1 (MQ=255)     aggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgc                                                  >  2:404222/1‑101 (MQ=255)      ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcg                                                 <  2:4015066/101‑1 (MQ=255)      ggaagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcg                                                 <  2:3371638/101‑1 (MQ=255)      ggaagTAGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcg                                                 <  1:1108747/101‑1 (MQ=255)        aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGgcgcc                                               >  1:1630644/1‑100 (MQ=255)        aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCt                                               <  2:4235542/101‑1 (MQ=255)        aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCt                                               <  1:2741890/101‑1 (MQ=255)        aagTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCt                                               <  1:889204/101‑1 (MQ=255)         agTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtg                                              <  2:2345220/101‑1 (MQ=255)         agTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtg                                              >  2:1633078/1‑101 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                           <  2:3943798/101‑4 (MQ=255)            cgacgaCTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccc                                           <  2:410178/101‑4 (MQ=255)             gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                          >  1:2461293/1‑97 (MQ=255)             gacgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccc                                          >  1:3314977/1‑97 (MQ=255)               cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                        <  1:3720408/101‑7 (MQ=255)               cgacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgcccccc                                        <  2:2111240/101‑7 (MQ=255)                gccTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGGGCCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       <  1:2188219/99‑8 (MQ=255)                gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCCCCGG.TT.TCGGTCCGCCCGGCGAGGGGGGGTGCtctcccccc                                       >  1:1409147/1‑93 (MQ=255)                gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       >  1:660020/1‑94 (MQ=255)                gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       >  1:2967542/1‑94 (MQ=255)                gacTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGGCCGCCCGAAGAGGGGGGGCGCtgccccccg                                       <  1:1791089/101‑8 (MQ=255)                 acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTCAGGTCC..CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccac                                      >  1:3766505/1‑93 (MQ=255)                 acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgt                                      >  1:914153/1‑93 (MQ=255)                 acTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgt                                      <  1:1149137/101‑9 (MQ=255)                  cTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtc                                     <  2:3154124/101‑10 (MQ=255)                  cTTCCTGTCGAAACTCTGACCAAGGGCTGCGGGACCGGAACCGGGGGCTTTC.GGGCCC.CCGCTT.TAGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtc                                     <  2:4461020/101‑10 (MQ=255)                   ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                    >  1:1421850/1‑91 (MQ=255)                   ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                    <  1:2454384/101‑11 (MQ=255)                   ttCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCtgccccccgtcc                                    <  1:3167746/101‑11 (MQ=255)                                   ggccTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCG.CCTTCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                  <  1:2810349/99‑1 (MQ=2)                                   gACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTcccccc                  <  2:2817969/101‑1 (MQ=9)                                    gccTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCg                 <  2:2998951/100‑1 (MQ=2)                                    aCCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCg                 <  2:262693/101‑1 (MQ=9)                                      cTAGGGGTGCCGGACCGGAACGGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTGCCCCCCGC..CCCTGCGGGACTCCCCCCggg               >  2:169957/1‑101 (MQ=2)                                       tAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  2:3177762/1‑101 (MQ=2)                                       tAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  2:3257410/1‑101 (MQ=2)                                       tAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGa              >  1:2740666/1‑101 (MQ=2)                                        aGGGGTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGat             <  2:3528564/101‑1 (MQ=2)                                         ggggTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCgggcaa            >  2:4286807/1‑98 (MQ=2)                                          gggTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGatat           <  2:2204250/101‑1 (MQ=2)                                            gTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt         <  2:4666616/101‑1 (MQ=0)                                            gTGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATAttt         <  1:2707785/101‑1 (MQ=0)                                             tGCCGGACCGGAACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTc        >  2:409650/1‑101 (MQ=2)                                               ccggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCg      <  1:1997034/101‑1 (MQ=2)                                                cggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGCTATTTCCGc     >  1:1594806/1‑101 (MQ=2)                                                cggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc     <  1:90527/101‑1 (MQ=2)                                                cggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGc     >  1:260502/1‑101 (MQ=2)                                                 ggaccggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGcc    <  1:2163655/101‑1 (MQ=2)                                                   accggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  2:3490024/101‑1 (MQ=2)                                                   accggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  1:2491559/101‑1 (MQ=2)                                                   accggaACCGGGGGCCTTC.GGGCCC.CCGCTT.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..CCCCCCGTCCCTGCGGGACTCCCCCCGGGATATTTCCGCCaa  <  2:2536190/101‑1 (MQ=2)                                                                                 |                                                                            AAGAAGGAAGTCGACGACTTCCTGTCGAAACTCTGACCTAGGGGTGCCGGACCGGAACCGGGGGCCTTC.GG.....CTG.TC.TCGGTCCGCCCGAAGAGGGGGGGCGCTG..TGGTAC..CCCTGCGGGACTCCCCCCGGGATATTTCCGCCAA  >  NC_014034/1255795‑1255938 Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG