Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,699,949 G→C 23.6% intergenic (+22/+21) RCAP_rcc01568 → / ← purF RES domain family protein/amidophosphoribosyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,699,9490G→C23.6% 23.2 ‑0.3 ‑1.043intergenic (+22/+21)RCAP_rcc01568/purFRES domain family protein/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  new base (5/5):  ref base (18/13):  total (25/18)
Fisher's exact test strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.04e-01

TCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699849‑1700046
                                                                                                     |                                                                                                  
tgtGGTGCTGTTTGGGTGGGGGGAGCCGCCGCTCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                                    <  1:3498239/99‑1 (MQ=255)
tcgggTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                                    <  2:1844320/98‑1 (MQ=255)
tCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                                    <  2:338217/101‑1 (MQ=255)
tCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                                    >  1:347680/1‑101 (MQ=255)
 cTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc                                                                                                   >  1:4846696/1‑101 (MQ=255)
   ggTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc.gc                                                                                                >  2:4770768/1‑101 (MQ=255)
   ggTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc.gc                                                                                                >  2:2560204/1‑101 (MQ=255)
    gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                               <  2:1360301/101‑1 (MQ=255)
    gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                               >  2:1747330/1‑101 (MQ=255)
    gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                               <  1:2597953/101‑1 (MQ=255)
     tgtttttgcGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                              <  1:4045289/93‑1 (MQ=255)
     tGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGGTTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                              >  1:2534349/1‑101 (MQ=255)
     tGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                              <  2:2705192/101‑1 (MQ=255)
      gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGa                                                                                             >  2:4173577/1‑101 (MQ=255)
        tGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                           >  2:3067103/1‑101 (MQ=255)
        tGTTTCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                           <  2:630569/101‑1 (MQ=255)
         ttttCGCTGGAGGGAGCCTCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                          <  1:1827850/100‑1 (MQ=255)
         gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                          >  1:2843410/1‑101 (MQ=255)
         gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                          >  2:1145843/1‑101 (MQ=255)
         gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                          >  2:1091932/1‑101 (MQ=255)
         gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGACGc                                                                                          >  2:819462/1‑101 (MQ=255)
           ttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGa                                                                                        >  2:2359826/1‑101 (MQ=255)
             cGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacc                                                                                      <  2:493345/101‑3 (MQ=255)
             cGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacc                                                                                      <  1:392726/101‑3 (MQ=255)
             cGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacc                                                                                      >  1:3883290/1‑99 (MQ=255)
               ggggggggggCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccg                                                                                    <  2:1518826/96‑5 (MQ=38)
                tGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgc                                                                                   >  1:2111922/1‑96 (MQ=255)
                tGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgc                                                                                   >  2:120113/1‑96 (MQ=255)
                 ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcc                                                                                  >  1:1349699/1‑95 (MQ=255)
                 ggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcc                                                                                  >  2:3441923/1‑95 (MQ=255)
                  gagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgccc                                                                                 >  2:2622421/1‑94 (MQ=255)
                   agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgccgg                                                                                >  2:4181632/1‑93 (MQ=255)
                     ggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcccggg                                                                              >  2:3828006/1‑91 (MQ=255)
                     ggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCt                                                                              >  2:374325/1‑101 (MQ=255)
                              gATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAAAACCCCCTCAGCCCTTc                                                                     >  2:33599/1‑101 (MQ=255)
                                 ctcGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGCGGCGAAATAACCCCCCAGCCCTTCATc                                                                  >  1:4406060/1‑101 (MQ=255)
                                            tGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCtt                                                       <  1:192407/101‑1 (MQ=255)
                                                 ggCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGa                                                  <  1:374378/101‑1 (MQ=255)
                                                      tctccAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGa                                             <  2:319271/98‑1 (MQ=255)
                                                         ccAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGaaga                                          >  1:319324/1‑98 (MQ=255)
                                                                gCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGc                                    >  1:846498/1‑101 (MQ=38)
                                                                            tGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAACCGCCCCCCCCGCCCTTCCTCTTGCCGCCCCTGCCGCTCTGCTCGGAGGGCTCCACAGGATaac                       >  1:4103892/1‑100 (MQ=19)
                                                                               gcggcccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                     <  2:1831146/91‑1 (MQ=25)
                                                                               ccccgccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                     >  2:2831891/11‑101 (MQ=25)
                                                                                  ccccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGccc                  <  2:1131170/94‑1 (MQ=37)
                                                                                   cccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCg                 <  2:3552265/95‑1 (MQ=37)
                                                                                    ccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCga                <  1:1441376/96‑1 (MQ=37)
                                                                                      gggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgcga              >  2:4291865/4‑101 (MQ=37)
                                                                                      gggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgaga              >  1:198195/4‑101 (MQ=38)
                                                                                         tCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc           >  2:2734351/1‑101 (MQ=38)
                                                                                         tCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc           <  2:3480562/101‑1 (MQ=38)
                                                                                         tCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc           >  1:2897567/1‑101 (MQ=38)
                                                                                          cGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa          >  2:681178/1‑101 (MQ=38)
                                                                                              tCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGt     <  1:2237897/101‑1 (MQ=255)
                                                                                                 gC.TGC.GCCGAGGCGAAATCCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGc  <  2:114850/101‑1 (MQ=255)
                                                                                                     |                                                                                                  
TCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699849‑1700046

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG