Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,330,295 +A 84.4% intergenic (+270/+42) 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,2954.→A84.4% 11.1 ‑4.0 ‑0.232intergenic (+270/+42)RCAP_rcc02161/mrcAhypothetical protein/penicillin‑binding protein 1A
Reads supporting (aligned to +/- strand):  new base (16/10):  ref base (0/5):  total (16/16)
Fisher's exact test strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.82e-01

TATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCCGGG..G....G.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCGCCTTCTTCG  >  NC_014034/2330195‑2330383
                                                                                                       |                                                                                                
tATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..g                                                                                                   >  2:2523131/1‑101 (MQ=255)
  agtgggtcAAGAAAAAGCCGGGGCCCGCATGCGGGGCCCCGGATCACCGCCGGCCCGGGCGGCCTTCACCTGGGGCCCAAGAtcccgggggtccgggg..gcg                                                                                                 >  2:3828701/8‑99 (MQ=255)
    tGCGCCAAGAAAANGCNGGGGCGCNCNNNNNGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcgga                                                                                               <  1:590716/101‑5 (MQ=38)
        ccAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcggagccc                                                                                           >  1:2137467/1‑93 (MQ=255)
         cAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcggagcccc                                                                                          >  2:3976145/1‑92 (MQ=38)
          aaGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGAGGCTTTCACCTGGGGCCAAAGATCCCGGGGGT.CCGGG..G....G.....cggagccccc                                                                               <  1:4136423/101‑10 (MQ=37)
           aGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGT.CCGGG..G....G.....cggagcccccg                                                                              <  1:2426598/101‑11 (MQ=38)
                       ggCGCGCAGGCGGGGCCCCGGATCCCCGGCGGCCCGGGCGGCTTTCCCCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGa                                                                       <  1:3371816/101‑1 (MQ=18)
                            gcATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCCCGGCAt                                                                  >  1:2255344/1‑101 (MQ=21)
                              aTGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCGAGGccacc                                                                >  2:2723659/1‑97 (MQ=21)
                                 gggggCCCCGGCCCCCGGGCGGCCCGGGCGGC.TCCCCCGGGGGGCAAAAGCCCCGGGGGGCGGGGG..G....GCGGCCCCCCCGGCGGGTCGCCGAGGCATCGCCt                                                             <  2:1905488/100‑1 (MQ=255)
                                  ggggCCCCGGATCCCCGGCGGCCCGGGCGGCTTTCCCCGGGGGCCAAAGATCCGGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTc                                                            <  1:2632290/101‑1 (MQ=18)
                                  ggggCCCCGGATCACCGGCGGCCGGGGCGGCTTTCCCCGGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTc                                                            <  1:1361394/101‑1 (MQ=18)
                                  ggggCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGGCGCGGCACCGcccc                                                            >  1:824036/1‑99 (MQ=18)
                                   gggCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGGCGGGGCAGCGccccc                                                           >  2:3098302/1‑98 (MQ=18)
                                   gggCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCGGGGCGACGCCCCg                                                           >  1:2811404/1‑101 (MQ=18)
                                     gCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGGCGCGGGCCCGCCTCggg                                                         >  2:1370207/1‑100 (MQ=16)
                                      ccccGGACCCCGGGCGGCCGGGGCGGCTTTCGCCGGGGGCCAAAGACCCGGGGGGTCGGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCg                                                        <  1:202039/101‑1 (MQ=12)
                                            aTCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGGCGCGGCCTCGCCCCGGCGGAAAcc                                                  >  2:3536347/1‑101 (MQ=12)
                                                    cggcCCGGGGGGCTTGCACGTGGGGCCAAGGTGCCGGGGGGTCGGGGG..GCGGAG.....CCCCGGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTaa                                          <  1:4313617/101‑1 (MQ=255)
                                                      gcCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGGGCCGAGGCATCGCCCCGGGGAAAAccccttttaaaa                                        >  2:4086862/1‑92 (MQ=12)
                                                         cGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGAcc                                     <  1:4092047/101‑1 (MQ=25)
                                                         cGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGGGCCCGGGGATCGCCCCGGCGACAAACCTTATTaaaaaac                                     >  2:68326/1‑97 (MQ=9)
                                                         cGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCGGGGGAGCGCCCCGGGGGAAACCCTTAGTaaaaaac                                     >  1:3467715/1‑97 (MQ=9)
                                                          gggCGGCTTTCCCCTGGGGCCAAAGACCCCGGGGGTCCGGGGGCG....GAC...CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGAccg                                    <  1:3098546/101‑1 (MQ=21)
                                                          gggCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGGGGCGAGGCGTCGCCTCCGCGGCATCCCCTTGGAAAAAccg                                    >  1:1609374/1‑101 (MQ=9)
                                                           ggcggcTTTCCCCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGAccgc                                   <  2:2971730/101‑1 (MQ=21)
                                                             cggcTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGGCGAGGCGGCGCCCCCGCGGAAACCCCTCGTACAAAccgccc                                 >  1:1623131/1‑101 (MQ=9)
                                                              gtcTTTCGCCGGGGGCCAAAGACCCGGGGGGTCGGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGAccgcccg                                <  2:2799186/99‑1 (MQ=16)
                                                              gtcTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGAccgcccg                                <  2:4542101/99‑1 (MQ=21)
                                                               gcTTTCCCCTGGGGCCAAAGATCCCGGGGGTCCGGGGGCG....GCC...CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGAccgcccgc                               <  1:3843164/101‑1 (MQ=25)
                                                                 tttCACCTGGGGCCAAAGATCCCGGGGGTCGGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCgg                             <  1:3074994/101‑1 (MQ=21)
                                                                 tttCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCGAGGGAACGCCCCGGCGAAAACCCCTTGTAAAAAACCNCCNCgg                             >  1:737228/1‑101 (MQ=9)
                                                                  ttCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCCGGGGAACGCCTCGGCGAAAACCCCTTGTAAAAACCACCCCCGGa                            >  2:1469658/1‑101 (MQ=2)
                                                                       cTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCGGACAggg                       <  2:3498500/101‑1 (MQ=21)
                                                                             ccAAAGATCCGGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGa                 <  2:2870567/101‑1 (MQ=21)
                                                                             ccAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGa                 <  2:651872/101‑1 (MQ=21)
                                                                              cAAAGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCGGGGGGTCGCCCCGGCGCAATACCCTAGTAAAAAACGCCCCCGGACCGGGAGCCGaa                >  1:199387/1‑101 (MQ=2)
                                                                               aaaGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCGAGGGATCGCCTCGGCGAAAACCCTTAGTAAAAACCCCCCGCGGGCAggggaccgcag               >  2:1128968/1‑93 (MQ=12)
                                                                                 aGATCCCGGGGGTCCGGGG..GCGGTG.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGt             <  1:4255922/101‑1 (MQ=25)
                                                                                 aGATCCCGGGGGTCCGGGG..GCGGAG.....CCCCCGGCGGGGCGCCGCGGGATCGCCTCGGCGCAAACCCTTAGTAAAAACCCCCCCCGGACAcggggccccaggg             >  2:1815489/1‑91 (MQ=2)
                                                                                        ggccgggggcgg..g....c.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCGCCTTCTTCg  <  2:1616190/91‑1 (MQ=38)
                                                                                                       |                                                                                                
TATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCCGGG..G....G.....CCCCCGGCGGGTCGCCGAGGCATCGCCTCGGCGAAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCGCCTTCTTCG  >  NC_014034/2330195‑2330383

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 23 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG