Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,165,277 A→C E106A (GAA→GCA)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,165,2770AC98.0% 206.0 / NA 52E106A (GAA→GCA) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base C (25/25);  minor base G (0/1);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATAT  >  REL606/4165204‑4165349
                                                                         |                                                                        
ggCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCAc                                                                        >  1:349780/1‑76 (MQ=255)
     cATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCt                                                                   >  1:656161/1‑76 (MQ=255)
     cATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCt                                                                   <  2:1370177/76‑1 (MQ=255)
         tgtgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTccc                                                               >  1:1080669/1‑76 (MQ=255)
          gtgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTcccc                                                              >  2:1404195/1‑76 (MQ=255)
           tgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCg                                                             <  1:1325220/76‑1 (MQ=255)
           tgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCg                                                             <  2:636316/76‑1 (MQ=255)
            gagaAGTGCGGCGTTGAAGTGACCCTGAGTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGa                                                            <  1:324024/76‑1 (MQ=255)
            gagaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGa                                                            <  1:215835/76‑1 (MQ=255)
             agaAGTGCGGCGTTNAAGTGACCCAGACTAAAGTACGCCGTGTGCGTATGGGCCACATCGGACTGGCTTCCCCGAc                                                           <  1:1750708/76‑1 (MQ=255)
             agaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGAc                                                           <  1:918885/76‑1 (MQ=255)
                  tGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTgcgc                                                      <  1:975369/76‑1 (MQ=255)
                   gcggcgTTGAAGTGACCCAGGCTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCATCCCCGACTGCGca                                                     >  1:98360/1‑76 (MQ=255)
                     ggcgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGcaca                                                   >  1:833147/1‑76 (MQ=255)
                         ttGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTg                                               >  1:905861/1‑76 (MQ=255)
                          tGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTgg                                              >  1:1449023/1‑76 (MQ=255)
                          tGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTgg                                              >  2:1336485/1‑76 (MQ=255)
                               tGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCt                                         <  1:438420/76‑1 (MQ=255)
                                gACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTg                                        <  1:976797/76‑1 (MQ=255)
                                 aCCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGa                                       >  1:1645793/1‑76 (MQ=255)
                                 aCCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGa                                       <  1:359296/76‑1 (MQ=255)
                                  cccAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGaa                                      <  2:1509730/76‑1 (MQ=255)
                                     aGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATc                                   <  1:979437/76‑1 (MQ=255)
                                      gACTAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCg                                  <  2:689643/76‑1 (MQ=255)
                                       aCTAAAGTACGCCGTGAGCGTATGGGCTACAACGNACTGACTNACNCANCTGCTCACCTCTGGTTCCTGAAATCGc                                 >  1:1759009/1‑76 (MQ=255)
                                         tAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGTGGAACTCTTGTTTCTGTCATCGcgt                               >  2:637416/1‑74 (MQ=255)
                                         tAAAGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGGAAACGCTg                               >  2:637442/1‑76 (MQ=255)
                                          aaaGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGc                              >  2:587231/1‑76 (MQ=255)
                                          aaaGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGc                              >  2:1104228/1‑76 (MQ=255)
                                          aaaGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGc                              <  2:1760277/76‑1 (MQ=255)
                                          aaaGTACGCCGTGAGCGTACGGGCCACATCGNCCTGACGTCACCGACTGCGCACATCTGGTTCGNGAAATCGCTGc                              >  1:245287/1‑76 (MQ=255)
                                           aaGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGcc                             <  1:1701237/76‑1 (MQ=255)
                                            aGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCg                            <  2:1037385/76‑1 (MQ=255)
                                            aGTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCg                            <  2:1538203/76‑1 (MQ=255)
                                             gTACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCCCTGCCGt                           >  1:949931/1‑76 (MQ=255)
                                              tACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTc                          <  1:1591069/76‑1 (MQ=255)
                                              tACGCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTc                          >  2:1238824/1‑76 (MQ=255)
                                                 gCCGTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCg                       >  1:724989/1‑76 (MQ=255)
                                                 gCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCg                       >  1:1461963/1‑76 (MQ=255)
                                                    gTGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTAt                    <  2:1771586/76‑1 (MQ=255)
                                                     tGAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATc                   >  1:550686/1‑76 (MQ=255)
                                                      gAGCGTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCg                  <  1:968257/76‑1 (MQ=255)
                                                          gTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTct              <  2:1301546/76‑1 (MQ=255)
                                                          gTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTct              >  2:1293231/1‑76 (MQ=255)
                                                          gTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTct              >  1:1685071/1‑76 (MQ=255)
                                                          gTATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTAGCTGTct              <  1:1914078/76‑1 (MQ=255)
                                                           tATGGGCCACATCGCACTGGCTTCCNCGACTGCGCACATCTGGTTCTNGAATTCGCTGCCGTCCCGTATCGGTctg             <  1:285663/76‑1 (MQ=255)
                                                           tATGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTctg             <  1:1628968/76‑1 (MQ=255)
                                                           tATGCGCCACATCGCACTGGCTTCACCGACTGCGCACATCTGGTTCCTGAAACCGCTGCCGTCCCGTATCGGTctg             <  2:1454120/76‑1 (MQ=255)
                                                            aTGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTAGCGGGctgc            >  2:1917815/1‑76 (MQ=255)
                                                             tGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTctgct           >  1:1806290/1‑76 (MQ=255)
                                                             tGGGCCACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGAATCGGTctgct           >  2:1824136/1‑76 (MQ=255)
                                                                gCCACATCGCGCTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTgctgct        <  2:1841043/76‑1 (MQ=255)
                                                                 ccACATCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTc       <  2:833146/76‑1 (MQ=255)
                                                                    caTCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGat    >  1:171005/1‑76 (MQ=255)
                                                                    caTCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCCGCCTGCTGCTCGat    >  2:499828/1‑76 (MQ=255)
                                                                      tCGCACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGatat  <  2:758197/76‑1 (MQ=255)
                                                                         |                                                                        
GGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATAT  >  REL606/4165204‑4165349

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