Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,699,950 +A 7.7% 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,9501.→A 7.7% 2.6 ‑1.8 ‑1.326intergenic (+23/+20)RCAP_rcc01568/purFRES domain family protein/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (9/15):  total (11/15)
Fisher's exact test strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 4.49e-02

CTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAA  >  NC_014034/1699850‑1700036
                                                                                                      |                                                                                      
gtGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc                                                                                         <  2:2239135/100‑1 (MQ=255)
cTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc                                                                                         >  1:1139436/1‑101 (MQ=255)
   gTGCTTTTTCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                     <  2:2549642/101‑1 (MQ=255)
   gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGG.Gcc                                                                                     >  1:2265883/1‑101 (MQ=255)
    tGCTGTTTCGCGGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                    <  1:1054701/101‑1 (MQ=255)
      ttGTTTCGCGGGAGGGACCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                  <  2:1643866/100‑1 (MQ=255)
      cTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCgcg                                                                                  >  1:3026027/1‑99 (MQ=255)
      cTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                  <  2:4342832/101‑1 (MQ=255)
       ttttccGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                 <  1:3880977/96‑1 (MQ=255)
       tGTTTCGGTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                 <  2:376805/101‑1 (MQ=255)
       tGTTTCGCTTGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                 >  2:1769871/1‑101 (MQ=255)
        ttttCGCGGGAGGGGGCCGCCGATCTCGGCGCATCTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                <  1:814848/100‑1 (MQ=255)
        gTTGCGGTGGGGGGGGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                <  2:1575955/101‑1 (MQ=255)
         tttctgtGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCAGCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGG.TGC.GCCGCGGCg                                                                               >  2:912924/7‑101 (MQ=255)
         tttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCg                                                                               <  2:3703957/101‑1 (MQ=255)
          ttCGCTGGAGGGAGCCGCCGATCTCGGCTCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGa                                                                              <  2:518860/101‑1 (MQ=255)
             gcgggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGaccc                                                                           <  1:3256440/98‑4 (MQ=255)
             gcgggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGaccc                                                                           <  2:2646898/98‑4 (MQ=255)
              cgggggggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccg                                                                          <  1:656376/96‑5 (MQ=255)
               tGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgc                                                                         >  2:1997626/1‑96 (MQ=255)
                  agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCgccccgcccg                                                                      >  2:4333289/1‑92 (MQ=255)
                    ggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcccggg                                                                    >  2:2663589/1‑91 (MQ=255)
                                 tcGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCt                                                       >  2:4069866/1‑101 (MQ=255)
                                                  ccATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATc                                      <  2:584049/101‑1 (MQ=255)
                                                         cAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGGCCTTCATCTTGAAGCCCTTGCCGATCTGATCgg                               >  2:2056507/1‑101 (MQ=255)
                                                            ccTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGAT.AGATCGGaaga                           >  2:3875073/1‑98 (MQ=255)
                                                                             aGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg           >  2:4386716/1‑101 (MQ=255)
                                                                              ccccgccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg           >  2:1395496/11‑101 (MQ=25)
                                                                                      tttCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaaa  <  2:1626246/101‑1 (MQ=255)
                                                                                      ggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaa   >  2:2325175/3‑101 (MQ=38)
                                                                                      ggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaa   >  2:4310254/3‑101 (MQ=38)
                                                                                      ggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaa   >  2:907216/3‑101 (MQ=38)
                                                                                                      |                                                                                      
CTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAA  >  NC_014034/1699850‑1700036

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