Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,732,721 T→C 6.5% I36T (ATC→ACC)  bcr_2 → Bicyclomycin resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,732,7210TC 6.5% 188.6 / 3.5 61I36T (ATC→ACC) bcr_2Bicyclomycin resistance protein
Reads supporting (aligned to +/- strand):  ref base T (38/19);  new base C (2/2);  total (40/21)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

CCGCCGCTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGTTCGCCTTCGCCCAATTACTGGT  >  ACB122/3732594‑3732853
                                                                                                                               |                                                                                                                                    
ccgccgCGGGGAGGCCAACGTAGGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAAc                                                                                                                               <  6:50080/135‑1 (MQ=255)
 cgccgcTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACCCGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCCATCGCTGCCATCAGCAACt                                                                                                                              >  2:136598/1‑135 (MQ=255)
 cgccgcTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACt                                                                                                                              >  1:33787/1‑135 (MQ=255)
 cgccgcTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACt                                                                                                                              >  4:324490/1‑135 (MQ=255)
 cgccgcTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCACCAGCAACt                                                                                                                              >  6:331430/1‑135 (MQ=255)
 cgccgcTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCACCAGCAACt                                                                                                                              >  6:325663/1‑135 (MQ=255)
  gccgcGGGGCGGCCAACGAATCACCCCCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACtt                                                                                                                             <  5:355999/135‑1 (MQ=255)
  gccgcGGGGAGGACAACGAATGACCCGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACtt                                                                                                                             <  1:109305/135‑1 (MQ=255)
     gcTGGGAGGCCACCTTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTc                                                                                                                          <  6:175572/135‑1 (MQ=255)
      cTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTcc                                                                                                                         >  4:119516/1‑135 (MQ=255)
         ggAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTccac                                                                                                                       >  3:321396/1‑134 (MQ=255)
           aGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTccacca                                                                                                                     >  8:341570/1‑134 (MQ=255)
                     aTGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTg                                                                                                          >  1:319204/1‑135 (MQ=255)
                                aCTCGAAACACACGCCAACGCAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTcc                                                                                               >  1:279083/1‑135 (MQ=255)
                                aCTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTcc                                                                                               >  5:124398/1‑135 (MQ=255)
                                aCTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTcc                                                                                               >  3:192939/1‑135 (MQ=255)
                                aCTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTcc                                                                                               >  1:8960/1‑135 (MQ=255)
                                  tCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCt                                                                                              >  4:210392/1‑134 (MQ=255)
                                  tCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGCGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTAGTCCACCAATATCATCCTGCCCGCTTTTCCt                                                                                              >  4:214908/1‑134 (MQ=255)
                                              ccAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCg                                                                                   >  3:370772/1‑133 (MQ=255)
                                              ccAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCg                                                                                   >  1:104776/1‑133 (MQ=255)
                                                aaCGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAg                                                                               <  5:329164/135‑1 (MQ=255)
                                                 aCGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGc                                                                              <  3:304977/135‑1 (MQ=255)
                                                 aCGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTATCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGc                                                                              <  8:316398/135‑1 (MQ=255)
                                                   gAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCt                                                                             >  1:139174/1‑134 (MQ=255)
                                                          aaGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTTGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATACTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGtc                                                                      >  7:197675/1‑134 (MQ=255)
                                                          aaGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACg                                                                        <  8:331675/132‑1 (MQ=255)
                                                          aaGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACg                                                                        >  7:331675/1‑132 (MQ=255)
                                                          aaGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGtct                                                                     >  1:79188/1‑135 (MQ=255)
                                                                   cTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATgg                                                                                       <  4:345237/108‑1 (MQ=255)
                                                                   cTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATgg                                                                                       >  3:345237/1‑108 (MQ=255)
                                                                    tCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           >  8:295490/1‑135 (MQ=255)
                                                                    tCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           >  3:250190/1‑135 (MQ=255)
                                                                    tCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           >  5:80905/1‑135 (MQ=255)
                                                                    tCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCGCGACAGCTCAACGTCTCATCCCAGGa                                                           <  5:32083/135‑1 (MQ=255)
                                                                    tCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCACCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           <  5:325663/135‑1 (MQ=255)
                                                                    tCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCACCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           <  5:331430/135‑1 (MQ=255)
                                                                          ggCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCgg                                                     >  7:4787/1‑135 (MQ=255)
                                                                             cAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCt                                                  <  4:218268/135‑1 (MQ=255)
                                                                             cAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGCAACTCGGGc                                                   >  8:225144/1‑134 (MQ=255)
                                                                                cGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGAc                                               >  2:219321/1‑135 (MQ=255)
                                                                                        cGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCa                                       >  4:28878/1‑135 (MQ=255)
                                                                                              aGGTCATCCTGTTTCTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCATGAAATGGCGCGACAGCTCAACGTTTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCtctt                                 <  5:39608/135‑1 (MQ=255)
                                                                                                gTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTc                               >  8:299389/1‑135 (MQ=255)
                                                                                                       tGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           >  6:220553/1‑100 (MQ=255)
                                                                                                       tGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           <  5:220553/100‑1 (MQ=255)
                                                                                                            cTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGttcgc                    <  1:281822/134‑1 (MQ=255)
                                                                                                              ggCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGttcgcctt                 >  8:87496/1‑135 (MQ=255)
                                                                                                                cagcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGttcgccttc                <  1:219321/134‑1 (MQ=255)
                                                                                                                cagcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGttcgccttc                <  3:119516/134‑1 (MQ=255)
                                                                                                                  gcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGt                                                                       <  7:218165/77‑1 (MQ=255)
                                                                                                                  gcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGt                                                                       >  8:218165/1‑77 (MQ=255)
                                                                                                                  gcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           >  3:342090/1‑89 (MQ=255)
                                                                                                                  gcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGa                                                           <  4:342090/89‑1 (MQ=255)
                                                                                                                  gcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGttcgccttcgc              >  7:212694/1‑134 (MQ=255)
                                                                                                                  gcaCTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTAATTACGttcgcggtggc              >  5:58894/1‑128 (MQ=255)
                                                                                                                     cTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGTTCGCCTTCGCCCaa          >  5:333025/1‑135 (MQ=255)
                                                                                                                      tCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCCGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGTTCGCCTTCGCCCaa          >  1:263619/1‑134 (MQ=255)
                                                                                                                      tCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGTTCGCCTTCGCCCaa          >  7:10995/1‑134 (MQ=255)
                                                                                                                          tGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCCTCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGTTCGCCTTCGCCCAATTACt     >  4:92373/1‑135 (MQ=255)
                                                                                                                             cATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTATTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGTTCGCCTTCGCCCAATTACTGGt  <  8:12550/135‑1 (MQ=255)
                                                                                                                               |                                                                                                                                    
CCGCCGCTGGGAGGCCAACGTATGAACAGCCGACTCGAAACACACGCCAACGAAGCTGAAGGTGATGCTCAGAGGGCCAACGTCCTGACGGCCAAGGTCATCCTGTTACTGGCAGCACTCGCTGCCATCAGCAACTTGTCCACCAATATCATCCTGCCCGCTTTTCCTGAAATGGCTCGACAGCTCAACGTCTCATCCCAGGAACTCGGGCTGACGCTCTCCAGCTTTTTCATTACGTTCGCCTTCGCCCAATTACTGGT  >  ACB122/3732594‑3732853

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