Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,164,745 A→C D1297A (GAC→GCC)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,164,7450AC100.0% 165.8 / NA 41D1297A (GAC→GCC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (20/21);  total (20/21)

GAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGA  >  REL606/4164675‑4164817
                                                                      |                                                                        
gAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTNATGCCGTGa                                                                     <  1:31484/76‑1 (MQ=255)
gAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGa                                                                     >  2:1129557/1‑76 (MQ=255)
  aGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATNCTCACCTTTAAGTCTGATGCCGTGAAc                                                                   <  1:1400102/76‑1 (MQ=255)
       gggNGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATNCTCACCGTTAAGTCTGATGCCGTGAACGGTCg                                                              <  2:1305088/76‑1 (MQ=255)
        ggCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGt                                                             <  1:540705/76‑1 (MQ=255)
         gcgcTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTa                                                            <  1:621032/76‑1 (MQ=255)
         gcgcTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTa                                                            <  2:1396515/76‑1 (MQ=255)
            cTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCg                                                              <  2:421514/71‑1 (MQ=255)
            cTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCg                                                              >  1:421514/1‑71 (MQ=255)
            cTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCg                                                              <  2:302618/71‑1 (MQ=255)
            cTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCg                                                              >  1:302618/1‑71 (MQ=255)
            cTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCa                                                         >  1:144310/1‑76 (MQ=255)
              ggAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGNCGTGAACGGTCGTACCAAg                                                       <  1:1302050/76‑1 (MQ=255)
                aaGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGAt                                                     >  2:952381/1‑76 (MQ=255)
                 aGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCCTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATg                                                    >  2:1029111/1‑76 (MQ=255)
                       cGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATaaa                                              <  2:15522/76‑1 (MQ=255)
                         gcgcAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATaaaaa                                            <  1:1221651/76‑1 (MQ=255)
                          cgcAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAAc                                           >  1:192140/1‑76 (MQ=255)
                            cagcaTACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACAt                                         >  1:779619/1‑76 (MQ=255)
                            cagcaTACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTCAACGGTCGTACCAAGATGTATAAAAACAt                                         >  2:376809/1‑76 (MQ=255)
                             agcaTACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATc                                        <  2:991008/76‑1 (MQ=255)
                             agcaTACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATc                                        <  1:1368301/76‑1 (MQ=255)
                                aTACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTg                                     >  2:1372622/1‑76 (MQ=255)
                                 tACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTgg                                    <  2:737524/76‑1 (MQ=255)
                                   cacCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGAc                                  <  2:160336/76‑1 (MQ=255)
                                    acCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACg                                 >  2:551281/1‑76 (MQ=255)
                                    acCCTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACg                                 >  2:763171/1‑76 (MQ=255)
                                     cccTGCAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACgg                                >  1:256999/1‑76 (MQ=255)
                                          cAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAAcc                           >  1:356376/1‑76 (MQ=255)
                                          cAGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTAGTACCAAGATGTATAAAAACATCGTGGACGGCAAcc                           <  1:871296/76‑1 (MQ=255)
                                           aGGAAATGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAAcct                          >  1:304731/1‑75 (MQ=255)
                                            ggAAATGCTCACCGTTAAGTCTGATTCCNTGAACGGTCCAAACAAGAGGTATAAAAACATCGTGGACGGCAACCAt                         <  1:122684/76‑1 (MQ=255)
                                            ggAAACGCCCACCGTTAAGTCTGATCCCGTGAACGGTCGTACCAAGATGTATAAAACCACCGTGGACGGCAACCAt                         <  2:1030131/76‑1 (MQ=255)
                                              aaaTGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCa                       >  1:383485/1‑76 (MQ=255)
                                               aaTGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAg                      >  2:808035/1‑76 (MQ=255)
                                               aaTGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAg                      <  1:969324/76‑1 (MQ=255)
                                               aaTGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAg                      >  2:808037/1‑76 (MQ=255)
                                                 tGCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGAt                    <  1:278979/76‑1 (MQ=255)
                                                  gCTCACCGTTAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATg                   <  1:413928/76‑1 (MQ=255)
                                                   cTCACCGTTAAGTCNGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATgn                  <  2:193693/76‑2 (MQ=255)
                                                           tAAGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGc          >  2:561683/1‑76 (MQ=255)
                                                            aaGTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCa         <  1:749770/76‑1 (MQ=255)
                                                              gTCTGATGCCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATg       >  1:1061266/1‑76 (MQ=255)
                                                                 tgatgcCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCa    <  2:658360/76‑1 (MQ=255)
                                                                 tgatgcCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCa    >  2:291587/1‑76 (MQ=255)
                                                                  gatgcCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCGTCAGATGGAGCCTGGCATGCCAg   >  2:1025348/1‑76 (MQ=255)
                                                                  gatgcCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAg   <  1:835859/76‑1 (MQ=255)
                                                                  gatgcCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAg   >  1:568506/1‑76 (MQ=255)
                                                                   atgcCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGa  <  1:251925/76‑1 (MQ=255)
                                                                   atgcCGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGa  >  1:818225/1‑76 (MQ=255)
                                                                      |                                                                        
GAAGTGTGGGCGCTGGAAGCATACGGCGCAGCATACACCCTGCAGGAAATGCTCACCGTTAAGTCTGATGACGTGAACGGTCGTACCAAGATGTATAAAAACATCGTGGACGGCAACCATCAGATGGAGCCGGGCATGCCAGA  >  REL606/4164675‑4164817

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