Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,168,967 C→T A1336V (GCC→GTC)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,168,9670CT95.9% 167.8 / ‑0.2 49A1336V (GCC→GTC) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (22/25);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCA  >  REL606/4168896‑4169037
                                                                       |                                                                      
gTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGtt                                                                    <  2:546209/76‑1 (MQ=255)
gTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGtt                                                                    <  2:1495487/76‑1 (MQ=255)
 tctcTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTg                                                                   <  2:313835/76‑1 (MQ=255)
 tctcTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTg                                                                   <  2:1399081/76‑1 (MQ=255)
  ctctGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGc                                                                  >  1:331920/1‑76 (MQ=255)
    ctGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCgg                                                                >  2:15198/1‑76 (MQ=255)
    ctGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCgg                                                                >  1:763528/1‑76 (MQ=255)
     tGGCAACCGTGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCggg                                                               >  2:1160747/1‑76 (MQ=255)
      ggCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGc                                                              >  2:811850/1‑76 (MQ=255)
      ggCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGc                                                              >  2:1429937/1‑76 (MQ=255)
        cAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCaa                                                            <  1:610068/76‑1 (MQ=255)
        cAACCGAGTCCTTCATCTACGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCaa                                                            >  1:254557/1‑76 (MQ=255)
         aaCCGAGTCCTTCATCTCCGCGGCATCGTTTCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCaaa                                                           <  1:30555/76‑1 (MQ=255)
         aaCCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCaaa                                                           >  1:1103940/1‑76 (MQ=255)
           ccGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAAcg                                                         <  2:925842/76‑1 (MQ=255)
           ccGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCGGACCGAAGCAGACGATGCGGGGAAAcg                                                         >  1:1681811/1‑76 (MQ=255)
              aGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGcga                                                      <  2:895459/76‑1 (MQ=255)
              aGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGcga                                                      <  2:574277/76‑1 (MQ=255)
                 ccTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGcgacga                                                   <  2:614355/76‑1 (MQ=255)
                  cTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGaa                                                  >  1:617021/1‑76 (MQ=255)
                    tCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACt                                                <  1:633582/76‑1 (MQ=255)
                    tCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACt                                                <  1:1460248/76‑1 (MQ=255)
                    tCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACt                                                <  1:9824/76‑1 (MQ=255)
                    tCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACt                                                <  1:186975/76‑1 (MQ=255)
                      aTCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTgc                                              >  1:402614/1‑76 (MQ=255)
                           cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGc                                         >  2:1361156/1‑76 (MQ=255)
                           cgcgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGc                                         >  2:1208243/1‑76 (MQ=255)
                             ctgCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCt                                       <  2:618004/74‑1 (MQ=255)
                             cgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCt                                       <  2:124615/76‑1 (MQ=255)
                             cgGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCt                                       <  2:380507/76‑1 (MQ=255)
                               gCATCGGTCCAGGAGACCCCCCGCGGGGTGGCCGAAGCCGGCGTTGCGGGCAAACGCGACGCACCGCGCGGCCTGa                                     >  2:241777/1‑76 (MQ=255)
                                cATCGTTCCAGGAGCCCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGaa                                    <  1:1307463/76‑1 (MQ=255)
                                cATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGaa                                    >  1:1445985/1‑76 (MQ=255)
                                     ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAgaga                               <  2:1134737/76‑1 (MQ=255)
                                     ttCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAgaga                               <  2:833461/76‑1 (MQ=255)
                                     tgccAGGAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAgaga                               <  2:1092604/74‑1 (MQ=255)
                                      tCCAGGAGACCACTCGCGTGCTGACCGAAGCAGTCTTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGaa                              <  2:1115787/76‑1 (MQ=255)
                                        cAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACg                            <  1:663636/76‑1 (MQ=255)
                                          ggAGACCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGtt                          <  1:643888/76‑1 (MQ=255)
                                           gagaCCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTa                         <  1:1281314/76‑1 (MQ=255)
                                            agaCCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTAt                        >  2:1014078/1‑76 (MQ=255)
                                            agaCCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTAt                        >  2:15103/1‑76 (MQ=255)
                                            agaCCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTAt                        <  2:1563690/76‑1 (MQ=255)
                                             gaCCACTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATc                       <  2:670502/76‑1 (MQ=255)
                                                  cTCGCGTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTggg                  >  1:86433/1‑76 (MQ=255)
                                                     gcgTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGgtcg               <  2:1192094/76‑1 (MQ=255)
                                                      cgTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGgtcgt              >  2:205006/1‑76 (MQ=255)
                                                       gTGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGgtcgtc             <  2:632176/76‑1 (MQ=255)
                                                        tGCTGACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCt            >  1:1439455/1‑76 (MQ=255)
                                                            gACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGGTc        >  1:981607/1‑76 (MQ=255)
                                                            gACCGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATc        >  1:515584/1‑76 (MQ=255)
                                                              ccGAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATccc      >  1:445854/1‑76 (MQ=255)
                                                               cGAAGCAGTCGGTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCg     <  1:14430/76‑1 (MQ=255)
                                                                gAAGCAGTCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCgg    >  1:1051038/1‑76 (MQ=255)
                                                                  agcagtCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCa  <  1:451771/76‑1 (MQ=255)
                                                                       |                                                                      
GTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCGGGCAAACGCGACGAACTGCGCGGCCTGAAAGAGAACGTTATCGTGGGTCGTCTGATCCCGGCA  >  REL606/4168896‑4169037

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