Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,161,105 G→A E84K (GAA→AAA)  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,1050GA97.7% 163.9 / ‑3.9 45E84K (GAA→AAA) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/23);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

TCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGC  >  REL606/4161032‑4161176
                                                                         |                                                                       
tCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGaaa                                                                       >  1:1052425/1‑76 (MQ=255)
tCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGaaa                                                                       <  1:40210/76‑1 (MQ=255)
  agagCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATg                                                                     >  1:222732/1‑76 (MQ=255)
    agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc                                                                   >  2:368864/1‑76 (MQ=255)
    agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc                                                                   <  1:713823/76‑1 (MQ=255)
    agCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTc                                                                   >  1:1807076/1‑76 (MQ=255)
       tACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGCAATGTCaaa                                                                >  1:1282670/1‑76 (MQ=255)
       tACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCaaa                                                                <  1:613074/76‑1 (MQ=255)
        aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt                                                               <  2:925780/76‑1 (MQ=255)
        aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt                                                               >  2:1690714/1‑76 (MQ=255)
        aCAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAAt                                                               <  2:825590/76‑1 (MQ=255)
               tAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGc                                                        >  2:524633/1‑76 (MQ=255)
                aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCg                                                       >  2:75262/1‑76 (MQ=255)
                aaTTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAACCCGTGGCg                                                       >  2:75245/1‑76 (MQ=255)
                  ttCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTg                                                     <  2:1468207/76‑1 (MQ=255)
                  ttCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTg                                                     <  2:1856855/76‑1 (MQ=255)
                   tCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGa                                                    <  2:583249/76‑1 (MQ=255)
                    ccGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGAc                                                   >  1:149582/1‑76 (MQ=255)
                      gAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCt                                                 >  2:573768/1‑76 (MQ=255)
                      gAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCt                                                 >  2:956096/1‑76 (MQ=255)
                       aGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTa                                                >  2:164592/1‑76 (MQ=255)
                         cTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTAtt                                              >  2:1508739/1‑76 (MQ=255)
                                   tCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCt                                    <  2:509168/76‑1 (MQ=255)
                                     aGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTgc                                  >  1:333893/1‑76 (MQ=255)
                                     aGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTgc                                  >  2:1469339/1‑76 (MQ=255)
                                          ccGCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGtt                             <  2:1117745/76‑1 (MQ=255)
                                            gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa                           >  1:1201179/1‑76 (MQ=255)
                                            gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa                           <  2:563303/76‑1 (MQ=255)
                                            gCCTTGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTaa                           >  2:376307/1‑76 (MQ=255)
                                               ttGGCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACt                        <  1:1000082/76‑1 (MQ=255)
                                                  gCGAACCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCg                     >  1:124138/1‑76 (MQ=255)
                                                     aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt                  <  1:791103/76‑1 (MQ=255)
                                                     aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt                  >  1:1497435/1‑76 (MQ=255)
                                                     aaCCGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCATGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCt                  >  1:1345740/1‑76 (MQ=255)
                                                       ccGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTgg                <  2:1151066/76‑1 (MQ=255)
                                                        cGGTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGt               <  2:1509742/76‑1 (MQ=255)
                                                         ggTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg              <  1:402648/76‑1 (MQ=255)
                                                         ggTGTTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTg              <  2:376228/76‑1 (MQ=255)
                                                           tgtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGAt            <  1:1615596/76‑1 (MQ=255)
                                                            gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc           <  1:376307/76‑1 (MQ=255)
                                                            gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc           <  1:1071087/76‑1 (MQ=255)
                                                            gtTTGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATc           >  1:1131320/1‑76 (MQ=255)
                                                               tGACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTAt        <  2:814515/76‑1 (MQ=255)
                                                                gACGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATg       <  2:963818/76‑1 (MQ=255)
                                                                  cGTGCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAg     <  1:387146/76‑1 (MQ=255)
                                                                  cGTCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAg     <  1:1500552/76‑1 (MQ=255)
                                                                    tCCAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcg   <  1:111121/76‑1 (MQ=255)
                                                                     ccAGAAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAgcgc  <  1:166771/76‑1 (MQ=255)
                                                                         |                                                                       
TCAGAGCTACAGCGGTAATTCCGAGCTGCAATACGTCAGCTACCGCCTTGGCGAACCGGTGTTTGACGTCCAGGAATGTCAAATCCGTGGCGTGACCTATTCCGCACCGCTGCGCGTTAAACTGCGTCTGGTGATCTATGAGCGC  >  REL606/4161032‑4161176

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