Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_014034	1,732,849	C→G	22.1%	intergenic (‑273/‑130)	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,849	0	C→G	22.1%	8.1 ‑1.2 0.0	43	intergenic (‑273/‑130)	RCAP_rcc01601/tgt	FeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/10):  ref base (6/27):  total (6/37)
Fisher's exact test strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GTC....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCG  >  NC_014034/1732749‑1732946                                                                                                     |                                                                                                             aGCGGGACCAGGCCCCCCTGCTGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                              >  2:573858/1‑101 (MQ=255) aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                              <  1:898027/101‑1 (MQ=255) aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                              >  2:2137016/1‑101 (MQ=255)    gggACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggg                                                                                                           >  1:520824/1‑101 (MQ=255)     ggACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgggg                                                                                                          <  1:527480/101‑1 (MQ=255)       aCCAGGCCTCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgggtgg                                                                                                        >  1:551052/1‑101 (MQ=255)         cAGGCCCCCCTGGAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgggggggg                                                                                                      >  1:2564523/1‑99 (MQ=255)         cAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCgggggggg                                                                                                      >  1:3486226/1‑99 (MQ=255)         cAGCCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTg                                                                                                      <  1:2875573/101‑1 (MQ=255)          aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                     >  1:3573260/1‑98 (MQ=255)          aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                     >  1:1114196/1‑98 (MQ=255)          aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                     >  2:4292290/1‑98 (MQ=255)          aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                     >  1:1497029/1‑98 (MQ=255)           ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGaa                                                                                                    <  2:1892978/101‑1 (MQ=255)            gCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggggg                                                                                                   >  1:967293/1‑96 (MQ=255)            gCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCggggggggggg                                                                                                   >  1:328606/1‑96 (MQ=255)             ccccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...c                                                                                               <  2:3988796/101‑1 (MQ=255)              cccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...cc                                                                                              >  1:341403/1‑101 (MQ=255)               ccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGCGCgggcggaggtggcc                                                                                                >  2:1576195/1‑93 (MQ=255)                cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...cccc                                                                                            <  2:4553171/101‑1 (MQ=255)                cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...cccc                                                                                            <  1:1765949/101‑1 (MQ=255)                 ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...ccccc                                                                                           <  1:458013/101‑1 (MQ=255)                 ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...ccccc                                                                                           >  2:3759142/1‑101 (MQ=255)                  cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...cccccc                                                                                          <  1:440951/101‑1 (MQ=255)                   tGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....g                                                                                     <  1:2245203/101‑1 (MQ=255)                    gCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...CCCCCC....gc                                                                                    >  2:1356922/1‑100 (MQ=255)                     cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....gcg                                                                                   <  1:1302915/101‑3 (MQ=255)                     cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....gcg                                                                                   >  2:378361/1‑99 (MQ=255)                     cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...ccccccccg                                                                                       >  2:768732/1‑98 (MQ=255)                      aGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gccccccccg                                                                                      >  1:1794362/1‑97 (MQ=255)                        ccGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGGGGGGGGTGAAG...CCCCCC....gcgggt                                                                                <  1:2226069/101‑6 (MQ=255)                        ccGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...ccgcccccgggg                                                                                    >  2:510491/1‑95 (MQ=255)                         cGCATGCCGTGCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAGC...ccccccccggcgg                                                                                   >  2:3154175/1‑94 (MQ=255)                          gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gcccccccgcgggg                                                                                  >  1:3378073/1‑93 (MQ=255)                          gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gccccccccggcgg                                                                                  >  2:427839/1‑93 (MQ=255)                           cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gcccccccggggggc                                                                                 >  2:4198391/1‑92 (MQ=255)                           cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAAG...gcccccccgggggcc                                                                                 >  1:3167740/1‑92 (MQ=255)                                         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)                                                                                                     |                                                                                                             AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GTC....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCG  >  NC_014034/1732749‑1732946 Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG