Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,608,042	G→C	14.3%	intergenic (+312/‑35)	RCAP_rcc01484 → / → murD	hypothetical protein/UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,608,042	0	G→C	14.3%	5.4 ‑1.7 ‑0.1	28	intergenic (+312/‑35)	RCAP_rcc01484/murD	hypothetical protein/UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (4/20):  total (6/22)
Fisher's exact test strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.27e-01

TTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607942‑1608068                                                                                                     |                           ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          <  2:2466092/100‑1 (MQ=38) ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:368559/2‑101 (MQ=38) ctCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcc                          >  1:4203524/2‑101 (MQ=38)  tCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGCGGCcca                         <  2:558129/101‑2 (MQ=38)    cGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAt                      <  1:2243143/101‑1 (MQ=37)    cGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCCTGGTCCGC.TGGCGC.GGCCCAAt                      >  2:1782260/1‑101 (MQ=255)     gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  1:1061500/101‑1 (MQ=37)     gACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAAtt                     <  2:11731/101‑1 (MQ=37)       ccTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCa                   <  1:828935/101‑1 (MQ=37)          tttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccc                <  2:2724014/101‑1 (MQ=38)          tttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccc                <  1:3708017/101‑1 (MQ=37)           ttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcccc               >  2:3088370/1‑101 (MQ=38)           ttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcccc               <  1:959689/101‑1 (MQ=38)           ttCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAcccc               <  1:677789/101‑1 (MQ=37)            tCCCAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCAccccc              <  1:3004073/101‑1 (MQ=38)             cccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCg             <  2:159319/101‑1 (MQ=38)              ccAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGc            <  1:1880069/101‑1 (MQ=38)               cAGCGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcg           >  1:1497020/1‑100 (MQ=38)                aGCGGATCCCGCGACGGGGTGCGACGCC.TGCAGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgg          <  1:4066694/101‑3 (MQ=37)                  cGGATCCCGCGACGGGGTGCGGCGCC.TGCGGCCCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcggga        <  1:2791024/101‑5 (MQ=25)                  cGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcggga        >  2:1869800/1‑97 (MQ=37)                  cGGATCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATCCACCCCCGcggga        >  2:1133022/1‑97 (MQ=37)                   agATCTCGAGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGC.CCGGCCCAATCCACCCCCGcgggat       <  2:3422457/100‑6 (MQ=19)                     gtCCCGCGGCGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatat     <  1:4066160/100‑8 (MQ=25)                     aTCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatat     <  1:10768/101‑8 (MQ=37)                      tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  1:4175435/101‑9 (MQ=37)                      tCCCGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatt    <  1:1048526/101‑9 (MQ=37)                        ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  1:2662624/101‑11 (MQ=25)                        ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  2:3109789/101‑11 (MQ=25)                        ccGCGACGGGGTGCGACGCC.TGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGTCCGC.TGGCGC.GGCCCAATTCACCCCCGcgggatatttc  <  2:439932/101‑11 (MQ=25)                                                                                                     |                           TTCCGACCTTTTCCCAGCGGATCCCGCGACGGGGTGCGACGCCTTGCGGCGCCGCCCGAAAGGGGGGCTGCGCGCCCCCCTCTTGGGCCGCGCGGCGCGGGCCCAATTCACCCCCGCAGCTTACGGA  >  NC_014034/1607942‑1608068 Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG