Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,163,133 A→C N760H (AAC→CAC)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,163,1330AC97.9% 1005.6 / 4.6 233N760H (AAC→CAC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base C (120/108);  total (124/109)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 236 total aligned reads displayed.

AGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGC  >  REL606/4163060‑4163206
                                                                         |                                                                         
aGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCAc                                                                         <  2:846561/76‑1 (MQ=255)
aGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCAc                                                                         >  2:3817007/1‑76 (MQ=255)
 gTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCAcc                                                                        <  2:765796/76‑1 (MQ=255)
  ttAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCa                                                                       >  1:2500829/1‑76 (MQ=255)
  ttAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCa                                                                       >  1:546794/1‑76 (MQ=255)
  ttAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCa                                                                       >  1:6781587/1‑76 (MQ=255)
  ttAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCa                                                                       >  1:2069650/1‑76 (MQ=255)
  ttAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCa                                                                       >  1:1291617/1‑76 (MQ=255)
   tAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAg                                                                      >  1:163911/1‑76 (MQ=255)
   tAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAg                                                                      <  1:3058280/76‑1 (MQ=255)
   tAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAg                                                                      >  1:7122658/1‑76 (MQ=255)
    aaCGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGa                                                                     >  2:2034417/1‑76 (MQ=255)
     aCGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGaa                                                                    >  1:800569/1‑76 (MQ=255)
      cGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAac                                                                   <  2:2442297/76‑1 (MQ=255)
       gAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAaca                                                                  <  2:1619737/76‑1 (MQ=255)
       gAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAaca                                                                  <  1:6867458/76‑1 (MQ=255)
       gAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAaca                                                                  >  1:4456644/1‑76 (MQ=255)
        aaGACGCGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAacac                                                                 <  2:6026984/76‑1 (MQ=255)
           aCGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTg                                                              <  1:7815189/76‑1 (MQ=255)
           aCGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTg                                                              >  2:1185261/1‑76 (MQ=255)
           aCGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTg                                                              >  1:7081074/1‑76 (MQ=255)
            cGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGt                                                             >  2:1662034/1‑76 (MQ=255)
            cGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGt                                                             <  2:7041310/76‑1 (MQ=255)
              agaTGTATGCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTAt                                                           <  2:6380826/76‑1 (MQ=255)
              agaTGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTAt                                                           <  1:1608994/76‑1 (MQ=255)
               gaTGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATc                                                          >  2:432894/1‑76 (MQ=255)
               gaTGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATc                                                          <  2:7740437/76‑1 (MQ=255)
                aTGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCa                                                         <  1:4609138/76‑1 (MQ=255)
                aTGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCa                                                         <  1:5177133/76‑1 (MQ=255)
                  gTATCCGGGTGAAGCAGGTATCGACATCTCCAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAAc                                                       >  1:6646533/1‑76 (MQ=255)
                  gTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAAc                                                       >  1:4020329/1‑76 (MQ=255)
                  gTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAAc                                                       >  1:6747219/1‑76 (MQ=255)
                  gTACCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAAc                                                       <  1:4931140/76‑1 (MQ=255)
                   tATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAAcc                                                      >  1:6543574/1‑76 (MQ=255)
                    aTCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCa                                                     <  1:2431667/76‑1 (MQ=255)
                    aTCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCa                                                     <  1:4390161/76‑1 (MQ=255)
                     tCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAg                                                    <  2:823940/76‑1 (MQ=255)
                      ccGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGa                                                   <  1:3217141/76‑1 (MQ=255)
                      ccGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGa                                                   <  1:5903386/76‑1 (MQ=255)
                       cGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGAt                                                  >  2:684148/1‑76 (MQ=255)
                         tgTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGc                                                >  2:2857039/2‑76 (MQ=255)
                           tGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCg                                              <  1:1132418/76‑1 (MQ=255)
                            gAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCgt                                             <  2:3752440/76‑1 (MQ=255)
                             aaGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCgtg                                            <  2:1510034/76‑1 (MQ=255)
                                cAGGTATCGACATCTACAACCTGACCAAATACTCCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCgtgtgt                                         <  1:4740676/76‑1 (MQ=255)
                                cAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCgtgtgt                                         >  1:7953550/1‑76 (MQ=255)
                                  ggTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCgtgtgtgt                                       <  2:6141784/76‑1 (MQ=255)
                                  ggTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCgtgtgtgt                                       >  1:6095155/1‑76 (MQ=255)
                                  ggTATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCgtgtgtgt                                       >  1:5300951/1‑76 (MQ=255)
                                    tATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtct                                     <  2:5104767/76‑1 (MQ=255)
                                    tATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtct                                     >  1:4544830/1‑76 (MQ=255)
                                    tATCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtct                                     <  2:6592290/76‑1 (MQ=255)
                                     aTCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtctc                                    >  1:1839237/1‑76 (MQ=255)
                                     aTCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGtctc                                    >  2:3279998/1‑76 (MQ=255)
                                      tCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTctct                                   >  1:7292852/1‑76 (MQ=255)
                                      tCGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTctct                                   >  2:5339621/1‑76 (MQ=255)
                                       cGACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTg                                  >  1:3512946/1‑76 (MQ=255)
                                        gACATCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTgg                                 >  1:813202/1‑76 (MQ=255)
                                          cATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGt                               >  2:263378/1‑76 (MQ=255)
                                           aTCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTg                              >  1:4160290/1‑76 (MQ=255)
                                           aTCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTg                              >  2:5032134/1‑76 (MQ=255)
                                            tCTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGa                             <  1:7220638/76‑1 (MQ=255)
                                             cTACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGaa                            >  2:4521220/1‑76 (MQ=255)
                                              tACAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAAc                           <  2:5260971/76‑1 (MQ=255)
                                               aCAACCTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAAcc                          <  1:2809483/76‑1 (MQ=255)
                                                   ccTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCGGATGCCGTGTGTGTCTCTGGGTGAACCGGtt                      <  1:3299544/76‑1 (MQ=255)
                                                   ccTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGtt                      <  1:2161034/76‑1 (MQ=255)
                                                   ccTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGtt                      >  1:3011195/1‑76 (MQ=255)
                                                    cTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTg                     >  2:6726771/1‑76 (MQ=255)
                                                    cTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTg                     <  2:2917643/76‑1 (MQ=255)
                                                    cTGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTCTGTGTCTCTGGGTGAACCGGTTg                     >  1:5970374/1‑76 (MQ=255)
                                                     tGACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGa                    >  2:3530902/1‑76 (MQ=255)
                                                      gACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGaa                   >  1:6314234/1‑76 (MQ=255)
                                                      gACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGaa                   <  2:4871750/76‑1 (MQ=255)
                                                      gACCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGCGAACCGGTTGaa                   >  2:6258502/1‑76 (MQ=255)
                                                       aCCAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAAc                  <  2:4412043/76‑1 (MQ=255)
                                                        ccAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACg                 <  2:6747219/76‑1 (MQ=255)
                                                         cAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGt                >  1:6454560/1‑76 (MQ=255)
                                                         cAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGt                <  1:1505361/76‑1 (MQ=255)
                                                         cAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGt                >  1:2624902/1‑76 (MQ=255)
                                                         cAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGt                >  1:7392323/1‑76 (MQ=255)
                                                         cAAATACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGt                >  1:1011705/1‑76 (MQ=255)
                                                          aaaTACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTg               >  1:5724443/1‑76 (MQ=255)
                                                          aaaTACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTg               >  1:3568647/1‑76 (MQ=255)
                                                            aTACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGc             <  2:1447300/76‑1 (MQ=255)
                                                            aTACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGc             >  2:2115051/1‑76 (MQ=255)
                                                             tACACCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCg            >  2:3478392/1‑76 (MQ=255)
                                                              acacCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGa           >  1:3412250/1‑76 (MQ=255)
                                                              acacCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGa           >  1:4825841/1‑76 (MQ=255)
                                                              acacCCGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGa           >  2:7074268/1‑76 (MQ=255)
                                                                 cccGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGt        <  2:7426123/76‑1 (MQ=255)
                                                                  ccGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTg       >  2:3165583/1‑76 (MQ=255)
                                                                  ccGTTCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTg       >  2:3465536/1‑76 (MQ=255)
                                                                     ttCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTGTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTg    >  2:5600816/1‑76 (MQ=255)
                                                                     ttCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTg    <  1:3765954/76‑1 (MQ=255)
                                                                     ttCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTg    >  2:2087345/1‑76 (MQ=255)
                                                                     ttCTCACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTg    <  2:291541/76‑1 (MQ=255)
                                                                      tctcACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTgg   >  1:4983733/1‑76 (MQ=255)
                                                                      tctcACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTgg   >  2:3013237/1‑76 (MQ=255)
                                                                       cTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGc  >  2:5867507/1‑76 (MQ=255)
                                                                         |                                                                         
AGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGCCGTGTGTGTCTCTGGGTGAACCGGTTGAACGTGGCGACGTGCTGGC  >  REL606/4163060‑4163206

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