Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,162,569 A→C I572L (ATC→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,162,5690AC97.5% 155.3 / NA 40I572L (ATC→CTC) rpoBDNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (21/18);  total (22/18)
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

CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  REL606/4162497‑4162641
                                                                        |                                                                        
cGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCa                                                                       <  1:933942/76‑1 (MQ=255)
      gTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActctc                                                                 >  1:1065505/1‑76 (MQ=255)
       tACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAActctct                                                                <  2:559966/76‑1 (MQ=255)
        ccacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTg                                                               <  2:1198465/75‑1 (MQ=255)
        acacCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTg                                                               >  2:1156246/1‑76 (MQ=255)
           cccGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTcc                                                            <  1:1042326/76‑1 (MQ=255)
              gACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCgtg                                                         >  1:1209196/1‑76 (MQ=255)
               aCTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCgtgt                                                        >  2:1274122/1‑76 (MQ=255)
                 tCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTAc                                                      <  1:867371/76‑1 (MQ=255)
                     tACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGcac                                                  <  1:653321/76‑1 (MQ=255)
                      aCGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGcaca                                                 <  2:626256/76‑1 (MQ=255)
                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTa                                            <  1:239990/76‑1 (MQ=255)
                           cgcgTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTa                                            >  1:261299/1‑76 (MQ=255)
                           cgcgAATGTACAATCGAAACCGCTGAATGTCCTAGCCTTGGCCTGGTCATCTCTCTGAGCGTGTACGCACAGACTa                                            <  1:1341115/76‑1 (MQ=255)
                            gcgTATTTCAAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTNTANGCACAGACTaa                                           >  2:210203/1‑76 (MQ=255)
                            gcgTATGTCCAATCGAAACCTCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGGCTGTATACGCACAGACTaa                                           >  1:60470/1‑76 (MQ=255)
                               tATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGa                                        <  1:803933/76‑1 (MQ=255)
                                aTGTCCAATCGAAAGCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACTGACTAACGaa                                       <  1:115427/76‑1 (MQ=255)
                                aTGTCCAATCGAAACCCCTGAAGGTCCGAACCTCNGTCTACTNTANTCNCTGTACGTGTACGCACAGACTAACGaa                                       >  1:1163535/1‑76 (MQ=255)
                                aTGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGaa                                       >  1:836970/1‑76 (MQ=255)
                                    ccAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACg                                   >  2:173867/1‑76 (MQ=255)
                                     cAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACgg                                  >  2:1058401/1‑76 (MQ=255)
                                      aaTCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGc                                 >  1:532164/1‑76 (MQ=255)
                                       aTCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCt                                <  2:250849/76‑1 (MQ=255)
                                       aTCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTACGTGTACGCACAGACTAACGAATACGGCt                                >  2:33173/1‑76 (MQ=255)
                                        tCGAAACCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                               <  2:539259/76‑1 (MQ=255)
                                           aaaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCt                            <  2:173018/76‑1 (MQ=255)
                                            aaCCCCTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                           >  2:122802/1‑76 (MQ=255)
                                                ccTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgaga                       >  2:513226/1‑76 (MQ=255)
                                                ccTGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgaga                       >  2:320146/1‑76 (MQ=255)
                                                  tGAAGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTTAGACt                     >  1:598800/1‑76 (MQ=255)
                                                     aGGTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCg                  <  1:1137040/76‑1 (MQ=255)
                                                       gTCCGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTGACTAATACGGCTTCCTTGAGACTCCGTa                >  1:214856/1‑76 (MQ=255)
                                                          cGAACATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTNTCg             <  1:1207304/76‑1 (MQ=255)
                                                            aaCATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa           >  1:45482/1‑76 (MQ=255)
                                                            aaCATCGGTCNGATCAACTCTNTGTCCGTGTACGAACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTa           >  2:1323172/1‑76 (MQ=255)
                                                             aCATCGGTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaa          <  1:1063686/76‑1 (MQ=255)
                                                                 cGGTCTGCTCAACTCTCTGTCCGTGNATGCACAGACGAACGAATACANCTTGCTTGAGACTCCGTATCGTAAAGTg      <  1:194674/76‑1 (MQ=255)
                                                                  ggTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa     >  2:602558/1‑76 (MQ=255)
                                                                  ggTCTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGa     >  2:706774/1‑76 (MQ=255)
                                                                     cTGCTCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCg  <  2:732735/76‑1 (MQ=255)
                                                                        |                                                                        
CGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCG  >  REL606/4162497‑4162641

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: