Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,126,699 Δ1 bp 100% coding (335/879 nt) pflC → pyruvate formate lyase II activase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6064,126,6990T→.100.0% 55.9 13coding (335/879 nt)pflCpyruvate formate lyase II activase
Reads supporting (aligned to +/- strand):  new base (5/8):  ref base (0/0):  total (5/8)

TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCAC  >  REL606/4126664‑4126729
                                   |                              
tGACACTCTGGAGCGGGAAGTGAGGAAAGATGACAt                                <  1:4304825/36‑1 (MQ=25)
 gACGCTCTGGAACGGGAAGTGATGAAAGATGACAtt                               >  1:4852031/1‑36 (MQ=255)
 gACGCTCTGGAACGGGAAGTGATGAAAGATGACAtt                               <  1:3627589/36‑1 (MQ=255)
 gACGCTCTGGAACGGGAAGTGATGAAAGATGACAtt                               <  1:3543729/36‑1 (MQ=255)
   cGCTCTGGAACGGGAAGTGATGAAAGATGACAtttt                             <  1:3233599/36‑1 (MQ=255)
    gCTCTGGAACGGGAAGTGATGAAAGATGACAttttt                            >  1:5353895/1‑36 (MQ=255)
    gCTCTGGAACGGGAAGTGATGAAAGATGACAttttt                            >  1:2660009/1‑36 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                         <  1:5450326/36‑2 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                         <  1:671992/36‑2 (MQ=255)
       atGGAACGGGAAGTGATGAAAGATGACAtttttttc                         <  1:1524866/35‑2 (MQ=255)
        tGGAACGGGAAGTGATGAAAGATGACAtttttttcg                        >  1:5187936/1‑34 (MQ=255)
        tGGAACGGGAAGTGATGAAAGATGACAtttttttcg                        <  1:5029554/36‑3 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                       >  1:2841610/1‑33 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                       <  1:5795175/36‑4 (MQ=255)
           aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                    <  1:3705891/36‑1 (MQ=37)
            aCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCACg                   <  1:880814/36‑1 (MQ=37)
             cGGGAAGTGATGAAAGATGACA.TTTTTTTCGCACGt                  <  1:1628722/36‑1 (MQ=37)
              gggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTc                 <  1:972522/36‑1 (MQ=37)
                  aGTGATGAAAGATGACA.TTTTTTTCGCACGTCcggc             >  1:4794768/1‑36 (MQ=38)
                  aGTGATGAAAGATGACA.TTTTTTTCGCACGTCcggc             <  1:5480687/36‑1 (MQ=38)
                     gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc          <  1:5936755/36‑1 (MQ=38)
                       tgaAAGATGACA.TTTTTTTCGCACGACcggcggcgg        >  1:2351055/1‑36 (MQ=37)
                        gaAAGATGACA.TTTTTTTCGCACGTCCggcggcggc       <  1:407742/36‑1 (MQ=37)
                        gaAAGATGACA.TTTTTTTCGCACGTCCggcggcggc       >  1:3554623/1‑36 (MQ=37)
                         aaaGATGACA.TTTTTTTCGCACGTCCGgcggcggcg      >  1:4493698/1‑36 (MQ=37)
                            gATGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCa   >  1:502687/1‑36 (MQ=37)
                             aTGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCAc  <  1:2506367/36‑1 (MQ=37)
                                   |                              
TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCAC  >  REL606/4126664‑4126729

Base quality scores:  ATCG  < 3 ≤  ATCG  < 12 ≤  ATCG  < 21 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG