Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,161,948 G→A E365K (GAA→AAA)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,161,9480GA99.6% 922.0 / ‑3.2 223coding (1093/4029 nt)rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (158/64);  minor base . (0/1);  total (158/65)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

Only 100 of 263 total aligned reads displayed.

TCTGGATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGT  >  REL606/4161875‑4162020
                                                                         |                                                                        
tCTGGATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTaaaa                                                                        <  2:4688028/76‑1 (MQ=255)
tCTGGATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTaaaa                                                                        <  2:2935881/76‑1 (MQ=255)
tCTGGATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTaaaa                                                                        >  2:7262305/1‑76 (MQ=255)
   ggATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATc                                                                     >  2:69921/1‑76 (MQ=255)
   ggATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATc                                                                     >  2:2076392/1‑76 (MQ=255)
   ggATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATc                                                                     >  2:4492112/1‑76 (MQ=255)
    gATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCt                                                                    <  2:4909260/76‑1 (MQ=255)
     aTCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTa                                                                   >  2:2131764/1‑76 (MQ=255)
       cACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTAcc                                                                 >  1:1078692/1‑76 (MQ=255)
       cACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTAcc                                                                 >  2:2820061/1‑76 (MQ=255)
        aCGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCg                                                                <  1:2521745/76‑1 (MQ=255)
          ggCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCa                                                              >  1:4289042/1‑76 (MQ=255)
           gCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCat                                                             >  2:1488432/1‑76 (MQ=255)
           gCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCat                                                             >  2:4254389/1‑76 (MQ=255)
            cccATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCatg                                                            >  1:5240012/1‑76 (MQ=255)
             ccATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCatga                                                           >  1:4057065/1‑76 (MQ=255)
              cATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCatgat                                                          >  1:2527485/1‑76 (MQ=255)
               atatatCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCatgatg                                                         >  2:7095371/1‑76 (MQ=255)
               atatatCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCatgatg                                                         <  1:2866679/76‑1 (MQ=255)
               atatatCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCatgatg                                                         >  2:1070554/1‑76 (MQ=255)
                tatatCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATgc                                                        <  1:3710416/76‑1 (MQ=255)
                tatatCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATgc                                                        >  1:2435700/1‑76 (MQ=255)
                  tatCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATgcgc                                                      >  1:786533/1‑76 (MQ=255)
                   atCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGcc                                                     >  1:1548067/1‑76 (MQ=255)
                    tctcTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGccc                                                    >  2:1895393/1‑76 (MQ=255)
                    tctcTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGccc                                                    >  2:1404932/1‑76 (MQ=255)
                      tctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTg                                                  >  1:1741675/1‑76 (MQ=255)
                      tctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTg                                                  >  1:129607/1‑76 (MQ=255)
                      cctGAAACCCTACGGGGCGACCCAACCAACCGCCGTCTGAGCGCACCGGTAAAAATCTACCGCATGATGCGCCCTg                                                  <  1:701111/75‑1 (MQ=255)
                       ctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTgg                                                 >  2:3703242/1‑76 (MQ=255)
                       ctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTgg                                                 >  2:5192490/1‑76 (MQ=255)
                       ctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTgg                                                 >  1:3090839/1‑76 (MQ=255)
                       ctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTgg                                                 >  2:3313206/1‑76 (MQ=255)
                       ctGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTgg                                                 >  2:2977746/1‑76 (MQ=255)
                        tGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGc                                                <  1:4126420/76‑1 (MQ=255)
                        tGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGc                                                <  2:2383991/76‑1 (MQ=255)
                          aaaCCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGa                                              >  2:6804424/1‑76 (MQ=255)
                          aaaCCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGa                                              <  1:1469018/76‑1 (MQ=255)
                          aaaCCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGa                                              <  1:1778584/76‑1 (MQ=255)
                           aaCCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAg                                             >  1:1518915/1‑76 (MQ=255)
                            aCCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAgc                                            >  2:6025647/1‑76 (MQ=255)
                             ccTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAgcc                                           >  2:2210049/1‑76 (MQ=255)
                             ccTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAgcc                                           >  1:4911538/1‑76 (MQ=255)
                              cTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAgccg                                          >  2:6591749/1‑76 (MQ=255)
                                tACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAgccgcc                                        >  1:6428302/1‑76 (MQ=255)
                                tACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAgccgcc                                        >  2:7639436/1‑76 (MQ=255)
                                tACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAgccgcc                                        >  2:7551758/1‑76 (MQ=255)
                                 ccGTGTCGCCCCAACTACCGCCCGTCGGAGCGCACTGGAAAAAATCTCCCGCATGATGCCCCCGGGCGAGccgccg                                       <  2:6301460/75‑1 (MQ=255)
                                     gtCGACCCAACTAACGACCGTCTGAGCGCACTGGTAACAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTc                                   >  2:3109896/1‑76 (MQ=255)
                                     gtCGACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTc                                   <  1:6317875/76‑1 (MQ=255)
                                      tCGACCAAACTAACGACCGTCGGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCg                                  <  2:573819/76‑1 (MQ=255)
                                        gACCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTg                                <  1:7569032/76‑1 (MQ=255)
                                         aCCCAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGa                               <  1:1209359/76‑1 (MQ=255)
                                          cccAACAAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGaa                              <  1:357929/76‑1 (MQ=255)
                                            cAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTCCCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagc                            <  2:7041517/76‑1 (MQ=255)
                                            cAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagc                            <  1:3872709/76‑1 (MQ=255)
                                            cAACTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagc                            <  1:2022339/76‑1 (MQ=255)
                                             aaCTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagca                           <  2:666165/76‑1 (MQ=255)
                                             aaCTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagca                           >  2:5421169/1‑76 (MQ=255)
                                              aCTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagcag                          >  1:1171105/1‑76 (MQ=255)
                                              aCTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagcag                          >  1:186694/1‑76 (MQ=255)
                                              aCTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagcag                          >  1:5889529/1‑76 (MQ=255)
                                               cTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagcagc                         >  1:4861344/1‑76 (MQ=255)
                                               cTAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAagcagc                         >  2:6161944/1‑76 (MQ=255)
                                               cTAACGACCGTCTGAGAGCACTGGTAAAAATGTACCGCATGATGCGCCCTGGGGAGCCGCCGGCTCGTGAagcagc                         <  2:7804977/76‑1 (MQ=255)
                                                tAACGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCg                        >  2:1347729/1‑76 (MQ=255)
                                                 aaCGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCgg                       >  2:363210/1‑76 (MQ=255)
                                                 aaCGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCgg                       >  2:4279566/1‑76 (MQ=255)
                                                 aaCGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCgg                       >  2:6266978/1‑76 (MQ=255)
                                                 aaCGACCGTCTGAGCACACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCgg                       >  1:6762871/1‑76 (MQ=255)
                                                  aCGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGa                      <  1:5723318/76‑1 (MQ=255)
                                                   cGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGaa                     <  2:1763287/76‑1 (MQ=255)
                                                   cGACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGaa                     >  1:322064/1‑76 (MQ=255)
                                                    gACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGaaa                    >  1:55021/1‑76 (MQ=255)
                                                    gACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGaaa                    >  2:494460/1‑76 (MQ=255)
                                                    gACCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGaaa                    >  1:5825602/1‑76 (MQ=255)
                                                    gACCGTCTGAGCGCACTGGAAAAAATCTCCCGCATGAGGCCCCCGGGCGAGCCGCCGACTCGTGAAGCAGCGGaaa                    <  2:1669902/76‑1 (MQ=255)
                                                     aCCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAg                   <  1:1579796/76‑1 (MQ=255)
                                                     aCCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAg                   >  2:3933983/1‑76 (MQ=255)
                                                     aCCGTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCGGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAg                   <  1:7084340/76‑1 (MQ=255)
                                                        gccTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCGGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCt                <  2:1115460/74‑1 (MQ=255)
                                                        gTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCt                <  1:6650871/76‑1 (MQ=255)
                                                        gTCTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCt                >  2:5119283/1‑76 (MQ=255)
                                                        gTCTGAGCGCACTGGTAAAAATCTACCGCATGAGGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCt                <  2:5924325/76‑1 (MQ=255)
                                                         tCTGAGCGCACTGGTAAAAATCTCCCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTg               <  1:7262305/76‑1 (MQ=255)
                                                          cTGAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGt              >  1:979704/1‑76 (MQ=255)
                                                           tGAGCGCACTGGGAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGtt             >  2:7505294/1‑76 (MQ=255)
                                                           tGAGCCCACTGGTAAAAATCTCCCGCATGATGCCCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGtt             <  1:7346440/76‑1 (MQ=255)
                                                            gAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTc            >  1:3785487/1‑76 (MQ=255)
                                                            gAGCGCACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTc            >  2:437251/1‑76 (MQ=255)
                                                              gcgcACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCga          <  2:7243792/76‑1 (MQ=255)
                                                              gcgcACTGGTAAAAATCTACCGCATGATGCCCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCga          <  2:5542445/76‑1 (MQ=255)
                                                               cgcACTGGTAAAAATCTACCGCATGATGCCCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCgag         <  1:883196/76‑1 (MQ=255)
                                                                gcACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCgaga        <  1:6914084/76‑1 (MQ=255)
                                                                 cACTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGaa       <  1:6111401/76‑1 (MQ=255)
                                                                   cTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAAcc     >  2:2604735/1‑76 (MQ=255)
                                                                   cTGGTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAAcc     >  1:347469/1‑76 (MQ=255)
                                                                    gggTAAAAATCTCCCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCt    <  1:445446/75‑1 (MQ=255)
                                                                     ggTAAAAATCTCCCGCATGATGCGCCCGGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTg   <  2:6005155/76‑1 (MQ=255)
                                                                     ggTAAAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTg   <  2:2525786/76‑1 (MQ=255)
                                                                         |                                                                        
TCTGGATCACGGCCCATATATCTCTGAAACCTTACGTGTCGACCCAACTAACGACCGTCTGAGCGCACTGGTAGAAATCTACCGCATGATGCGCCCTGGCGAGCCGCCGACTCGTGAAGCAGCGGAAAGCCTGTTCGAGAACCTGT  >  REL606/4161875‑4162020

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