Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,330,295 +A 90.9% 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.→A90.9% 7.0 ‑3.1 0.044intergenic (+270/+42)RCAP_rcc02161/mrcAhypothetical protein/penicillin‑binding protein 1A
Reads supporting (aligned to +/- strand):  new base (23/17):  ref base (0/4):  total (23/21)
Fisher's exact test strand distribution p-value = 4.41e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

TATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCCGGG..G....G..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCGCCTTCTTCGGGAT  >  NC_014034/2330195‑2330387
                                                                                                       |                                                                                                    
tATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..g                                                                                                       >  1:333753/1‑101 (MQ=255)
tATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGAtcccgggggtccgggg..g                                                                                                       >  1:3680493/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                                                 <  1:3532975/101‑1 (MQ=25)
                                                 cggcggCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCGGGG..GCGGAG..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAg                                                 <  2:3680265/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)
                                                                                            gggggcgg..g....c..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCGCCTTCTTCGGGAt  <  2:469013/95‑1 (MQ=255)
                                                                                                       |                                                                                                    
TATTTGCGCCAAGAAAAAGCCGGGGCGCGCATGCGGGGCCCCGGATCACCGGCGGCCCGGGCGGCTTTCACCTGGGGCCAAAGATCCCGGGGGTCCCGGG..G....G..CCCCCGGCGGGTCGCCGAGGCATCGCCT.C.GGCG.AAATCCCTTAGTAAAGACCGCCCGCGGACAGGGATCCGAAGGTCGCCTTCTTCGGGAT  >  NC_014034/2330195‑2330387

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