Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,798,711 A→C 18.1% N215T (AAC→ACC)  mrcA → Penicillin‑binding protein 1A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,798,7110AC18.1% 90.5 / 2.4 56N215T (AAC→ACC) mrcAPenicillin‑binding protein 1A
Reads supporting (aligned to +/- strand):  ref base A (38/6);  new base C (4/7);  total (43/13)
Fisher's exact test for biased strand distribution p-value = 1.73e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGTGCCGACGCC  >  ACB122/5798580‑5798840
                                                                                                                                   |                                                                                                                                   
cAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAAccc                                                                                                                                  <  8:375447/135‑1 (MQ=255)
 aaCCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCg                                                                                                                                 >  4:7996/1‑135 (MQ=255)
 aaCCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCg                                                                                                                                 >  1:303259/1‑135 (MQ=255)
  aCCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCCAGGCACCGTCGCGCTTCAACCCGCTGGCCACCCCgg                                                                                                                                >  2:382780/1‑135 (MQ=255)
      cgcCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTgcg                                                                                                                           >  3:228956/1‑135 (MQ=255)
        ccTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTgcgcg                                                                                                                         <  3:152451/135‑1 (MQ=255)
                aTCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAagag                                                                                                                 >  3:202934/1‑135 (MQ=255)
                    aTGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCACCCCGGTGCGCGCCAAAGAGcgcc                                                                                                            >  4:154740/1‑135 (MQ=255)
                    aTGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGcgcc                                                                                                            >  1:263116/1‑135 (MQ=255)
                       cagcagcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCTCGCTTCAACCCGCTGGCCAACCCCGTGCGCGCCAAAGAGCGCcgcg                                                                                                         >  5:194478/1‑135 (MQ=255)
                          cagcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACACGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACt                                                                                                      >  4:333743/1‑135 (MQ=255)
                            gcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTgg                                                                                                    >  1:43358/1‑135 (MQ=255)
                            gcaCAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTgg                                                                                                    >  5:79189/1‑135 (MQ=255)
                                  gtgtATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTc                                                                                              >  1:258271/1‑135 (MQ=255)
                                  gtgtATTACGGCAAGTCCATCCGTGACGTGAGACTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCCACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTc                                                                                              >  2:2894/1‑135 (MQ=255)
                                    ttATTACGGCAAGGCCATCCGGGAGGGGCGCCTGGCGCTGTGGGCGATGAGCGCGGGCCTGCCCACGGCACCGGCGCGCTTCAACCCGCTGGCCACCCCGGGGCGCGCAAAAGCGCCCCGCGACTGGATCCTCgg                                                                                            <  2:303259/134‑1 (MQ=255)
                                     tATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCCCCCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCggg                                                                                           >  4:259333/1‑135 (MQ=255)
                                     tATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCCACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCGCggg                                                                                           >  6:361641/1‑135 (MQ=255)
                                     tATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCggg                                                                                           >  2:399567/1‑135 (MQ=255)
                                          cGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCCACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGCCGGGATCCTCGGGCGCAt                                                                                      >  2:220805/1‑135 (MQ=255)
                                          cGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAGCCCCGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCAt                                                                                      >  3:430945/1‑135 (MQ=255)
                                          cGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCAt                                                                                      >  3:383548/1‑135 (MQ=255)
                                          cGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCAt                                                                                      >  7:449665/1‑135 (MQ=255)
                                          cGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCAt                                                                                      >  5:218995/1‑135 (MQ=255)
                                             cAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCCACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGCTCCTCGGGCGCATgga                                                                                   >  5:280072/1‑133 (MQ=255)
                                             cAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGt                                                                                    >  2:337900/1‑134 (MQ=255)
                                               aGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCCACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGCTCCTCGGGCGCATGTa                                                                                   >  7:288964/1‑133 (MQ=255)
                                               aGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGCTCCTCGGGCGCATGTa                                                                                   >  8:291221/1‑133 (MQ=255)
                                               aGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTa                                                                                   >  3:99273/1‑133 (MQ=255)
                                               aGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTa                                                                                   >  5:449731/1‑133 (MQ=255)
                                               aGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTa                                                                                   >  3:200548/1‑133 (MQ=255)
                                                  ccATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGCCTGGATCCTCGGGCGCATGTACAAGc                                                                              >  5:452520/1‑135 (MQ=255)
                                                  ccATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGc                                                                              >  6:316369/1‑135 (MQ=255)
                                                  ccATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACAGGAACCAACGGCGCATGTACAAGc                                                                              >  6:358530/1‑135 (MQ=255)
                                                  ccATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCAGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCAAGTACAAga                                                                              >  7:161465/1‑134 (MQ=255)
                                                    aTCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTg                                                                            >  6:62892/1‑135 (MQ=255)
                                                        gtgacgtgaGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTggg                                                                          >  1:2870/1‑133 (MQ=255)
                                                             gtgaGCCTGGCGAAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTACAACCCGCTGGCCACCCCCGAGCGCGCCAAAGACCGCCCCCACTGGATCCTCGCGCGCCTGTACAAGCTGGGCAAGAt                                                                    >  7:2041/1‑134 (MQ=255)
                                                              tgaGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCg                                                                  >  8:276117/1‑135 (MQ=255)
                                                                   cTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAg                                                             >  1:273714/1‑135 (MQ=255)
                                                                             tGGCGATCAGCGCCGGCCCGCCCCAGCCACCGGCGCGCCTCACCCCGCGGCCCCCCCCGGGCCGCGCCAAAGCGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACg                                                   <  4:86749/135‑1 (MQ=255)
                                                                                      gccccGGCCGCCCCACGCCCCCGGCGCGCCCCACCCCGCGGCCCACCCGGGGCCGCCCCAACCCCCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGccc                                          <  7:104790/132‑1 (MQ=255)
                                                                                                    aaGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTc                            >  3:173262/1‑135 (MQ=255)
                                                                                                      gccACCGGGGCGCTTCCCCCCGCGGCCCACCCCGGGGCGCCCAACCGCGCCCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCaa                          <  3:120677/133‑1 (MQ=255)
                                                                                                         cccGGCGCGCTCCCCCCCGCGGGCCACCCCGGGCCGCGCCAAAGCGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGc                       <  7:7910/134‑1 (MQ=255)
                                                                                                              cgcgCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGCCTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGCGCCGCTCACCGCCAGCt                  >  1:295207/1‑135 (MQ=255)
                                                                                                              cgcgCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGCGCCGCTCAACGCCAGCt                  >  5:433912/1‑135 (MQ=255)
                                                                                                               gcgcTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTa                 >  3:287021/1‑135 (MQ=255)
                                                                                                                     cAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGt           >  6:392102/1‑135 (MQ=255)
                                                                                                                      accccGCGGCCCACCCGGGTGCGCGCCAACGCGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGt           <  1:283241/132‑1 (MQ=255)
                                                                                                                      aaCCCGGTGCCCACCCCGGGGCGCGCCACCCCCCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGt           <  3:454344/134‑1 (MQ=255)
                                                                                                                      aaCCCGCGGCCCAACCCGGGGCGCGCCAACCCCCCCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGt           <  3:446985/134‑1 (MQ=255)
                                                                                                                       ccccGCTGCCCCCCCCGGGGCGCCCCAACGCGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGTGc         <  7:336882/134‑1 (MQ=255)
                                                                                                                        cccgggggCCCACCCGGGGCGCCCCAAACCGCGCCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGTGc         <  2:258271/129‑1 (MQ=255)
                                                                                                                          cGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGTGCCgc      >  1:154351/1‑134 (MQ=255)
                                                                                                                              gcccAACCCGGGGCGCCCAAAAGCGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGTGCCGACGcc  <  6:187524/133‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                   
CAACCGCGCCTACGGTATCGATGCAGCAGCACAGGTGTATTACGGCAAGTCCATCCGTGACGTGAGCCTGGCGCAGATGGCGATGATCGCCGGCCTGCCCAAGGCACCGTCGCGCTTCAACCCGCTGGCCAACCCGGTGCGCGCCAAAGAGCGCCGCGACTGGATCCTCGGGCGCATGTACAAGCTGGGCAAGATCGACCAGGCCAGCTACGAGGCCGCCCTGGCCGAGCCGCTCAACGCCAGCTATCACGTGCCGACGCC  >  ACB122/5798580‑5798840

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: