Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,424,370 G→T 22.6% D68Y (GAT→TAT)  ECB_01344 → conserved hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,424,3700G→T22.6% 4.6 ‑4.2 ‑0.017202ECB_01344conserved hypothetical protein
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (3/10):  total (7/10)
Fisher's exact test strand distribution p-value = 1.47e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.95e-01

AGCGAATGAATATGCCGTGGTGGGGCGGTTTTGATCCGCGCGCGTGGGATGTGGCAAGCCTGATGGA  >  REL606/1424338‑1424404
                                |                                  
aGCGAATGAATATGCCGTGGTGGGGCGGTTTTTATc                                 >  1:4084546/1‑36 (MQ=12)
aGCGAATGAATATGCCGTGGTGGGGCGGTTTTGATc                                 <  1:2565478/36‑1 (MQ=32)
      tgaatATGCCGTGGTGGGGCGGTTTTGATCcgcgcg                           <  1:4269357/36‑1 (MQ=255)
        aatATGCCGTGGTGGGGCGGTTTTGATCcgcgcgcg                         <  1:1779718/36‑1 (MQ=255)
              ccGTGGTGGGGCGGTTTTGATCCGCGCGCGTGGGAt                   <  1:2111839/36‑1 (MQ=255)
               cGTGGTGGGGCGGTTTTTATCCCCGCGCGGGGGGtg                  >  1:2228668/1‑36 (MQ=12)
                gtggtgGGGCGGTTTTGATCCGCGCGCGTGGGAtgt                 <  1:3890016/36‑1 (MQ=255)
                 tggtggGGCGGTTTTTCTCCTCGCGCGTGGGGtgtg                >  1:4563197/1‑36 (MQ=16)
                    tggGGCGGTTTTTTTCCCCGCCAGTGGGGTGTgggg             >  1:2474731/1‑34 (MQ=2)
                    tggGGCGGTTTTGATCCGCGCGCGTGGGATGTGGCa             <  1:3337927/36‑1 (MQ=255)
                    tggGGCGGTTTTGATCCGCGCGCGTGGGATGTGGCa             <  1:1596716/36‑1 (MQ=255)
                     ggggCGGTTTTGATCCGCGCGCGTGGGATGTGGCaa            <  1:2279540/36‑1 (MQ=255)
                        gCGGTTTTGATCCGCGCGCGTGGGATGTGGCAAGcc         <  1:2974516/36‑1 (MQ=255)
                        gCGGTTTTGATCCGCGCGCGTGGGATGTGGCAAGcc         <  1:449120/36‑1 (MQ=255)
                            ttttGATCCGCGCGCGTGGGATGTGGCAAGCCTGAt     >  1:3712081/1‑36 (MQ=255)
                              ttGATCCGCGCGCGTGGGATGTGGCAAGCCTGATgg   >  1:3079162/1‑36 (MQ=255)
                               tGATCCGCGCGCGTCGGATGTGGCAAGCCTGATGGa  >  1:2670679/1‑36 (MQ=255)
                                |                                  
AGCGAATGAATATGCCGTGGTGGGGCGGTTTTGATCCGCGCGCGTGGGATGTGGCAAGCCTGATGGA  >  REL606/1424338‑1424404

Base quality scores:  ATCG  < 3 ≤  ATCG  < 7 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 37 ≤  ATCG