Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,732,869	C→G	12.5%	intergenic (‑293/‑110)	RCAP_rcc01601 ← / → tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	NC_014034	1,732,869	0	T→G	12.5%	6.8 0.0 ‑1.4	35	intergenic (‑293/‑110)	RCAP_rcc01601/tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/4):  ref base (6/23):  total (6/29)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.84e-02

CAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GTC....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGCG  >  NC_014034/1732789‑1732949                                                                                        |                                                                                     cAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTTTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCg                                                                 <  1:2328063/101‑1 (MQ=21)       gCAAAGGCCGTGTTTTCCCTGGGTAAATTAGCCCCTTCTTTCTTTTGGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGa                                                           <  2:1260467/101‑1 (MQ=2)            ggCCCCTCTTTTCCTAGGCAAAATACGCCTTTCTTTCTTTTCGGGGGGGGGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCt                                                      <  1:796689/101‑1 (MQ=18)              ccccctttttCCAAAGGAAAAATCCCCTTTTTTTTTTTTGGGGGGGGGGGGGGGGGAAC...CCCCCC....CGCGGGTCGCCTCGGGCTTGCCCGAGGCGAAAGCtt                                                     <  2:1029882/95‑1 (MQ=2)              ccccGCTTTTCCTAGGCAAAATACGCCTTTCTTTCTTTTGGGGGGGGGGGGGGGTGAAC...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTc                                                    <  1:2957995/101‑1 (MQ=18)              ccccGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTc                                                    <  2:433158/101‑1 (MQ=21)                 cGCTTTTCCTAGGCAAAATCCCCCCTTCTTTCTTTTGGGGGGGGCGGGGGGTGAAG...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTcccc                                                 <  2:4653146/101‑1 (MQ=17)                 cGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTcccc                                                 <  2:4695634/101‑1 (MQ=21)                   tttttCCCAGGAAAAACCCCCCTTTTTTTTTTTGGGGGGGGGGGGGGGGGGAAC...CCCCCCCGCGGGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTccccc                                                <  1:3492410/100‑1 (MQ=2)                     tttCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGt                                             <  2:431960/101‑1 (MQ=21)                      ttCCTAGTCAAAATCCGCCTTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:2829377/101‑1 (MQ=17)                      ttCCTAGGCAAATTGCGCCCTTGTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  1:1582472/101‑1 (MQ=17)                      ttCCTAGGCAAAATACGTCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:3091519/101‑1 (MQ=18)                      ttCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:1272338/101‑1 (MQ=21)                      ttCCGAGGCAAAATACGCCCTTCTTGCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:341369/101‑1 (MQ=18)                        ccaaGGCAAAATACGCCGTTATTTCTTTTGGGGGGGGCTGGGGTGGAAA...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaa                                          <  2:551001/98‑1 (MQ=9)                          tAGGCAAAATCCGCCCTTCTTTCTTTGCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGA....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaaaa                                        <  2:2422101/101‑1 (MQ=18)                          gaGGAAAAACCCGCCCTTTTTTCTTTGGGGGGGGGCGGGGGGGGAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaaaa                                        <  2:1662226/100‑1 (MQ=2)                               aaaaTACGCCCTTCTTTCTTTTGGGGGGGGGGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTg                                   <  1:3571116/101‑1 (MQ=18)                                   tACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGGCCCCCCCCG...GGG....GGCCCCGGGCTTCCCCCGGGGGGAAGCTTCCCCCCCCGGAAACCTTCccc                                >  2:2996899/1‑100 (MQ=255)                                    ccGCCTTTTTTTTTTTGGGGGGGGGGGGGGGGTGAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgc                              <  1:2659567/100‑1 (MQ=9)                                    aCGCCCTTCTTTCTTTTCGGGGGGGGGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgc                              <  2:1382433/101‑1 (MQ=21)                                     cGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                             <  1:666771/101‑1 (MQ=21)                                       cctttcttgctTTTGGGGGGGCCGGGGGGTGAGG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcgcg                           <  2:4616448/98‑1 (MQ=17)                                             tttttttgggggggggGGGGGGGTGACC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGc                     <  1:4134160/92‑1 (MQ=12)                                                tttttGGGGGGGGGGGGGGGGAAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                  <  1:761539/100‑1 (MQ=18)                                                ttttgggggggggCGGGGGGTAAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                  <  2:1821016/95‑1 (MQ=17)                                                     ggggggggCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCtt             <  2:293113/100‑1 (MQ=18)                                                      gggggggCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCttt            <  2:4647188/101‑1 (MQ=21)                                                           ggCGGGGGGTGATG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcc       <  1:768806/101‑1 (MQ=21)                                                           ggCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcc       <  2:1496928/101‑1 (MQ=21)                                                             cGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcccc     >  2:3334962/1‑100 (MQ=21)                                                              ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGGGGCGAAAGCTTCCCCCCGTGACAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg    >  1:1832212/1‑101 (MQ=21)                                                              ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTTCCTTTCGGCCCgg    >  2:3642628/1‑101 (MQ=21)                                                              ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGccccg    >  2:621535/1‑99 (MQ=21)                                                              ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg    >  2:196821/1‑101 (MQ=21)                                                                ggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCTCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGgcc  >  1:1860244/1‑100 (MQ=21)                                                                                        |                                                                                     CAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GTC....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGCG  >  NC_014034/1732789‑1732949 Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 20 ≤  ATCG  < 31 ≤  ATCG  < 39 ≤  ATCG