Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,166,077 G→A A373T (GCA→ACA)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,166,0770GA94.9% 213.3 / 0.5 59A373T (GCA→ACA) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (31/25);  total (33/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

CTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAA  >  REL606/4166004‑4166145
                                                                         |                                                                    
ctactCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACa                                                                    >  1:1361366/1‑76 (MQ=255)
     ccGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGa                                                               >  2:1013257/1‑76 (MQ=255)
     ccGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGa                                                               >  1:1313565/1‑76 (MQ=255)
     ccGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGa                                                               >  2:914935/1‑76 (MQ=255)
     ccGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGa                                                               >  2:913947/1‑76 (MQ=255)
      cGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAg                                                              <  1:1323678/76‑1 (MQ=255)
       ggTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGc                                                             <  1:773060/76‑1 (MQ=255)
          cGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGt                                                          >  2:1252968/1‑76 (MQ=255)
          cGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGt                                                          >  1:855457/1‑76 (MQ=255)
            ttCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTc                                                        <  1:156541/76‑1 (MQ=255)
            ttCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTc                                                        >  1:552247/1‑76 (MQ=255)
                gTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAAc                                                    >  2:925165/1‑76 (MQ=255)
                gTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAAc                                                    <  1:975134/76‑1 (MQ=255)
                   aTCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGt                                                 >  1:942880/1‑76 (MQ=255)
                         gTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCt                                           >  1:281446/1‑76 (MQ=255)
                          tAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTa                                          <  1:596345/76‑1 (MQ=255)
                          tAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTa                                          <  2:1145656/76‑1 (MQ=255)
                          tAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTa                                          <  2:1156239/76‑1 (MQ=255)
                          tAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTa                                          <  1:301930/76‑1 (MQ=255)
                          tAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTa                                          >  1:313839/1‑76 (MQ=255)
                             gTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACgg                                       <  2:383594/76‑1 (MQ=255)
                               ggaTACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCa                                     >  2:1374399/3‑76 (MQ=255)
                               ccATTCCCGCGTCTGCAGCATAGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCa                                     <  2:1196988/76‑1 (MQ=255)
                               ccATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCa                                     >  2:1399437/1‑76 (MQ=255)
                                cATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                                    <  2:165509/76‑1 (MQ=255)
                                cATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                                    >  1:678446/1‑76 (MQ=255)
                                cATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                                    <  2:42095/76‑1 (MQ=255)
                                cATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                                    <  1:994069/76‑1 (MQ=255)
                                cATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCaa                                    <  2:118432/76‑1 (MQ=255)
                                 aTACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAg                                   >  2:753365/1‑76 (MQ=255)
                                    ccTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTg                                >  1:557856/1‑76 (MQ=255)
                                    ccTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTg                                >  1:1550737/1‑76 (MQ=255)
                                       gCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGaa                             <  1:151567/76‑1 (MQ=255)
                                         gTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACt                           >  2:554945/1‑76 (MQ=255)
                                         gTCTGCACCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACt                           >  1:1105629/1‑76 (MQ=255)
                                          tCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTg                          <  1:1662875/76‑1 (MQ=255)
                                            tGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCg                        >  1:888917/1‑76 (MQ=255)
                                            tGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCg                        >  1:459133/1‑76 (MQ=255)
                                             gCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGt                       >  2:973327/1‑76 (MQ=255)
                                             gCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGt                       >  2:2070/1‑76 (MQ=255)
                                              cATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTg                      >  2:987771/1‑76 (MQ=255)
                                              cATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTg                      >  2:424281/1‑76 (MQ=255)
                                               aTCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTgg                     <  2:124545/76‑1 (MQ=255)
                                               aTCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTgg                     <  2:466956/76‑1 (MQ=255)
                                                 cAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTc                   >  1:1020459/1‑76 (MQ=255)
                                                   gTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCtt                 <  1:586315/76‑1 (MQ=255)
                                                   gTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCtt                 >  1:865924/1‑76 (MQ=255)
                                                       ggTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTa             <  2:181805/76‑1 (MQ=255)
                                                       ggTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTa             <  2:448095/76‑1 (MQ=255)
                                                       ggTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTa             >  2:722912/1‑76 (MQ=255)
                                                         tCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTacc           >  1:601117/1‑76 (MQ=255)
                                                           tGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTaccac         >  1:691825/1‑76 (MQ=255)
                                                            gCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTaccacc        >  2:168052/1‑76 (MQ=255)
                                                            gCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTaccacc        >  2:913443/1‑76 (MQ=255)
                                                             ccGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTAccacca       <  2:1204406/76‑1 (MQ=255)
                                                             ccGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACGGCGTGGTCTTGCTAccacca       <  1:143959/76‑1 (MQ=255)
                                                               gaagaaAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCAtt     <  1:956288/76‑1 (MQ=255)
                                                                 agaaAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTaa   <  2:976612/76‑1 (MQ=255)
                                                                  gaaAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTaaa  <  1:29976/76‑1 (MQ=255)
                                                                         |                                                                    
CTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAA  >  REL606/4166004‑4166145

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