Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,126,699 Δ1 bp 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% 63.7 16coding (335/879 nt)pflCpyruvate formate lyase II activase
Reads supporting (aligned to +/- strand):  new base (6/10):  ref base (0/0):  total (6/10)

TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCACG  >  REL606/4126664‑4126730
                                   |                               
tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 <  2:1072133/36‑1 (MQ=255)
tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 <  1:75775/36‑1 (MQ=255)
tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 <  1:2693199/36‑1 (MQ=255)
  aCGCTCTGGAACGGGAAGTGATGAAAGATGACAttt                               <  1:755321/36‑1 (MQ=255)
   cGCTCTGGAACGGGAAGGGATGAAAGATGACAtttt                              <  2:3364878/36‑1 (MQ=255)
    gCTCTGGAACGGGAAGTGATGAAAGATGACAttttt                             <  1:313383/36‑1 (MQ=255)
    gCTCTAGAACGGGAAGTGATGAAAGATGACAttttt                             <  1:3315030/36‑1 (MQ=255)
     ctctGGAACGGGAAGTGATGAAAGATGACAtttttt                            <  2:2432221/36‑1 (MQ=255)
      tctGGAACGGGAAGTGATGAAAGATGACAttttttt                           >  2:2159331/1‑36 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                          <  2:392631/36‑2 (MQ=255)
        tGGAACGGGAAGTGATGAAAGATGACAtttttttcg                         <  2:1120800/36‑3 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        >  2:262501/1‑33 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  2:2286450/36‑4 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  1:1745955/36‑4 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  2:1595583/36‑4 (MQ=255)
           aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                     <  2:1483611/36‑1 (MQ=37)
             cGGGAAGTGATGAAAGATGACA.TTTTTTTCGCACGt                   <  1:3082321/36‑1 (MQ=37)
              gggAAGTGATGATAGATGACA.TTTTTTTCGCACGTc                  <  1:3090296/36‑1 (MQ=25)
               ggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTcc                 <  1:1567857/36‑1 (MQ=37)
                 aaGTGATGAAAGATGACA.TTTTTTTCGCACGTCcgg               <  2:3077957/36‑1 (MQ=37)
                 aaGTGATGAAAGATGACA.TTTTTTTCGCACGTCcgg               >  2:88739/1‑36 (MQ=37)
                 aaGTGATGAAAGATGACA.TTTTTATCGCACCTCcgg               >  1:2994927/1‑36 (MQ=18)
                  aGTGATGAAAGATGACA.TTTTTTTCGCACGTCcggc              >  2:2425229/1‑36 (MQ=38)
                   gTGATGAAAGATGACA.TTTTTTTCGCACGTCcggcg             <  1:3557922/36‑1 (MQ=38)
                    tgatgaAAGATGACA.TTTTTTTCGCACGTCcggcgg            <  2:2816324/36‑1 (MQ=38)
                     gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc           <  2:3751841/36‑1 (MQ=38)
                      atgaAAGATGACA.TTTTTTTCGCACGTCcgccggcg          >  1:3437022/1‑36 (MQ=37)
                           aGATGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTc     <  2:2524207/36‑1 (MQ=37)
                            gATGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCa    >  1:3676875/1‑36 (MQ=37)
                             aTGACA.TTTTTTTCGCACGTCCGGCGGAGGCGTCa  >  1:550706/1‑35 (MQ=25)
                              tGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCACg  <  2:1645332/36‑1 (MQ=37)
                                   |                               
TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCACG  >  REL606/4126664‑4126730

Base quality scores:  ATCG  < 3 ≤  ATCG  < 10 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG