Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,868 T→C 21.7% intergenic (‑292/‑111) 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,8680G→C21.7% 8.7 ‑1.3 ‑0.523intergenic (‑292/‑111)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (6/12):  total (6/17)
Fisher's exact test strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.27e-01

AACAGGCAAGGGCCCCG.CTTTTCCTAGGC.AAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC...GT...C....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGC  >  NC_014034/1732790‑1732948
                                                                                      |                                                                                           
aaCAGGCAAGGCCCCCG.TTTTTCCTAGGC.AAAATCCGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAC...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGa                                                                 <  2:2401650/101‑1 (MQ=18)
   ggggcaGGGCCCCG.TTTTTCCCTAGC.AAAAACCCCCCTTTTTTTTTTTGGGGGGGGGGGGGGGGGAAG...CCCCCC...CG...CGGGTCG..CCTCGGGCTTGCCCGAgg                                                               <  1:3802626/96‑1 (MQ=2)
   aggccgggcccccG.CTTTCCATGGAA.AAAATGC.CCCTTTTTTTTTTGCGGGGGGGGGGGGGGATAAA...CCCCCC...GG...CGGGTCG..CCTCGGGCTTGCCCGAGGc                                                              <  2:2330100/92‑1 (MQ=2)
       aaGGGCCCCG.CTTTTCCTAGGC.AAAATACGCCTTTTTTTCTTTGGGGGGGGGCGGGGGGGGAAC...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.Gaaa                                                         <  2:970752/101‑1 (MQ=9)
       aaGGGCCCCG.CTTTTCCTAGGC.AAAATACGCCCTTCTTTCTTTTCGGGGGGGTGGGGGGTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.Gaaa                                                         <  1:3505702/101‑1 (MQ=18)
           gCCCCGTCTTTCCCTGGGA.AAAATCCCCCTTTCTTTTTTTGGGGGGGGGGGGGGGGGGACG...CCCCCC...GCGGGT....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCt                                                      <  2:368073/101‑1 (MQ=2)
           gCCCCG.CTTTTCCTAGGC.AAAATACGCCCTTCTTTCTTTGCGGGGGGGGGGGGGGTGAAG...CCCCCC...GCGGGT....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCtt                                                     <  2:3313039/101‑1 (MQ=18)
              ccG.CTTTCCCAAGGA.AAAATCCGCCTTTCTTTCTTTGGGGGGGGGCGGGGGGTGAAG...CCCCCC...GCGGGT....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.ccc                                                 <  1:1283984/101‑1 (MQ=2)
               cccgttttCCTAAGCAAAAATCCCCCTTTTTTTTTTTTGGGGGGGGCGGGGGGTGAAC...CCCCCC...GC...GGGT.CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.cc                                                  <  2:359166/97‑1 (MQ=2)
               cG.CTTTTCCTAGGC.AAAATACCCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC...GC...GGGT.CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.cccc                                                <  2:1617701/101‑1 (MQ=18)
                   ttttCCTAGGC.AAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCg                                             <  2:2785093/101‑1 (MQ=21)
                      tCCTAGGC.AAAATCCGCCCTTCTTTCTTTTCGGGGGGGGGGGGGGTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGa                                          <  1:3103737/101‑1 (MQ=18)
                        cTAGGA.AAAATCCGCCCTTTTTTCTTTCGGGGGGGGCGGGGGGTGAAC...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGaaa                                        <  2:329583/101‑1 (MQ=9)
                           tac.AAAAAGCGCCGTTTTGTCTTTGGGGGGGGGGGGGGGGGGAAC...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAAcc                                     <  1:1310269/99‑1 (MQ=2)
                                  aTACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAAGGCCCCCCCCG..GG...C....GG..CCCCGGGCTTCCCCGGGGG.GGAACTTT.CCCCCCCGCAAAACCTTCcc                                >  2:3522820/1‑101 (MQ=255)
                                   tACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAAGGC.CCCCCCG..GG...G....CG..CCCCCGGCGTGGCCCAGGCGGAAAGCTTCCCCCCCGTGaaacctgccct                                >  2:3430024/1‑93 (MQ=255)
                                              ttctttTCGGGGGGGCGGGGGGGGGAGAC.CCCCCC...GG...G....GGGCCCCCGGGCTGGCCCGGGGG.GAAAGCTC.CCCCCCCGGAAAAAC.TGCCCGCGCGCGGCGcg                   >  2:2540037/1‑100 (MQ=255)
                                               tttttgggggggggGGGGGGGTGACC...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCtt                  <  1:3167391/94‑1 (MQ=17)
                                                  tttCGGGGGGGCGGGGGGTGAGG...CCCCCCGCGGG...T....CG..GCCCAGGCCTGCCCCGGGC.GAAAGAGT.CCCCCCCCGAAAACCCCCCCCCCGGGCGTCGCCTTcc               >  2:91638/1‑101 (MQ=255)
                                                       ggggggCGGGGGGGGACC...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTc          <  1:341857/101‑1 (MQ=18)
                                                         ggggGGGGGGGTGAAC...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCgg        <  1:3523073/101‑1 (MQ=21)
                                                            gCGGGGGGTGAAG...CCCCCCGCGGG...G....CG..CCTCGGGCCTGCCCGAGGC.GAAAGCCT.TCCCCCCTGCAAAACTTGCCCGCGCGCGTCGCTTTTCCTTTCGGccc     >  2:2788887/1‑101 (MQ=12)
                                                             cGGGGGGTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGGGAAAACCTTCCCTCCCCCCGCCCCTTTCCCTTTCGCCCCg    >  1:4034232/1‑101 (MQ=9)
                                                             cGGGGGGTGAAG...CCCCCCGCGGG...G....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCCTCGCTTTCCCTTTCGGCCCg    >  2:4157530/1‑101 (MQ=21)
                                                              ggggggTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg   >  2:3048376/1‑101 (MQ=21)
                                                              ggggggTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg   >  1:3420310/1‑101 (MQ=21)
                                                              ggggggTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAACTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGccccg   >  1:3548913/1‑99 (MQ=21)
                                                               gggggTGAAG...CCCCCCGCGGG...T....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGgc  <  2:3477174/101‑1 (MQ=21)
                                                                                      |                                                                                           
AACAGGCAAGGGCCCCG.CTTTTCCTAGGC.AAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC...GT...C....CG..CCTCGGGCTTGCCCGAGGC.GAAAGCTT.CCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGC  >  NC_014034/1732790‑1732948

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 15 ≤  ATCG  < 30 ≤  ATCG  < 40 ≤  ATCG