Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,733,649 A→C 3.9% N229H (AAC→CAC)  pykF → pyruvate kinase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,733,6490A→C 3.9% 3.5 ‑0.5 ‑0.1102685pykFpyruvate kinase
Reads supporting (aligned to +/- strand):  new base (1/3):  ref base (46/52):  total (47/55)
Fisher's exact test strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.45e-01

Only 100 of 110 total aligned reads displayed.

TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684
                                   |                                   
catttcatcgtccaacaaaatgggtgcagtacaaca                                     <  1:1006445‑M2/3‑1 (MQ=255)
catttcatcgtccaacaaaatgggtgcagtacaaca                                     <  1:7774606‑M2/3‑1 (MQ=255)
catttcatcgtccaacaaaatgggtgcagtacaaca                                     >  1:1327762‑M2/34‑36 (MQ=255)
 catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                    >  1:6862688‑M2/1‑35 (MQ=255)
 catcTCCAAAATCGAAAACCAGGAAGGCCTCAACAt                                    >  1:2248579‑M2/1‑35 (MQ=255)
 catcTCCAAAATAGAAAAACAGGAAAGCCTcaacat                                    >  1:7854007/1‑35 (MQ=21)
 atttcatcgtccaacaaaatgggtgcagtacaacaa                                    <  1:7612709‑M2/4‑1 (MQ=255)
 atttcatcgtcaaacaaaatgggtgcagtacaacaa                                    >  1:1667201‑M2/33‑36 (MQ=255)
   tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                  >  1:7193711‑M2/1‑33 (MQ=255)
   tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                  >  1:4958776‑M2/1‑33 (MQ=255)
   tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                  >  1:3962631‑M2/1‑33 (MQ=255)
   tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                  <  1:2265101‑M2/36‑4 (MQ=255)
   tctcCAAAATCGAAAACCAGGAAGGCCTCAACAtgt                                  >  1:1442553‑M2/1‑33 (MQ=255)
    tcatcgtccaacaaaatgggtgcagtacaACAACtt                                 <  1:4953085‑M2/7‑1 (MQ=255)
     tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                                <  1:7767779‑M2/36‑6 (MQ=255)
     tcCAAAATCGAAAACCAGGAAGGCCTCAACAtgtac                                >  1:4231114‑M2/1‑31 (MQ=255)
     catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:3330933‑M2/8‑1 (MQ=255)
     catcgtccaacaaaatgggtgcagtacaACAACTTc                                <  1:5314379‑M2/8‑1 (MQ=255)
      ccAAAATCGAAAACCAGGAAGGCCTCAACAtgtact                               <  1:7081885‑M2/36‑7 (MQ=255)
      atcgtccaacaaaatgggtgcagtacaACAACTTcg                               <  1:5830176‑M2/9‑1 (MQ=255)
       tcgtccaacaaaatgggtgcagtacaACAACTTcga                              >  1:1720457‑M2/27‑36 (MQ=255)
        cgtccaacaaaatgggtgcagtacaACAACTTcgac                             >  1:6665799‑M2/26‑36 (MQ=255)
        aaaaTCGAAAACCAGGAAGGCCTCAACAtgtactga                             <  1:214213‑M2/36‑9 (MQ=255)
        aaaaTCGAAAACCAGGAAGGCCTCAACAtgtactga                             <  1:6198348‑M2/36‑9 (MQ=255)
        aaaaTCGAAAACCAGGAAGGACTCAACAACTTcgac                             >  1:6366001/1‑36 (MQ=255)
         aaaTCGAAAACCAGGAAGGCCTCAACAtgtactgac                            <  1:886220‑M2/36‑10 (MQ=255)
          tccaacaaaatgggtgcagtacaACAACTTcgacga                           >  1:3592864‑M2/24‑36 (MQ=255)
           ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          <  1:3833232‑M2/14‑1 (MQ=255)
           ccaacaaaatgggtgcagtacaACAACTTCGACGaa                          <  1:5823721‑M2/14‑1 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                          >  1:848609‑M2/1‑25 (MQ=255)
           aTCGAAAACCAGGAAGGCCTCAACAtgtactgaccc                          >  1:598156‑M2/1‑25 (MQ=255)
            tCGAAAACCAGGAAGGCCTCAACAtgtactgacccc                         >  1:4552755‑M2/1‑24 (MQ=255)
            caacaaaatgggtgcagtacaACAACTTCGACGaaa                         <  1:7461775‑M2/15‑1 (MQ=255)
             cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                        <  1:7462026‑M2/36‑14 (MQ=255)
             cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                        <  1:1713749‑M2/36‑14 (MQ=255)
             cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                        <  1:1480102‑M2/36‑14 (MQ=255)
             cGAAAACCAGGAAGGCCTCAACAtgtactgacccca                        <  1:7154245‑M2/36‑14 (MQ=255)
             aacaaaatgggtgcagtacaACAACTTCGACGAAAt                        >  1:2695741‑M2/21‑36 (MQ=255)
             aacaaaatgggtgcagtacaACAACTTCGACGAAAt                        >  1:4214091‑M2/21‑36 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                       >  1:4310454‑M2/1‑22 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                       <  1:5025707‑M2/36‑15 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAtgtactgaccccaa                       >  1:6836706‑M2/1‑22 (MQ=255)
              gAAAACCAGGAAGGCCTCAACAtgcactgaccccaa                       <  1:1149210‑M2/36‑15 (MQ=255)
              acaaaatgggtgcagtacaACAACTTCGACGAAATc                       >  1:665546‑M2/20‑36 (MQ=255)
               caaaatgggtgcagtacaACAACTTCGACGAAATcc                      <  1:2564286‑M2/18‑1 (MQ=255)
                aaaatgggtgcagtacaACAACTTCGACGAAATCCt                     >  1:1282515‑M2/18‑36 (MQ=255)
                aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                     >  1:3335562‑M2/1‑20 (MQ=255)
                aaaCCAGGAAGGCCTCAACAtgtactgaccccaaaa                     >  1:5081604‑M2/1‑20 (MQ=255)
                 aaatgggtgcagtacaACAACTTCGACGAAatccac                    >  1:663926‑M2/17‑34 (MQ=255)
                 aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    <  1:5380533‑M2/20‑1 (MQ=255)
                 aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    >  1:7104671‑M2/17‑36 (MQ=255)
                 aaatgggtgcagtacaACAACTTCGACGAAATCCTc                    <  1:3173666‑M2/20‑1 (MQ=255)
                 aaCCAGGAAGGCCTCAACAtgtactgaccccaaaaa                    <  1:377793‑M2/36‑18 (MQ=255)
                  aatgggtgcagtacaACAACTTCGACGAAATCCTCg                   <  1:6507697‑M2/21‑1 (MQ=255)
                  aatgggtgcagtacaACAACTTCGACGAAATCCTCg                   >  1:1677531‑M2/16‑36 (MQ=255)
                   ccAGGAAGGCCTCAACAtgtactgaccccaaaaagt                  >  1:5067228‑M2/1‑17 (MQ=255)
                   ccAGGAAGGCCTCAACAtgtactgaccccaaaaagt                  >  1:3465141‑M2/1‑17 (MQ=255)
                    tgggtgaagtacaACAACTTCGACGAAATCCTCGaa                 <  1:6265446‑M2/23‑1 (MQ=255)
                    cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt                 <  1:2208025‑M2/36‑21 (MQ=255)
                    cAGGAAGGCCTCAACAtgtactgaccccaaaaagtt                 <  1:124213‑M2/36‑21 (MQ=255)
                     gggtgaagtacaACAACTTCGACGAAATCCTCGAAg                <  1:123758‑M2/24‑1 (MQ=255)
                     aGGAAGGCCTCAACAtgtactgaccccaaaaagttg                <  1:1740802‑M2/36‑22 (MQ=255)
                      ggtgcagtacaACAACTTCGACGAAATCCTCGAAGc               >  1:308886‑M2/12‑36 (MQ=255)
                      ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:1767742‑M2/36‑23 (MQ=255)
                      ggAAGGCCTCAACAtgtactgaccccaaaaagttgg               <  1:5638919‑M2/36‑23 (MQ=255)
                       gtgcagtacaACAACTTCGACGAAATCCTCGGAGcc              >  1:6310895‑M2/11‑36 (MQ=255)
                       gtgcagtacaACAACTTCGACGAAATCCTCGAAGcc              >  1:2864339‑M2/11‑36 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             <  1:363499‑M2/36‑25 (MQ=255)
                        aaGGCCTCAACAtgtactgaccccaaaaagttggac             <  1:7013053‑M2/36‑25 (MQ=255)
                         aGGCCTCAACAtgtactgaccccaaaaagttggaca            <  1:5791549‑M2/36‑26 (MQ=255)
                          ggCCTCAACAtgtactgaccccaaaaagttggacag           <  1:7355592‑M2/36‑27 (MQ=255)
                          ggCCTCAACAtgtactgaccccaaaaagttggacag           <  1:618903‑M2/36‑27 (MQ=255)
                          cagtacaACAACTTCGACGAAATCCTCGAAGCctct           >  1:5921742‑M2/8‑36 (MQ=255)
                           gCCTCAACAtgtactgaccccaaaaagttggacagt          >  1:4369694‑M2/1‑9 (MQ=255)
                           agtacaACAACTTCGACGAAATCCTCGAAGCCTCTg          >  1:3803752‑M2/7‑36 (MQ=255)
                            gtacaACAACTTCGACGAAATCCTCGAAGCCTCTGa         >  1:908203‑M2/6‑36 (MQ=255)
                            gtacaACAACTTCGACGAAAGCCTCGAAGCCTCTGa         <  1:2359024‑M2/31‑1 (MQ=255)
                            ccTCAACAtgtactgaccccaaaaagttggacagtt         >  1:5863688‑M2/1‑8 (MQ=255)
                             tacaACAACTTCGACGAAATCCTCGAAGCCTCTGAc        <  1:6711578‑M2/32‑1 (MQ=255)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        <  1:1024346‑M2/36‑30 (MQ=255)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        <  1:1498665‑M2/36‑30 (MQ=255)
                             cTCAACAtgtactgaccccaaaaagttggacagtta        <  1:5503595‑M2/36‑30 (MQ=255)
                              tCAACAtgtactgaccccaaaaagttggacagttaa       >  1:6819818‑M2/1‑6 (MQ=255)
                              tCAACAtgtactgaccccaaaaagttggacagttaa       <  1:1397773‑M2/36‑31 (MQ=255)
                              tCAACAtgtactgaccccaaaaagttggacagttaa       >  1:4499705‑M2/1‑6 (MQ=255)
                              gcaACCtgtactgaccccaaaaagttggacagttaa       <  1:2563011‑M2/35‑31 (MQ=255)
                              acaACAACTTCGACGAAATCCTCGAAGCCTCTGACg       <  1:7077003‑M2/33‑1 (MQ=255)
                                aaCCtgtactgaccccaaaaagttggacagttaaac     <  1:4863894‑M2/36‑33 (MQ=255)
                                aaCCtgtactgaccccaaaaagttggacagttaaac     <  1:2534997‑M2/36‑33 (MQ=255)
                                aaCAtgtactgaccccaaaaagttggacagttaaac     >  1:5722040‑M2/1‑4 (MQ=255)
                                aaCAtgtactgaccccaaaaagttggacagttaaac     >  1:3099488‑M2/1‑4 (MQ=255)
                                aaCAtgtactgaccccaaaaagttggacagttaaac     >  1:2237493‑M2/1‑4 (MQ=255)
                                aaCAtgtactgaccccaaaaagttggacagttaaac     >  1:912950‑M2/1‑4 (MQ=255)
                                aaCAtgtactgaccccaaaaaggtggacagttaaac     <  1:1195056‑M2/36‑33 (MQ=255)
                                aACAACTTCAACGAAATCCTCGAAGCCTCTGACGGc     <  1:2098290‑M2/35‑1 (MQ=255)
                                 ACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCa    >  1:5783622‑M2/1‑36 (MQ=255)
                                  cAtgtactgaccccaaaaagttggacagttaaacac   <  1:5175378‑M2/36‑35 (MQ=255)
                                  CAACTTCGACGAAATCCTCGAAGCCTCTGACGGCAt   <  1:1907612‑M2/36‑1 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  <  1:5958491‑M2/36‑36 (MQ=255)
                                   atgtactgaccccaaaaagttggacagttaaacacg  <  1:3796057‑M2/36‑36 (MQ=255)
                                   |                                   
TCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATC  >  REL606/1733614‑1733684

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG