Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,575,638	5 bp→GGGGC	100%	intergenic (+40/+129)	corA → / ← RCAP_rcc01455	magnesium/cobalt transport protein CorA/hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,575,638	0	C→G	100.0%	94.8	28	intergenic (+40/+133)	corA/RCAP_rcc01455	magnesium/cobalt transport protein CorA/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/28):  ref base (0/0):  total (0/28)
*	NC_014034	1,575,639	0	C→G	100.0%	96.0	28	intergenic (+41/+132)	corA/RCAP_rcc01455	magnesium/cobalt transport protein CorA/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/28):  ref base (0/0):  total (0/28)
*	NC_014034	1,575,640	0	A→G	100.0%	93.5	28	intergenic (+42/+131)	corA/RCAP_rcc01455	magnesium/cobalt transport protein CorA/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/28):  ref base (0/0):  total (0/28)
*	NC_014034	1,575,641	0	T→G	100.0%	82.2	28	intergenic (+43/+130)	corA/RCAP_rcc01455	magnesium/cobalt transport protein CorA/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/28):  ref base (0/0):  total (0/28)
*	NC_014034	1,575,642	0	T→C	100.0%	92.9	28	intergenic (+44/+129)	corA/RCAP_rcc01455	magnesium/cobalt transport protein CorA/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (0/28):  ref base (0/0):  total (0/28)

CCTGGGCCTTCTTCAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGCCATTTG.CGGCGCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAAGAAGAAACCCC  >  NC_014034/1575564‑1575720                                                                           |||||                                                                                ccTGGGCCTTCTTCAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGc                                                           <  2:1059869/101‑1 (MQ=37)  cTGGGCCTTCTTCAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGcc                                                          <  1:2782932/101‑1 (MQ=37)    gggCCTTCTTCAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTg                                                        <  2:58884/101‑1 (MQ=37)        cttcttCAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGACCGCGGGGGGCGGCCGCCCCCCGCACCCCCCGCCTGCggg                                                    <  2:4621167/101‑1 (MQ=25)        cttcttAAAATGGAAGGACTGGTTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCggg                                                    <  1:3392113/101‑1 (MQ=25)         ttcttcAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCgggg                                                   <  2:3311460/101‑1 (MQ=37)         ttcttcAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCgggg                                                   <  2:3268341/101‑1 (MQ=37)         ttcttcAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCgggg                                                   <  2:5409676/101‑1 (MQ=37)             tcAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACg                                               <  1:1479805/101‑1 (MQ=37)                 aTGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGt                                           <  1:5413323/101‑1 (MQ=37)                     aaGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTcccc                                       <  2:1867958/101‑1 (MQ=37)                         gcTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCcgcg                                   <  2:30312/100‑1 (MQ=25)                         aCTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCcgcg                                   <  2:4768684/101‑1 (MQ=37)                               tttAGGGCTTGCGCCCCCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCt                             <  2:750789/101‑1 (MQ=25)                                ttAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCtt                            <  2:5142916/101‑1 (MQ=37)                                 tAGGGCTGGCGCCCTCGCGCGGCTCGGGCGCCCGCAGCGCGGGGGTCTGCCCCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTg                           <  1:4967697/101‑1 (MQ=21)                                  agcgCTGGCGCCCTCGCGGGGCGCGGGCGCCCGGACCGCGGGGGTCTGCCTCCCCCCGTCCCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGc                          <  2:2878575/98‑1 (MQ=255)                                   gggCTTGCGCCCACGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCTCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCg                         <  1:2261406/101‑1 (MQ=25)                                    ggCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGt                        <  2:3877375/101‑1 (MQ=37)                                      cTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCCCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTAt                      <  1:1512405/101‑1 (MQ=25)                                            cccTCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTgcgc                <  2:4217412/101‑1 (MQ=37)                                               tCGCGCGGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCaa             <  2:1790579/101‑1 (MQ=37)                                                cgcgcgGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCaag            <  2:1802107/101‑1 (MQ=37)                                                   gcgGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCaagaag         <  1:2856676/101‑1 (MQ=37)                                                   gcgGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCaagaag         <  1:1393234/101‑1 (MQ=37)                                                    cgGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAagaaga        <  1:4452369/101‑1 (MQ=37)                                                    cgGCTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAagaaga        <  1:3752272/101‑1 (MQ=37)                                                       cTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAAGAAGAAAc     <  2:1623086/101‑1 (MQ=37)                                                       cTCGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAAGAAGAAAc     <  2:1250852/101‑1 (MQ=37)                                                         cGGGCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAAGAAGAAAccc   <  1:2743955/101‑1 (MQ=37)                                                          gggCGCCCGGAACGCGGGGGGCTGCCGCCCCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAAGAAGAAAcccc  <  1:4305162/101‑1 (MQ=37)                                                                           |||||                                                                                CCTGGGCCTTCTTCAAATGGAAGGACTGGCTTTAGGGCTTGCGCCCTCGCGCGGCTCGGGCGCCCGGAACGCGGCCATTTG.CGGCGCCCGCACCCCCCGCCTGCGGGGGACGATGTCCCCCGCGCCCCCTTGCGTATTTGCGCCAAGAAGAAACCCC  >  NC_014034/1575564‑1575720 Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 22 ≤  ATCG  < 31 ≤  ATCG  < 39 ≤  ATCG