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% 69.0 51intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (30/1):  ref base (14/6):  total (44/7)
Fisher's exact test strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.25e-03
*NC_0140341,253,8902.→G100.0% 69.0 51intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (30/1):  ref base (14/6):  total (44/7)
Fisher's exact test strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.38e-04
*NC_0140341,253,8903.→G100.0% 69.0 51intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (30/1):  ref base (14/6):  total (44/7)
Fisher's exact test strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.19e-04
*NC_0140341,253,8904.→G100.0% 69.0 51intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (30/1):  ref base (14/6):  total (44/7)
Fisher's exact test strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.25e-03
*NC_0140341,253,8905.→C100.0% 69.4 51intergenic (+48/+43)RCAP_rcc01181/hisHhypothetical protein/imidazole glycerol phosphate synthase subunit HisH
Reads supporting (aligned to +/- strand):  new base (30/1):  ref base (14/6):  total (44/7)
Fisher's exact test strand distribution p-value = 1.10e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.87e-03

TGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAA  >  NC_014034/1253790‑1253954
                                                                                                     |||||                                                                    
tGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg                                                                           >  2:265434/1‑101 (MQ=255)
  gTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCCTCTTGCCCTAAATACCTCGGGGGTGCg.....gg                                                                    >  2:773657/1‑101 (MQ=255)
   tGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....ggg                                                                   >  2:3794108/1‑100 (MQ=255)
   tGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCg.....ggg                                                                   >  1:2640140/1‑100 (MQ=255)
    gCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCACATCTTGCCCTAAATACCTCGGGGGTGCg.....gggg                                                                  >  1:3624506/1‑99 (MQ=255)
      tGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                            >  1:2337399/1‑101 (MQ=255)
      tGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                            >  1:2931106/1‑101 (MQ=255)
      tGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cg                                                            >  2:3186399/1‑101 (MQ=255)
       gggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGGGCG.....GGGg....cgg                                                           >  1:3080398/1‑100 (MQ=255)
        gtcAAGCTCGGCTTCCAAGCCGGTTCCGCCTTCAGGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                          <  2:2697613/99‑3 (MQ=37)
        ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggg                                                          >  2:1758546/1‑99 (MQ=255)
        ggCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cgga                                                          <  1:2679440/101‑3 (MQ=255)
          cAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagc                                                        >  2:1534348/1‑97 (MQ=255)
          cAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagc                                                        >  1:212548/1‑97 (MQ=255)
           aagggcGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                       >  2:3582048/6‑96 (MQ=38)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTTATCTTGCCCTAAATACCTCGGGGGGGCG.....GGGg....cggagcc                                                       >  2:779777/1‑96 (MQ=255)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                       >  1:4049554/1‑96 (MQ=255)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                       <  2:3216606/101‑6 (MQ=255)
           aaGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTGGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcc                                                       >  1:3634017/1‑96 (MQ=255)
             gCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcccc                                                     >  1:50021/1‑94 (MQ=255)
              cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                    <  2:1127441/101‑9 (MQ=255)
              cTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                    <  2:2159290/101‑9 (MQ=255)
              cTCGGCTCCCAGACCGGCTGCTCCTTCGAGCCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagccccc                                                    <  2:1664232/101‑9 (MQ=37)
                cGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagtccccgc                                                  >  2:3856971/1‑91 (MQ=38)
                cGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGGg....cggagcccccgc                                                  >  1:1317538/1‑91 (MQ=38)
                                           ggCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCg                                >  2:1675891/1‑101 (MQ=12)
                                              gAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGGGaaa                             >  1:964550/1‑101 (MQ=9)
                                                aaTAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATc                           >  2:398008/1‑101 (MQ=9)
                                                  tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATccc                         >  1:2539774/1‑100 (MQ=9)
                                                  tAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCct                         >  1:3904326/1‑101 (MQ=9)
                                                   aaaCCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCctc                        >  1:207399/1‑101 (MQ=9)
                                                      ccGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCa                     >  2:3899640/1‑101 (MQ=9)
                                                       cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCAc                    >  2:3007818/1‑101 (MQ=9)
                                                       cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCAc                    >  2:864274/1‑101 (MQ=9)
                                                       cGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCAc                    >  2:1050253/1‑101 (MQ=9)
                                                        gATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACg                   >  2:2903907/1‑101 (MQ=9)
                                                          tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                 >  2:241577/1‑101 (MQ=9)
                                                          tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                 >  1:4194087/1‑101 (MQ=9)
                                                          tCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGc                 >  2:3765258/1‑101 (MQ=9)
                                                             cgccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTc              >  2:1150350/1‑101 (MQ=9)
                                                              gccccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTcc             >  2:2713197/1‑101 (MQ=9)
                                                               ccccgcTCTCATCTTGCCCTAAATCCCTGGGGGGGGCGGGGGCGGAGTCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCa            <  1:3857230/101‑1 (MQ=9)
                                                                cccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTGGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAg           >  2:1099618/1‑101 (MQ=9)
                                                                 ccgcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGt          >  1:1004273/1‑101 (MQ=9)
                                                                   gcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTc        >  2:3790670/1‑101 (MQ=9)
                                                                   gcTCTCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTc        >  1:1845532/1‑101 (MQ=9)
                                                                    ctctCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCa       >  1:881059/1‑101 (MQ=9)
                                                                    ctctCATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCa       >  2:2714141/1‑101 (MQ=9)
                                                                     tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCCCGGCTTCCAGTTCAg      >  1:3260606/1‑101 (MQ=9)
                                                                     tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg      >  2:187343/1‑101 (MQ=9)
                                                                     tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg      >  2:418273/1‑101 (MQ=9)
                                                                     tctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAg      >  2:2498087/1‑101 (MQ=9)
                                                                      ctcATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCACTCACGGCTTCCAGTTCAgg     >  2:611908/1‑101 (MQ=2)
                                                                        cATCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGGCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaa   >  2:293927/1‑101 (MQ=9)
                                                                         aTCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaaa  >  2:1496405/1‑101 (MQ=9)
                                                                         aTCTTGCCCTAAATACCTCGGGGGTGCGGGGGCGGAGCCCCCGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGaaa  >  2:2224928/1‑101 (MQ=9)
                                                                                                     |||||                                                                    
TGGTGCTGGGCAAGCTCGGCTTCCAGACCGGCTGCGCCTTCGAGGCGAAATAAACCGATCCCGCCCCGCTCTCATCTTGCCCTAAATACCTCGGGGGTGCG.....GGAG....CGCCGGGTCGCCTCGGCGAAAGCCGAGGCGAAATCCTCTCACGGCTTCCAGTTCAGGAAA  >  NC_014034/1253790‑1253954

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