Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,849 C→G 14.6% 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→G14.6% 9.2 ‑3.3 ‑0.048intergenic (‑273/‑130)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/7):  ref base (18/23):  total (18/30)
Fisher's exact test strand distribution p-value = 3.56e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCC....GT...C....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTT  >  NC_014034/1732749‑1732939
                                                                                                     |                                                                                                        
aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                          <  2:3650977/101‑1 (MQ=255)
aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                          >  2:2941807/1‑101 (MQ=255)
aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                          >  2:2410506/1‑101 (MQ=255)
aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGc                                                                                                          >  1:3344194/1‑101 (MQ=255)
 gCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCg                                                                                                         <  2:886095/101‑1 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                        >  2:3845033/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                        >  2:2841205/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                        >  1:4515584/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                        >  1:551836/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgg                                                                                                        >  2:1411593/1‑101 (MQ=255)
  cGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGAgg                                                                                                        >  2:3931340/1‑101 (MQ=255)
   gggACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCggg                                                                                                       >  2:2471737/1‑101 (MQ=255)
   gggACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCggg                                                                                                       >  2:234130/1‑101 (MQ=255)
     gACCGGGCCCCCCGGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCggggg                                                                                                     <  2:3836981/101‑1 (MQ=255)
     gACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCggggg                                                                                                     >  1:4639247/1‑101 (MQ=255)
     gACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCggggg                                                                                                     >  2:126881/1‑101 (MQ=255)
      aCCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgggggg                                                                                                    >  2:213456/1‑101 (MQ=255)
      aCCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgggggg                                                                                                    >  1:2753521/1‑101 (MQ=255)
      aCCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCggcggg                                                                                                    >  1:3060671/1‑101 (MQ=255)
        cAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgggggggg                                                                                                  >  2:17083/1‑99 (MQ=255)
        cAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTg                                                                                                  >  2:1285305/1‑101 (MQ=255)
        cAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTg                                                                                                  >  2:1108267/1‑101 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCggggggggg                                                                                                 >  1:1021472/1‑98 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGa                                                                                                 >  2:2131525/1‑101 (MQ=255)
         aGGCCCCCCCGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGa                                                                                                 <  1:3845318/101‑1 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCgggggggggg                                                                                                >  2:4383692/1‑97 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGaa                                                                                                <  1:595074/101‑1 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGaa                                                                                                <  1:1917170/101‑1 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGaa                                                                                                <  1:801820/101‑1 (MQ=255)
            gcccccCGCAGCCGCATGCCGTCCGCGTCAACATGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..c                                                                                            <  2:2290234/100‑1 (MQ=255)
             cccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..cc                                                                                           <  2:3474848/101‑1 (MQ=255)
             cccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG..cg                                                                                           >  2:1747849/1‑100 (MQ=255)
                 cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..cccccc                                                                                       >  1:2041480/1‑101 (MQ=255)
                 cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..cccccc                                                                                       <  2:852805/101‑1 (MQ=255)
                   acAGCCGCATGCCGTCCGCGTGAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG..GCGCCC....gc                                                                                 >  1:2639940/2‑100 (MQ=255)
                    cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGACCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCC....gcg                                                                                >  2:3573503/1‑99 (MQ=255)
                       ccGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCC....gcgggt                                                                             <  1:788266/101‑6 (MQ=255)
                        cGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGC.GGGGGTGGAG..CCCCCC....gcggggcg                                                                           >  2:1249094/1‑94 (MQ=255)
                          cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG..gcccccccggggggg                                                                              >  2:4429586/1‑92 (MQ=255)
                          cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG..gccccccccgggggc                                                                              >  2:562252/1‑92 (MQ=255)
                          cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAAG..CCCCCC....gcggggcgc                                                                          >  2:1663018/1‑93 (MQ=255)
                                         aaCAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTGCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGa                                                       <  1:3070274/101‑1 (MQ=21)
                                              gCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.Ga                                                 <  1:3828337/101‑1 (MQ=21)
                                              gCAAGGCCCCTGTTTTTCCCAGGGAAAATT.CCCCTTTTTTTTTTTGGGGGGGGGGGGGGGGGGCCC..CCCCCC....GG...CGGGTCGCCGCGGGCTTGCCCGAGGC.g                                                  <  2:2473163/101‑1 (MQ=2)
                                                aaGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCC....GCGGGT....CGCCTCGGGCTTGCCCGAGGC.Gaaa                                               <  2:625536/101‑1 (MQ=18)
                                                 agacccctGCTTTTCCTAGGAAAATGC.CGCCCTTTTTTC.TTTCGGGGGGGGCGGGGGTAAAG..CCCCCCCGCGGG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAg                                              <  2:2041340/97‑1 (MQ=2)
                                                      cccGTTTTTCCCAGGCAAAATG.CGCCCTTCTTTCTTTTGGGGGGGGGGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTcc                                         <  1:2131675/101‑1 (MQ=17)
                                                      cccGCTTTTCCTAGGCAAAATC.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTcc                                         <  1:2819242/101‑1 (MQ=18)
                                                         gCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAC..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTccccc                                      <  1:612871/101‑1 (MQ=21)
                                                          cTTTTCCTAGGCAAAATC.CGCCCTGCTTTCTTTGCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTcccccc                                     <  2:2216636/101‑1 (MQ=17)
                                                               ccTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGaa                                <  2:3038981/101‑1 (MQ=21)
                                                                 tAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGaaaa                              <  1:4631361/101‑1 (MQ=21)
                                                                 tAGCCAAAATC.CGCCTTTTTTTCTTTGGGGGGGGGCGGGGGGTAAAC..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGaaaa                              <  1:1388593/101‑1 (MQ=9)
                                                                 caagCAAAATCCCGCCTTTTTTTCTTTGGGGGGGGGGGGGGGGTAACC..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGaaa                               <  2:4139376/98‑1 (MQ=2)
                                                                   ggCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAAcc                            <  1:2228877/101‑1 (MQ=21)
                                                                    gaaaagTA.CCCCCTTCTTTCTTTCGGGGGGGGCGGGGGGTGAAC..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCt                           <  2:2568441/99‑1 (MQ=17)
                                                                        aaTA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAAGGCCCCCCCG...GG...G....GGCCTCGGGGTTGGCCGCGGCGGGAGGCTCCCCCCCCTGAAAActcgcc                        >  2:3650040/1‑96 (MQ=255)
                                                                         aTA.CGCCCTTCTTTCTTTTGGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCt                      <  2:4638946/101‑1 (MQ=18)
                                                                          atccGCCTTGCTTTCTTTGGGGGGGGGGGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCt                      <  2:4515289/98‑1 (MQ=9)
                                                                           tccGCCTTTTTTTCTTTTGGGGGGGGGGGGGGGTGAGC..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTg                     <  1:1265839/99‑1 (MQ=12)
                                                                                ccTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTgcgcgc                <  2:2866880/101‑1 (MQ=21)
                                                                                  ttctttttTTTGGGGGGGGGGGGGGGTAAAC..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGt              <  1:2787169/101‑1 (MQ=18)
                                                                                            tggggggggGGGGGGGGGAAC..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCt    <  1:1178798/99‑1 (MQ=18)
                                                                                            gggggggggGGGGGGGTGACG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCt    <  1:4502018/99‑1 (MQ=18)
                                                                                             ggggggggGGGGGGGTGAAG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCtt   <  1:2841421/100‑1 (MQ=18)
                                                                                              gggggggGGGGGGGGGACG..CCCCCCGCG.GG...T....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCttt  <  1:331732/101‑1 (MQ=18)
                                                                                                     |                                                                                                        
AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATA.CGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG..CCCCCC....GT...C....CGCCTCGGGCTTGCCCGAGGC.GAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTT  >  NC_014034/1732749‑1732939

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG