Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,995,643 +C 6.2% intergenic (+183/+37) RCAP_rcc02802 → / ← RCAP_rcc02803 membrane dipeptidase/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,995,6431.→C 6.2% 3.5 0.0 ‑0.148intergenic (+183/+37)RCAP_rcc02802/RCAP_rcc02803membrane dipeptidase/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (3/0):  ref base (40/5):  total (43/5)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.99e-01

ACCCCTTCGGTTGCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCC  >  NC_014034/2995548‑2995743
                                                                                                 |                                                                                                      
aCCCCTTCGGTTGCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.Cgggc                                                                                                   >  1:2832347/1‑100 (MQ=255)
    cTTCGGTTGCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGcc                                                                                               >  1:2252209/1‑101 (MQ=255)
     ttgggTTGCCGCCCCGCACCCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CGGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGccc                                                                                              <  2:586805/98‑1 (MQ=255)
                    gCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCCGGGCg                                                                             >  1:94853/1‑101 (MQ=255)
                     cAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCggggggg                                                                            >  1:4006701/1‑98 (MQ=255)
                      aaCCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGaa                                                                           >  1:3498010/1‑101 (MQ=255)
                           gTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATccc                                                                      >  1:3112198/1‑101 (MQ=255)
                             cggcggCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGGAATccccc                                                                    >  2:1978130/1‑101 (MQ=255)
                                       gcgcCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTccc                                                          >  2:432538/1‑101 (MQ=255)
                                            ttgtttttCCTTGGAAAAATACGCCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttc                                                        >  2:1285997/5‑98 (MQ=0)
                                                ttttCCTTGGAAAAATACGCCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatct                                                    >  2:4195359/1‑94 (MQ=2)
                                                  ttCCTTGGAAAAATACGCCCGGGGGGTGCGGGGGGCGGAGCCCCCCGCCGGGTCGCCTCGGCGTCAGCCGAGGCGAAATACCCC.CTCAGCccttcatcttg                                                  >  1:3319971/1‑92 (MQ=0)
                                                   tCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGa                                              >  1:1538854/1‑101 (MQ=255)
                                                      ttGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCt                                           >  1:3327980/1‑101 (MQ=255)
                                                            aaaatcccgcGGGGGGTTGGGGGGGCGGCCCCCCCCG.GGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAAc                                      <  2:2242409/92‑1 (MQ=255)
                                                                       gggggTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..TCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGCCAGAcg                           >  1:3500114/1‑101 (MQ=255)
                                                                       gggggTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCCCAGTCCCACAGACTTTCGACCCCGAACCCGGGCCGAcg                           >  1:2880052/1‑101 (MQ=37)
                                                                         gggTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGAcgcg                         >  2:3208128/1‑101 (MQ=255)
                                                                          ggTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCCGGTCAGAcgccg                        >  1:1845/1‑99 (MQ=255)
                                                                          ggTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCCGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGCCCTCGAACCGGGTCAGACGCgg                        >  1:1341200/1‑101 (MQ=255)
                                                                          ggTGCGCGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGGCAGACGGgg                        >  2:122288/1‑101 (MQ=255)
                                                                           gTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGa                       >  2:2393342/1‑101 (MQ=255)
                                                                            tGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGag                      >  2:487477/1‑100 (MQ=255)
                                                                            tGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGAt                      >  1:738041/1‑101 (MQ=255)
                                                                               ggggggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGGTgcg                   >  2:2348761/1‑101 (MQ=255)
                                                                                tctggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgcgc                  <  2:173741/98‑1 (MQ=255)
                                                                                gggggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgcgc                  >  2:3238310/1‑101 (MQ=255)
                                                                                gggggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgcgc                  >  1:266601/1‑101 (MQ=255)
                                                                                 gggggggccccccccG.CGGGTCGCCCCGGCTT..GCCGGTGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCa                 <  1:622405/92‑1 (MQ=39)
                                                                                 ggggCGGAGCCCCCCG.CGGGTGGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCa                 >  1:2410128/1‑101 (MQ=255)
                                                                                 ggggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATgcgcc                 >  2:1454700/1‑100 (MQ=255)
                                                                                 ggggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGGCAGACGCGGATGCGCa                 >  1:2978527/1‑101 (MQ=255)
                                                                                  gggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAg                >  1:1990706/1‑101 (MQ=255)
                                                                                   ggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAgg               >  2:3759665/1‑101 (MQ=255)
                                                                                   ggCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAga               >  1:173779/1‑100 (MQ=255)
                                                                                   ggCCGAGCCCCCCG.CGGGGCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGGCCCAGACATTTCTGACCCCTAACCGGGCTGGACCCGGATGGGCCgg               >  2:4122086/1‑101 (MQ=255)
                                                                                     cGGAGCCCCCCG.CGGGTCGCCCCGGCTT..CCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCTACCTCGAACCGGGCCAGACGCGGATGCGCACGCg             >  2:1085213/1‑101 (MQ=39)
                                                                                      ggAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCgg            >  2:4169143/1‑101 (MQ=255)
                                                                                      ggAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCgg            >  1:3585038/1‑101 (MQ=255)
                                                                                       gAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCggg           >  1:3963514/1‑101 (MQ=255)
                                                                                       gAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCggg           >  2:688071/1‑101 (MQ=255)
                                                                                        aGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGc          <  1:4230058/101‑1 (MQ=255)
                                                                                         gCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGcc         >  2:2848319/1‑101 (MQ=255)
                                                                                          ccccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGGCg        >  2:3940515/1‑101 (MQ=255)
                                                                                          ccccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCg        >  2:2225828/1‑101 (MQ=255)
                                                                                          ccccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCg        >  1:1644648/1‑101 (MQ=255)
                                                                                           cccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGa       >  1:3733842/1‑101 (MQ=255)
                                                                                           cccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCCCCCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGAGGGGGCGGCGGGCCGa       >  1:4098024/1‑101 (MQ=38)
                                                                                            ccccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGaa      <  1:3759922/101‑1 (MQ=255)
                                                                                             cccG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGaaa     >  2:807859/1‑101 (MQ=255)
                                                                                                g.cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGcc  >  2:289353/1‑101 (MQ=255)
                                                                                                g.cgGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGcc  >  1:166179/1‑101 (MQ=255)
                                                                                                 |                                                                                                      
ACCCCTTCGGTTGCCGCCCCGCAACCCGTCGGCGGCCCAGCGCCCTGCTTTTCCTTGGCAAAATACG.CCGGGGGGTGCGGGGGGCGGAGCCCCCCG.CGGGTCGCCCCGGCTT..GCCGGGGCGAAATCCCCCGCTCAGTCCCAGAGACTTTCGACCTCGAACCGGGTCAGACGCGGATGCGCAGGCGGGCCGAAAGCC  >  NC_014034/2995548‑2995743

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 20 ≤  ATCG  < 30 ≤  ATCG  < 39 ≤  ATCG