Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,162,513 C→T T553I (ACT→ATT)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,162,5130CT99.4% 682.3 / ‑4.7 169T553I (ACT→ATT) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (92/76);  total (93/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 179 total aligned reads displayed.

CGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGT  >  REL606/4162440‑4162585
                                                                         |                                                                        
cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTc                                                                        >  2:3296833/1‑76 (MQ=255)
cgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTc                                                                        <  2:2439057/76‑1 (MQ=255)
 gtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCa                                                                       <  2:2625654/76‑1 (MQ=255)
  tATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCAc                                                                      >  1:3888197/1‑76 (MQ=255)
   aTCTCCTCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACt                                                                     <  1:4969950/76‑1 (MQ=255)
   aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACCCCCGATTCACt                                                                     >  1:173158/1‑76 (MQ=255)
   aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACt                                                                     >  2:4885104/1‑76 (MQ=255)
   aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACt                                                                     >  1:3509988/1‑76 (MQ=255)
   aTCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACt                                                                     >  2:2323561/1‑76 (MQ=255)
    tctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTa                                                                    >  1:1507348/1‑76 (MQ=255)
    tctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTa                                                                    <  2:201570/76‑1 (MQ=255)
    tctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTa                                                                    <  1:3662987/76‑1 (MQ=255)
    tctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTa                                                                    >  2:3135477/1‑76 (MQ=255)
     ctcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTAc                                                                   <  2:2496935/76‑1 (MQ=255)
      tcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACg                                                                  >  2:4936255/1‑76 (MQ=255)
       ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACgg                                                                 >  2:200691/1‑76 (MQ=255)
       ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACgg                                                                 >  2:4360530/1‑76 (MQ=255)
       ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACgg                                                                 <  1:4698947/76‑1 (MQ=255)
       ccGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTAAgg                                                                 >  2:2908744/1‑76 (MQ=255)
          cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTcg                                                              <  1:397704/76‑1 (MQ=255)
          cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTcg                                                              <  1:4271596/76‑1 (MQ=255)
          cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTcg                                                              <  2:221096/76‑1 (MQ=255)
          cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTcg                                                              >  2:386823/1‑76 (MQ=255)
          cACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTcg                                                              <  1:2079802/76‑1 (MQ=255)
             tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGt                                                           >  2:748842/1‑76 (MQ=255)
               ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTAt                                                         >  2:4106994/1‑76 (MQ=255)
               ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTAt                                                         >  2:2182735/1‑76 (MQ=255)
               ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTAt                                                         <  1:3913251/76‑1 (MQ=255)
               ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTAt                                                         >  2:4643680/1‑76 (MQ=255)
               ggCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTAt                                                         <  2:3943558/76‑1 (MQ=255)
                gCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATg                                                        >  1:926431/1‑76 (MQ=255)
                gCCAAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTCCGAGACGTACACCCGATTCACTACGGTCGCGTATg                                                        <  1:293228/76‑1 (MQ=255)
                  ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTc                                                      <  2:708800/76‑1 (MQ=255)
                  ccAGGCGGTCTGACCCGCGAACGTGCAGGCTTCGGCGTTCGCGACGTACACCCGATTCACTACGCTCGCGAATGTc                                                      >  2:3524365/1‑76 (MQ=255)
                   cAGGCTGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTcc                                                     <  2:4752380/76‑1 (MQ=255)
                   cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTcc                                                     >  2:2939333/1‑76 (MQ=255)
                   cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTcc                                                     <  2:1882296/76‑1 (MQ=255)
                   cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTcc                                                     >  2:241282/1‑76 (MQ=255)
                   cAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTcc                                                     <  1:2110310/76‑1 (MQ=255)
                    aGGCGGTCTGTCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCa                                                    <  1:3402581/76‑1 (MQ=255)
                    aGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCa                                                    <  2:3066665/76‑1 (MQ=255)
                      gCGGTCTGCCCCGTGAACGTGCAGGCTTCGAGGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAAt                                                  <  1:4518246/76‑1 (MQ=255)
                      gCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAAt                                                  <  1:2052046/76‑1 (MQ=255)
                       cGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATc                                                 <  2:811130/76‑1 (MQ=255)
                       cGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATc                                                 <  2:3570719/76‑1 (MQ=255)
                        ggTCTGACCCGTGAGCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCg                                                >  1:4134038/1‑76 (MQ=255)
                        ggTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCg                                                >  2:3439359/1‑76 (MQ=255)
                         gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGa                                               >  1:1254598/1‑76 (MQ=255)
                         gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGa                                               >  2:4455706/1‑76 (MQ=255)
                           cTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGaaa                                             >  2:3179141/1‑76 (MQ=255)
                           cTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACCACGGTGGCGTATGTCCGATGGaaa                                             >  1:4412811/1‑76 (MQ=255)
                            tGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAAc                                            <  2:3836387/76‑1 (MQ=255)
                            tGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAAc                                            >  1:84490/1‑76 (MQ=255)
                             gACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAAcc                                           >  1:1812615/1‑76 (MQ=255)
                             gACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAAcc                                           >  2:4808882/1‑76 (MQ=255)
                             gACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAAcc                                           >  1:1157054/1‑76 (MQ=255)
                                ccGTGAACGTGCAGGCTTCGAAGTTTGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                        <  1:1272523/76‑1 (MQ=255)
                                ccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCt                                        <  1:1827474/76‑1 (MQ=255)
                                  gTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa                                      <  2:4200073/76‑1 (MQ=255)
                                  gTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGa                                      <  1:3315876/76‑1 (MQ=255)
                                   tGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGac                                     >  2:1251481/1‑75 (MQ=255)
                                   tGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                                     <  2:4556910/76‑1 (MQ=255)
                                   tGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGaa                                     <  1:4582608/76‑1 (MQ=255)
                                    gAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAg                                    >  1:2731235/1‑76 (MQ=255)
                                     aaCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAgg                                   >  1:4186324/1‑76 (MQ=255)
                                        gTGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTcc                                >  2:3633573/1‑76 (MQ=255)
                                         tGCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCg                               >  1:3247503/1‑76 (MQ=255)
                                          gCAGGTTTCGAAGTCCGAGCCGTCCACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGa                              <  1:2759496/76‑1 (MQ=255)
                                          gCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGa                              >  1:4632475/1‑76 (MQ=255)
                                          gCAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGa                              >  1:656226/1‑76 (MQ=255)
                                           cAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGaa                             <  1:3688427/76‑1 (MQ=255)
                                           cAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGaa                             >  2:4627804/1‑76 (MQ=255)
                                           cAGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGaa                             >  1:286550/1‑76 (MQ=255)
                                            aGGCTTCGAAGTTCGAGACGTACGCCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                            <  1:1035966/76‑1 (MQ=255)
                                            aGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                            <  1:2401405/76‑1 (MQ=255)
                                            aGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                            <  2:2539168/76‑1 (MQ=255)
                                            aGGCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                            <  1:1185407/76‑1 (MQ=255)
                                            aGGCTTCAAAGTGCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAAc                            <  1:4798099/76‑1 (MQ=255)
                                              gCTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACAt                          >  2:2183924/1‑76 (MQ=255)
                                               cTTCGAAGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATc                         >  1:2171685/1‑76 (MQ=255)
                                                    aaGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCt                    >  2:1481523/1‑76 (MQ=255)
                                                    aaGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCt                    <  1:106115/76‑1 (MQ=255)
                                                    aaGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCt                    >  2:2847811/1‑76 (MQ=255)
                                                     aGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTg                   <  2:4208988/76‑1 (MQ=255)
                                                     aGTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTg                   >  2:1579340/1‑76 (MQ=255)
                                                     aGTGCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTg                   <  1:160297/76‑1 (MQ=255)
                                                      gTTCGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGa                  >  1:1951387/1‑76 (MQ=255)
                                                         cGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa               >  1:4839070/1‑76 (MQ=255)
                                                         cGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa               <  1:3244917/76‑1 (MQ=255)
                                                         cGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa               <  2:1861052/76‑1 (MQ=255)
                                                         cGAGACGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa               <  2:2748249/76‑1 (MQ=255)
                                                         cGAGACGGACACCGGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCa               <  1:4092924/76‑1 (MQ=255)
                                                           agaCGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAAc             >  2:4316244/1‑76 (MQ=255)
                                                             aCGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActc           >  2:2131754/1‑76 (MQ=255)
                                                             aCGTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActc           >  1:1981344/1‑76 (MQ=255)
                                                               gTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctc         <  2:3740171/76‑1 (MQ=255)
                                                               gTACACCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAActctc         <  2:4930497/76‑1 (MQ=255)
                                                                 acacCCGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTg       >  1:3172219/1‑76 (MQ=255)
                                                                    cccGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc    <  1:3933278/76‑1 (MQ=255)
                                                                    cccGATTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTcc    >  2:3526111/1‑76 (MQ=255)
                                                                         |                                                                        
CGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGT  >  REL606/4162440‑4162585

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: