Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,163,089 A→T E745V (GAA→GTA)  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,0890AT93.2% 272.8 / 9.7 74E745V (GAA→GTA) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base T (34/35);  total (36/38)
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.32e-01

GTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGC  >  REL606/4163017‑4163160
                                                                        |                                                                       
gtggtgTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGc                                                                      >  1:861975/1‑76 (MQ=255)
gtggtgTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGc                                                                      <  2:1402695/76‑1 (MQ=255)
 tggtgTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCa                                                                     >  2:236908/1‑76 (MQ=255)
  ggtgTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACAAGAAGTATCCGGGTGTAGCAg                                                                    >  1:1119174/1‑76 (MQ=255)
   gtgtCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAgg                                                                   <  1:1330465/76‑1 (MQ=255)
    tgtCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGt                                                                  <  1:619240/76‑1 (MQ=255)
      tCGTTCAGTACGTGGATGCTTCGCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTAt                                                                <  1:121718/76‑1 (MQ=255)
      tCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTAt                                                                <  1:1861625/76‑1 (MQ=255)
      tCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTAt                                                                <  1:670427/76‑1 (MQ=255)
       cGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATc                                                               >  2:805282/1‑76 (MQ=255)
        gTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCg                                                              <  1:912844/76‑1 (MQ=255)
        gTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCg                                                              <  1:986388/76‑1 (MQ=255)
         ttCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGa                                                             <  2:21022/76‑1 (MQ=255)
          tCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGAc                                                            <  1:1505530/76‑1 (MQ=255)
           cAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACa                                                           >  1:564924/1‑76 (MQ=255)
           cAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACa                                                           >  2:102786/1‑76 (MQ=255)
           cAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACa                                                           >  1:1720762/1‑76 (MQ=255)
            aGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACAt                                                          <  2:148106/76‑1 (MQ=255)
            aGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACAt                                                          <  2:867917/76‑1 (MQ=255)
            aGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACAt                                                          <  1:215255/76‑1 (MQ=255)
              tACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCt                                                        <  1:210433/76‑1 (MQ=255)
              tACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCt                                                        <  1:1467242/76‑1 (MQ=255)
               aCGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTa                                                       >  1:763998/1‑76 (MQ=255)
               aCGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTa                                                       <  2:1184547/76‑1 (MQ=255)
                 gTGGATGCTTCCCGTATCGTTATCAAAGTTAACGACGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACa                                                     >  2:191425/1‑76 (MQ=255)
                 gTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACa                                                     >  2:1586510/1‑76 (MQ=255)
                 gTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACa                                                     >  1:1525435/1‑76 (MQ=255)
                  tGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACaa                                                    >  1:1661232/1‑76 (MQ=255)
                    gATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAAcc                                                  >  2:879137/1‑76 (MQ=255)
                    gATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAAcc                                                  >  1:1043668/1‑76 (MQ=255)
                        cTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGAc                                              >  2:1667912/1‑76 (MQ=255)
                          tCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCa                                            >  2:98855/1‑76 (MQ=255)
                            ccGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCaaa                                          <  2:455564/76‑1 (MQ=255)
                            ccGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCAGGTGTAGCAGGTATCGACATCTACAACCTGACCaaa                                          <  2:455548/76‑1 (MQ=255)
                             cGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAAt                                         >  1:104544/1‑76 (MQ=255)
                             cGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAAt                                         <  2:460944/76‑1 (MQ=255)
                               tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATac                                       >  1:582380/1‑76 (MQ=255)
                               tATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATac                                       >  1:47643/1‑76 (MQ=255)
                                aTCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATaca                                      <  1:596579/76‑1 (MQ=255)
                                 tCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATacac                                     <  1:706537/76‑1 (MQ=255)
                                  cGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACAcc                                    <  1:399695/76‑1 (MQ=255)
                                   gTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACAccc                                   <  1:677794/76‑1 (MQ=255)
                                    ttATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCg                                  >  2:222993/1‑76 (MQ=255)
                                    ttATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCg                                  >  1:1070700/1‑76 (MQ=255)
                                    ttATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCg                                  >  2:286825/1‑76 (MQ=255)
                                    ttATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCg                                  >  1:285434/1‑76 (MQ=255)
                                     tATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGt                                 >  1:158833/1‑76 (MQ=255)
                                     tATCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGt                                 >  1:1663799/1‑76 (MQ=255)
                                      aTCAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGtt                                <  2:1084357/76‑1 (MQ=255)
                                        cAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCt                              <  1:928641/76‑1 (MQ=255)
                                        cAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCt                              <  2:704022/76‑1 (MQ=255)
                                        cAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCt                              <  1:1213603/76‑1 (MQ=255)
                                        cAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCt                              >  2:1403975/1‑76 (MQ=255)
                                        cAAAGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCt                              >  1:562292/1‑76 (MQ=255)
                                         aaaGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTa                             >  1:1872712/1‑76 (MQ=255)
                                          aaGTTAACGAAGACGATATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTaa                            <  1:48076/76‑1 (MQ=255)
                                          aaGTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAACTACACCCGTTCTaa                            >  1:1655636/1‑76 (MQ=255)
                                            gTTAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAAcc                          <  1:57266/76‑1 (MQ=255)
                                             ttAACGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCa                         >  1:714743/1‑76 (MQ=255)
                                                aCGAAGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGaa                      >  1:1629447/1‑76 (MQ=255)
                                                  gAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAaca                    <  1:1175022/76‑1 (MQ=255)
                                                   aaGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGATCTAACCAGAacac                   >  1:773460/1‑76 (MQ=255)
                                                    aGACGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACAcc                  >  2:745024/1‑76 (MQ=255)
                                                      aCGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTg                <  2:493421/76‑1 (MQ=255)
                                                      aCGAGATGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTg                <  2:1019946/76‑1 (MQ=255)
                                                        gagaTGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTa              <  1:131566/76‑1 (MQ=255)
                                                         agaTGTAGCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTAt             <  1:16965/76‑1 (MQ=255)
                                                          gaTGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATc            <  1:1724968/76‑1 (MQ=255)
                                                            tGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCaa          <  2:707725/76‑1 (MQ=255)
                                                            tGTATCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCaa          >  1:1031712/1‑76 (MQ=255)
                                                               aTCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCa       >  1:911562/1‑76 (MQ=255)
                                                               aTCCGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCa       >  2:465052/1‑76 (MQ=255)
                                                                  cGGGTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGAt    <  1:396641/76‑1 (MQ=255)
                                                                   gggTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATg   <  2:762804/76‑1 (MQ=255)
                                                                   gggTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATg   <  2:303953/76‑1 (MQ=255)
                                                                    ggTGTAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGc  >  1:222645/1‑76 (MQ=255)
                                                                        |                                                                       
GTGGTGTCGTTCAGTACGTGGATGCTTCCCGTATCGTTATCAAAGTTAACGAAGACGAGATGTATCCGGGTGAAGCAGGTATCGACATCTACAACCTGACCAAATACACCCGTTCTAACCAGAACACCTGTATCAACCAGATGC  >  REL606/4163017‑4163160

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