Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,602,297 C→A 21.3% D68Y (GAT→TAT)  ECB_01505 ← conserved hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,602,2970C→A21.3% 4.9 ‑4.6 ‑0.032202ECB_01505conserved hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/6):  ref base (18/7):  total (19/13)
Fisher's exact test strand distribution p-value = 1.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.83e-01

CCATCAGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCGCCCCACCACGGCATGTTAATCCGCTTGC  >  REL606/1602264‑1602332
                                  |                                   
ccATCAGGCTTGCCACATCCCAAACA.CGCGGATCaa                                   <  1:3785803/36‑1 (MQ=255)
ccATCAGGCTTGCCACATCCCAAACA.CGCGGATCaa                                   <  1:7004702/36‑1 (MQ=255)
    cAGGCTTGCCACATCCCAAACA.CGCGGATCAAAAcc                               >  1:5452131/1‑36 (MQ=255)
    cAGGCTTGCCACATCCCAAACA.CGCGGATAAAAAcc                               >  1:5225868/1‑36 (MQ=255)
     aGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCg                              >  1:4332580/1‑36 (MQ=255)
     aGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCg                              <  1:3894868/36‑1 (MQ=255)
     aGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCg                              <  1:1577844/36‑1 (MQ=255)
      ggCTTGCCACATCCCAAACA.CGCGGATCAAAACCGc                             <  1:4938078/36‑1 (MQ=255)
        cTTGCCACATCCCAAACA.CGCGGATCAAAACCGccc                           >  1:6999892/1‑36 (MQ=255)
         ttGCCACATCCCAAACA.CGCGGATCAAAACCGcccc                          >  1:7068346/1‑36 (MQ=255)
          tGCCACATCCCAAACA.CGCGGATCAAAACCGCCcca                         >  1:3708941/1‑36 (MQ=255)
           gCCACATCCCAAACA.CGCGGATCAAAACCGCCacac                        >  1:7750290/1‑36 (MQ=255)
           gCCACACCCCAAACA.CGCGGATCAAAACCGCCccac                        <  1:2654754/36‑1 (MQ=255)
           gCCACACCCCAAACA.CGCGAATCAAAACCGCCccac                        <  1:397183/36‑1 (MQ=18)
             cacaTCCCAAACA.CGCGGATCAAAACCGCCccacca                      >  1:7931161/1‑36 (MQ=255)
             cacaTCCCAAACA.CGCGGATCAAAACCCCCccacca                      >  1:2122270/1‑36 (MQ=255)
               caTCCCAAACA.CGCGGATCAAAACCGCCCCACCACg                    >  1:4714562/1‑36 (MQ=255)
               caTCCCAAACA.CGCGGATCAAAACCGCCCAACCACg                    >  1:3237960/1‑36 (MQ=255)
                cacccAAAAACCGCGGAAAAAAACCGCCCCACCACg                    <  1:3567183/34‑1 (MQ=0)
                aTCCCAAACA.CGCGGATCAAAACCGCCCCACCACgg                   >  1:7496004/1‑36 (MQ=255)
                 tCCCAAACA.CGCGGATCAAAACCTTCCCACCACGGc                  >  1:731508/1‑36 (MQ=18)
                 tCCCAAACA.CGCGGATCAAAACCGCCCCACCACGGc                  >  1:7539352/1‑36 (MQ=255)
                  ccaaaaCA.CCCGAATAAAAACCGCCCCACCACGGCa                 <  1:6356746/33‑1 (MQ=9)
                   caaaaCA.CGCGAATAAAAACCGCCCCACCACGGCAt                <  1:2714389/34‑1 (MQ=17)
                      aaCA.CGCGGATCAAAACCGCCACAACACGGCATGtt             >  1:5573608/1‑36 (MQ=18)
                       aca.cGCGGATAAAAACCGCCCCACCACGGCATGTTa            <  1:6687436/36‑1 (MQ=255)
                        ca.cGCGAATAAAAACCGCCCCACCACGGCATGTTaa           <  1:4806137/36‑1 (MQ=18)
                           cgcgGATCAAAACCGCCCCACCACGGCATGTTAATc         >  1:5384374/1‑36 (MQ=255)
                           cgcgGATAAAAACCGCCCCACCACGGCATGTTAATc         <  1:4064734/36‑1 (MQ=255)
                             cgGATCAAAACCGCCCCACCACGGCATGTTAATCCg       >  1:3085316/1‑36 (MQ=255)
                               gATCAAAACCGCCCCACCACGGCATGTTAATCCGCt     >  1:240397/1‑36 (MQ=255)
                                aTCAAAACCGCCCCACCACGGCATGTTAATCCGCtt    >  1:225556/1‑36 (MQ=255)
                                 acAAACCCGCCCCACCACGGCATGTTAATCCGCTTg   <  1:3159940/35‑1 (MQ=255)
                                  cAAAACCGCCCCACCACGGCATGTTAATCCGCTTGc  >  1:1344541/1‑36 (MQ=255)
                                  |                                   
CCATCAGGCTTGCCACATCCCAAACA.CGCGGATCAAAACCGCCCCACCACGGCATGTTAATCCGCTTGC  >  REL606/1602264‑1602332

Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG