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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,699,9500C→A19.5% 19.7 ‑0.8 ‑0.741intergenic (+23/+20)RCAP_rcc01568/purFRES domain family protein/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  new base (5/3):  ref base (15/18):  total (20/21)
Fisher's exact test strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.02e-01

TGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699851‑1700046
                                                                                                     |                                                                                                 
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                                                                                     >  1:235163/1‑98 (MQ=255)
            gCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGaccc                                                                                     >  2:2367339/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)
                                                                                                     |                                                                                                 
TGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTC.GCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTCGC  >  NC_014034/1699851‑1700046

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