Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,330,175 +C 83.3% intergenic (+150/+162) RCAP_rcc02161 → / ← mrcA hypothetical protein/penicillin‑binding protein 1A

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,330,1757.→C83.3% 22.0 ‑1.7 ‑0.148intergenic (+150/+162)RCAP_rcc02161/mrcAhypothetical protein/penicillin‑binding protein 1A
Reads supporting (aligned to +/- strand):  new base (3/37):  ref base (2/6):  total (5/43)
Fisher's exact test strand distribution p-value = 1.89e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.41e-01

GCTTGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......CGGCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACC  >  NC_014034/2330077‑2330244
                                                                                                    |                                                                           
gCTTGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGgc.ggg                                                                            <  1:3054466/101‑3 (MQ=255)
   tGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggg                                                                  >  2:783732/1‑99 (MQ=255)
    ggCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggct                                                                 <  1:2023306/101‑4 (MQ=255)
     gCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctc                                                                <  2:2589968/101‑5 (MQ=255)
     gCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctc                                                                <  2:3030/101‑5 (MQ=255)
      ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctcc                                                               <  1:3210951/101‑6 (MQ=255)
      ccTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctcc                                                               <  1:1738679/101‑6 (MQ=255)
         ccGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggggcgcgc                                                            >  1:2374538/1‑93 (MQ=255)
         ccGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......Ggggctccgcc                                                            <  1:2267797/101‑9 (MQ=255)
             cAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGCGG.......GGGCGCCGCCcccc                                                         >  2:1376707/1‑101 (MQ=38)
                  gcggcAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCCGGC.G.......GGGGCGCCCCccccccaccc                                                   >  1:261050/1‑96 (MQ=255)
                               gcAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgc                                             <  1:2508570/101‑1 (MQ=18)
                                cAGAGCCCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTgcg                                            <  1:143234/101‑1 (MQ=18)
                                   agCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCa                                         <  2:137965/101‑1 (MQ=18)
                                      gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGa                                      <  1:254394/101‑1 (MQ=18)
                                      gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGa                                      <  1:791891/101‑1 (MQ=18)
                                      gccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGa                                      <  1:3550770/101‑1 (MQ=18)
                                       ccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaa                                     <  1:2771336/101‑1 (MQ=18)
                                       ccgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaa                                     <  2:3300611/101‑1 (MQ=18)
                                        cgcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaa                                    <  2:158467/101‑1 (MQ=18)
                                         gcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaaa                                   <  2:174791/101‑1 (MQ=18)
                                         gcACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTCGCGCCAAGaaaa                                   <  1:2287543/101‑1 (MQ=18)
                                          cACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGaaaaa                                  <  1:1565773/101‑1 (MQ=17)
                                           aCGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCACGGCGTATTTGCGCCAAGAAAAAg                                 <  1:1158062/101‑1 (MQ=18)
                                            cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  1:4231225/101‑1 (MQ=18)
                                            cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  2:470202/101‑1 (MQ=18)
                                            cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  1:4109795/101‑1 (MQ=18)
                                            cGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGc                                <  1:2804222/101‑1 (MQ=18)
                                             ggATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCTCCCCCGGCGTATTTGCGCCAAGAAAAAGcc                               <  2:3682415/101‑1 (MQ=17)
                                              gATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCg                              >  2:1947682/1‑101 (MQ=17)
                                               aTCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCgg                             <  2:2525313/101‑1 (MQ=18)
                                               aTCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCgg                             <  2:2950817/101‑1 (MQ=18)
                                                 cGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCgggg                           <  2:2374362/101‑1 (MQ=17)
                                                  ggCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgc                          <  2:2248284/101‑1 (MQ=17)
                                                   gCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcg                         <  2:1650494/101‑1 (MQ=17)
                                                    ccGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgc                        <  1:68240/101‑1 (MQ=17)
                                                     cGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGgcgcg                       <  2:2013409/101‑1 (MQ=17)
                                                       aTGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCCCACCCCCGGCGGATTTGCGCCAAGAAAAAGCCGGGGGGCGCa                     >  2:633029/1‑101 (MQ=12)
                                                        tGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCAt                    <  1:3119717/101‑1 (MQ=17)
                                                         ggCGGAGAGTAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                   <  2:433213/101‑1 (MQ=17)
                                                         ggCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                   <  2:1124815/101‑1 (MQ=17)
                                                         ggCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATg                   <  2:2385734/101‑1 (MQ=17)
                                                          gCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGc                  <  1:60966/101‑1 (MQ=18)
                                                           cGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCg                 <  2:1152945/101‑1 (MQ=17)
                                                           cGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGA.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGGGACAAGAAAAAGACGCGGCGCGGATGCg                 >  2:1102766/1‑101 (MQ=9)
                                                                 ggAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGccc           <  1:1376807/101‑1 (MQ=17)
                                                                  gAATTTCGCCGCGCAGGCGCGGCGCCCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGcccc          <  1:938465/101‑1 (MQ=17)
                                                                   aaTTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCg         <  2:2083351/101‑1 (MQ=17)
                                                                    aTTTCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCgg        <  2:2754206/101‑1 (MQ=17)
                                                                      ttCGCCGCGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGAt      >  1:2454994/1‑101 (MQ=17)
                                                                        cgcagcGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCCCCTCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCa    <  2:969230/101‑1 (MQ=12)
                                                                         gccgcGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAc   <  1:606938/101‑1 (MQ=18)
                                                                          ccgcGCAGGCGCGGCGACCCGGGC.GGGGGCTCCGCCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCAcc  <  2:3432700/101‑1 (MQ=18)
                                                                                                    |                                                                           
GCTTGGCCTCCGCCAAGGGCGGCAACGGCGCGCAGAGCGCCGCACGGATCGGCCGATGGCGGAGAGGAATTTCGCCGCGCAGGCGCGGCGACCCGGGC.G.......CGGCCCCGCACCCCCGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACC  >  NC_014034/2330077‑2330244

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG