Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,699,949 G→C 17.5% 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→C17.5% 10.6 ‑1.1 ‑0.036intergenic (+22/+21)RCAP_rcc01568/purFRES domain family protein/amidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  new base (3/3):  ref base (7/21):  total (12/24)
Fisher's exact test strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.51e-01

TCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTC  >  NC_014034/1699849‑1700044
                                                                                                     |                                                                                                
tCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                                  <  1:247379/101‑1 (MQ=255)
tCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGACGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tg                                                                                                  <  1:3982237/101‑1 (MQ=255)
 cgggtttggTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc                                                                                                 <  2:2864539/93‑1 (MQ=255)
   ggTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tggtc                                                                                               >  2:264820/1‑98 (MQ=255)
   ggTGCTGTTTCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.Tgc.gc                                                                                              <  2:2835555/101‑1 (MQ=255)
    gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                             <  2:386146/101‑1 (MQ=255)
    gTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.Gcc                                                                                             >  1:784508/1‑101 (MQ=255)
     tGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                            >  1:1229122/1‑101 (MQ=255)
     tGCTGTTTCGCGGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCg                                                                                            <  1:1293317/101‑1 (MQ=255)
       ttGTTTCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                          <  1:994359/100‑1 (MQ=255)
       ctttttCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                          <  1:2226586/98‑1 (MQ=255)
       cTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                          <  2:2798326/101‑1 (MQ=255)
       cTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAg                                                                                          <  1:2250557/101‑1 (MQ=255)
        tGTTTCGCGGGAGGGACCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAgg                                                                                         <  2:238734/101‑1 (MQ=255)
         ttttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                        <  2:1005453/100‑1 (MQ=255)
         ttttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGc                                                                                        <  1:2324222/100‑1 (MQ=255)
          tttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCg                                                                                       <  1:4342146/101‑1 (MQ=255)
          tttCGCGGGAGGGACCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCg                                                                                       <  2:4360138/101‑1 (MQ=255)
           ttCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGGGa                                                                                      >  2:4104224/1‑101 (MQ=255)
           tgcgtgggggggcgCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGa                                                                                      <  2:2977235/92‑1 (MQ=255)
            tCGCGGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGac                                                                                     <  1:3167548/101‑2 (MQ=255)
             cGCGGGGGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacc                                                                                    <  1:1640000/101‑3 (MQ=255)
              gCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGaccc                                                                                   >  1:4328922/1‑98 (MQ=255)
               cgggagggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccg                                                                                  <  1:3371097/99‑5 (MQ=255)
                tGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGGgggggcccgc                                                                                 >  2:3957670/1‑93 (MQ=255)
                tGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgc                                                                                 <  2:1920282/101‑6 (MQ=255)
                   agggagCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGacccgcccc                                                                              >  1:312732/1‑93 (MQ=255)
                            ccgATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCt                                                                     >  1:3717993/1‑101 (MQ=255)
                                  tcGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCt                                                               <  1:445212/101‑1 (MQ=255)
                                    ggCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTg                                                             <  1:3914188/101‑1 (MQ=255)
                                         aaCTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCCAGGCGAAATACCCCCTCCGCCCTTCATCTTGAAgac                                                        >  2:3405546/1‑99 (MQ=255)
                                              aCCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGc                                                   <  1:3782175/101‑1 (MQ=255)
                                                                         ggTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGAt                        <  2:3717718/101‑1 (MQ=255)
                                                                               cacgcccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCg                   <  2:540674/91‑1 (MQ=25)
                                                                                cccgccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCGCCGGATAAGCGc                  >  2:145965/10‑101 (MQ=25)
                                                                                  gcccgggtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGccc                <  1:361603/94‑1 (MQ=37)
                                                                                  aaGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCg               >  1:2482818/1‑101 (MQ=255)
                                                                                     cgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgag             >  2:4561554/5‑101 (MQ=37)
                                                                                     cgggtCGCCTCGGC.GTCAGCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCgag             >  1:2832272/5‑101 (MQ=37)
                                                                                        gtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaaa          <  2:4432783/100‑1 (MQ=38)
                                                                                        gtCGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGaaa          >  1:2412615/2‑101 (MQ=38)
                                                                                         tCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa        <  1:4507670/101‑1 (MQ=255)
                                                                                          cGCCTCGGCGTCA.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACa        >  2:1153614/1‑101 (MQ=38)
                                                                                             cTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAggag    >  2:1892201/1‑99 (MQ=255)
                                                                                               cGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTc  >  2:795182/1‑101 (MQ=255)
                                                                                                     |                                                                                                
TCTGGTGCTGTTTCGCTGGAGGGAGCCGCCGATCTCGGCGCAACTGACCGGCCATCGCCAGCCTGCGGATTTCGGTTGAGGGAAGAGTTTCGCCTCGGC.TGC.GCCGAGGCGAAATACCCCCTCAGCCCTTCATCTTGAAGCCCTTGCCGATCTGATCGGAGGGCTTCACCGGATAATCGCCCGAGAAACAGGCGTC  >  NC_014034/1699849‑1700044

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG