Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,164,097 C→T P1081L (CCG→CTG)  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,0970CT95.9% 174.7 / ‑0.4 49P1081L (CCG→CTG) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (18/29);  total (20/29)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

CGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGT  >  REL606/4164027‑4164168
                                                                      |                                                                       
cGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCg                                                                    <  2:1443329/76‑1 (MQ=255)
 gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGa                                                                   >  1:340037/1‑76 (MQ=255)
 gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGa                                                                   >  1:1501965/1‑76 (MQ=255)
 gCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGa                                                                   >  1:1400278/1‑76 (MQ=255)
  ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGaa                                                                  >  1:894070/1‑76 (MQ=255)
  ccGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGCGTAATTTCTAAGATCAACCTGATCGaa                                                                  >  1:448402/1‑76 (MQ=255)
    gTATCCAGCCTGGCGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGa                                                                >  1:582774/1‑76 (MQ=255)
     tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGat                                                               <  1:1002568/76‑1 (MQ=255)
     tATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGat                                                               <  1:1750118/76‑1 (MQ=255)
       tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGatat                                                             <  1:1516595/76‑1 (MQ=255)
       tCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGCAATTTCTAAGATCAACCTGATCGAAGatat                                                             <  1:619264/76‑1 (MQ=255)
        ccAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATg                                                            <  2:205450/76‑1 (MQ=255)
         cAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGc                                                           >  2:199263/1‑76 (MQ=255)
            ccTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCtt                                                        <  1:1701718/76‑1 (MQ=255)
             cTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTa                                                       <  2:202164/76‑1 (MQ=255)
              tGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTAc                                                      <  2:1745554/76‑1 (MQ=255)
               ggTGACAAGATGGCAGGTCGTCGCGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACg                                                     <  2:1046364/76‑1 (MQ=255)
                gTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGa                                                    <  1:146905/76‑1 (MQ=255)
                  gACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATg                                                  <  1:1419700/76‑1 (MQ=255)
                   aCAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGa                                                 <  2:1624178/76‑1 (MQ=255)
                    cAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGaa                                                <  1:1155825/76‑1 (MQ=255)
                      aGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGGAGATATGCCTTATGATGaaac                                              >  2:1236493/1‑75 (MQ=255)
                       gATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAAc                                             <  1:1691738/76‑1 (MQ=255)
                          ggCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGt                                          <  1:52176/76‑1 (MQ=255)
                               gtcgtcACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGcc                                     >  1:1009928/1‑76 (MQ=255)
                                    cACGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAg                                >  1:1136619/1‑76 (MQ=255)
                                      cGGTAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGAc                              <  2:975211/76‑1 (MQ=255)
                                         tAACAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATc                           <  2:1084645/76‑1 (MQ=255)
                                          aaCAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCg                          >  2:379941/1‑76 (MQ=255)
                                           aCAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGt                         >  2:1277796/1‑76 (MQ=255)
                                            cAAGGGTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTa                        <  2:76218/76‑1 (MQ=255)
                                                ggTGTAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGa                    <  1:871095/76‑1 (MQ=255)
                                                ggTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGa                    >  1:1076185/1‑76 (MQ=255)
                                                 gtgtAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGaa                   <  2:366059/76‑1 (MQ=255)
                                                 gtgtAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGaa                   >  2:19321/1‑76 (MQ=255)
                                                  tgtAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAc                  <  1:983029/76‑1 (MQ=255)
                                                  tgtAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAAc                  <  1:620533/76‑1 (MQ=255)
                                                    tAATTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTAAGCCGGTAGACATCGTACTGAAccc                <  2:877655/76‑1 (MQ=255)
                                                      aTTTCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGc              >  2:316864/1‑76 (MQ=255)
                                                         tCTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTgg           <  1:1487057/76‑1 (MQ=255)
                                                          cTAAGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTggg          >  2:1006440/1‑76 (MQ=255)
                                                           tAAGATCAACCCAATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGc         >  1:1250259/1‑76 (MQ=255)
                                                             aGATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGt       <  1:107037/76‑1 (MQ=255)
                                                              gATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTa      >  1:1415949/1‑76 (MQ=255)
                                                              gATCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTa      >  1:751742/1‑76 (MQ=255)
                                                               aTCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAANCCGCTGGGCGTAc     <  1:1194560/76‑1 (MQ=255)
                                                                tCAACCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTAcc    >  2:1067997/1‑76 (MQ=255)
                                                                  aaCCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGt  <  2:1684122/76‑1 (MQ=255)
                                                                  aaCCTGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGt  <  2:561670/76‑1 (MQ=255)
                                                                      |                                                                       
CGCCGTATCCAGCCTGGTGACAAGATGGCAGGTCGTCACGGTAACAAGGGTGTAATTTCTAAGATCAACCCGATCGAAGATATGCCTTACGATGAAAACGGTACGCCGGTAGACATCGTACTGAACCCGCTGGGCGTACCGT  >  REL606/4164027‑4164168

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: