Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,169,030 T→A I1357N (ATC→AAC)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,169,0300TA97.6% 153.4 / ‑3.5 41I1357N (ATC→AAC) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (27/13);  total (28/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

ACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTG  >  REL606/4168957‑4169101
                                                                         |                                                                       
aCCGAAGAAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAAcc                                                                       >  1:655263/1‑76 (MQ=255)
 ccGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAAccc                                                                      >  1:150428/1‑76 (MQ=255)
      gcagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAg                                                                 <  2:5293/76‑1 (MQ=255)
       cagcCTGTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAgg                                                                <  2:1804626/76‑1 (MQ=255)
       cagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAgg                                                                >  1:1627925/1‑76 (MQ=255)
       cagcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGCCTGAACCCGGCAgg                                                                >  1:663646/1‑76 (MQ=255)
        agcCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAggg                                                               >  2:854447/1‑75 (MQ=255)
          ccGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTAc                                                             <  2:891651/76‑1 (MQ=255)
          ccGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTAc                                                             <  2:796769/76‑1 (MQ=255)
            gTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCg                                                           <  1:696464/76‑1 (MQ=255)
               gCGGGCAAACTCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGtt                                                        <  2:1268635/76‑1 (MQ=255)
                cGGGCAAGCGCGGCGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTa                                                       <  2:79967/76‑1 (MQ=255)
                cGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTa                                                       >  2:579283/1‑76 (MQ=255)
                cGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTa                                                       <  2:1132770/76‑1 (MQ=255)
                cGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGCACCGGTTa                                                       <  2:25560/76‑1 (MQ=255)
                    cAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTAcgcg                                                   >  1:1476614/1‑76 (MQ=255)
                     aaaCGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGt                                                  <  2:791816/76‑1 (MQ=255)
                     aaaCGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGt                                                  >  1:673448/1‑76 (MQ=255)
                         gcgACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTacca                                              <  1:818312/76‑1 (MQ=255)
                           nacgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTaccacc                                            >  1:117491/2‑76 (MQ=255)
                            acgaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTAccacca                                           <  1:485403/76‑1 (MQ=255)
                              gaACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAgg                                         >  1:201725/1‑76 (MQ=255)
                               aaCTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGAACCGGTTACGCGAAACAACAGGa                                        >  2:676866/1‑76 (MQ=255)
                                aCTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGAt                                       >  2:1152842/1‑76 (MQ=255)
                                         ccTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGt                              >  2:1038227/1‑76 (MQ=255)
                                          cTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTc                             >  1:1623609/1‑76 (MQ=255)
                                          cTGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTc                             >  2:360566/1‑76 (MQ=255)
                                           tGAAAGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCg                            >  1:1319496/1‑76 (MQ=255)
                                             aaaGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGcc                          >  1:737840/1‑76 (MQ=255)
                                             aaaGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGcc                          >  1:899451/1‑76 (MQ=255)
                                             aaaGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGcc                          >  2:156113/1‑76 (MQ=255)
                                              aaGAGAACGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCg                         >  2:205779/1‑76 (MQ=255)
                                                   aaCGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGtgctg                    >  2:465562/1‑76 (MQ=255)
                                                   aaCGTTATCGGGGGTCGTCTGACCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGtgctg                    <  1:645072/76‑1 (MQ=255)
                                                    aCGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGtgctgc                   >  2:1038900/1‑76 (MQ=255)
                                                     cGTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCg                  <  2:1181508/76‑1 (MQ=255)
                                                      gTTATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCgg                 >  1:1191774/1‑76 (MQ=255)
                                                        tATCGTGGGTCGTCTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCgggg               >  1:134467/1‑75 (MQ=255)
                                                        tATCGTGGGGCGTCTGCACCCGGCAGGTACCGGTTACGCGCACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGt               <  2:1620066/76‑1 (MQ=255)
                                                           cGTGGGTCGTCTGCACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCGGCGGGGGaa            >  1:1214060/1‑76 (MQ=255)
                                                                  cgtcTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCgg     >  1:1685767/1‑76 (MQ=255)
                                                                  cgtcTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCgg     >  1:1650868/1‑76 (MQ=255)
                                                                  cgtcTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGGGCTGCGGGTGAAGCGCCgg     >  1:111061/1‑76 (MQ=255)
                                                                    tcTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGGGAAGCTCCGGCt   >  1:716238/1‑76 (MQ=255)
                                                                    tcTGAACCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGGGAAGCTCCGGCt   >  1:278743/1‑76 (MQ=255)
                                                                     cTGATCCCGGCCGGAACCGGGTACGCGTACCACCAGGATCGAATGCGTCGCCGTGCTGCGGGGGAAGCTCCGGcgg  >  2:623131/1‑74 (MQ=255)
                                                                         |                                                                       
ACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCAGGTACCGGTTACGCGTACCACCAGGATCGTATGCGTCGCCGTGCTGCGGGTGAAGCTCCGGCTG  >  REL606/4168957‑4169101

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