Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,162,551 G→A G566S (GGT→AGT)  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,5510GA97.6% 1008.8 / 7.1 246G566S (GGT→AGT) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base A (119/121);  total (121/125)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 248 total aligned reads displayed.

ACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACT  >  REL606/4162478‑4162622
                                                                         |                                                                       
aCGTGCAGGCTTCGAAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGt                                                                       <  1:1876787/76‑1 (MQ=255)
  gTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTcc                                                                     >  2:4676381/1‑76 (MQ=255)
  gTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTcc                                                                     >  2:5433444/1‑76 (MQ=255)
    gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGa                                                                   <  2:2623077/76‑1 (MQ=255)
    gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGa                                                                   <  1:6990231/76‑1 (MQ=255)
    gCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAATCCGa                                                                   >  2:5037957/1‑76 (MQ=255)
     cAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGaa                                                                  >  2:7057449/1‑76 (MQ=255)
          ttCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGACGGTCCGAACATCg                                                             >  1:6917722/1‑76 (MQ=255)
           tCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCgg                                                            >  1:1196932/1‑76 (MQ=255)
            cGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGt                                                           <  2:3194043/76‑1 (MQ=255)
            cGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGt                                                           >  1:2736025/1‑76 (MQ=255)
            cGAAGTGCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGt                                                           <  2:4857457/76‑1 (MQ=255)
             gAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTc                                                          >  2:2502376/1‑76 (MQ=255)
              aaGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCt                                                         >  2:6186142/1‑76 (MQ=255)
              aaGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCt                                                         <  2:3862835/76‑1 (MQ=255)
              aaGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCt                                                         >  1:4495506/1‑76 (MQ=255)
              aaGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCt                                                         >  1:229431/1‑76 (MQ=255)
                 ttCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACNCNTGAAAGTCCGAACATCGGTCTGAt                                                      <  1:173432/76‑1 (MQ=255)
                 ttCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGAt                                                      <  2:4163699/76‑1 (MQ=255)
                   cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCa                                                    >  1:3025383/1‑76 (MQ=255)
                   cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCa                                                    >  1:2419446/1‑76 (MQ=255)
                    gagaCGTGCCCCCGACTCACGACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCaa                                                   <  2:2731509/76‑1 (MQ=255)
                    gagaCGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCaa                                                   <  1:647393/76‑1 (MQ=255)
                    gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCaa                                                   <  1:2810879/76‑1 (MQ=255)
                    gagaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCaa                                                   <  1:5015668/76‑1 (MQ=255)
                      gaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAAct                                                 >  2:165147/1‑76 (MQ=255)
                        cGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAActct                                               <  2:5290670/76‑1 (MQ=255)
                        cGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAActct                                               >  1:4237525/1‑76 (MQ=255)
                         gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAActctc                                              >  1:5537366/1‑76 (MQ=255)
                          tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAActctct                                             >  2:173607/1‑76 (MQ=255)
                               ccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCg                                        >  1:518067/1‑76 (MQ=255)
                               ccGACGCACGACGGTCGCGTATGTCCAATCGAAACCCCGGAAAGGCCCAACATCGGGCGGAGCAACGCTCTGTCCg                                        <  1:1303802/76‑1 (MQ=255)
                                cGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt                                       >  2:1143219/1‑76 (MQ=255)
                                cGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt                                       >  1:4292848/1‑76 (MQ=255)
                                cGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCgt                                       <  2:5678738/76‑1 (MQ=255)
                                 gACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtg                                      <  2:6834012/76‑1 (MQ=255)
                                  aCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                     <  2:4421377/76‑1 (MQ=255)
                                  aCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCgtgt                                     <  2:3905385/76‑1 (MQ=255)
                                   cTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTa                                    >  2:5631342/1‑76 (MQ=255)
                                   cTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTa                                    >  2:3417670/1‑76 (MQ=255)
                                   cTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTa                                    >  1:1503780/1‑76 (MQ=255)
                                    tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAc                                   >  1:5907647/1‑76 (MQ=255)
                                    tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAc                                   <  2:629640/76‑1 (MQ=255)
                                    tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTAc                                   >  2:3480119/1‑76 (MQ=255)
                                     cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                  <  2:195718/76‑1 (MQ=255)
                                     cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                  >  1:652124/1‑76 (MQ=255)
                                     cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                  >  1:1536384/1‑76 (MQ=255)
                                     cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACg                                  >  1:3660433/1‑76 (MQ=255)
                                      aCTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGc                                 >  2:1815971/1‑76 (MQ=255)
                                       cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                <  1:6239798/76‑1 (MQ=255)
                                       cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                >  2:6692774/1‑76 (MQ=255)
                                       cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                <  2:4020156/76‑1 (MQ=255)
                                       cTACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCCCTGTCCGTGTACGca                                <  2:4001861/76‑1 (MQ=255)
                                        tACGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcac                               >  1:3522368/1‑76 (MQ=255)
                                          cgatCGCGTATGTCCAATCGAAACCCCTGGAAGTCCGAACATCTGTCTGATCAACTCTCTGTCCGTGTACGCACAg                             <  2:2219513/73‑1 (MQ=255)
                                          cGGTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAg                             >  2:5113757/1‑76 (MQ=255)
                                           ggTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGa                            <  2:2316183/76‑1 (MQ=255)
                                            gTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                           <  2:554503/76‑1 (MQ=255)
                                            gTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                           <  1:107063/76‑1 (MQ=255)
                                            gTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                           <  1:5046354/76‑1 (MQ=255)
                                            gTCGCGTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                           <  1:2633864/76‑1 (MQ=255)
                                              cgcgTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                         >  1:2098011/1‑76 (MQ=255)
                                              cgcgTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                         <  1:806499/76‑1 (MQ=255)
                                              cgcgTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                         <  2:6531122/76‑1 (MQ=255)
                                              cgcgTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTa                         <  2:6886163/76‑1 (MQ=255)
                                                cgTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                       <  1:1386018/76‑1 (MQ=255)
                                                cgTATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                       >  2:4724047/1‑76 (MQ=255)
                                                  tATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGa                     >  2:1358114/1‑76 (MQ=255)
                                                  tATGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGa                     >  1:4825814/1‑76 (MQ=255)
                                                   aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                    <  2:4286086/76‑1 (MQ=255)
                                                   aTGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                    >  1:4107086/1‑76 (MQ=255)
                                                    tGTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAt                   >  1:1263384/1‑76 (MQ=255)
                                                     gTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                  <  1:6570290/76‑1 (MQ=255)
                                                     gTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                  <  2:3763186/76‑1 (MQ=255)
                                                     gTCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                  >  2:6787945/1‑76 (MQ=255)
                                                      tCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                 <  2:835008/76‑1 (MQ=255)
                                                      tCCAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                 <  1:195184/76‑1 (MQ=255)
                                                       ccAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                >  2:3672284/1‑76 (MQ=255)
                                                       ccAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                <  2:6096790/76‑1 (MQ=255)
                                                        cAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg               <  1:6732691/76‑1 (MQ=255)
                                                        cAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg               <  2:6856074/76‑1 (MQ=255)
                                                        cAATCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg               <  1:3760545/76‑1 (MQ=255)
                                                           tCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt            <  2:2948120/76‑1 (MQ=255)
                                                           tCGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt            <  2:1536982/76‑1 (MQ=255)
                                                            cGAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTc           >  1:3253552/1‑76 (MQ=255)
                                                             gAAACCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc          >  2:3371124/1‑76 (MQ=255)
                                                              aaaCCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTTCt         >  2:4984215/1‑76 (MQ=255)
                                                              aaaCCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt         <  2:4226790/76‑1 (MQ=255)
                                                               aaCCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt        <  2:6154468/76‑1 (MQ=255)
                                                               aaCCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt        <  2:6113203/76‑1 (MQ=255)
                                                               aaCCCCTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt        >  2:779583/1‑76 (MQ=255)
                                                                 ccccTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTTCtttc      >  1:5621077/1‑74 (MQ=255)
                                                                  cccTGAAAGGCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag     <  1:3470460/76‑1 (MQ=255)
                                                                   ccTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTTCTTgaga    >  1:6347365/1‑76 (MQ=255)
                                                                    cTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc   <  1:4343348/76‑1 (MQ=255)
                                                                    cTGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGAc   >  1:2904863/1‑76 (MQ=255)
                                                                     tGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt  <  1:5113757/76‑1 (MQ=255)
                                                                     tGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt  <  1:4167236/76‑1 (MQ=255)
                                                                     tGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt  <  1:6094685/76‑1 (MQ=255)
                                                                     tGAAAGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACt  >  1:6823537/1‑76 (MQ=255)
                                                                         |                                                                       
ACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACT  >  REL606/4162478‑4162622

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: