Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,699,947 +G 17.4% intergenic (+20/+23) RCAP_rcc01568 → / ← purF RES domain family protein/amidophosphoribosyltransferase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,699,9471.→G17.4% 24.0 ‑1.3 ‑0.546intergenic (+20/+23)RCAP_rcc01568/purFRES domain family protein/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  new base (5/3):  ref base (15/23):  total (20/26)
Fisher's exact test strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.14e-01

CATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699847‑1700046
                                                                                                      |                                                                                                    
cATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGc                                                                                                       <  1:3867497/101‑1 (MQ=255)
 ctCTGGTGCTTTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.t                                                                                                     <  1:1569355/100‑1 (MQ=255)
 aTCTGGTGCTGTTTGGCTGGAGGGGGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.t                                                                                                     <  1:1564147/101‑1 (MQ=255)
 aTCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.t                                                                                                     <  1:1329409/101‑1 (MQ=255)
 aTCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.t                                                                                                     >  2:3633019/1‑101 (MQ=255)
 aTCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.t                                                                                                     >  1:758322/1‑101 (MQ=255)
 aTCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.t                                                                                                     >  2:1442001/1‑101 (MQ=255)
 aTCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.t                                                                                                     >  1:3511279/1‑101 (MQ=255)
  tcgggTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.Tg                                                                                                    <  2:2383615/98‑1 (MQ=255)
  gcTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.Tg                                                                                                    <  2:1448750/100‑1 (MQ=255)
    tGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.Tgc.g                                                                                                 <  1:1014237/101‑1 (MQ=255)
    tGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.Tgc.g                                                                                                 <  1:5476970/101‑1 (MQ=255)
     ggTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.Tgc.gc                                                                                                <  1:3110039/101‑1 (MQ=255)
      gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.Gcc                                                                                               >  1:4714908/1‑101 (MQ=255)
      gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.Gcc                                                                                               <  1:2277342/101‑1 (MQ=255)
       tGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCg                                                                                              >  2:2297486/1‑101 (MQ=255)
       tGCTGTTTCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCg                                                                                              <  1:1166513/101‑1 (MQ=255)
       ggCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCg                                                                                              <  1:247631/100‑1 (MQ=255)
        gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGa                                                                                             <  2:733549/101‑1 (MQ=255)
        gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGa                                                                                             >  1:733617/1‑101 (MQ=255)
        gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGa                                                                                             >  2:939068/1‑101 (MQ=255)
        gCTGTTTCGCTGGAGGGACCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGa                                                                                             <  2:2730547/101‑1 (MQ=255)
         cTGTTTCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAg                                                                                            <  1:1456328/101‑1 (MQ=255)
         cTGTTTCGCGGGAGGGACCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAg                                                                                            <  2:5129541/101‑1 (MQ=255)
          tGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAgg                                                                                           >  2:3262592/1‑101 (MQ=255)
           gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGCCCGGACATCGACAGCCTGCGGATTTCGGTTGCGGGCAGAGTTTC.GCCTCGGC.TGC.GCCGAgga                                                                                          >  2:869832/1‑100 (MQ=255)
           gTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGc                                                                                          <  1:725148/101‑1 (MQ=255)
            tttCGCTGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCg                                                                                         <  2:2499335/101‑1 (MQ=255)
            tttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCg                                                                                         <  1:1861199/101‑1 (MQ=255)
            tttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCg                                                                                         >  1:3656739/1‑101 (MQ=255)
            tttCGCGGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCg                                                                                         <  1:3168377/101‑1 (MQ=255)
              tcgggggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGac                                                                                       <  2:1926820/96‑2 (MQ=255)
              tCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCCAGGCGac                                                                                       >  2:2068950/1‑100 (MQ=255)
               cGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGacc                                                                                      >  2:749953/1‑99 (MQ=255)
                gggggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGaccc                                                                                     <  1:672034/98‑4 (MQ=255)
                gcgggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGaccc                                                                                     <  2:1108272/98‑4 (MQ=255)
                gCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGaccc                                                                                     >  2:2367339/1‑98 (MQ=255)
                gCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGaccc                                                                                     >  1:235163/1‑98 (MQ=255)
                 cTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGacccg                                                                                    >  2:487565/1‑97 (MQ=255)
                     agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGGCCTCGCC.GCC.GCCGAggggccccgccc                                                                                 >  1:4011638/1‑91 (MQ=255)
                     agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATAcccc                                                                                <  1:687105/101‑1 (MQ=255)
                      gggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGacccgcccgg                                                                               >  2:275994/1‑92 (MQ=255)
                      gggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATAccccc                                                                               >  2:399637/1‑101 (MQ=255)
                                 aTCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCa                                                                    <  1:399681/101‑1 (MQ=255)
                                       gcgcAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGG.TGC.GCCGAGGCGAAAAACCCCCCCCCCCCTTCATCttta                                                              >  2:3961700/1‑99 (MQ=255)
                                         gcAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAg                                                            >  1:2128785/1‑101 (MQ=255)
                                             cTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCt                                                        <  1:2287970/101‑1 (MQ=255)
                                              tGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCtt                                                       <  2:2128668/101‑1 (MQ=255)
                                                             aGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGa                                        >  2:342918/1‑101 (MQ=255)
                                                                                 ccccgccgggtC.GCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                     >  1:3755147/11‑101 (MQ=25)
                                                                                    ggccgggtC.GCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGccc                  <  1:3906481/94‑1 (MQ=37)
                                                                                      ccgggtC.GCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCga                >  1:4158472/6‑101 (MQ=37)
                                                                                       cgggtC.GCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgag               <  2:3125218/97‑1 (MQ=37)
                                                                                       cgggtC.GCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgag               >  2:5043845/5‑101 (MQ=37)
                                                                                       cgggtC.GCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgag               >  2:3388180/5‑101 (MQ=37)
                                                                                         ggtC.GCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaa             >  1:4985811/3‑101 (MQ=38)
                                                                                           tC.GCCTCGGCGTCA.GCCGAGGCGAAATCCCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAAc           <  1:746008/101‑1 (MQ=255)
                                                                                            c.GCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATAGGAGGGCTTCACCGGATAATCGCCCGAGAAACa          >  1:4962728/1‑101 (MQ=38)
                                                                                                   ggC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCg   >  2:5421828/1‑101 (MQ=255)
                                                                                                    gC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGc  >  2:1626663/1‑101 (MQ=255)
                                                                                                      |                                                                                                    
CATCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699847‑1700046

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG