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

TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCACG  >  REL606/4126664‑4126730
                                   |                               
tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 >  1:4300279/1‑36 (MQ=255)
tGACGCTCTGGAACGGGAAGTGATGAAAGATGACAt                                 <  1:2518060/36‑1 (MQ=255)
 gACGCTCTGGAACGGGAAGTGATGAAAGATGACAtt                                <  1:431566/36‑1 (MQ=255)
 gACGCTCTGGAACGGGAAGTGATGAAAGATGACAtt                                <  1:671461/36‑1 (MQ=255)
  aCGCTCTGGAACGGGAAGTGATGAAAGATGACAttt                               <  1:2180024/36‑1 (MQ=255)
  aCGCTCTGGAACGGGAAGTGATGAAAGATGACAttt                               >  1:5022615/1‑36 (MQ=255)
    gCTCTGGAACGGGAAGTGATGAAAGATGACAttttt                             <  1:7999624/36‑1 (MQ=255)
    gCTCTGGAACGGGAAGTGATGAAAGATGACAttttt                             <  1:4413845/36‑1 (MQ=255)
    gCTCTGGAACGGGAAGTGATGAAAGATGACAttttt                             >  1:2773244/1‑36 (MQ=255)
    gCTCGGGAACGGGAAGTGATGAAAGATGACAttttt                             <  1:2876117/36‑1 (MQ=255)
     ctctGGAACGGGAAGTGATGAAAGATGACAtttttt                            <  1:4064679/36‑1 (MQ=255)
     ctctGGAACGGGAAGTGATGAAAGATGACAtttttt                            <  1:2003858/36‑1 (MQ=255)
     ctctGGAACGGGAAGTGATGAAAGATGACAtttttt                            >  1:6545742/1‑36 (MQ=255)
     ctctGGAACGGGAAGTGATGAAAGATGACAtttttt                            >  1:1138854/1‑36 (MQ=255)
      tctGGAACGGGAAGTGATGAAAGATGACAttttttt                           >  1:4461363/1‑36 (MQ=255)
      tctGGAACGGGAAGTGATGAAAGATGACAttttttt                           <  1:983301/36‑1 (MQ=255)
      tctGGAACGGGAAGTGATGAAAGATGACAttttttt                           >  1:3782973/1‑36 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                          <  1:6059872/36‑2 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                          <  1:6316085/36‑2 (MQ=255)
       ctGGAACGGGAAGTGATGAAAGATGACAtttttttc                          >  1:2852575/1‑35 (MQ=255)
        tGGAACGGGAAGTGATGAAAGATGACAtttttttcg                         <  1:1146410/36‑3 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  1:4535230/36‑4 (MQ=255)
         ggAACGGGAAGTGATGAAAGATGACAtttttttcgc                        <  1:1177227/36‑4 (MQ=255)
           aaCGGGAAGTGATGCAAGATGACATTTTTTTTccac                      >  1:6140654/1‑33 (MQ=37)
           aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                     <  1:8175680/36‑1 (MQ=37)
           aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                     <  1:2882330/36‑1 (MQ=37)
           aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                     <  1:5881278/36‑1 (MQ=37)
           aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                     <  1:8269365/36‑1 (MQ=37)
           aaCGGGAAGTGATGAAAGATGACA.TTTTTTTCGCAc                     <  1:2536451/36‑1 (MQ=37)
             cGGGAAGTGATGAAAGATGACA.TTTTTTTCGCACGt                   <  1:6253046/36‑1 (MQ=37)
              gggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTc                  <  1:8019875/36‑1 (MQ=37)
               ggAAGTGATGAAAGATGACATTTTTTTTCGCACGTc                  >  1:4305580/1‑36 (MQ=255)
               ggAAGTGATGAAAGATGACA.TTTTTTTCGCACGTcc                 >  1:7682918/1‑36 (MQ=37)
                gAAGTGATGAAAGATGACATTTTTTTTCGCACGTcc                 <  1:3112103/36‑1 (MQ=255)
                gAAGTGATGAAAGATGACA.TTTTTTTCGCACGTCcg                >  1:8304976/1‑36 (MQ=37)
                gAAGTGATGAAAGATGACA.TTTTTTTCGCACGTCcg                >  1:2638955/1‑36 (MQ=37)
                 aaGTGATGAAAGATGACA.TTTTTTTCGCACGTCcgg               >  1:5965014/1‑36 (MQ=37)
                 aaGTGATGAAAGATGACA.TTTTTTTCGCACGTCcgg               >  1:2017323/1‑36 (MQ=37)
                  aGTGATGAAAGATGACA.TTTTTTTCGCACGTCcggc              >  1:5745003/1‑36 (MQ=38)
                   gTGATGAAAGATGACA.TTTTTTTCGCACGTCcggcg             >  1:854216/1‑36 (MQ=38)
                   gTGATGAAAGATGACA.TTTTTTTCGCACGTCcggcg             >  1:6574584/1‑36 (MQ=38)
                    tgatgaAAGATGACATTTTTTTTCGCACGTCcggcg             <  1:3953754/36‑1 (MQ=255)
                    tgatgaAAGATGACA.TTTTTTTCGCACGTCcggcgg            >  1:4976665/1‑36 (MQ=38)
                     gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc           <  1:6173813/36‑1 (MQ=38)
                     gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc           <  1:3028369/36‑1 (MQ=38)
                     gatgaAAGATGACA.TTTTTTTCGCACGTCcggcggc           >  1:299400/1‑36 (MQ=38)
                      atgaAAGATGACA.TTTTTTTCGCACGTCcggcggcg          <  1:8305460/36‑1 (MQ=37)
                      atgaAAGATGACA.TGTTTTTCGCACGTcc.gcggcgg         >  1:1752049/1‑36 (MQ=16)
                       tgaAAGATGACA.TTTTTTTCGCACGTCcggcggcgg         >  1:5080155/1‑36 (MQ=37)
                       tgaAAGATGACA.TTTTTTTCGCACGTCcggcggcgg         <  1:4541853/36‑1 (MQ=37)
                        gaAAGATGACA.TTTTTTTCGCACGTCCggcggcggc        >  1:8011827/1‑36 (MQ=37)
                        gaAAGATGACA.TTTTTTTCGCACGTCCggcggcggc        <  1:5313490/36‑1 (MQ=37)
                        gaAAGATGACA.TTTTTTTCGCACGTCCggcggcggc        <  1:3489340/36‑1 (MQ=37)
                         aaaGATGACA.TTTTTTTCGCACGTCCGgcggcggcg       >  1:3576688/1‑36 (MQ=37)
                         aaaGATGACA.TTTTTTTCGCACGTCCGgcggcggcg       >  1:993065/1‑36 (MQ=37)
                          aaGATGACA.TTTTTTTCGCACGTCCGGCGGCGGCGt      >  1:2096946/1‑36 (MQ=37)
                           aGATGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTc     <  1:1875312/36‑1 (MQ=37)
                             aTGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCAc   <  1:7672837/36‑1 (MQ=37)
                             aTGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCAc   <  1:6851444/36‑1 (MQ=37)
                             aTGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCAc   <  1:7936145/36‑1 (MQ=37)
                             aTGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCAc   <  1:6580625/36‑1 (MQ=37)
                              tGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCACg  <  1:1570603/36‑1 (MQ=37)
                              tGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCACg  <  1:229000/36‑1 (MQ=37)
                              tGACA.TTTTTTTCGCACGTCCGGCGGCGGCGTCACg  >  1:2843861/1‑36 (MQ=37)
                                   |                               
TGACGCTCTGGAACGGGAAGTGATGAAAGATGACATTTTTTTTCGCACGTCCGGCGGCGGCGTCACG  >  REL606/4126664‑4126730

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 37 ≤  ATCG