Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,699,949 G→C 20.2% 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→C20.2% 17.1 ‑0.8 ‑0.941intergenic (+22/+21)RCAP_rcc01568/purFRES domain family protein/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  new base (5/3):  ref base (14/18):  total (20/21)
Fisher's exact test strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.21e-01

TCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699849‑1700046
                                                                                                      |                                                                                                  
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:2277342/101‑1 (MQ=255)
    gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.Gcc                                                                                               >  1:4714908/1‑101 (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:939068/1‑101 (MQ=255)
      gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGa                                                                                             <  2:733549/101‑1 (MQ=255)
      gCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGa                                                                                             >  1:733617/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)
                                                                                                      |                                                                                                  
TCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699849‑1700046

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