Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,849 C→G 19.3% 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→G19.3% 5.8 ‑3.5 0.035intergenic (‑273/‑130)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/7):  ref base (13/15):  total (13/22)
Fisher's exact test strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....GT...C...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCC  >  NC_014034/1732749‑1732944
                                                                                                      |                                                                                                                  
aGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGc                                                                                                                    >  1:2451077/1‑101 (MQ=255)
 gCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACGGGCAAGGGCCCCGCTTTTCAT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCg                                                                                                                   >  1:3806407/1‑101 (MQ=255)
  cGGGATCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCACGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCgg                                                                                                                  >  2:210593/1‑101 (MQ=255)
     gACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCggggg                                                                                                               <  1:2143278/101‑1 (MQ=255)
      aCCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGTCggaggg                                                                                                              >  2:1475999/1‑101 (MQ=255)
       ccAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCggggggg                                                                                                             >  2:3316413/1‑100 (MQ=255)
       ccAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCggggggg                                                                                                             >  1:470503/1‑100 (MQ=255)
        cAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTg                                                                                                            >  2:216577/1‑101 (MQ=255)
         aGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCggggggggg                                                                                                           >  1:846945/1‑98 (MQ=255)
          ggCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGaa                                                                                                          <  1:2082998/101‑1 (MQ=255)
           gCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGagg                                                                                                         >  2:439808/1‑99 (MQ=255)
              ccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCggggggggggggcc                                                                                                      >  1:439026/1‑93 (MQ=255)
              ccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGCGGGCgggggggggggggg                                                                                                      >  1:47186/1‑93 (MQ=255)
                ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCgggggggggggggccc                                                                                                    >  1:3082352/1‑91 (MQ=255)
                 ctccatccGCATGCCTGCCGCGTCACCAGGCAAGGCCCCCGTTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGGGGGGGCTGAAG...tccccc                                                                                                <  2:911180/95‑1 (MQ=255)
                 cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGGGGAG...cgcccc                                                                                                >  2:4210060/1‑101 (MQ=255)
                 cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGAGatcggaaga                                                                                                   >  1:572735/1‑93 (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.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGc   <  2:2950318/101‑1 (MQ=21)
                                                                                                     gCGGGGGGTGAAG...CCCCCCGCG..GG...G...GG.CCCC.GGGCTTGCCCC.GGGC.GCAAAGTTTCCCCCCGGACAAACCTTGCCCCCCCCCGCGCCTTTCCCTTCCGGcc  >  1:333941/1‑101 (MQ=2)
                                                                                                      |                                                                                                                  
AGCGGGACCAGGCCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCT.AGGCAAAATACGCCCTTCTTTCTTTT.CGGGGGGGCGGGGGGTGAAG...CCCCCC.....GT...C...CG.CCTC.GGGCTTGCCCG.AGGC.G.AAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCC  >  NC_014034/1732749‑1732944

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG