Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,866 +G 87.1% intergenic (‑290/‑113) 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_0140341,732,8663.→G87.1% 7.3 0.0 ‑0.231intergenic (‑290/‑113)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (5/22):  ref base (0/4):  total (5/26)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.47e-01

CTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...GT...C...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGC  >  NC_014034/1732766‑1732948
                                                                                                        |                                                                                           
cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGccccc..c                                                                                              <  2:422046/101‑1 (MQ=255)
cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGGGAAGccccc..c                                                                                              >  1:1609247/1‑101 (MQ=255)
cTGCAGCCGCATGCCGGCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGccccc..c                                                                                              <  2:3685906/101‑1 (MQ=255)
 tGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...g                                                                                          <  2:1246823/101‑1 (MQ=255)
  gCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...gc                                                                                         >  2:3477643/1‑100 (MQ=255)
   cggCCTCAGGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...gcg                                                                                        <  2:1540772/99‑3 (MQ=255)
    aGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGGGAAGacccc..cccgg                                                                                          >  2:229621/1‑97 (MQ=255)
        gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGGAGgcccc..cccgggggg                                                                                      >  1:4175091/1‑93 (MQ=255)
        gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGGGGAGgcccc..ccccggggg                                                                                      >  2:2807130/1‑93 (MQ=255)
           tGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...gcggggcgcgc                                                                                >  2:1639461/1‑91 (MQ=255)
                      ccAACAGGCAAGGGCCCCGCTTTTCCTAGGAAAAATACCCCTTTCTTTCTTTGG.GGGGGGGGGGGGGGGGAAGCCCCC..C...GCGGGT...CG.CCTCGGGCTTGccc                                                                 <  2:4476187/100‑1 (MQ=12)
                       cAACGGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCg                                                                <  2:4201378/101‑1 (MQ=21)
                           aGGCAGGGGCCCCGCTGTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGc                                                            <  2:3335678/101‑1 (MQ=18)
                            ggcgggggCCCCGTTTTTCCTGGGCAAAACCCCCCTTTTTTTCTTTGG.GGGGGGGGGGGGGGTGAACCCCCC..C...GCGGGT...CG.CCTCGGGCTTGCCCGAGGCg                                                           <  1:1639593/96‑1 (MQ=2)
                               aggcacccgtcTTTTCCGAGGCAAAAAACGCCCTTCTTTCTTTGC.GGGGGGGGGGGGTGGAAATCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGaaa                                                        <  2:1609115/91‑1 (MQ=9)
                                aGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTGTTC.GGGGGGGCGGGGGGGGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAg                                                       <  1:2118920/101‑1 (MQ=18)
                                 gggCCCCGCTTTTCCTAGCCAAAATCCGCCCTTCTTTCTTTTG.GGGGGGGCGGGGGGTGAAGCCCCC..C...GCGGGT...CG.CCTCGGGCTTGCCCGAGGCGAAAGc                                                      <  2:389012/101‑1 (MQ=17)
                                    ccccGTTTTCCTTGGGAAAAATCCGCCTTTTTTTCTTTCG.GGGGGGGGGGGGGGTGAGCCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTc                                                   <  2:4907067/101‑1 (MQ=0)
                                       cgttttgcttAGGAAAAATCCGCCTTTTTTTTTTTGG.GGGGGGGGGGGGGGGGAACCCCCC..C...GC...GGGTCG.CCTCGGGCTTGCCCGAGGCGAAAGCTTcccc                                                <  1:3700412/92‑1 (MQ=2)
                                                  ggCAAAATCCGCCTTTCTTTCTTTGG.GGGGGGGCGGGGGGGGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAAcc                                     <  2:1946435/101‑1 (MQ=14)
                                                     aaaaTACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTTAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGCGGAAACCTTg                                  >  1:403538/1‑101 (MQ=18)
                                                      aagtACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCTCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGc                                 <  1:219711/98‑1 (MQ=18)
                                                      aaaTACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGGGGAGCCCCC..CCCGGG...G...CGCCCCCGGGCTTGCCCGGGGCGAAGGGTCCCCCCCCCGAAAAcccct                                  >  2:4911729/1‑98 (MQ=255)
                                                         tccGCCCTTTTTTGTTTGG.GGGGGGGCGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTg                              <  1:2235457/99‑1 (MQ=14)
                                                           cccccTTTTTTTTTTTG.GGGGGGGGGGGGGGTAACCCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                            <  1:1294113/99‑1 (MQ=12)
                                                           cGCCTTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                            <  1:783118/101‑1 (MQ=18)
                                                           cGCCCTTTTTTCTTTTG.GGGGGGGCGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                            <  2:4808478/101‑1 (MQ=18)
                                                                 tttttctttCG.GGGGGGGGGGGGGGTGAACCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTc                      <  2:2186733/99‑1 (MQ=17)
                                                                   tttgcggtC.GGGGGGGAGGGGGGCGAAGCCNCCGGCG..GG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGc                    <  1:2459183/94‑1 (MQ=17)
                                                                      ttttgC.GGGGGGGCGGGGGGTGACCCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                 <  2:1353154/100‑1 (MQ=18)
                                                                       ttttCGGGGGGGGCGGGGGGTAAAACCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                 <  1:806752/101‑1 (MQ=18)
                                                                        tttG.GGGGGGGGGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTcc               <  1:744353/101‑1 (MQ=18)
                                                                        tttC.GGGGGGGGGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTcc               <  1:2322561/101‑1 (MQ=21)
                                                                           tggggggggGGGGGGGTGCACCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCt             <  2:2424321/99‑1 (MQ=18)
                                                                             gggggggGGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCttt           <  2:3314660/101‑1 (MQ=21)
                                                                             gggggggCGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCttt           <  2:696480/101‑1 (MQ=21)
                                                                              ggggggGGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc          <  2:384356/101‑1 (MQ=21)
                                                                               gggggCGGGGGGTGAACCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCg         <  1:1676318/101‑1 (MQ=21)
                                                                                   gCGGGGGGTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGccc     <  1:4199792/101‑1 (MQ=21)
                                                                                     ggggggTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGGGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCCGCCCgg   >  1:1217163/1‑101 (MQ=18)
                                                                                     ggggggTGAAGCCCCC..CGCGGG...G...CG.CCTCGGGCTTGCCCGAGGGGAAAGCTTCCCCCCCTGAAAACCTTGCCTGCGCGCGGCGCTTTCCCTTTTGGccccc   >  1:4268330/1‑99 (MQ=18)
                                                                                      gggggTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGGAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGGCGCCTTCCCTTTCGGCCCCgc  >  2:2069423/1‑101 (MQ=18)
                                                                                      gggggTGAAGCCCCC..CGCGGG...T...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGgc  >  1:1100806/1‑101 (MQ=21)
                                                                                                        |                                                                                           
CTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...GT...C...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGC  >  NC_014034/1732766‑1732948

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 20 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG