Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,732,860 G→C 17.3% 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→C17.3% 3.0 ‑1.2 0.042intergenic (‑284/‑119)RCAP_rcc01601/tgtFeS assembly SUF system protein/queuine tRNA‑ribosyltransferase
Reads supporting (aligned to +/- strand):  new base (0/8):  ref base (7/27):  total (7/35)
Fisher's exact test strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

CCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GTC....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGCG  >  NC_014034/1732761‑1732949
                                                                                                   |                                                                                                    
ccccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...c                                                                                                  <  2:3988796/101‑1 (MQ=255)
 cccccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...cc                                                                                                 >  1:341403/1‑101 (MQ=255)
   cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...cccc                                                                                               <  1:1765949/101‑1 (MQ=255)
   cccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...cccc                                                                                               <  2:4553171/101‑1 (MQ=255)
    ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...ccccc                                                                                              <  1:458013/101‑1 (MQ=255)
    ccTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...ccccc                                                                                              >  2:3759142/1‑101 (MQ=255)
     cTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...cccccc                                                                                             <  1:440951/101‑1 (MQ=255)
      tGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....g                                                                                        <  1:2245203/101‑1 (MQ=255)
       gCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...CCCCCC....gc                                                                                       >  2:1356922/1‑100 (MQ=255)
        cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....gcg                                                                                      <  1:1302915/101‑3 (MQ=255)
        cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....gcg                                                                                      >  2:378361/1‑99 (MQ=255)
        cAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...ccccccccg                                                                                          >  2:768732/1‑98 (MQ=255)
         aGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gccccccccg                                                                                         >  1:1794362/1‑97 (MQ=255)
           ccGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGGGGGGGGTGAAG...CCCCCC....gcgggt                                                                                   <  1:2226069/101‑6 (MQ=255)
           ccGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...ccgcccccgggg                                                                                       >  2:510491/1‑95 (MQ=255)
            cGCATGCCGTGCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAGC...ccccccccggcgg                                                                                      >  2:3154175/1‑94 (MQ=255)
             gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gcccccccgcgggg                                                                                     >  1:3378073/1‑93 (MQ=255)
             gCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gccccccccggcgg                                                                                     >  2:427839/1‑93 (MQ=255)
              cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAG...gcccccccggggggc                                                                                    >  2:4198391/1‑92 (MQ=255)
              cATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGAAG...gcccccccgggggcc                                                                                    >  1:3167740/1‑92 (MQ=255)
                            cAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTTTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCg                                                                 <  1:2328063/101‑1 (MQ=21)
                                  gCAAAGGCCGTGTTTTCCCTGGGTAAATTAGCCCCTTCTTTCTTTTGGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGa                                                           <  2:1260467/101‑1 (MQ=2)
                                       ggCCCCTCTTTTCCTAGGCAAAATACGCCTTTCTTTCTTTTCGGGGGGGGGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCt                                                      <  1:796689/101‑1 (MQ=18)
                                         ccccctttttCCAAAGGAAAAATCCCCTTTTTTTTTTTTGGGGGGGGGGGGGGGGGAAC...CCCCCC....CGCGGGTCGCCTCGGGCTTGCCCGAGGCGAAAGCtt                                                     <  2:1029882/95‑1 (MQ=2)
                                         ccccGCTTTTCCTAGGCAAAATACGCCTTTCTTTCTTTTGGGGGGGGGGGGGGGTGAAC...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTc                                                    <  1:2957995/101‑1 (MQ=18)
                                         ccccGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTc                                                    <  2:433158/101‑1 (MQ=21)
                                            cGCTTTTCCTAGGCAAAATCCCCCCTTCTTTCTTTTGGGGGGGGCGGGGGGTGAAG...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTcccc                                                 <  2:4653146/101‑1 (MQ=17)
                                            cGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GCGGGT.CGCCTCGGGCTTGCCCGAGGCGAAAGCTTcccc                                                 <  2:4695634/101‑1 (MQ=21)
                                              tttttCCCAGGAAAAACCCCCCTTTTTTTTTTTGGGGGGGGGGGGGGGGGGAAC...CCCCCCCGCGGGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTccccc                                                <  1:3492410/100‑1 (MQ=2)
                                                tttCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGt                                             <  2:431960/101‑1 (MQ=21)
                                                 ttCCTAGTCAAAATCCGCCTTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:2829377/101‑1 (MQ=17)
                                                 ttCCTAGGCAAATTGCGCCCTTGTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  1:1582472/101‑1 (MQ=17)
                                                 ttCCTAGGCAAAATACGTCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:3091519/101‑1 (MQ=18)
                                                 ttCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:1272338/101‑1 (MQ=21)
                                                 ttCCGAGGCAAAATACGCCCTTCTTGCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTg                                            <  2:341369/101‑1 (MQ=18)
                                                   ccaaGGCAAAATACGCCGTTATTTCTTTTGGGGGGGGCTGGGGTGGAAA...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaa                                          <  2:551001/98‑1 (MQ=9)
                                                     tAGGCAAAATCCGCCCTTCTTTCTTTGCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGA....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaaaa                                        <  2:2422101/101‑1 (MQ=18)
                                                     gaGGAAAAACCCGCCCTTTTTTCTTTGGGGGGGGGCGGGGGGGGAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGaaaa                                        <  2:1662226/100‑1 (MQ=2)
                                                          aaaaTACGCCCTTCTTTCTTTTGGGGGGGGGGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTg                                   <  1:3571116/101‑1 (MQ=18)
                                                              tACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGGGGAGGCCCCCCCCG...GGG....GGCCCCGGGCTTCCCCCGGGGGGAAGCTTCCCCCCCCGGAAACCTTCccc                                >  2:2996899/1‑100 (MQ=255)
                                                               ccGCCTTTTTTTTTTTGGGGGGGGGGGGGGGGTGAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgc                              <  1:2659567/100‑1 (MQ=9)
                                                               aCGCCCTTCTTTCTTTTCGGGGGGGGGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgc                              <  2:1382433/101‑1 (MQ=21)
                                                                cGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcg                             <  1:666771/101‑1 (MQ=21)
                                                                  cctttcttgctTTTGGGGGGGCCGGGGGGTGAGG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTgcgcg                           <  2:4616448/98‑1 (MQ=17)
                                                                        tttttttgggggggggGGGGGGGTGACC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGc                     <  1:4134160/92‑1 (MQ=12)
                                                                           tttttGGGGGGGGGGGGGGGGAAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                  <  1:761539/100‑1 (MQ=18)
                                                                           ttttgggggggggCGGGGGGTAAAC...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCttt                  <  2:1821016/95‑1 (MQ=17)
                                                                                ggggggggCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCtt             <  2:293113/100‑1 (MQ=18)
                                                                                 gggggggCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCttt            <  2:4647188/101‑1 (MQ=21)
                                                                                      ggCGGGGGGTGATG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcc       <  1:768806/101‑1 (MQ=21)
                                                                                      ggCGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcc       <  2:1496928/101‑1 (MQ=21)
                                                                                        cGGGGGGTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGcccc     >  2:3334962/1‑100 (MQ=21)
                                                                                         ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGGGGCGAAAGCTTCCCCCCGTGACAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg    >  1:1832212/1‑101 (MQ=21)
                                                                                         ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTTCCTTTCGGCCCgg    >  2:3642628/1‑101 (MQ=21)
                                                                                         ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGccccg    >  2:621535/1‑99 (MQ=21)
                                                                                         ggggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCgg    >  2:196821/1‑101 (MQ=21)
                                                                                           ggggTGAAG...CCCCCCGCG.GGT....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCTCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGgcc  >  1:1860244/1‑100 (MQ=21)
                                                                                                   |                                                                                                    
CCCCCCTGCAGCCGCATGCCGTCCGCGTCAACAGGCAAGGGCCCCGCTTTTCCTAGGCAAAATACGCCCTTCTTTCTTTTCGGGGGGGCGGGGGGTGAAG...CCCCCC....GTC....CGCCTCGGGCTTGCCCGAGGCGAAAGCTTCCCCCCGTGAAAACCTTGCCTGCGCGCGTCGCTTTCCCTTTCGGCCCGGCG  >  NC_014034/1732761‑1732949

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 21 ≤  ATCG  < 32 ≤  ATCG  < 40 ≤  ATCG