Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,930 +C 88.0% intergenic (+203/‑302) mtlK → / → RCAP_rcc02028 mannitol 2‑dehydrogenase/FeoA family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,185,9301.→C88.0% 14.4 ‑0.4 ‑0.350intergenic (+203/‑302)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (22/21):  ref base (2/4):  total (24/26)
Fisher's exact test strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.02e-01

TGCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCGCCGGGTC  >  NC_014034/2185829‑2185986
                                                                                                        |                                                        
tGCGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCg                                                           <  2:1605340/101‑1 (MQ=37)
 gCGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCCAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.c                                                         >  2:3681964/1‑101 (MQ=25)
  cGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cc                                                        >  1:1677143/1‑101 (MQ=37)
    ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACCAGAAACGCCGAA.CCC.cccc                                                      >  1:1312641/1‑101 (MQ=25)
    ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cccc                                                      <  1:2813353/101‑1 (MQ=37)
       gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCctt                                                   <  1:2590385/101‑1 (MQ=37)
       gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCctt                                                   <  1:1843519/101‑1 (MQ=37)
       gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCctt                                                   <  2:502260/101‑1 (MQ=37)
        cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCcttt                                                  <  1:1488086/101‑2 (MQ=37)
         tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCctttc                                                 >  1:2330710/6‑99 (MQ=37)
          ccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCa                                                 <  2:2213473/97‑1 (MQ=37)
          ccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCctttcc                                                >  1:1797228/5‑98 (MQ=38)
            gcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCAcc                                               <  1:4220803/99‑1 (MQ=38)
             cccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCt                                              <  2:4541455/100‑1 (MQ=38)
             cccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCt                                              <  2:3119675/100‑1 (MQ=38)
                   gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCt                                         >  1:4280145/4‑101 (MQ=39)
                   gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCt                                         >  2:3813176/4‑101 (MQ=255)
                    cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTc                                        <  1:213107/99‑1 (MQ=39)
                    cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGACCCCTTTCACCTTGGCTc                                        <  1:940393/99‑1 (MQ=39)
                       ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCaaa                                     >  2:4566712/1‑101 (MQ=255)
                        cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAat                                    <  2:3170923/101‑1 (MQ=255)
                        cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAat                                    <  1:503206/101‑1 (MQ=255)
                        cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAat                                    <  1:3778738/101‑1 (MQ=255)
                         ccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAata                                   <  2:2330548/101‑1 (MQ=255)
                         ccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTT.TTTTTGATGGCACAAGAAACGCCGAACCCC.CCCCTTTCACCTTGGCTCAAata                                   >  1:2166800/1‑101 (MQ=255)
                           ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAACTATc                                 >  1:2229733/1‑101 (MQ=255)
                           ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATc                                 <  2:2453051/101‑1 (MQ=255)
                           ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATc                                 <  1:1089342/101‑1 (MQ=255)
                             cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAATAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATccc                               <  1:2859911/101‑1 (MQ=255)
                             cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATccc                               <  2:2846664/101‑1 (MQ=255)
                             cGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATccc                               >  1:838381/1‑101 (MQ=255)
                              gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCg                              >  2:2907970/1‑101 (MQ=255)
                              gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCg                              >  1:702103/1‑101 (MQ=255)
                              gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATACCg                              >  1:3275870/1‑101 (MQ=255)
                               tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCgg                             >  2:2570277/1‑101 (MQ=255)
                                 tttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCgggg                           <  1:3197724/101‑1 (MQ=255)
                                 tttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCgggg                           >  2:1593539/1‑101 (MQ=255)
                                    aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTc                        <  1:3813461/101‑1 (MQ=255)
                                     tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                       >  2:2368700/1‑101 (MQ=255)
                                     tCCAAAGAAAAAGCAGCGGCTTGCCCCGAAGGACTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTcc                       <  1:2685057/101‑1 (MQ=38)
                                      cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCg                      <  2:1296528/101‑1 (MQ=255)
                                          aGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCggggg                  >  1:2772505/1‑101 (MQ=255)
                                              aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa              >  2:3707582/1‑101 (MQ=255)
                                                agcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGc            <  1:2119360/101‑1 (MQ=255)
                                                 gcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcc           <  1:1732596/101‑1 (MQ=255)
                                                  cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGccc          >  2:2041401/1‑101 (MQ=255)
                                                  cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGccc          >  1:828243/1‑101 (MQ=255)
                                                  cagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGccc          <  2:2675189/101‑1 (MQ=255)
                                                    gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGccccc        >  1:1375967/1‑99 (MQ=37)
                                                    gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGccccc        <  1:3682239/101‑3 (MQ=37)
                                                     cGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccg       <  1:314406/101‑4 (MQ=37)
                                                     cGGCTTCCCCCGAAGGCC.TGTTTTTGATGCCAAAAGAAACGCGGAA.CCGCCCCCTTTCACCTTGGCCCAAATTCCCCGGGGGTCCGGGGGCGGAGcccccg       <  2:1797309/101‑4 (MQ=18)
                                                      ggCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccgg      >  1:4653522/1‑97 (MQ=25)
                                                        cTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcg    >  2:1270081/1‑95 (MQ=25)
                                                         ttGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcgg   >  1:2738581/1‑94 (MQ=25)
                                                          tGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcggg  >  2:1105762/1‑93 (MQ=21)
                                                          tGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcggg  >  1:3001391/1‑93 (MQ=21)
                                                          tGCCCCGAAGGCCTT.ATTTTGATGGCACAAGAAACGCCGAA.CCGCCCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGcccccggcggg  >  1:3373897/1‑93 (MQ=21)
                                                                                                        |                                                        
TGCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAA.CCG.CCCCTTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCGCCGGGTC  >  NC_014034/2185829‑2185986

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