Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,696 +C 94.3% intergenic (‑120/‑283) 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,6961.→C94.3% 3.7 0.0 ‑0.153intergenic (‑120/‑283)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (6/44):  ref base (0/3):  total (6/47)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.97e-01

CTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.CCGGGC..GGGCGG..G..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATG  >  NC_014034/1732596‑1732778
                                                                                                        |                                                                                        
cTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCCC.GG.CG.a                                                                                           >  2:3289220/1‑101 (MQ=255)
   gagTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.ccc                                                                                       <  1:496254/101‑2 (MQ=255)
    agTTTATATAGGGTATCGGGGGGCGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCCTGCCAAGCAGGGCTTCGCCGCGGCGGGGCCGC.GG.CG.aaccg                                                                                       >  2:1349825/1‑97 (MQ=255)
    agTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.cccg                                                                                      >  2:331700/1‑99 (MQ=255)
     gTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.Cccgg                                                                                     <  2:3503900/101‑1 (MQ=255)
      tttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGGCGC.GG.CG.A.ccccgg                                                                                    >  2:4121300/1‑97 (MQ=255)
       ttATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.Cccgggc                                                                                   <  1:2132761/101‑3 (MQ=255)
        tatataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..g                                                                                 >  2:3693152/1‑101 (MQ=255)
          tataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..ggg                                                                               <  1:3327517/101‑1 (MQ=255)
          tataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.cccgg...gggcg                                                                             >  1:1077418/1‑101 (MQ=255)
           ataTGGTATAGGAGGGATAGGTAAAGGCCGTCGGTGCGGGGCGGAGAGCCGCGTTCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.AG.ACCcgggc..gggg                                                                              <  2:3624362/101‑2 (MQ=255)
           ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGGGGGGGGGCCGC.GG.CG.G.ccccggggggg                                                                               >  1:1111452/1‑92 (MQ=255)
           ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggg                                                                              >  2:673887/1‑100 (MQ=255)
           ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggg                                                                              <  1:244695/101‑2 (MQ=255)
           ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggg                                                                              <  2:4265780/101‑2 (MQ=255)
           ataGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CC.A.Cccggggcggg                                                                               >  1:1190021/1‑95 (MQ=255)
           ataGGGTATCGGGGGGAGAGGCAAAGGCCGTCGGTGCGGGGCGGAAACCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.ccccgtggggg                                                                               >  2:2161596/1‑92 (MQ=255)
            taGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGGCGC.GG.CG.G.Cccggggccggg                                                                              >  2:930105/1‑94 (MQ=255)
              gggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggct                                                                           <  1:217549/101‑5 (MQ=255)
               ggTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAGGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggctc                                                                          <  2:2742811/101‑6 (MQ=255)
                gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggctcc                                                                         <  2:4083487/101‑7 (MQ=255)
                gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggctcc                                                                         <  1:3203335/101‑7 (MQ=255)
                gTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggctcc                                                                         <  2:2487314/101‑7 (MQ=255)
                 tATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGGCGCGGG.GG.C.CCCGGC..GGGGgg..g..ggg                                                                    >  2:1651099/1‑99 (MQ=255)
                  aTCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGC..CGGGgg..g..gggc                                                                   >  2:4325521/1‑98 (MQ=255)
                   tCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggctccgcc                                                                      <  1:2091978/101‑10 (MQ=38)
                   gcGGGGGGAGAGGTAAAGGCCGTCGGGGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggctccgcc                                                                      <  2:3858938/100‑10 (MQ=38)
                    cGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCcgggc..gggggctccgccc                                                                     <  2:76431/101‑11 (MQ=38)
                      gggggAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.cccgg...gggcgg..gcgggggcgcg                                                               >  2:76043/1‑91 (MQ=255)
                              ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCAccccc                                                         <  1:711268/101‑1 (MQ=25)
                              ggTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCAccccc                                                         <  2:2075254/101‑1 (MQ=25)
                                 aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGc                                                      <  2:1973429/101‑1 (MQ=25)
                                 aaaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGc                                                      <  2:1150798/101‑1 (MQ=25)
                                  aacgCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGC..GGGGGG..G..CCGCCCCCCccccccccacc                                                   >  2:2993324/4‑95 (MQ=255)
                                  aaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCg                                                     <  1:1001759/101‑1 (MQ=25)
                                  aaGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCg                                                     <  1:187209/101‑1 (MQ=25)
                                    atctgTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTa                                                   <  2:3350387/97‑1 (MQ=21)
                                     gCCGTCGGTGCGGGGCGGAAAGCGGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTAt                                                  <  1:1614284/101‑1 (MQ=21)
                                     gCCGTCAGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTAt                                                  <  1:770119/101‑1 (MQ=21)
                                         tCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTgc                                              <  2:3544459/101‑1 (MQ=21)
                                         tCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTgc                                              <  2:1664309/101‑1 (MQ=21)
                                               cGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaag                                        <  2:2079693/101‑1 (MQ=21)
                                               cGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaag                                        <  1:3009094/101‑1 (MQ=21)
                                               cGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaag                                        <  1:1035101/101‑1 (MQ=21)
                                                  ggCGGAATGCCGGGTGCCAAGCAGTGTGTCGCTTCGGCGGGTCCGC.GG.CGCC.CCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCaagaag                                     <  2:3965149/101‑1 (MQ=9)
                                                     ggAAAGCGGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGaaa                                  <  1:2146661/101‑1 (MQ=18)
                                                       aaaGCCGCGTGCCAAGCAGGCTTGCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGTGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACa                                <  2:3652841/101‑1 (MQ=18)
                                                        aaGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACaa                               <  1:2602226/101‑1 (MQ=18)
                                                         aGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAg                              <  1:930185/101‑1 (MQ=18)
                                                              cgggCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGa                         <  1:4327914/98‑1 (MQ=18)
                                                               gTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGG..GGGGGG..GCCGCCCCCCCCCCCCCCGCGGATTTGGGGCCAGAAAAAACAAACggggc                        >  2:2228715/1‑99 (MQ=255)
                                                                tGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAcc                       <  2:3612074/101‑1 (MQ=18)
                                                                tGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGAcc                       <  1:3693427/101‑1 (MQ=18)
                                                                  ccAAGCAGGGTTTCGCTGCGGCGGGGCGGC.GGAAG.C.ACGAGC..GGGGAGTAC..GCCACCGAACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAg                     <  2:2927147/101‑1 (MQ=12)
                                                                       cAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGcccc                <  1:4410554/101‑1 (MQ=18)
                                                                          ggTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCt             <  2:840507/101‑1 (MQ=18)
                                                                           gTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTg            <  1:4121586/101‑1 (MQ=18)
                                                                           gTTTCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTg            <  2:40717/101‑1 (MQ=18)
                                                                              tCGCCGCGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAg         <  2:4664067/101‑1 (MQ=18)
                                                                                gccgcGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGcc       <  1:673949/101‑1 (MQ=18)
                                                                                gccgcGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGcc       <  1:3223963/101‑1 (MQ=18)
                                                                                 ccgcGGCGGGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCg      <  1:619161/101‑1 (MQ=18)
                                                                                  cgcgGCGGGGCCGC.GG.CG.ACCCGGGC..GGGGGGTCC..GCCCCCCCACCCCCCGCGTATTTGCGCCAAGAAGAAACAAGCGGGGGCAAGCCCCCCCGCCGCCGc     >  2:3070030/1‑101 (MQ=19)
                                                                                    cggcggGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCa  <  2:2915139/101‑2 (MQ=18)
                                                                                    cggcggGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCAt   <  1:4251077/101‑1 (MQ=18)
                                                                                     ggcggGGCCGC.GG.CG.ACCCGGGCGGGGGCTC..C..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATg  <  1:3502259/101‑1 (MQ=18)
                                                                                                        |                                                                                        
CTAGAGTTTATATAGGGTATCGGGGGGAGAGGTAAAGGCCGTCGGTGCGGGGCGGAAAGCCGCGTGCCAAGCAGGGTTTCGCCGCGGCGGGGCCGC.GG.CG.A.CCGGGC..GGGCGG..G..GCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAGAAACAAGCGGGACCAGGCCCCCCTGCAGCCGCATG  >  NC_014034/1732596‑1732778

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 38 ≤  ATCG