Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 4,126,699 Δ1 bp 50.0% 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→.50.0% 68.4 0.0 ‑0.132coding (335/879 nt)pflCpyruvate formate lyase II activase
Reads supporting (aligned to +/- strand):  new base (7/9):  ref base (6/10):  total (13/19)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.79e-01

TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCAC  >  REL606/4126664‑4126729
                                   |                              
tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                <  1:1746603/36‑1 (MQ=255)
 gACGCTCTGGAACGGGAAGTGATGAAAGATGACAtt                               <  1:2065987/36‑1 (MQ=255)
  aCGCTCTGGAACGGGAAGTGATGAAAGATGACAttt                              >  1:520035/1‑36 (MQ=255)
  aCGCTCTGGAACGGGAAGTGATGAAAGATGACAttt                              <  1:6454639/36‑1 (MQ=255)
  aCGCTCTGGAACGGGAAGTGATGAAAGATGACAttt                              <  1:3261712/36‑1 (MQ=255)
   ctctctGGAACGTGAAGTGATGAAAGATGACAtttt                             <  1:3595884/34‑1 (MQ=37)
   cGCTCTGGAACGGGAAGTGATGAAAGATGACAtttt                             <  1:2938953/36‑1 (MQ=255)
     ctctGGAACGGGAAGTGATGAAAGATGACAtttttt                           <  1:599275/36‑1 (MQ=255)
      tctGGAACGGGAAGTGATGAAAGATGACAttttttt                          >  1:1961950/1‑36 (MQ=255)
      tctGGAACGGGAAGTGATGAAAGATGACAttttttt                          <  1:6550066/36‑1 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttt                         <  1:2230789/36‑1 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                         <  1:3595215/36‑2 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                         <  1:7496953/36‑2 (MQ=255)
       atGGAACGGGAAGTGATGAAAGATGACAtttttttc                         <  1:5828018/35‑2 (MQ=255)
        tGGAACGGGAAGTGATGAAAGATGACAtttttttcg                        >  1:4536883/1‑34 (MQ=255)
        tGGAACGGGAAGTGATGAAAGATGACATTTTTTTTc                        >  1:3461036/1‑36 (MQ=255)
        tGGAACGGGAAGTGATGAAAGATGACATTTTTTTTc                        <  1:5332976/36‑1 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                       >  1:5023483/1‑33 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACATTTTTTTTCg                       >  1:3444797/1‑36 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACATTTTTTTTCg                       >  1:2495689/1‑36 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACATTTTTTTTCg                       <  1:6828263/36‑1 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACATTTTTTTTCg                       <  1:3869547/36‑1 (MQ=255)
           aaCGGGAAGTGATGAAAGATGACATTTTTTTTCGCa                     >  1:1131082/1‑36 (MQ=255)
            aCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCACg                   >  1:743772/1‑36 (MQ=37)
              gggAAGTGATGAAAGATGACATTTTTTTTCGCACGt                  <  1:5916686/36‑1 (MQ=255)
              gggAAGTGATGAAAGATGACATTTTTTTTCGCACGt                  >  1:559970/1‑36 (MQ=255)
              gggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTc                 >  1:7443733/1‑36 (MQ=37)
              gggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTc                 <  1:3865682/36‑1 (MQ=37)
               ggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTcc                <  1:1076321/36‑1 (MQ=37)
               ggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTcc                <  1:387996/36‑1 (MQ=37)
               ggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTcc                <  1:1378364/36‑1 (MQ=37)
                gAAGTGATGAAAGATGACATTTTTTTTCGCACGTcc                <  1:3613413/36‑1 (MQ=255)
                gAAGTGATGAAAGATGACATTTGTTTTCGCACGTcc                >  1:7882095/1‑36 (MQ=255)
                gAAGTGATGAAAGATGACA.TTTTTTTCGCACGTCcg               >  1:3758303/1‑36 (MQ=37)
                gAAGTGATGAAAGATGACA.TTTTTTTCGCACGTCcg               <  1:1402691/36‑1 (MQ=37)
                 aaGTGATGAAAGATGACA.TTTTTTTCGCACGTCcgg              <  1:3549102/36‑1 (MQ=37)
                  aGTGATGAAAGATGACA.TTTTTTTCGCACGTCcggc             <  1:6173274/36‑1 (MQ=38)
                  aGTGATGAAAGATGACA.TTTTTTTCGCACGTCcggc             >  1:7350348/1‑36 (MQ=38)
                   gTGATGAAAGATGACATTTTTTTTCGCACGTCcggc             <  1:6178784/36‑1 (MQ=255)
                    tgatgaAAGATGACA.TTTTTTTCGCACGTCcggcgg           >  1:7452722/1‑36 (MQ=38)
                     gatgaAAGATGCCATTTTTTTTCGCACGTCcggcgg           <  1:7154527/36‑1 (MQ=255)
                     gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc          <  1:5033960/36‑1 (MQ=38)
                     gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc          >  1:6172801/1‑36 (MQ=38)
                       tgaAAGATGACATTTTTTTTCGCACGTCcggcggcg         <  1:6527770/36‑1 (MQ=255)
                       tgaAAGATGACA.TTTTTTTCGCACGTCcggcggcgg        >  1:4590821/1‑36 (MQ=37)
                         aaaGATGACATTTTTTTTCGCTCGTCCggcggcggc       <  1:1772622/36‑1 (MQ=255)
                            gATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTc    <  1:6011346/36‑1 (MQ=255)
                             aTGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTTAc  <  1:5275054/36‑1 (MQ=21)
                                   |                              
TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCAC  >  REL606/4126664‑4126729

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG