Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 6,001,307 A→C 28.0% S313R (AGC→CGC)  mdtK → Multidrug resistance protein MdtK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1226,001,3070AC28.0% 41.5 / 2.5 52S313R (AGC→CGC) mdtKMultidrug resistance protein MdtK
Reads supporting (aligned to +/- strand):  ref base A (23/13);  new base C (10/5);  total (33/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAACTGGCGGCGCACCAGATTGCCCTGCAGATCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGCCGCGATCTTCC  >  ACB122/6001176‑6001439
                                                                                                                                   |                                                                                                                                    
ccGAACTGGCGGCGCACCAGATGGCCCTGCAGAGCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATCCGGGGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCGGCGGGCGGCGCGCAgcg                                                                                                                                   <  8:447083/135‑1 (MQ=255)
      tGGCGGCGCACCAGATTGCCCTGCAGATCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCcgggc                                                                                                                             >  4:42989/1‑135 (MQ=255)
      gggcggcGCACCAGATGGCCCTGCAGATCGGGTCCACGGCATTCATGGTGCCGACTGGCTTGCCCTATGCGGGCACCATGCGCGGCGGCCTGTACTCCGGCGCGGGCAACCTGCGGGCGGCGCGCAGCGCcgggc                                                                                                                             <  6:422734/134‑1 (MQ=255)
           gcgcACCAGATTGCCCTGCAGATCGTGTCAACGGCATTCATGGTGCCGACTGGCTTGTCCTAAGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGGCGGGCGCAACCCcgcgcggggg                                                                                                                        >  5:423178/1‑133 (MQ=255)
            cgcACCTGTTTGCCCTGCTTAGCGTGCCCAGGGCATTCATGGGGCCTACTGGGTTGTCCTATGGGGTGACCATGCGCGGGGGCCTGTACTCCGGGGCTGGCAACCTGCGGGGGGCGCGCAGCGCCGGGCGGGTgg                                                                                                                       <  8:144825/135‑1 (MQ=255)
                            gCAGATCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCGACCTGCTGGCGGCCCGCCGCGCCGGGCGGGTGGGGATCGGGGTTTGGGc                                                                                                       >  7:340955/1‑135 (MQ=255)
                            gCAGATCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCCGGCGGGTGGGGATCGGGTTTgggg                                                                                                        >  6:445431/1‑134 (MQ=255)
                                aTCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCCGCGCCGGGCGGGTGGGGCTCGGGGTTggggggggg                                                                                                   >  7:315114/1‑130 (MQ=255)
                                aTCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCCGCGCCGGGCGGGTGGGGATCCGGTTTGGGGGCagg                                                                                                   >  3:461468/1‑133 (MQ=255)
                                     gtCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCCCACCGCCGGGCGGGTGGGGATCGGGTTTGGGGCCATGATCAt                                                                                              >  7:96398/1‑135 (MQ=255)
                                     gtCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGAGCTGGCAACCTGCTGGCCGCCCGCAGCGCCGGGCCGGTGGGCATCCGGTTTGGGGCGATGATCAt                                                                                              >  1:76935/1‑135 (MQ=255)
                                        cagggCATTCTTGGGGCCGAGTGGCTGGTCATCGGGGGTCACCAGGCGCGTCGGCCTGTACACCGGCGCGGGCCACGTGCGGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGtt                                                                                           <  8:340955/132‑1 (MQ=255)
                                                   tgggcccGACTGGCTTGGCCTATGCGGGGACCATGCGCGGCGGCCTGTCCTACGGCGCGGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCa                                                                                <  2:48027/130‑1 (MQ=255)
                                                    ggggCCGACTGGCTTGTCCTATGCGGGGACCATGCGCGGCGGCCTGTACTCCGGCGCTGGCAACCGGCGGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCAc                                                                               <  2:3949/132‑1 (MQ=255)
                                                          gACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCCGCGCCGGGCGGGTGGGCATCGGGGTTGGGGGCCTGATCATTTTTGGCTTCGCCCCGCGGtt                                                                         >  7:334543/1‑135 (MQ=255)
                                                          gACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCcgggc                                                                                                                             >  7:100458/1‑83 (MQ=255)
                                                          gACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCcgggc                                                                                                                             <  8:100458/83‑1 (MQ=255)
                                                                     ccTATGCGGGGACCCGGGGCGGCGGCCTGTACTCCGGCGCGGGCACCCGGCGGGCGGCGCGCCGCGCCGGGCGGGTGGGCATCGGGTTTGGGTCCATGATCATGTTTGCCTTCGCCACGCTGTTCCtgctgctgc                                                              <  7:215088/135‑1 (MQ=255)
                                                                     ccTAGGCGGTGACCATGCGCGGCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCctgctgctg                                                               <  4:303706/134‑1 (MQ=255)
                                                                            ggTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTTCCTTCGCCACGCTGTTCCTGCTGCTGCCCGCAGc                                                       >  3:142804/1‑135 (MQ=255)
                                                                            ggTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGc                                                       >  7:39972/1‑135 (MQ=255)
                                                                            ggTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCCCCGCGCCGGGGGGGTGGGCATCGGGTGTTGGGCCCTTAACATTTTTGCCGTCGCCACGCTGTTGCTGCTTCTTCCCGaatc                                                       >  8:95217/1‑133 (MQ=255)
                                                                              gccccAGGCGGGGCGGCCTGACCTCGGGCGGGGGCACCCGGCGGGGGGCGCGCCGGGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCAc                                                     <  2:409575/132‑1 (MQ=255)
                                                                                aCCTTGCGGGGCGGCCGGTACTCCGGGGCTGCCAACGGGGGGGGGGCGCGCCGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTg                                                   <  8:334543/135‑1 (MQ=255)
                                                                                 ccATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCCGCGCCGGGCGGGGGGGCATCGGGGTTGGGGCCATGATCATGTTTGCCTTCGCCCCGCTGTTCCTGCTGCTGCCCGAAGCACggg                                                  >  8:419648/1‑135 (MQ=255)
                                                                                 ccATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGGCTCGGGTTTTGGGCCATGGTCATGGTTTCCTTAGGCCCGCAGTTCCTGCTGCTGCCCGAAGCACTg                                                   >  1:78341/1‑134 (MQ=255)
                                                                                 ccATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTTGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTg                                                   >  6:407760/1‑134 (MQ=255)
                                                                                 ccATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCCCAGCGCCGGGCGGGGGGGCATCGGGTTTGGGGCCCTGCTCATGTTTGCCTTCGCCACGCTGTTCCTGCAGCGGCCCCCAGCACTgg                                                  >  4:106197/1‑135 (MQ=255)
                                                                                   aTGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCCGCGCCGGGCGGGTGGGCATCCGGTTTGGGGCCAAGATGATGTTTTCCTTCGCCCCCCGGTTCCAGCAGCGGCCCGAACGACTGGTa                                                >  7:355323/1‑135 (MQ=255)
                                                                                              ccTGTAGTAGGGCGCTGGAAACCTGCTGGGGGGGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCAt                                     <  4:204068/135‑1 (MQ=255)
                                                                                              ccTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCAt                                     <  8:39972/135‑1 (MQ=255)
                                                                                              ccTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGGTGCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTAAt                                     >  4:159246/1‑135 (MQ=255)
                                                                                                   actacGTGGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGAcc                                <  8:68166/135‑1 (MQ=255)
                                                                                                     tacGGCGCTGGCAACCTGCTGGCGGCGCGCCGCGCCGGGCGGGGGGGGCTCGGGGTTTGGGGCATGATCAAGTTTTGCGTCGCCCCGCGGGTGCGGCGGCTGGCCGCAGCCCAGGTTGGCCGGTTGATCGACCg                               >  2:106281/1‑134 (MQ=255)
                                                                                                          cgcTGGCAACCTGCTGGCGGCGCGCCGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGCTCCTGTTTGCCTTCGCCACGCTGGTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAAtta                         >  8:197067/1‑133 (MQ=255)
                                                                                                          cgcTGGCAACCTGCTGGCGGCGCGCCCCGCCGGGGGGGGGGGGCTCGGGTTTGGGGGCCTGATCATGTTTTGCTTGGCCACGCTGGTGCTGCTGCTGCCCGGAGCCCTGGTAGGGCTGTTCATCGCCCGGAATGa                         >  4:239041/1‑135 (MQ=255)
                                                                                                          cgcTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGGGGGTGGGCATCGGGGTTGGGGCCATGATCATGTTTGCCTTAGCCCCGCTGTTCCGGCTGCTGCCCGAAGCACTGGTAGGCCTGTTGATCGGCCGTAATGa                         >  6:215329/1‑135 (MQ=255)
                                                                                                          cgcTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGGATCCGGTTTGGGGGCATGATCATGTTTGGCTTCGCCACCCTGTTCCTGCTGCTGCCCGCAGAACAGGTAGGGCTTTTGATCGACCGGAAAGa                         >  6:76622/1‑135 (MQ=255)
                                                                                                          cgcTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATAGGGTTTGGGGGCCTGATCATGTTTGCCTTCGCCACGGGGGTCCTGCTGCTGCCCGGAGCACTGGTAGGGCAGGTGATCGCCCGTCATgc                         >  8:235685/1‑134 (MQ=255)
                                                                                                          cgcTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGGGGGGCTCGGGTTTGGGGCCATGATCATGTTTTCCTTCGCCCCCCAGTTGCTGCGGCTGCCCCCAGCACAGGTAAGCCTGTTAATCGACCCTaaata                         >  2:348881/1‑132 (MQ=255)
                                                                                                          cgcTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGGGGGGATCGGGTTTTGGGCCATGCTGATGTTTTCCTTAGCCACGCGGTTCCAGCGGCTTCCCGAATCACTGGTAGGCCTGTTCATCGCCCGGAATGa                         >  6:326313/1‑135 (MQ=255)
                                                                                                           gcTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGAc                        >  7:378871/1‑135 (MQ=255)
                                                                                                             tGGCAACCTGCTGGCGGCGCGCAGCCCCGGGCGGGGGGGGCTCGGGGTTGGGGGCAGGGTCATTTTTGCGGTAGCCCCGCTGTTGCTGCTGCTGCCCGCACCCCTGGTAGGCCGGTTGCTAGACCGAGATGaaac                      >  6:26231/1‑132 (MQ=255)
                                                                                                               gCAACCTGCTGGCGGCGCGCAGCGCCGGGGGGGTGGGCCTCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGGTGCTGCCCGAAGCACGGGTAGGCCGGGTCATCGCCCGTAATGACCCgg                    >  1:283514/1‑135 (MQ=255)
                                                                                                               gCAACCTGCTGGCGGCGCGCAGCGCCGGGCCGGTGGGCATAGGGGTTGGGGCCATGATAATGTTTGCCTTCGCCACGCGGTTGCTGCTGCGGCCCGCAGCACGGGTAGGCCTGTTCATCGACCGTAATGACCCgg                    >  6:128796/1‑135 (MQ=255)
                                                                                                                      ggggggggcgCCCCGCGCCGGGCGGGTGGGCATCGGGGTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGc             <  2:213060/129‑1 (MQ=255)
                                                                                                                      gggggcggcGCGCGGCGCCGGGCGGGTGGGCTTCGGGTTTGGGGACATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGc             <  6:423178/132‑1 (MQ=255)
                                                                                                                      gggggcggcGCGCCGCCCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGc             <  1:106281/132‑1 (MQ=255)
                                                                                                                       ctgGCGGCGGGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCGACGCTGTTGCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGcc            <  2:76935/135‑1 (MQ=255)
                                                                                                                                 gcAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGCCGCGATCtt    >  7:164101/1‑133 (MQ=255)
                                                                                                                                 gcAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGCCGCGATCTTcc  <  2:7745/135‑1 (MQ=255)
                                                                                                                                 gcAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGAACATGTTTGCCTTCGCCACCCTGTTCCTGCTGCTGCCCCCAGCACAGGTAGGCCTGTTCATCGACCGAAAAGACCCGGCCTTTGCCGCGAACtt    >  6:119033/1‑133 (MQ=255)
                                                                                                                                   |                                                                                                                                    
CCGAACTGGCGGCGCACCAGATTGCCCTGCAGATCGTGTCCACGGCATTCATGGTGCCGACTGGCTTGTCCTATGCGGTGACCATGCGCGTCGGCCTGTACTACGGCGCTGGCAACCTGCTGGCGGCGCGCAGCGCCGGGCGGGTGGGCATCGGGTTTGGGGCCATGATCATGTTTGCCTTCGCCACGCTGTTCCTGCTGCTGCCCGAAGCACTGGTAGGCCTGTTCATCGACCGTAATGACCCGGCCTTTGCCGCGATCTTCC  >  ACB122/6001176‑6001439

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