Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,172,748 C→T 76.9% 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,7480CT76.9% 1340.1 / inf 615H526Y (CAC→TAC) rpoBRNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (85/57);  new base T (301/171);  total (386/229)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 666 total aligned reads displayed.

AGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTC  >  CP009273/4172681‑4172814
                                                                   |                                                                  
agagTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGc                                                                    >  1:25649834/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                  >  1:10264887/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                  >  1:22360257/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                  >  1:31268724/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                  >  1:37896832/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                  >  1:40455851/1‑68 (MQ=255)
  agTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtac                                                                  >  1:8795641/1‑68 (MQ=255)
   gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGtaca                                                                 >  1:29805968/1‑68 (MQ=255)
   gTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGca                                                                   <  1:27748237/66‑1 (MQ=255)
    ttcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACaa                                                                >  1:38150079/1‑68 (MQ=255)
     tcttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACaaa                                                               <  1:14667198/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAAGAACCCGCTGTCTGAGATTACGTACAAAcg                                                             <  1:39416171/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                             <  1:41309666/68‑1 (MQ=255)
       ttcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                             <  1:1421348/68‑1 (MQ=255)
        tcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                             >  1:43567326/1‑67 (MQ=255)
        tcGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAAcg                                                             >  1:7422055/1‑67 (MQ=255)
         cGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAAcgtc                                                           <  1:9530174/68‑1 (MQ=255)
          ggTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAAcgtcg                                                          >  1:19320660/1‑68 (MQ=255)
             tCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTa                                                        >  1:18914558/1‑67 (MQ=255)
             tCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTa                                                        >  1:23546530/1‑67 (MQ=255)
             tCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTa                                                        >  1:35315978/1‑67 (MQ=255)
              ccagccagCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTAtc                                                      <  1:11311170/68‑1 (MQ=255)
                  ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTcc                                                   >  1:320180/1‑67 (MQ=255)
                  ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTcc                                                   >  1:391469/1‑67 (MQ=255)
                  ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTcc                                                   >  1:13881192/1‑67 (MQ=255)
                  ccagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTcc                                                   >  1:22803255/1‑67 (MQ=255)
                   cagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtctc                                                    <  1:36259915/65‑1 (MQ=255)
                   cagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtctc                                                    <  1:13210615/65‑1 (MQ=255)
                   cagcTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTAtctc                                                    <  1:26710721/65‑1 (MQ=255)
                       tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCAc                                               <  1:37430393/66‑1 (MQ=255)
                       tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCAc                                               <  1:36953293/66‑1 (MQ=255)
                       tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTc                                             >  1:1218166/1‑68 (MQ=255)
                       tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTc                                             >  1:2530913/1‑68 (MQ=255)
                       tGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTc                                             >  1:36478242/1‑68 (MQ=255)
                        gTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCg                                            >  1:7910560/1‑68 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                           <  1:3478974/67‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                           <  1:8382627/67‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                           <  1:40191759/67‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                           <  1:18169879/67‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                           <  1:43344145/67‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCgg                                           <  1:10943022/67‑1 (MQ=255)
                          ctcAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGc                                          <  1:27030893/68‑1 (MQ=255)
                            cAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGccc                                        >  1:4291088/1‑68 (MQ=255)
                                ttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCGGCACTCGGCCCAGGc                                    >  1:18589155/1‑68 (MQ=255)
                                ttATGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCa                                       <  1:2833990/65‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  <  1:41712755/68‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  <  1:24993600/68‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  <  1:20239938/68‑1 (MQ=255)
                                  aTGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  <  1:40919206/68‑1 (MQ=255)
                                   tGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  >  1:27391583/1‑67 (MQ=255)
                                   tGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  >  1:34388251/1‑67 (MQ=255)
                                   tGGACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGt                                 <  1:16719384/68‑1 (MQ=255)
                                    ggACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTc                                <  1:33729036/68‑1 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  >  1:17079735/1‑65 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  >  1:16049688/1‑65 (MQ=255)
                                     gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCgg                                  >  1:41985748/1‑65 (MQ=255)
                                       ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCCGGCGGTCTg                              >  1:4976/1‑67 (MQ=255)
                                       ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                              >  1:4835420/1‑67 (MQ=255)
                                       ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTg                              >  1:19975036/1‑67 (MQ=255)
                                       ccAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGa                             <  1:39856888/68‑1 (MQ=255)
                                        cAGAACAACCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAc                            <  1:6780273/68‑1 (MQ=255)
                                        cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAc                            <  1:25879054/68‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAACCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGAcc                           <  1:10061386/65‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGt                        >  1:7568651/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                        >  1:9265972/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                        >  1:5785908/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                        >  1:10519314/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                        >  1:5167/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                        >  1:6588020/1‑68 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCg                         <  1:19659478/67‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                        <  1:29596694/68‑1 (MQ=255)
                                            acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGt                        >  1:20789203/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACg                   >  1:42391799/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                   >  1:33169530/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                   >  1:29107069/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                   >  1:41089219/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                   >  1:1646881/1‑68 (MQ=255)
                                                 ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACg                   >  1:36706191/1‑68 (MQ=255)
                                                     tGTCTGAGATTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCa               >  1:27927928/1‑68 (MQ=255)
                                                        cTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGc            <  1:21188132/68‑1 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:1762997/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:21896061/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:5588474/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:6097929/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:21604739/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:24103931/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:24246343/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:26146099/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:11882900/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:30347140/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:12799361/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:15158901/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:38647991/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:33456145/1‑68 (MQ=255)
                                                             aTTACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:35019394/1‑68 (MQ=255)
                                                             aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:3177635/1‑68 (MQ=255)
                                                             aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGa       >  1:35822833/1‑68 (MQ=255)
                                                              ttACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGaa      <  1:32459230/68‑1 (MQ=255)
                                                               tACGTACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAg     >  1:15698249/1‑68 (MQ=255)
                                                                  gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTc  >  1:24656870/1‑68 (MQ=255)
                                                                   |                                                                  
AGAGTTCTTCGGTTCCAGCCAGCTGTCTCAGTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTC  >  CP009273/4172681‑4172814

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: