Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,161,145 G→T R97L (CGC→CTC)  rpoB → DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,161,1450GT97.5% 147.4 / ‑2.4 40R97L (CGC→CTC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (18/21);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGA  >  REL606/4161072‑4161217
                                                                         |                                                                        
tACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGctcg                                                                        >  1:1008409/1‑76 (MQ=255)
 aCCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGt                                                                       <  1:1139044/76‑1 (MQ=255)
 aCCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGt                                                                       <  1:1591299/76‑1 (MQ=255)
  ccGCCTTGGCGAGCCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGtt                                                                      <  1:1045978/76‑1 (MQ=255)
  ccGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGtt                                                                      <  2:1089496/76‑1 (MQ=255)
   cGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTa                                                                     <  1:228438/76‑1 (MQ=255)
     ccTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTaaa                                                                   <  2:432921/76‑1 (MQ=255)
       ttGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACt                                                                 <  1:2170395/76‑1 (MQ=255)
        tGGCGAACCGGTGTTTGACGTCCAGGGATGTAAAATCCGTGGCGTGACCTGTTCCGCACCGCTGCTCGTTAAACTg                                                                >  1:1897964/1‑76 (MQ=255)
        tGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTg                                                                >  2:1907583/1‑76 (MQ=255)
           cGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGt                                                             <  1:900818/76‑1 (MQ=255)
            gAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTc                                                            >  2:1060654/1‑76 (MQ=255)
             aaCCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCGGCTCGTTAAACAGCgcgg                                                           >  2:76438/1‑73 (MQ=255)
                   tgtTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCCCCGCTGCTCGTTAAACTGCGTCTGGTGAt                                                     <  1:2018652/76‑1 (MQ=255)
                    gtTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc                                                    >  2:8491/1‑76 (MQ=255)
                      ttGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTa                                                  <  1:1880021/76‑1 (MQ=255)
                       tGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTAt                                                 >  1:968757/1‑76 (MQ=255)
                       tGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTAt                                                 >  1:2166009/1‑76 (MQ=255)
                         aCGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGa                                               >  1:243378/1‑76 (MQ=255)
                         aCGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGa                                               <  1:1701882/76‑1 (MQ=255)
                             ccAGGAATGTCAAATCCGTCGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAgcgc                                           <  2:1562849/76‑1 (MQ=255)
                              cAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGcgcg                                          <  1:1983223/76‑1 (MQ=255)
                               aGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGa                                         >  1:448142/1‑76 (MQ=255)
                               aGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGa                                         >  2:123066/1‑76 (MQ=255)
                                ggAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGaa                                        <  1:272082/76‑1 (MQ=255)
                                      tCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGATGCGCCg                                  <  2:1771686/76‑1 (MQ=255)
                                      tCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCg                                  <  2:743466/76‑1 (MQ=255)
                                      tCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCg                                  <  2:1004013/76‑1 (MQ=255)
                                        aaaTCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGa                                <  2:467137/76‑1 (MQ=255)
                                        aaaTCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGa                                >  1:2098951/1‑76 (MQ=255)
                                         aaTCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGaa                               >  1:533647/1‑76 (MQ=255)
                                          aTCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATTTATGAGCGCGAAGCGCCGGAAg                              >  2:91964/1‑76 (MQ=255)
                                          aTCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAg                              >  1:164837/1‑76 (MQ=255)
                                           tCCGTGGCGTGCCCTATTCCGCCCCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAgg                             <  2:1555553/76‑1 (MQ=255)
                                           tCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAgg                             >  2:1901982/1‑76 (MQ=255)
                                           tCCGTGGCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAgg                             <  2:261674/76‑1 (MQ=255)
                                                ggCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCg                        <  2:140379/76‑1 (MQ=255)
                                                 gCGTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGt                       >  1:1437366/1‑76 (MQ=255)
                                                  ggTGACCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTAGGAGCGCGAAGCGCCGGAAGGCACCGTa                      <  2:171567/75‑1 (MQ=255)
                                                  cGTGCCCTATTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTa                      <  2:698811/76‑1 (MQ=255)
                                                      aCCTATTCCGCGCCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAg                  <  2:1850206/76‑1 (MQ=255)
                                                          aTTCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACAt              >  1:780540/1‑76 (MQ=255)
                                                           ttCCGCACCGCTGCTCGTTAAACTGCGTCTGGTAATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACAtt             <  2:2192739/76‑1 (MQ=255)
                                                            tCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTa            >  2:1209776/1‑76 (MQ=255)
                                                            tCCGCACCGCTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTa            >  1:1461041/1‑76 (MQ=255)
                                                             ccGCCCCGCTGCTCTTTAAACTGCGTTTGGTGATCTAAGAGCGCGAGGCGCCGGAAGGCACCGTAAAAGACATTaa           <  1:2104163/76‑1 (MQ=255)
                                                                  ccGCTGCTCTTAAACCGGCGTGGGGTGATTTGGGGGCGGGAGGCCCCGGAAGGCACCATAAAAGACATTAAAGAAc      <  1:338717/76‑1 (MQ=255)
                                                                     cTGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACaag   >  2:1959666/1‑76 (MQ=255)
                                                                      tGCTCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACaaga  >  1:1570307/1‑76 (MQ=255)
                                                                      tGCTCGTTAAACGGCGTCTGGTGATCTAGGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACaaga  <  2:65999/76‑1 (MQ=255)
                                                                         |                                                                        
TACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGCGAAGCGCCGGAAGGCACCGTAAAAGACATTAAAGAACAAGA  >  REL606/4161072‑4161217

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