Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,166,065 C→T P369S (CCG→TCG)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,166,0650CT98.4% 243.4 / ‑4.6 64P369S (CCG→TCG) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (36/27);  total (37/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AGCGTGTTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCA  >  REL606/4165994‑4166137
                                                                       |                                                                        
aGCGTGTTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCgaa                                                                      <  2:639047/76‑1 (MQ=255)
 gCGTGTTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCgaag                                                                     <  1:1488874/76‑1 (MQ=255)
 gCGTGTTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCgaag                                                                     <  2:1609827/76‑1 (MQ=255)
   gtgtTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCgaagaa                                                                   <  1:1217960/76‑1 (MQ=255)
    tgtTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGaaa                                                                  <  1:1372852/76‑1 (MQ=255)
    tgtTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGaaa                                                                  >  2:125423/1‑76 (MQ=255)
     gtTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGaaaa                                                                 >  1:508733/1‑76 (MQ=255)
          ctactCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCa                                                            >  1:1799887/1‑76 (MQ=255)
           tactCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCAc                                                           >  2:1306944/1‑76 (MQ=255)
             ctCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTg                                                         >  2:1122826/1‑76 (MQ=255)
             ctCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTg                                                         >  1:572050/1‑76 (MQ=255)
              tCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTgg                                                        <  1:1959979/76‑1 (MQ=255)
              tCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTgg                                                        >  1:174763/1‑76 (MQ=255)
                    cGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGt                                                  >  2:750788/1‑76 (MQ=255)
                    cGNTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGt                                                  >  2:1422753/1‑76 (MQ=255)
                     gTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGtt                                                 >  1:611435/1‑76 (MQ=255)
                     gTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGtt                                                 >  2:135987/1‑76 (MQ=255)
                      ttCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTc                                                >  1:598870/1‑76 (MQ=255)
                      ttCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTc                                                >  2:335321/1‑76 (MQ=255)
                       tCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCa                                               >  1:1833462/1‑76 (MQ=255)
                       tCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCa                                               >  1:39327/1‑76 (MQ=255)
                         tGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCaaa                                             >  1:243273/1‑76 (MQ=255)
                          gTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAAc                                            >  1:805506/1‑76 (MQ=255)
                                  cGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATc                                    <  1:86900/76‑1 (MQ=255)
                                   gTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCt                                   <  2:294937/76‑1 (MQ=255)
                                     aGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTAc                                 <  2:1304305/76‑1 (MQ=255)
                                     aGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTAc                                 <  1:450061/76‑1 (MQ=255)
                                     aGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTAc                                 <  1:1576932/76‑1 (MQ=255)
                                      ggTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACg                                <  1:363564/76‑1 (MQ=255)
                                      ggTCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACg                                <  1:750788/76‑1 (MQ=255)
                                        tCCATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGc                              >  1:1768316/1‑76 (MQ=255)
                                         ccATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCa                             >  1:1866600/1‑76 (MQ=255)
                                         ccATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCa                             <  2:1670073/76‑1 (MQ=255)
                                          cATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                            >  1:495635/1‑76 (MQ=255)
                                          cATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                            >  2:26633/1‑76 (MQ=255)
                                          cATACCTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                            <  2:516178/76‑1 (MQ=255)
                                              ccTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTg                        <  1:1190789/76‑1 (MQ=255)
                                              ccTGCGTCTCCAACGGGGCGGTTTGACGAAGAAAATTGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTg                        <  2:1797215/76‑1 (MQ=255)
                                               cTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTgg                       >  1:164621/1‑76 (MQ=255)
                                               cTGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTgg                       <  2:1309394/76‑1 (MQ=255)
                                                tGCGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGa                      <  2:988384/76‑1 (MQ=255)
                                                 gCGTCTGCATCAGTGCGGTCTTTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGaa                     <  2:912395/76‑1 (MQ=255)
                                                 gCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGaa                     >  1:1129356/1‑76 (MQ=255)
                                                  cGTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAAc                    >  2:770503/1‑76 (MQ=255)
                                                  cGTCTGCAACAGTGCGGTCCGCGGAAGAAAATGGGTCTGGTGCCGGTCAGAGCCTTCATCTCCGGCAAGCTGGAAc                    <  2:1224199/76‑1 (MQ=255)
                                                   gTCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACt                   >  2:10306/1‑76 (MQ=255)
                                                    tCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTg                  >  1:1439626/1‑76 (MQ=255)
                                                    tCTGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTg                  <  2:1062957/76‑1 (MQ=255)
                                                      tGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCg                >  1:481475/1‑76 (MQ=255)
                                                      tGCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCg                <  2:1704480/76‑1 (MQ=255)
                                                       gCATCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGt               >  1:1222966/1‑76 (MQ=255)
                                                         aTCAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTgg             >  1:1543676/1‑76 (MQ=255)
                                                           cAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTc           <  1:1456593/76‑1 (MQ=255)
                                                           cAGTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTc           >  2:976132/1‑76 (MQ=255)
                                                             gTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCtt         <  2:753737/76‑1 (MQ=255)
                                                             gTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCtt         <  2:1918295/76‑1 (MQ=255)
                                                             gTGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCtt         >  1:1334881/1‑76 (MQ=255)
                                                              tGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACTGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTg        >  1:1288949/1‑76 (MQ=255)
                                                              tGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTg        >  1:430135/1‑76 (MQ=255)
                                                              tGCGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTg        >  1:430116/1‑76 (MQ=255)
                                                                cGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCt      >  2:1715112/1‑76 (MQ=255)
                                                                cGGTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCTCCACGGCAAGCTGGAACTGCGTGGTCTTGGTa     >  1:1864036/1‑76 (MQ=255)
                                                                 ggTCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTa     <  2:2015224/76‑1 (MQ=255)
                                                                   tCTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTacc   >  2:962056/1‑76 (MQ=255)
                                                                    cTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTacca  <  1:810302/76‑1 (MQ=255)
                                                                    cTGTCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTacca  <  1:212686/76‑1 (MQ=255)
                                                                       |                                                                        
AGCGTGTTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCA  >  REL606/4165994‑4166137

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