Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,172,748 C→T 91.8% H526Y (CAC→TAC)  rpoB → RNA polymerase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,172,7480CT91.8% 2799.3 / 151.8 832H526Y (CAC→TAC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (41/27);  new base T (444/319);  total (485/347)
Fisher's exact test for biased strand distribution p-value = 7.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 888 total aligned reads displayed.

AGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGT  >  CP009273/4172681‑4172812
                                                                   |                                                                
agagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGc                                                                  >  1:17461609/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAAGCCGCTGTCTGAGATTACGtac                                                                >  1:53585205/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                >  1:40581754/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                >  1:52972729/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                >  1:36955234/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                >  1:35951059/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                >  1:13879006/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                >  1:16664953/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATAACGtac                                                                >  1:50262921/1‑68 (MQ=255)
     tcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACaaa                                                             <  1:8478179/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATTGACCAGATCAACCCGCTGTCTGAGATTACGTACAAAcg                                                           <  1:6258847/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                           <  1:16001955/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                           <  1:52571084/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                           <  1:47277931/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                           <  1:51852863/68‑1 (MQ=255)
         cGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcgtc                                                         <  1:14630853/68‑1 (MQ=255)
          ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcgtcg                                                        >  1:19451195/1‑68 (MQ=255)
             tCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTa                                                      <  1:6661688/67‑1 (MQ=255)
             tCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTa                                                      >  1:20818852/1‑67 (MQ=255)
               cagccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGGAtct                                                   <  1:23401520/68‑1 (MQ=255)
                  ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTcc                                                 >  1:49912128/1‑67 (MQ=255)
                  ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTcc                                                 >  1:30142722/1‑67 (MQ=255)
                   cagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtctc                                                  <  1:47954423/65‑1 (MQ=255)
                     gcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGc                                               >  1:18335191/1‑66 (MQ=255)
                     gcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGc                                               >  1:2606655/1‑66 (MQ=255)
                       tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTc                                           >  1:43751026/1‑68 (MQ=255)
                       tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTc                                           >  1:1254206/1‑68 (MQ=255)
                        gTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCg                                          <  1:50866256/68‑1 (MQ=255)
                         tctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCg                                          <  1:7778437/67‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                        <  1:5922476/68‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                        <  1:41510245/68‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                        <  1:42602455/68‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCGGGCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                         <  1:20694510/67‑1 (MQ=255)
                           tcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGcc                                       <  1:7495402/68‑1 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:24712566/1‑67 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:21210323/1‑67 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:54649840/1‑67 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:38010936/1‑67 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:38369650/1‑67 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:44840751/1‑67 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:26954996/1‑67 (MQ=255)
                               tttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAg                                    >  1:47517389/1‑67 (MQ=255)
                                 tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCg                                 >  1:15339026/1‑68 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                <  1:38719210/68‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                <  1:33636042/68‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                <  1:3757057/68‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                <  1:57653336/68‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCACGCgg                                <  1:31188052/68‑1 (MQ=255)
                                   tGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:23055791/1‑67 (MQ=255)
                                   tGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:57193351/1‑67 (MQ=255)
                                   tGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:453957/1‑67 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:56073805/1‑65 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:51543180/1‑65 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:19053186/1‑65 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:2723680/1‑65 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                >  1:46958367/1‑65 (MQ=255)
                                       ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                            >  1:9229289/1‑67 (MQ=255)
                                       ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                            >  1:32369227/1‑67 (MQ=255)
                                       ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGa                           <  1:56857862/68‑1 (MQ=255)
                                        cAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAc                          <  1:40115210/68‑1 (MQ=255)
                                        cAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAc                          <  1:16912623/68‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                         <  1:54074929/65‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                         <  1:29568323/65‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                         <  1:28871037/65‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                      >  1:3184111/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                      >  1:28092516/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                      >  1:23991073/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                      >  1:22634827/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                      >  1:20235411/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGCTCTGCCCCGt                      >  1:11391381/1‑68 (MQ=255)
                                              aacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGa                    >  1:35411142/1‑68 (MQ=255)
                                              aacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGa                    >  1:24341471/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:9651186/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:45978885/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:43347232/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:41936710/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:41645609/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:36495542/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:21779671/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                 >  1:5502777/1‑68 (MQ=255)
                                                   gCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGt                >  1:25601127/1‑67 (MQ=255)
                                                    cTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGc              >  1:39159404/1‑68 (MQ=255)
                                                        cTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAgg           <  1:2769807/67‑1 (MQ=255)
                                                        cTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGc          <  1:11882533/68‑1 (MQ=255)
                                                        cTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGc          <  1:1294645/68‑1 (MQ=255)
                                                             aTTACGTACAAACGTCGTATTTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa     >  1:10343696/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:871116/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:11045533/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:14458514/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:33195214/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:34602728/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:34896810/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:55419411/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:39723784/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:41115033/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:44332368/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa     >  1:48861870/1‑68 (MQ=255)
                                                                aCGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGt  <  1:10049789/68‑1 (MQ=255)
                                                                aCGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGt  <  1:46343057/68‑1 (MQ=255)
                                                                aCGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGt  <  1:49493386/68‑1 (MQ=255)
                                                                   |                                                                
AGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGT  >  CP009273/4172681‑4172812

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