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

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,602,2970C→A13.5% 2.8 ‑4.4 0.045202ECB_01505conserved hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/6):  ref base (26/12):  total (27/18)
Fisher's exact test strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

ATCCATCAGGCTTGCCACATCCCAAACACGCGGATCAAAACCGCCCCACCACGGCATGTTAATCCGCTT  >  REL606/1602262‑1602330
                                   |                                 
aTCCATCAGGCTTGCCACATCCCAAACACGCGGATc                                   <  1:3480607/36‑1 (MQ=255)
aTCCATCAGGCTTGCCACATCCCAAACACGCGGATc                                   >  1:1431634/1‑36 (MQ=255)
aTCCATCAGGCTTGCCACATCCCAAACACGCGGATc                                   <  1:1519319/36‑1 (MQ=255)
 tCCATCATGCTTGCCACATCCCAAACACGCGGATCa                                  <  1:4197194/36‑1 (MQ=255)
 tCCATCAGGCTTGCCACATCCCAAACACGCGGATCa                                  >  1:3324830/1‑36 (MQ=255)
 tCCATCAGGCTTGCCACATCCCAAACACGCGGATCa                                  <  1:2038642/36‑1 (MQ=255)
  ccATCAGGCTTGCCACATCCCAAACACGCGGATCaa                                 <  1:6338566/36‑1 (MQ=255)
   cATCAGGCTTGCCACATCCCAAACACGCGGATCaaa                                <  1:2410101/36‑1 (MQ=255)
   cATCAGGCTTGCCACATCCCAAACACGCGGATCaaa                                <  1:1551963/36‑1 (MQ=255)
    aTCAGGCTTGCCACATCCCAAACACGCGGATCaaaa                               >  1:1676085/1‑36 (MQ=255)
    aTCAGGCTTGCCACATCCCAAACACGCGGATCaaaa                               >  1:8242632/1‑36 (MQ=255)
     tCAGGCTTGCCACATCCCAAACACGCGGATCAAAAc                              >  1:2915891/1‑36 (MQ=255)
      cAGGCTTGCCACATCCCAAACACGCGGATCAAAAcc                             <  1:701448/36‑1 (MQ=255)
      cAGGCTTGCCACATCCCAAACACGCGGATCAAAAcc                             >  1:7713391/1‑36 (MQ=255)
       aGGCTTGCCACATCCCAAACACGCTGATCAAAACCg                            <  1:1231864/36‑1 (MQ=255)
       aGGCTTGCCACATCCCAAACACGCGGATCAAAACCg                            <  1:3755665/36‑1 (MQ=255)
        ggCTTGCCACATCCCAAACACGCGGATCAAAACCGc                           <  1:3332768/36‑1 (MQ=255)
          cTTGCCACATCCCAAACACGCGGATCAAAACCGccc                         <  1:3684828/36‑1 (MQ=255)
          cTTGCCACATCCCAAACACGCGGATCAAAACCGccc                         >  1:6275261/1‑36 (MQ=255)
          cTTGCCACATCCCAAACACGCGGATCAAAACCGccc                         <  1:2478735/36‑1 (MQ=255)
          cTTGCCACATCCCAAACACGCGGATCAAAACCGccc                         >  1:7220255/1‑36 (MQ=255)
           ttGCCACATCCCAAACACGCGGATCAAAACCTcccc                        >  1:6433790/1‑36 (MQ=255)
           ttGCCACATCCCAAACACGCGGATCAAAACCGcccc                        >  1:6789362/1‑36 (MQ=255)
             gCCACATCCCAAACACGCGGATCAAAACCGCCccac                      >  1:3250562/1‑36 (MQ=255)
             gCCACATCCCAAACACGCGGATCAAAACCGCCccac                      >  1:1695121/1‑36 (MQ=255)
              ccACATCCCAAACACGCGGATCAAAACCGCCccacc                     >  1:4701041/1‑36 (MQ=255)
              ccACATCCCAAACACGCGGATCAAAACCGCCccacc                     >  1:2767496/1‑36 (MQ=255)
               cacaTCCCAAACACGCGGACCAAAACCGCCacacca                    >  1:7804933/1‑36 (MQ=18)
                acaTCCCAAACACGCGGATCAAAACCGCCCcaccac                   >  1:1673526/1‑36 (MQ=255)
                  cccccAAACACGCGGATAAAAACCGCCCCACCACgg                 <  1:4633022/34‑1 (MQ=18)
                  aTCCCAAACACGCGGATCAAAACCGCCCCACCACgg                 >  1:1148628/1‑36 (MQ=255)
                    ccaaaaCACGCGGAAAAAAACCGCCCCACCACGGCa               <  1:3951786/33‑1 (MQ=17)
                     ccAAACCCGCGGATCAAAACCGCCCCACCACGGCAt              <  1:1420795/36‑1 (MQ=255)
                     caaaaCACGCGGATAAAAACCGCCCCACCACGGCAt              <  1:8036184/34‑1 (MQ=17)
                      cAAACCCGCGTATAAAAACCGCCCCACCACGGCATg             <  1:5684927/36‑1 (MQ=14)
                       aaaCACGCGGAAAAAAACCGCCCCACCACGGCATGt            <  1:586978/36‑1 (MQ=18)
                         acacGCGGATCAAAACCGCCCCACCACGGCATGTTa          >  1:6319901/1‑36 (MQ=255)
                         acacGCGGATCAAAACCGCCCCACCACGGCATGTTa          >  1:5704858/1‑36 (MQ=255)
                          cacGCGGATCAAAACCGCCCCACCACGGCATGTTaa         >  1:6187329/1‑36 (MQ=255)
                          cacGCGGATCAAAACCGCCCCACCACGGCATGTTaa         >  1:3181411/1‑36 (MQ=255)
                          cacGCGAATAAAAACCGCCCCACCACGGCATGTTaa         <  1:3516665/36‑1 (MQ=18)
                            cgcgGATCAAAACCGCCCCACCACGGCATGTTAATc       >  1:6541641/1‑36 (MQ=255)
                             gcgGATCAAAACCGCCCCACCACGGCATGTTAATcc      >  1:3499224/1‑36 (MQ=255)
                             gcgGATCAAAACCGCCCCACCACGGCATGTTAATcc      >  1:2387598/1‑36 (MQ=255)
                             gcgGATCAAAACCGCCCCACCACGGCATGTTAATcc      >  1:128341/1‑36 (MQ=255)
                             gcgGATAAAAACCGCCCCAACACGGCATGTTAATcc      >  1:1451802/1‑36 (MQ=18)
                              cgGATCAAAACCGCCCCACCACGGCATGTTAATCCg     >  1:6020175/1‑36 (MQ=255)
                                 aacAAAACCGCCCCACCACGGCATGTTAATCCGCtt  <  1:6478277/34‑1 (MQ=255)
                                 aTCAAAACCGCCCCACCACGGCATGTTAATCCGCtt  >  1:2961897/1‑36 (MQ=255)
                                   |                                 
ATCCATCAGGCTTGCCACATCCCAAACACGCGGATCAAAACCGCCCCACCACGGCATGTTAATCCGCTT  >  REL606/1602262‑1602330

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG