Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,162,470 A→C T539P (ACC→CCC)  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,4700AC96.9% 257.9 / ‑1.8 65T539P (ACC→CCC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (29/34);  total (30/35)
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

TATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162397‑4162542
                                                                         |                                                                        
tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGccc                                                                        <  2:1065663/76‑1 (MQ=255)
tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGccc                                                                        >  1:1069451/1‑76 (MQ=255)
tATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGccc                                                                        >  1:1363552/1‑76 (MQ=255)
  tGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCg                                                                      <  2:1131630/76‑1 (MQ=255)
  tGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCg                                                                      <  1:1791782/76‑1 (MQ=255)
     aCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGa                                                                   >  2:473841/1‑76 (MQ=255)
      ccAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGaa                                                                  >  1:458615/1‑76 (MQ=255)
         gAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGt                                                               <  2:1160405/76‑1 (MQ=255)
          aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTg                                                              >  2:339497/1‑76 (MQ=255)
           acaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGc                                                             <  1:77014/76‑1 (MQ=255)
            caacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCa                                                            >  2:1246604/1‑76 (MQ=255)
             aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCan                                                           <  1:472763/76‑2 (MQ=255)
             aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAg                                                           >  2:1919884/1‑76 (MQ=255)
               cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGc                                                         <  1:956506/76‑1 (MQ=255)
                ccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCt                                                        >  1:113455/1‑76 (MQ=255)
                 cGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGAACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCtt                                                       <  1:770955/76‑1 (MQ=255)
                    tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa                                                    >  2:185116/1‑76 (MQ=255)
                    tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa                                                    >  2:136414/1‑76 (MQ=255)
                     gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGaa                                                   <  2:1589251/76‑1 (MQ=255)
                      tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAg                                                  <  1:639839/76‑1 (MQ=255)
                       cTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGCGCAGGCTTCGAAGt                                                 >  1:559577/1‑76 (MQ=255)
                        tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGtt                                                <  2:206602/76‑1 (MQ=255)
                        tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGtt                                                >  1:1021588/1‑76 (MQ=255)
                         gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTc                                               >  1:1393243/1‑76 (MQ=255)
                          agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              <  2:815298/76‑1 (MQ=255)
                          agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCg                                              <  2:1945026/76‑1 (MQ=255)
                            aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag                                            <  2:997069/76‑1 (MQ=255)
                            aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag                                            <  1:1189235/76‑1 (MQ=255)
                            aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag                                            <  1:1448253/76‑1 (MQ=255)
                             ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaga                                           <  2:1759479/76‑1 (MQ=255)
                             ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaga                                           >  1:1041907/1‑76 (MQ=255)
                             ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgaga                                           <  1:488063/76‑1 (MQ=255)
                             ttACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGGTCgaga                                           >  1:601761/1‑76 (MQ=255)
                              tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc                                          >  1:1164345/1‑76 (MQ=255)
                               aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                         <  2:247888/76‑1 (MQ=255)
                               aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                         >  1:352645/1‑76 (MQ=255)
                                cGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGt                                        <  1:502227/76‑1 (MQ=255)
                                 gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                       <  2:1072322/76‑1 (MQ=255)
                                  cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac                                      >  2:487282/1‑76 (MQ=255)
                                    caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac                                    <  1:330853/76‑1 (MQ=255)
                                    caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac                                    <  1:474624/76‑1 (MQ=255)
                                     aaaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAcc                                   >  1:862877/1‑76 (MQ=255)
                                            gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                            <  1:1066264/76‑1 (MQ=255)
                                            gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                            <  2:311399/76‑1 (MQ=255)
                                            gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                            >  2:1674918/1‑76 (MQ=255)
                                             tATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAc                           >  1:1457528/1‑76 (MQ=255)
                                              aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          >  2:1965647/1‑76 (MQ=255)
                                              aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          <  2:514945/76‑1 (MQ=255)
                                              aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          >  1:1514872/1‑76 (MQ=255)
                                              aTCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          <  2:1457528/76‑1 (MQ=255)
                                              aTCTCCGCACTCGGCCCAGGCGGTCTGCCACGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACt                          >  2:764013/1‑76 (MQ=255)
                                                 tcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                       >  2:213096/1‑76 (MQ=255)
                                                 tcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                       >  2:74884/1‑76 (MQ=255)
                                                   cgcgCTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGt                     <  1:642087/76‑1 (MQ=255)
                                                     cACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcg                   <  2:1867104/76‑1 (MQ=255)
                                                      aCTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc                  >  2:993607/1‑76 (MQ=255)
                                                      aCTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTcgc                  >  1:1548414/1‑76 (MQ=255)
                                                         cGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTa               <  2:1598799/76‑1 (MQ=255)
                                                          ggCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTAt              <  2:842642/76‑1 (MQ=255)
                                                           gCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg             >  2:1266268/1‑76 (MQ=255)
                                                           gCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg             <  1:137097/76‑1 (MQ=255)
                                                                 gCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt       <  1:551511/76‑1 (MQ=255)
                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCga                                             <  1:970943/36‑1 (MQ=255)
                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCga                                             >  2:970943/1‑36 (MQ=255)
                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg     <  2:779662/76‑1 (MQ=255)
                                                                   ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg     >  1:26704/1‑76 (MQ=255)
                                                                    gTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGa    >  1:1604476/1‑76 (MQ=255)
                                                                     tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa   >  1:536939/1‑76 (MQ=255)
                                                                     tCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaa   <  2:1671039/76‑1 (MQ=255)
                                                                      cTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGaaa  <  1:878548/76‑1 (MQ=255)
                                                                         |                                                                        
TATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAA  >  REL606/4162397‑4162542

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