Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,849 C→G 21.7% 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_0140341,732,8490C→G21.7% 11.1 ‑1.2 ‑0.141intergenic (‑273/‑130)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/9):  ref base (6/26):  total (6/35)
Fisher's exact test strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.59e-01

AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...GT...C...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCC  >  NC_014034/1732749‑1732945
                                                                                                     |                                                                                                            
aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGc                                                                                                              >  2:3067514/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCgg                                                                                                            >  1:384409/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCgg                                                                                                            >  1:2572735/1‑101 (MQ=255)
     gACCAGGCCCCCCTGCAGCCGCATGCCGCCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCggggg                                                                                                         <  1:2554209/101‑1 (MQ=255)
       ccAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGt                                                                                                       <  1:3934840/101‑1 (MQ=255)
       ccAGGCCCCCCTGCAGCCGCATGCCCTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATTCCCCCTTCTTTCTTTTC.GGGGGGGCggggggg                                                                                                       <  2:1104357/101‑2 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGa                                                                                                     <  1:1567838/101‑1 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGa                                                                                                     <  1:4421487/101‑1 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGaa                                                                                                    <  2:3273589/101‑1 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGaa                                                                                                    <  1:3866241/101‑1 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGaa                                                                                                    <  1:3089479/101‑1 (MQ=255)
           gCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCgggggggggcg                                                                                                   >  2:4974108/1‑96 (MQ=255)
           accccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCgggggggggag                                                                                                   >  2:2541007/2‑96 (MQ=255)
            ccccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGTGCGGGGGGTGAAGc                                                                                                  <  1:5006427/101‑1 (MQ=255)
               cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGGGGAGgccc                                                                                               >  1:896026/1‑101 (MQ=255)
                 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)
                                                                                                     |                                                                                                            
AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTC.GGGGGGGCGGGGGGTGAAGCCCCC..C...GT...C...CG.CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCC  >  NC_014034/1732749‑1732945

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG