Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,733,865 G→T 100% A301S (GCC→TCC)  pykF → pyruvate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,733,8650GT98.4% 268.7 / ‑3.1 80A301S (GCC→TCC) pykFpyruvate kinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (34/44);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCA  >  REL606/1733818‑1733912
                                               |                                               
tCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCa                                              >  1:6826482/1‑51 (MQ=255)
 cAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCAt                                             <  1:6400361/51‑1 (MQ=255)
  aaaaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATc                                            <  1:9489829/51‑1 (MQ=255)
  aaaaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATc                                            >  1:3815594/1‑51 (MQ=255)
   aaacCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATcc                                           <  1:5657942/51‑1 (MQ=255)
    aaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCt                                          <  1:9716894/51‑1 (MQ=255)
    aaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCt                                          >  1:2247746/1‑51 (MQ=255)
    aaaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCt                                          >  1:6287512/1‑51 (MQ=255)
     aaCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTc                                         <  1:8030950/51‑1 (MQ=255)
      aCCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCg                                        >  1:5110575/1‑51 (MQ=255)
       cccACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATcc                                           <  1:4295110/47‑1 (MQ=255)
       cccACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCg                                        <  1:5240969/50‑1 (MQ=255)
       cccACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCg                                        <  1:5826351/50‑1 (MQ=255)
         cACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACg                                     >  1:607228/1‑51 (MQ=255)
         cACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACg                                     >  1:1744650/1‑51 (MQ=255)
           cGCCTGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACgg                                    >  1:650149/1‑50 (MQ=255)
           cGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGt                                   >  1:8288424/1‑51 (MQ=255)
             cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATcc                                           <  1:299822/41‑1 (MQ=255)
             cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTc                                         <  1:6702808/43‑1 (MQ=255)
             cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGAc                                      <  1:7851910/46‑1 (MQ=255)
             cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGAc                                      <  1:929578/46‑1 (MQ=255)
             cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTAc                                 <  1:2264026/51‑1 (MQ=255)
             cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTAc                                 >  1:7339306/1‑51 (MQ=255)
             cccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTAc                                 >  1:8831426/1‑51 (MQ=255)
              ccGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACt                                >  1:9427028/1‑51 (MQ=255)
               cGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTg                               >  1:6960133/1‑51 (MQ=255)
               cGACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTg                               <  1:5971219/51‑1 (MQ=255)
                gACTCGCGCAGAAGCCGGTGACGTTGTAAACTCCATCCTCGACGGTACTg                               <  1:342501/50‑1 (MQ=255)
                gACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGa                              >  1:2733556/1‑51 (MQ=255)
                gACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGa                              >  1:533388/1‑51 (MQ=255)
                gACTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGa                              >  1:4515818/1‑51 (MQ=255)
                  cTCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTc                                         <  1:5266819/38‑1 (MQ=255)
                   tCGCGCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGc                           >  1:4384470/1‑51 (MQ=255)
                   tCGCGCAGAAGCCCGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGc                           >  1:3471668/1‑51 (MQ=38)
                    cgcgCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTg                               <  1:2704829/46‑1 (MQ=255)
                    cgcgCAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCa                          >  1:5645318/1‑51 (MQ=255)
                      cgcAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACg                                     <  1:5698053/38‑1 (MQ=255)
                      cgcAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAg                         <  1:5561716/50‑1 (MQ=255)
                       gcAGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTg                       <  1:392966/51‑1 (MQ=255)
                         aGAAGCCGGTGACGTTGCAAACTCCATCCTCGAc                                      <  1:5512781/34‑1 (MQ=255)
                         aGAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGAt                     >  1:4574225/1‑51 (MQ=255)
                          gAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGAc                             >  1:1127742/1‑42 (MQ=255)
                          gAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCGGTGATg                    <  1:44611/51‑1 (MQ=38)
                          gAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATg                    <  1:6648000/51‑1 (MQ=255)
                          gAAGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATg                    >  1:6336905/1‑51 (MQ=255)
                           aaGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGc                   >  1:6829710/1‑51 (MQ=255)
                           aaGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGc                   <  1:10716874/51‑1 (MQ=255)
                           aaGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGc                   >  1:2334082/1‑51 (MQ=255)
                           aaGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGAc                             <  1:3412989/41‑1 (MQ=255)
                            aGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCt                  <  1:4647055/51‑1 (MQ=255)
                            aGCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCt                  <  1:7920288/51‑1 (MQ=255)
                             gCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTg                 <  1:6401127/51‑1 (MQ=255)
                             gCCGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTg                 <  1:5157878/51‑1 (MQ=255)
                              ccGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGt                <  1:1228623/51‑1 (MQ=255)
                              ccGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGt                >  1:8269054/1‑51 (MQ=255)
                               cGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTc               >  1:2824997/1‑51 (MQ=255)
                               cGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTc               <  1:1185869/51‑1 (MQ=255)
                               cGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTc               <  1:4027422/51‑1 (MQ=255)
                               cGGTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTc               <  1:7905834/51‑1 (MQ=255)
                                ggTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt              <  1:2844923/51‑1 (MQ=255)
                                ggTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCt              >  1:4524827/1‑51 (MQ=255)
                                 gTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg             <  1:3371246/51‑1 (MQ=255)
                                 gTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg             <  1:4233493/51‑1 (MQ=255)
                                 gTGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTg             <  1:8156808/51‑1 (MQ=255)
                                  tGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg            >  1:405623/1‑51 (MQ=255)
                                  tGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg            >  1:823198/1‑51 (MQ=255)
                                  tGACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg            <  1:8697669/51‑1 (MQ=255)
                                   gACGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTgg            <  1:10571045/50‑1 (MQ=255)
                                     cGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGa         <  1:6335154/51‑1 (MQ=255)
                                     cGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGa         >  1:8619886/1‑51 (MQ=255)
                                     cGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGa         >  1:10797706/1‑51 (MQ=255)
                                     cGTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGa         >  1:792395/1‑51 (MQ=255)
                                      gTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTg          <  1:8588231/49‑1 (MQ=255)
                                      gTTGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGaa        <  1:2598317/51‑1 (MQ=255)
                                       ttGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAAt       >  1:9685009/1‑51 (MQ=255)
                                        tGCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATc      <  1:7038910/51‑1 (MQ=255)
                                         gCAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATcc     <  1:3264512/51‑1 (MQ=255)
                                          cAAACTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCg    <  1:9591359/51‑1 (MQ=255)
                                          cAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCg    >  1:9091936/1‑51 (MQ=255)
                                            aaCTCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCa  >  1:3605733/1‑51 (MQ=255)
                                               |                                               
TCAAAAACCCACGCCCGACTCGCGCAGAAGCCGGTGACGTTGCAAACGCCATCCTCGACGGTACTGACGCAGTGATGCTGTCTGGTGAATCCGCA  >  REL606/1733818‑1733912

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: