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,4700AC97.6% 165.0 / ‑4.0 41T539P (ACC→CCC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (20/20);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

ATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCG  >  REL606/4162398‑4162539
                                                                        |                                                                     
aTGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGcccc                                                                    >  1:1462001/1‑76 (MQ=255)
 tGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCg                                                                   >  2:1390191/1‑76 (MQ=255)
     ccAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGaa                                                               <  1:46967/76‑1 (MQ=255)
      cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAAc                                                              >  1:897898/1‑76 (MQ=255)
         aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCTGCCCAGGCGGTCTGCCCCGTGAACGTg                                                           <  2:887288/76‑1 (MQ=255)
         aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTg                                                           <  1:1118017/76‑1 (MQ=255)
            accccGCGGTGTGAGAGTACGCACAAGCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAg                                                        <  2:150282/74‑1 (MQ=255)
            aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAg                                                        >  2:264725/1‑76 (MQ=255)
            aacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAg                                                        >  1:1059559/1‑76 (MQ=255)
             acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAgg                                                       <  2:374901/76‑1 (MQ=255)
             acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAgg                                                       <  1:400746/76‑1 (MQ=255)
             acCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAgg                                                       <  1:526155/76‑1 (MQ=255)
              cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGc                                                      >  2:90169/1‑76 (MQ=255)
                 gCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTc                                                   <  1:721815/76‑1 (MQ=255)
                   tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa                                                 <  1:109740/76‑1 (MQ=255)
                   tGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGa                                                 >  1:729367/1‑76 (MQ=255)
                    gTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGaa                                                >  2:1487328/1‑76 (MQ=255)
                     tCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAg                                               <  1:1301011/76‑1 (MQ=255)
                      cTGAGATTACGCACAAACGTCGTATCTCCGCAGTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGt                                              <  1:581737/76‑1 (MQ=255)
                           aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag                                         >  2:1119491/1‑76 (MQ=255)
                           aTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCgag                                         <  1:1326337/76‑1 (MQ=255)
                             tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGAc                                       <  2:842415/76‑1 (MQ=255)
                              aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                      <  1:33469/76‑1 (MQ=255)
                              aCGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACg                                      >  1:235125/1‑76 (MQ=255)
                                gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                    <  1:428811/76‑1 (MQ=255)
                                gCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTa                                    <  1:208919/76‑1 (MQ=255)
                                 cacaAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTac                                   >  1:517055/1‑76 (MQ=255)
                                   caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTacac                                 >  2:1401496/1‑76 (MQ=255)
                                     aaCGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAccc                               <  2:480584/76‑1 (MQ=255)
                                          cgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTc                          <  1:450709/76‑1 (MQ=255)
                                          cgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTc                          >  1:1308098/1‑76 (MQ=255)
                                          cgtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTc                          >  2:62534/1‑76 (MQ=255)
                                           gtATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCa                         >  2:978925/1‑76 (MQ=255)
                                            tATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAc                        >  2:500361/1‑76 (MQ=255)
                                            tATCTCCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCAc                        >  2:902865/1‑76 (MQ=255)
                                               ctcCGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACAACCGACTCACTAc                     >  1:516502/1‑76 (MQ=255)
                                                cccGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACg                    <  1:1429339/75‑1 (MQ=255)
                                                    cACTCGGCCCAGTCGGGCTGCCCCCCAAACGAGCAGGCTTCGAAGCTCGAGACGCCCACCCCACTCACTACGGTcg                >  2:615970/1‑76 (MQ=255)
                                                          gCCCGGTCGGTCTGCCCGGTTGGCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATg          <  2:465039/76‑1 (MQ=255)
                                                             cAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTcc       <  1:1013417/76‑1 (MQ=255)
                                                               ggCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa     >  2:642550/1‑76 (MQ=255)
                                                               ggCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCaa     >  1:343843/1‑76 (MQ=255)
                                                                gCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAAt    <  1:1271051/76‑1 (MQ=255)
                                                                 cGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACTCCCGACTCACTACGGTCGCGTATGTCCAATc   <  1:214276/76‑1 (MQ=255)
                                                                  ggTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCg  >  2:148567/1‑76 (MQ=255)
                                                                        |                                                                     
ATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCG  >  REL606/4162398‑4162539

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