Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,162,569 A→T I572F (ATC→TTC)  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,5690AT98.7% 294.3 / ‑3.7 74I572F (ATC→TTC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (35/38);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  REL606/4162497‑4162641
                                                                        |                                                                        
cGAGACGTACACCCGACTCACTATGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCa                                                                       >  2:688751/1‑76 (MQ=255)
cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCa                                                                       >  1:682146/1‑76 (MQ=255)
cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCa                                                                       >  2:225638/1‑76 (MQ=255)
cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCa                                                                       >  2:1161305/1‑76 (MQ=255)
  agaCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAAc                                                                     <  2:778928/76‑1 (MQ=255)
       tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAActctct                                                                <  2:349393/76‑1 (MQ=255)
        acccccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTg                                                               <  2:1632125/73‑1 (MQ=255)
        acacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTg                                                               <  2:1690448/76‑1 (MQ=255)
         cccccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGt                                                              <  2:657926/74‑1 (MQ=255)
         cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGCACATCGGGCTGTTCAACTCTCTGt                                                              >  1:968339/1‑76 (MQ=255)
         cacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGt                                                              <  1:173247/76‑1 (MQ=255)
            ccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCg                                                           <  2:1918022/76‑1 (MQ=255)
               aCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCgtgt                                                        <  2:778927/76‑1 (MQ=255)
               aCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCgtgt                                                        <  2:876531/76‑1 (MQ=255)
                cTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTa                                                       <  2:842658/76‑1 (MQ=255)
                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACg                                                     >  2:1204204/1‑76 (MQ=255)
                  cACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACg                                                     <  1:542811/76‑1 (MQ=255)
                      aCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGcaca                                                 >  2:469205/1‑76 (MQ=255)
                      aCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGcaca                                                 <  1:1181281/76‑1 (MQ=255)
                       cGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAg                                                <  1:260619/76‑1 (MQ=255)
                       cGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAg                                                <  1:69483/76‑1 (MQ=255)
                        ggTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGa                                               >  2:962088/1‑76 (MQ=255)
                         gTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGAc                                              <  2:696018/76‑1 (MQ=255)
                          tCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACt                                             >  2:247176/1‑76 (MQ=255)
                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGCCTa                                            <  2:1121226/76‑1 (MQ=255)
                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTa                                            >  1:323404/1‑76 (MQ=255)
                            gcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTaa                                           <  2:1875481/76‑1 (MQ=255)
                            gcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTaa                                           <  1:74970/76‑1 (MQ=255)
                              gTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACg                                         >  1:1054739/1‑76 (MQ=255)
                               tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                        >  2:1338906/1‑76 (MQ=255)
                               tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                        >  2:1519865/1‑76 (MQ=255)
                               tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                        <  1:1526975/76‑1 (MQ=255)
                               tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                        >  1:1012789/1‑76 (MQ=255)
                                aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                       >  2:997391/1‑76 (MQ=255)
                                  gTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                     >  2:178747/1‑76 (MQ=255)
                                  gTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                     <  1:513744/76‑1 (MQ=255)
                                   tCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATAc                                    >  2:34406/1‑76 (MQ=255)
                                    ctaaTCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                   >  2:494841/3‑76 (MQ=255)
                                     cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                  <  2:1196839/76‑1 (MQ=255)
                                      aaTCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                                 <  2:294615/76‑1 (MQ=255)
                                       aTCGAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCt                                <  2:1604464/76‑1 (MQ=255)
                                          gAAACCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTcc                             >  2:1097730/1‑76 (MQ=255)
                                            aaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           >  1:1310502/1‑76 (MQ=255)
                                            aaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           <  2:876062/76‑1 (MQ=255)
                                            aaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           <  2:1062955/76‑1 (MQ=255)
                                            aaCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           <  1:340779/76‑1 (MQ=255)
                                             aCCCCTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTg                          <  2:450911/76‑1 (MQ=255)
                                               cccTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag                        <  1:1451833/76‑1 (MQ=255)
                                               cccTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag                        <  2:624781/76‑1 (MQ=255)
                                               cccTGAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgag                        <  1:214618/76‑1 (MQ=255)
                                                   gAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                    >  2:866837/1‑76 (MQ=255)
                                                   gAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                    <  2:1851100/76‑1 (MQ=255)
                                                   gAAGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                    <  1:676054/76‑1 (MQ=255)
                                                    naGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTcc                   >  1:468519/2‑76 (MQ=255)
                                                     aGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                  >  2:1082964/1‑76 (MQ=255)
                                                     aGGTCCGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                  >  2:1906926/1‑76 (MQ=255)
                                                          cGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg             >  1:387903/1‑76 (MQ=255)
                                                          cGAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCg             <  2:1449581/76‑1 (MQ=255)
                                                           gAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt            >  1:1711300/1‑76 (MQ=255)
                                                           gAACATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt            >  1:489360/1‑76 (MQ=255)
                                                            aaCATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa           >  1:1738463/1‑76 (MQ=255)
                                                             aCATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa          >  1:1200629/1‑76 (MQ=255)
                                                             aCATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa          <  2:1568698/76‑1 (MQ=255)
                                                              cATCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa         >  2:551427/1‑76 (MQ=255)
                                                               aTCGGTCTGTTCAACTCTCTGTCCGTGTACGGACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        <  1:1131874/76‑1 (MQ=255)
                                                               aTCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        <  2:450360/76‑1 (MQ=255)
                                                               aTCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        >  2:1338507/1‑76 (MQ=255)
                                                               aTCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        >  2:723670/1‑76 (MQ=255)
                                                               aTCGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg        <  1:361532/76‑1 (MQ=255)
                                                                 cGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg      >  1:106608/1‑76 (MQ=255)
                                                                 cGGTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTg      <  1:1723056/76‑1 (MQ=255)
                                                                   gTCTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCTTAAAGTGAc    >  1:563179/1‑76 (MQ=255)
                                                                     cTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCg  <  1:1277145/76‑1 (MQ=255)
                                                                     cTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCg  >  1:961671/1‑76 (MQ=255)
                                                                     cTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCg  <  1:684812/76‑1 (MQ=255)
                                                                     cTGTTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCCTGAGACTCCGTATCGTAAAGTGACCg  >  1:1255337/1‑76 (MQ=255)
                                                                        |                                                                        
CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  REL606/4162497‑4162641

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