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

CTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....GT...C...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTC.GGCCCGG  >  NC_014034/1732766‑1732947
                                                                                                |                                                                                                           
ctccatccGCATGCCTGCCGCGTCACCAGGCAAGGCCCCCGTTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGGGGGGGCTGAAG...tccccc                                                                                                    <  2:911180/95‑1 (MQ=255)
cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGGAG...cgcccc                                                                                                    >  2:4210060/1‑101 (MQ=255)
 tGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGGAG...gcccccc                                                                                                   >  2:1011195/1‑100 (MQ=255)
  gCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....gc                                                                                             <  1:216604/101‑2 (MQ=255)
  gCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....gc                                                                                             >  1:782739/1‑100 (MQ=255)
         cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....gcgggtcgc                                                                                      <  2:2528483/101‑9 (MQ=255)
          aTGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGGAG...gccccccccggggcgc                                                                                          >  2:3906379/1‑91 (MQ=255)
           tGCCGTCCGCGTCAACAGGCAAGGCCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....gcgggtcgcct                                                                                    <  2:2047324/101‑11 (MQ=37)
                       acgcgggCCAGGGCCCCCTTTTTCCTGAGGCAAAATCCCCCTTTTTTTTTTTG.GGGGGGGGGGGGGGGTGAAC...CCCCCC.....GC...GGGTCG.CCTC.GGGCTTGccc                                                                    <  2:225427/96‑1 (MQ=2)
                         aCAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATCCGCCCTTCTTTCTTTC.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....GCGGGT...CG.CCTC.GGGCTTGCCCG.Ag                                                                <  2:3309517/101‑1 (MQ=18)
                          cAGGCCAGGGCCCCGCTTTTCCT.AGCCAAAATCCTCCCTTCTTTTCTTG.CGGGTGGGCGGGAGGTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.Agg                                                               <  2:509640/101‑1 (MQ=9)
                               cgcggCCCCTCTTTTCCA.AGACAAAA.CCCCCCTTTTTTTTTTT.GGGGGGGGGGGGGGGGGCAC...CCCCCCCCGCGGG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.aa                                                         <  2:947871/98‑1 (MQ=2)
                                         ctgttCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTcccccc                                              <  1:2116335/98‑1 (MQ=18)
                                            ttCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTg                                           <  1:2103321/101‑1 (MQ=21)
                                               tagaGGAAAAATCCGCCTTTTTTTCTTTG.GGGGGGGGCGGGGGGTGAAC...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTTCCCG.AGGC.G.AAAGCTTCCCCCCGTGaa                                         <  2:4334514/98‑1 (MQ=2)
                                               ctttggCAAAATTCCCCTTTTTTTTTTTG.CGGGGGGGGGGGGGGTGACC...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGaa                                         <  2:3565551/97‑1 (MQ=9)
                                                    gCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCt                                    <  2:2080889/101‑1 (MQ=21)
                                                      aaaaTACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGGAGGC.CCCCCCG....GG...G...CG.CCCCGGGGCTGGCCCC.GGGG.G.AAAGATTCCCCCCCGCAAAACCtt                                   >  2:622029/1‑101 (MQ=255)
                                                       aaaTACTCCCTTCTTTGTTTT.GGGGGGGGGGGGGGTTTAAC...CCCGCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGc                                 <  1:1401343/101‑1 (MQ=12)
                                                       aaaTACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGc                                 <  1:1049144/101‑1 (MQ=21)
                                                         atccGCCCTTCTTTCTTTT.GGGGGGGGCGGGGGGGGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCt                               <  2:379879/98‑1 (MQ=17)
                                                           ccccccTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTgc                             <  2:675348/98‑1 (MQ=18)
                                                            cccccTTTTTTTTTTTGGGGGGGGGGGGGGGGTGACC...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTgc                             <  1:291304/99‑1 (MQ=9)
                                                            cccccTTTTTTTTTTT.GGGGGGGGGGGGGGGGGAAC...CCCCCCCGCG.GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTgc                             <  2:1057991/99‑1 (MQ=12)
                                                            cccccTTCTTTCTTTG.CGGGGGGGGGGGGGGTGAAC...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                            <  1:1281378/99‑1 (MQ=18)
                                                              cctttctttctTTG.GGGGGGGGCGGGGGGTGAAC...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTgcgcg                          <  2:2383413/98‑1 (MQ=17)
                                                                 ttctttctTTT.CGGGGGGGCGGGGGGGGGAGGCCCCCCCC.....GG...G...CGCCCCC.GGCCTCGCCCC.GGGG.G.AAGGCTCCCCCCCCGGAAAACCTTCCCtcccccgg                        >  2:2593836/1‑94 (MQ=255)
                                                                     ttcttgT.CGGGGGGGCGGGGGGTGAAC...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCt                   <  2:1771566/101‑1 (MQ=18)
                                                                         ttt.CGGGGGGGCGGGGGGGGAAG...CCCCCCCCG..GG...G...CG.CCCC.GGGCGTGGCCGCGGGG.G.AAACATCCCCCCCCGGAAAACCTT.CCCGCGCGCGGCGCctttc               >  2:1344531/1‑99 (MQ=255)
                                                                               ggggggCGGGGGGTGAAG...CCCCCCGCG..GG...G...CG.CCTC.CGGCTTGCCCG.GGGCGG.AAACTTTCCCCCCGCGAAAACC.TGCCCTCGCGCGGCGGTTTCCCCTTc          >  2:2010236/1‑101 (MQ=255)
                                                                                  gggCGGGGGGTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTC.GGc      <  2:2950318/101‑1 (MQ=21)
                                                                                    gCGGGGGGTGAAG...CCCCCCGCG..GG...G...GG.CCCC.GGGCTTGCCCC.GGGC.GCAAAGTTTCCCCCCGGACAAACCTTGCCCCCCCCCGCGCCTTTCCCTTCC.GGcc     >  1:333941/1‑101 (MQ=2)
                                                                                      ggggggTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTC.GGCCCgg  >  1:2504441/1‑101 (MQ=21)
                                                                                      ggggggTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTC.GGCCCgg  >  2:132787/1‑101 (MQ=21)
                                                                                      ggggggTGAAG...CCCCCCGCG..GG...T...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTC.GGCCCgc  >  2:3057275/1‑100 (MQ=21)
                                                                                       gggggTGAAG...CCCCCCGCG..GG...G...CG.CCTC.CTGCTTGCCCG.AGGC.G.ACACCTTCCCCCCTTGAGACTCTTGCGTGCGCGGGTCCCTTTGTCTCTCTGGCCCgg  >  2:3685466/1‑101 (MQ=0)
                                                                                                |                                                                                                           
CTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....GT...C...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTC.GGCCCGG  >  NC_014034/1732766‑1732947

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 20 ≤  ATCG  < 32 ≤  ATCG  < 39 ≤  ATCG