Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,253,890 +GGGGC 100% intergenic (+48/+43) RCAP_rcc01181 → / ← hisH hypothetical protein/imidazole glycerol phosphate synthase subunit HisH

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,253,8901.→G100.0% 71.1 55intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (35/0):  ref base (8/12):  total (43/12)
Fisher's exact test strand distribution p-value = 2.87e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.88e-01
*NC_0140341,253,8902.→G100.0% 71.1 55intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (35/0):  ref base (8/12):  total (43/12)
Fisher's exact test strand distribution p-value = 2.87e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00
*NC_0140341,253,8903.→G100.0% 71.1 55intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (35/0):  ref base (8/12):  total (43/12)
Fisher's exact test strand distribution p-value = 2.87e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.99e-01
*NC_0140341,253,8904.→G100.0% 71.1 55intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (35/0):  ref base (8/12):  total (43/12)
Fisher's exact test strand distribution p-value = 2.87e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00
*NC_0140341,253,8905.→C100.0% 70.9 55intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (34/0):  ref base (8/12):  total (43/12)
Fisher's exact test strand distribution p-value = 3.67e-07
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAT  >  NC_014034/1253790‑1253955
                                                                                                     |||||                                                                     
tGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg                                                                            >  2:3133959/1‑101 (MQ=255)
tGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg                                                                            >  2:3120174/1‑101 (MQ=255)
 ggTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....g                                                                      <  2:3211487/101‑1 (MQ=255)
 ggTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....g                                                                      <  2:2149643/101‑1 (MQ=255)
  gTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gg                                                                     <  2:1009168/101‑1 (MQ=255)
  gTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gg                                                                     <  1:1477081/101‑1 (MQ=255)
    gCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  1:560487/101‑3 (MQ=255)
    gCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  2:3357661/101‑3 (MQ=255)
    gCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  1:2491560/101‑3 (MQ=255)
    gCTGGGCAAGCCCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                   <  1:4101229/101‑3 (MQ=255)
     cTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....c                                                              <  1:1894408/101‑1 (MQ=255)
     cTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....c                                                              <  1:4472112/101‑1 (MQ=255)
      tGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                             >  2:826320/1‑101 (MQ=255)
      cgggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                             <  2:3855967/100‑1 (MQ=255)
       gggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgg                                                            >  1:3535510/1‑100 (MQ=255)
        ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                           >  2:1550662/1‑99 (MQ=255)
        ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                           <  2:4385871/101‑3 (MQ=255)
        ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                           <  1:3904558/101‑3 (MQ=255)
         gCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggag                                                          <  2:391647/101‑4 (MQ=255)
          taaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagc                                                         <  1:2121340/100‑5 (MQ=255)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                        <  2:4255640/101‑6 (MQ=255)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                        >  2:3299303/1‑96 (MQ=255)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                        <  2:2834656/101‑6 (MQ=255)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagca                                                        >  2:2121190/1‑96 (MQ=255)
            aGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccc                                                       >  2:1028516/1‑95 (MQ=255)
            aGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccc                                                       <  2:2022593/101‑7 (MQ=255)
              cTCGGCTTCCAGCCCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                     <  2:3792155/101‑9 (MQ=38)
              cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                     >  1:2540634/1‑93 (MQ=255)
              cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                     <  2:4254239/101‑9 (MQ=255)
               tCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcccccg                                                    >  2:3379154/1‑92 (MQ=38)
                   cTTCCAGCCCGGCTGCGCCTTGGAGGCGAAATAACCCGATCCCGCCCCGCTCCCATCTTCCCCAAAATCCCCCGGGGGTGCG.....GGGG....CGGCgcccccgccgg                                                <  2:2902753/101‑11 (MQ=19)
                                           ggCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCg                                 >  1:783322/1‑101 (MQ=12)
                                              gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                              >  2:3503396/1‑101 (MQ=12)
                                              gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                              >  1:2234569/1‑101 (MQ=9)
                                              gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGaaa                              >  2:2113155/1‑101 (MQ=12)
                                              gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGGGCGGGGGCGGAGCCCCCGCCGGGTAGCCTCGGCGAAAGCCGAGGCGaaa                              >  2:2855155/1‑101 (MQ=9)
                                               aaaTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAAt                             >  1:3902527/1‑101 (MQ=9)
                                                 aTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATcc                           >  2:3645734/1‑101 (MQ=9)
                                                  tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                          >  1:4383723/1‑101 (MQ=9)
                                                  tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                          >  1:4270036/1‑101 (MQ=9)
                                                  tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                          >  1:3380305/1‑101 (MQ=9)
                                                    aaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctct                        >  2:2745360/1‑101 (MQ=9)
                                                      ccGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCa                      >  2:128056/1‑101 (MQ=9)
                                                       cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCAc                     >  2:2252298/1‑101 (MQ=9)
                                                       cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGAGGAGCCCCCGCCGGGTCGCCGCGGCGAAAGCAGCGGCGAAATCCTCTCAc                     >  1:1160051/1‑101 (MQ=12)
                                                        gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACg                    >  1:2968515/1‑101 (MQ=9)
                                                        gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGGCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCg                    >  2:2784583/1‑101 (MQ=2)
                                                         aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                   >  2:2940972/1‑101 (MQ=9)
                                                         aTCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACgg                   >  2:3473969/1‑101 (MQ=9)
                                                          tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                  >  1:4090564/1‑101 (MQ=9)
                                                          tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                  >  1:776976/1‑101 (MQ=9)
                                                           cccgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCt                 >  2:2918898/1‑101 (MQ=9)
                                                              gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc              >  2:621781/1‑101 (MQ=9)
                                                              gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc              >  2:1351082/1‑101 (MQ=9)
                                                               ccccgcTCTCATCTTCCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCGGCTTCCa             >  1:1748467/1‑101 (MQ=9)
                                                               ccacgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAGGCCGAGGCGTAATCCTCTCACGGCTTccc             >  1:2658816/4‑100 (MQ=2)
                                                                cccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAg            >  1:13246/1‑101 (MQ=9)
                                                                  cgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGtt          >  2:4562557/1‑101 (MQ=9)
                                                                  cgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGtt          >  1:589026/1‑101 (MQ=9)
                                                                   gcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTc         >  1:1072279/1‑101 (MQ=9)
                                                                    ctctCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCGGCTTCCAGTTCa        >  2:1374257/1‑101 (MQ=9)
                                                                     tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg       >  2:3098563/1‑101 (MQ=9)
                                                                     tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg       >  1:3239096/1‑101 (MQ=9)
                                                                       tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa     >  1:1852214/1‑101 (MQ=9)
                                                                       tcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGa     >  2:1546584/1‑101 (MQ=9)
                                                                         aTCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaaa   >  2:2896212/1‑101 (MQ=9)
                                                                          tCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAt  >  2:2188324/1‑101 (MQ=9)
                                                                                                     |||||                                                                     
TGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAAT  >  NC_014034/1253790‑1253955

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG