Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,860 G→C 14.8% intergenic (‑284/‑119) 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,8600G→C14.8% 3.5 ‑2.4 0.047intergenic (‑284/‑119)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (12/27):  total (12/35)
Fisher's exact test strand distribution p-value = 9.32e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....GT...C...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGCG  >  NC_014034/1732761‑1732949
                                                                                                   |                                                                                                           
ccccccTGCAGCCGCATGCCGTCCGCGTCAACCGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAGGc                                                                                                            >  1:2718732/1‑101 (MQ=255)
  ccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGccc                                                                                                          >  1:1361725/1‑101 (MQ=255)
  ccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGccc                                                                                                          >  2:508623/1‑101 (MQ=255)
   cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc                                                                                                         <  1:2040799/101‑1 (MQ=255)
   cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc                                                                                                         <  1:4138062/101‑1 (MQ=255)
    ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c                                                                                                     >  1:1349763/1‑101 (MQ=255)
    ccTGCAGCCGCATCCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc...c                                                                                                     <  1:4972910/101‑1 (MQ=255)
     cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGcccc...c.c                                                                                                   <  2:3586581/101‑1 (MQ=255)
     cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGcgcc...c.c                                                                                                   >  1:689170/1‑101 (MQ=255)
      ggCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....g                                                                                              <  1:2636303/100‑1 (MQ=255)
        cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....gcg                                                                                            <  2:2782920/101‑3 (MQ=255)
        cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGGGcccc...c.cccg                                                                                                >  2:3538509/1‑98 (MQ=255)
         aGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c.cccgg                                                                                               >  2:4642426/1‑97 (MQ=255)
         aGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGcccc...c.cccgg                                                                                               >  1:5503465/1‑97 (MQ=255)
          gCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c.ccccgg                                                                                              >  1:5588846/1‑96 (MQ=255)
             gCATGCCGTCCGCGTCAACGGGCAAGGCCCCCTCTTTTCCTAGGCAAAATACCCCCTTTTTTCTTTGCGGGGGGGCGGGGGGTGAAGCCCC...C.C....gcgggtcg                                                                                       <  1:1370606/101‑8 (MQ=25)
             gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGgccc...c.cccggggcg                                                                                           >  1:3769685/1‑93 (MQ=255)
              cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGTGGGGGGTGAAGCCCC...C.C....gcggggcgc                                                                                      >  1:4863383/1‑93 (MQ=255)
                             aaCAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGa                                                                <  1:3845554/101‑1 (MQ=21)
                                aGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGc                                                             <  2:991574/101‑1 (MQ=21)
                                          ccggtttttCCTAGGAAAAATACGCCTTTCTTTCTTTTCGGGGGGGCGGGGGGTGAACCCCC...C.C....GC...GGGTCG....CCTCGGGCTTGCCCGAGGCGAAAGCTTcc                                                   <  2:70871/96‑1 (MQ=17)
                                           ccGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....GCGGGT...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTccc                                                  <  1:3610843/101‑1 (MQ=21)
                                              cTTTTCCTAGGCAAAATCCGCCCTTCTTTCTTTGCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTcccccc                                               <  2:3336542/101‑1 (MQ=18)
                                              cTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTcccccc                                               <  2:5781576/101‑1 (MQ=21)
                                                 tcccTGGGCAAAATCCGCCCTTCTTTCTTTGCGGGGGGGCGGGGGGTGAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:3229874/99‑1 (MQ=17)
                                                   ccgaGGAAAAATTCGCCCTTTTTGCTTTGGGGGGGGGCGGGGGGGCAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaa                                          <  1:599072/98‑1 (MQ=2)
                                                   ccaaGGAAAAATCCGCCTTTTTTTCTTTTGGGGGGGGCGGGGGGGAAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaa                                          <  2:335494/98‑1 (MQ=9)
                                                    cTAGGCAAAATCCCCCCTTTTTTCTTTGGGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaaa                                         <  1:2496240/101‑1 (MQ=14)
                                                      agtcAAAATCCTCCTTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAAc                                       <  2:2537918/98‑1 (MQ=17)
                                                        gCAAAGTACCCCCTTTTTTTTTTTCGGGGGGGGGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCt                                     <  2:5385751/101‑1 (MQ=17)
                                                        gCAAAATACGCCCTTCTTTCTTGTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCt                                     <  2:1652602/101‑1 (MQ=18)
                                                         taaaaTCCGCCCTCGTTTCTTTTGGGGGGGGCGGGAGGTGAAGCCCCTGGC.G....GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                    <  2:2685915/100‑1 (MQ=12)
                                                         cAAAATCCCCCTTTTTTTTTTTGGGGGGGGGGGGGGGGTGACCCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                    <  1:2407497/101‑1 (MQ=2)
                                                         cAAAATCCCCCTTTCTTTTTTTGGGGGGGGGGGGGGGGTGAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                    <  2:5068540/101‑1 (MQ=9)
                                                         aaaattccGCCCTTCTTTCTTTGGGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCtt                                    <  2:1954077/94‑1 (MQ=14)
                                                         aaaaactcccccTTTTTTTTTTGGGGGGGGGGGGGGGGGGAACCCCC...C.CCGCGGG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCt                                     <  2:3011927/92‑1 (MQ=2)
                                                            aaTACGCCCTTCTTTCTTTTGGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGcc                                 <  2:1669152/101‑1 (MQ=18)
                                                             aTACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCt                                <  2:868409/101‑1 (MQ=21)
                                                                cGCTCTTCTTTCTTTGCGGGGGGGCGGGGGGGGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                             <  2:805220/101‑1 (MQ=17)
                                                                   ccTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcgcgc                          <  2:5079872/101‑1 (MQ=21)
                                                                        tttctttTCGGGGGGGCGGGGGGGGGAGCCCC...C.CCCG.GG...G...GG....CCCCCGGCGTGCCCCAGGCGGAAACGTCCCCCCCCGGAAACCCTGCCCCCGCGCGCCGc                     >  2:4423967/1‑101 (MQ=255)
                                                                           tttttGGGGGGGGGGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                  <  2:5312555/100‑1 (MQ=18)
                                                                            ttttCGGGGGGGGGGGGGGTGAAACCCC...C.CCGCGGG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                  <  2:1349675/101‑1 (MQ=21)
                                                                            ttttCGGGGGGGCGGGGGGGGAAGCCCC...C.C....CC...C...GGGGGCCCCCGGGCTTGCCCCGGGGGAAAGCTTCCCCCCCGGAAACCCTGCCCGCcgcgcgcgccttc                  >  2:2409782/1‑94 (MQ=255)
                                                                            tttgggggggggGGGGGGGTGACCCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTc                 <  1:4642702/96‑1 (MQ=18)
                                                                                ggggggggCGGGGGGTGACCCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCtt             <  2:1427717/100‑1 (MQ=18)
                                                                                 gggggggGGGGGGGGGAACCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCttt            <  2:5614901/101‑1 (MQ=18)
                                                                                  ggggggCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc           <  1:3473759/101‑1 (MQ=21)
                                                                                  ggggggCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc           <  2:3769442/101‑1 (MQ=21)
                                                                                  ggggggCGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc           <  2:5463259/101‑1 (MQ=21)
                                                                                   gggggGGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCg          <  2:4416929/101‑1 (MQ=21)
                                                                                    ggggCGGGGGGGGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCgg         <  1:3502339/101‑1 (MQ=21)
                                                                                      ggCGGGGGGTCAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcc       <  2:2413746/101‑1 (MQ=21)
                                                                                        cGGGGGGTGAAGCCCG...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTTCCTGCGCGCGTCGCTTTCCCTTTTGGCCCg     >  1:4881740/1‑101 (MQ=18)
                                                                                        cGGGGGGTGAAGCCCCT..CGC....GTG..T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCg     <  1:3548126/101‑1 (MQ=25)
                                                                                        cGGGGGGTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCg     >  2:2458531/1‑101 (MQ=21)
                                                                                         ggggggTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGGCTCGCTTTCCCTTTCGGCCCgg    >  2:296370/1‑101 (MQ=18)
                                                                                         ggggggTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGccccg    >  1:964044/1‑99 (MQ=21)
                                                                                         ggggggTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg    >  2:3760503/1‑101 (MQ=21)
                                                                                         ggggggTGAAGCCCC...C.CGCG.GG...T...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg    >  2:3026747/1‑101 (MQ=21)
                                                                                          gggggTGAAGCCCC...C.CGCG.GG...G...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGCCGCTTTCCCTTTCTGCCCGgc   >  2:2652351/1‑101 (MQ=21)
                                                                                           ggggTGAAGCCCC...C.C....GCGGGG...CG....CCTCGGGCTTGCCCGAGGCGTAAGCTTCCCCCCGTGGAAACCTTGCCTGCGCGCGTCGCTTTCCCTTGCGGCCCGgcg  >  2:1681885/1‑101 (MQ=21)
                                                                                                   |                                                                                                           
CCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAGCCCC...C.C....GT...C...CG....CCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGCG  >  NC_014034/1732761‑1732949

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