Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,330,289 Δ1 bp 6.8% intergenic (+264/+48) 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,2890C→. 6.8% 3.6 ‑2.3 ‑0.744intergenic (+264/+48)RCAP_rcc02161/mrcAhypothetical protein/penicillin‑binding protein 1A
Reads supporting (aligned to +/- strand):  new base (0/3):  ref base (23/18):  total (23/21)
Fisher's exact test strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.00e-01

CGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCCGGG..G....G..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCG  >  NC_014034/2330190‑2330374
                                                                                                   |                                                                                                
cGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtcc                                                                                                 >  1:1653761/1‑101 (MQ=255)
cGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtcc                                                                                                 >  1:657968/1‑101 (MQ=255)
  gCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgg                                                                                               >  1:1499983/1‑99 (MQ=255)
  gCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgg                                                                                               >  1:2181652/1‑99 (MQ=255)
   cGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGatcccggggg.cccggg                                                                                             >  1:3223310/1‑101 (MQ=255)
   cGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccggggggccggg                                                                                              >  1:2038793/1‑101 (MQ=255)
    gTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg                                                                                             >  1:2207657/1‑101 (MQ=255)
    gTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg                                                                                             >  2:3116455/1‑101 (MQ=255)
    gTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg                                                                                             >  1:5526742/1‑101 (MQ=255)
     tATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..g                                                                                          >  1:3680493/1‑101 (MQ=255)
     tATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..g                                                                                          >  1:333753/1‑101 (MQ=255)
     tATGTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..g                                                                                          <  2:5421633/101‑1 (MQ=255)
      aTTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gc                                                                                         >  1:2806085/1‑100 (MQ=255)
       tttGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcg                                                                                        <  2:1221700/101‑3 (MQ=255)
        ttGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcgg                                                                                       <  1:5424391/101‑4 (MQ=255)
         tGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcgga                                                                                      <  1:2344184/101‑5 (MQ=255)
           cgcCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcggagc                                                                                    <  2:1540063/101‑7 (MQ=255)
           cgcCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCAGCTTTCACCTGTGGCCAAAGAtcccgcgggtccgggg..gcggagc                                                                                    >  2:5460692/1‑95 (MQ=38)
             ccAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcggagccc                                                                                  <  2:1653667/101‑9 (MQ=255)
              cAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..gcggagcccc                                                                                 <  1:3681417/101‑10 (MQ=38)
              cAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccggggggccgggg..gcggagcccc                                                                                 >  2:1111262/1‑92 (MQ=255)
              cAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGT.CCGGG..G....G..cggagcccc                                                                          <  2:3702726/101‑9 (MQ=38)
               aaGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGT.CCGGG..G....G..cggagccccc                                                                         <  1:5354894/101‑10 (MQ=38)
                aGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGGCAAAGAtccctggggcccgggg..g....gggagcccccg                                                                           >  2:5087575/1‑91 (MQ=255)
                aGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGT.CCGGG..G....G..cggagcccccg                                                                        <  1:4511114/101‑11 (MQ=38)
                               gcgcATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGc                                                              <  2:711806/101‑1 (MQ=25)
                                cgcATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGGGCCGAgggc                                                             >  1:3725784/1‑99 (MQ=21)
                                 gcATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCGGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCAt                                                            <  1:2703923/101‑1 (MQ=21)
                                  cATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATc                                                           <  1:780240/101‑1 (MQ=25)
                                   aTGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGTGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCg                                                          <  2:1878473/101‑1 (MQ=25)
                                         ggCCCGGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.gg                                                  <  2:4815746/101‑1 (MQ=21)
                                          gCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGCGGCAttcgcccgcg                                                   >  2:3615164/1‑92 (MQ=18)
                                             ccGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGCGGCAGCGCCC.C.GCCgca                                              >  1:5723720/1‑99 (MQ=16)
                                               ggATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGGCGCGGGATCGCCC.C.CGCG.aaaa                                           >  2:1500891/1‑100 (MQ=16)
                                                 aTCGCCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATcc                                         <  1:4610341/101‑1 (MQ=21)
                                                     ccGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTa                                     <  1:470124/101‑1 (MQ=25)
                                                      cggcggCCGGGGCGGCTTTCACCTGGGGCCAAAGATCCGGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAg                                    <  1:4571700/101‑1 (MQ=21)
                                                      cggcggCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAg                                    <  2:3680265/101‑1 (MQ=25)
                                                      cggcggCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAg                                    <  1:3532975/101‑1 (MQ=25)
                                                       gggggcCGGGGCGGCTTCCCCCGGGGCCAAAAGATCCGGGGGGCCGGGGGGCG....GACCCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGt                                   <  2:3976790/98‑1 (MQ=12)
                                                          ggcCGGGGCGGCTTTCACCGGGGGCCAAAGACCCGGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTaaa                                <  2:2901894/101‑1 (MQ=21)
                                                          ggcCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGGGGCCGCGCCC.C.CGCG.GAAACCCCTCGTaaa                                >  2:4840662/1‑101 (MQ=9)
                                                              cGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGGGCCGGGGCGTCGCCC.C.CGCG.GAATCCCCTTGTaaaaaac                            >  1:3660852/1‑97 (MQ=9)
                                                                ggcggcTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGAccgc                          <  2:1718757/101‑1 (MQ=25)
                                                                 gcggcTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCAGCGCCC.C.CGCG.GAATCCCTTTGTAAAGAccccc                         >  2:1689011/1‑101 (MQ=12)
                                                                    gcTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGGGCCGGGGCATCGCCC.CTGGGC.AAACCCCTTTTTAAAAAccgcccc                       >  2:5756380/1‑100 (MQ=2)
                                                                    gcTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGAGGGCTCGCCC.C.GGCG.AAAACCCCTAATaaaaaacaccccc                      >  2:343808/1‑91 (MQ=9)
                                                                     cTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGGGGGGTCGCCT.C.CGCG.GAAACCCTTTTTAAAGAACGCCcccg                     >  1:1290974/1‑101 (MQ=9)
                                                                      tttCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGGGGCAAAGCCC.C.CGCG.GAAACCCCTTGTAAAAAACGCCCCCgg                    >  2:3985151/1‑101 (MQ=2)
                                                                         cACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACa                 <  1:2258589/101‑1 (MQ=21)
                                                                         cACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGGGGCATCGCCC.C.GGCG.GAAAACCTTAGTAAAAACCGCCCCCGGAcc                 >  2:1223338/1‑100 (MQ=12)
                                                                         cACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGCGGCATCGCCT.C.CGCG.AAAACCCCTCGTAAAGAACGCCCGCgggcc                 >  2:4954994/1‑98 (MQ=12)
                                                                         cACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCCGCGCGGGATCGCCTCC.GGCA.AACCCCCTTGTTAAAAACCCCCcccgcca                  >  1:2878178/1‑98 (MQ=2)
                                                                          aCCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGGGCCGGGGCAACGCCC.C.GGCG.GAATCCCCTAGTAAAAAACGCCCGCCGGCAg                >  1:5677056/1‑101 (MQ=2)
                                                                             tGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGa             <  2:750643/101‑1 (MQ=21)
                                                                             tGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGGGGGATCGCCT.C.CGCG.AAAACCCCTAGTAAAAAACGCCCCCGGGCCgggg             >  2:628769/1‑100 (MQ=2)
                                                                             tGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGCGGGAAGGCCC.C.GGCG.CAAAACCTTACTAAAGAACGCCCCCgggacgggg             >  1:469063/1‑94 (MQ=0)
                                                                              ggggCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGAt            <  1:3049425/101‑1 (MQ=21)
                                                                               gggCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGGCCCCAGGCACCGCCC.C.GGGG.AAAACCCCTAGTAAAAACAGCCCCCGGGCCggggcc           >  1:5226542/1‑98 (MQ=2)
                                                                                ggCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATcc          <  1:876749/101‑1 (MQ=21)
                                                                                 gCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGGGGCATCTCCT.C.GGCG.AAATCCCCTAGTAAAGACCCCCCCCGGGCCGGGAGCCg         >  1:745034/1‑101 (MQ=9)
                                                                                  ccAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGAGGCATCGCCC.C.CGCGCAAAACCGTAATAAAAAACCGCCC.CCGACAGGGGACCGa        >  2:5256341/1‑101 (MQ=2)
                                                                                    aaaGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGGCGGGGCATTCGCT.C.GGCG.AAAAACCTTAGTAAAGCACGCCCCCCGGCGAGGGTCCGAAg      >  2:665765/1‑101 (MQ=0)
                                                                                     aaGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGGGGGACCGCCC.C.GGCG.CAATCCCCTAGTAAAAACCCCCCGCGGGCAGGGAGCCGCAgg     >  1:3211471/1‑101 (MQ=9)
                                                                                     aaGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGGCGCCGAGGGATCGCCT.C.GGCG.AAAACCCTTAGTAAAAACAGCCCGCGGACAGGGATCCGAAgg     >  2:5014365/1‑101 (MQ=18)
                                                                                        aTCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCg  <  2:5591475/101‑1 (MQ=25)
                                                                                                   |                                                                                                
CGGCGTATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCCGGG..G....G..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCG  >  NC_014034/2330190‑2330374

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG