Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,916,109 A→G 100% intergenic (+12/‑116) yebK → / → pykA predicted DNA‑binding transcriptional regulator/pyruvate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,916,1090AG99.0% 358.7 / ‑4.7 95intergenic (+12/‑116)yebK/pykApredicted DNA‑binding transcriptional regulator/pyruvate kinase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (42/52);  total (43/52)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01
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.

Only 100 of 110 total aligned reads displayed.

GATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCACTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCAT  >  REL606/1916062‑1916156
                                               |                                               
gATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctg                                              >  1:1668581/1‑51 (MQ=255)
gATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctg                                              <  1:4138029/51‑1 (MQ=255)
 aTAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             >  1:3403427/1‑51 (MQ=255)
 aTAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             >  1:6237519/1‑51 (MQ=255)
 aTAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             <  1:6952239/51‑1 (MQ=255)
  tAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgtc                                            <  1:6773472/51‑1 (MQ=255)
   aaGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGactgtcgctgt                                             <  1:3949991/49‑1 (MQ=255)
   aaGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTcc                                           >  1:6569683/1‑51 (MQ=255)
    aGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCt                                          >  1:9575563/1‑51 (MQ=255)
    aGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCt                                          <  1:3716047/51‑1 (MQ=255)
    aGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCt                                          <  1:3633398/51‑1 (MQ=255)
    aGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCt                                          <  1:7561244/51‑1 (MQ=255)
     gCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCt                                          >  1:9821735/1‑50 (MQ=255)
     gCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTa                                         <  1:9673708/51‑1 (MQ=255)
      cAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTaa                                        <  1:10324113/51‑1 (MQ=255)
      cAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTaa                                        <  1:8581902/51‑1 (MQ=255)
      cAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTaa                                        >  1:4610596/1‑51 (MQ=255)
      cAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTaa                                        >  1:6706043/1‑51 (MQ=255)
       aGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAAt                                       >  1:1184595/1‑51 (MQ=255)
         ttACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCt                                     >  1:1093415/1‑51 (MQ=255)
         ttACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCt                                     <  1:6449959/51‑1 (MQ=255)
         ttACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCt                                     >  1:9656207/1‑51 (MQ=255)
         ttACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCACTGTCCTAATCt                                     >  1:9272300/1‑51 (MQ=255)
           aCTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTta                                   >  1:817644/1‑51 (MQ=255)
           aCTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTta                                   <  1:3329374/51‑1 (MQ=255)
           aCTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTta                                   >  1:4568166/1‑51 (MQ=255)
           aCTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTta                                   <  1:7688372/51‑1 (MQ=255)
            cTTAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTtat                                  >  1:1572731/1‑51 (MQ=255)
             ttAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTtata                                 <  1:6886091/51‑1 (MQ=255)
             ttAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTtata                                 >  1:1197804/1‑51 (MQ=255)
              tAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATAc                                <  1:2873792/51‑1 (MQ=255)
              tAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATAc                                <  1:6764802/51‑1 (MQ=255)
              tAATTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATAc                                >  1:1680300/1‑51 (MQ=255)
               aaTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACg                               <  1:1905398/51‑1 (MQ=255)
               aaTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACg                               >  1:5229148/1‑51 (MQ=255)
               aaTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACg                               <  1:10802365/51‑1 (MQ=255)
                aTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGa                              >  1:3720037/1‑51 (MQ=255)
                aTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGa                              <  1:10393991/51‑1 (MQ=255)
                aTTTAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGa                              <  1:5640485/51‑1 (MQ=255)
                 tttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGAc                             <  1:4038471/51‑1 (MQ=255)
                 tttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGAc                             >  1:69981/1‑51 (MQ=255)
                 tttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGAc                             >  1:8548682/1‑51 (MQ=255)
                 tttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGAc                             <  1:9656730/51‑1 (MQ=255)
                  ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            <  1:5759295/51‑1 (MQ=255)
                  ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            >  1:6804774/1‑51 (MQ=255)
                  ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            <  1:91188/51‑1 (MQ=255)
                  ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            <  1:1185371/51‑1 (MQ=255)
                  ttAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACa                            >  1:6575930/1‑51 (MQ=255)
                   tAAGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACAt                           >  1:3379292/1‑51 (MQ=255)
                    aaGTGACGATCGCTAAAAACGACTGTCGCTGCCCTAATCTTATACGACATc                          >  1:9730264/1‑51 (MQ=255)
                     aGTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATcc                         >  1:8390246/1‑51 (MQ=255)
                      gTGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCg                        <  1:204315/51‑1 (MQ=255)
                       tGACGATCGCTAAAAACGGCTGTCGCTGTCCTAATCTTATACGACATCCg                        <  1:8506535/50‑1 (MQ=255)
                       tGACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGa                       <  1:9464222/51‑1 (MQ=255)
                        gACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGaa                      <  1:1611855/51‑1 (MQ=255)
                        gACGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGaa                      <  1:2781167/51‑1 (MQ=255)
                         aCGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAAt                     <  1:6731845/51‑1 (MQ=255)
                          cGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATg                    <  1:6856472/51‑1 (MQ=255)
                          cGATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATg                    >  1:7841369/1‑51 (MQ=255)
                           gATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATga                   <  1:192138/51‑1 (MQ=255)
                           gATCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATga                   <  1:6484999/51‑1 (MQ=255)
                             tCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATgaga                 >  1:2568732/1‑51 (MQ=255)
                             tCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATgaga                 >  1:5793092/1‑51 (MQ=255)
                             tCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATgaga                 >  1:7484050/1‑51 (MQ=255)
                             tCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATgaga                 >  1:9558783/1‑51 (MQ=255)
                             tCGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATgaga                 <  1:7735215/51‑1 (MQ=255)
                              cGCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGAt                <  1:8798696/51‑1 (MQ=255)
                               gCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGAtt               >  1:9913430/1‑51 (MQ=255)
                               gCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGAtt               <  1:2890752/51‑1 (MQ=255)
                               gCTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGAtt               <  1:2487856/51‑1 (MQ=255)
                                cTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTa              <  1:4269308/51‑1 (MQ=255)
                                cTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTa              <  1:10465932/51‑1 (MQ=255)
                                cTAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTa              <  1:9309927/51‑1 (MQ=255)
                                 tAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTaa             >  1:7563955/1‑51 (MQ=255)
                                 tAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTaa             >  1:2948718/1‑51 (MQ=255)
                                 tAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTaa             >  1:1830470/1‑51 (MQ=255)
                                 tAAAAACGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTaa             <  1:2078019/51‑1 (MQ=255)
                                  aaaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAt            >  1:7965322/1‑51 (MQ=255)
                                  aaaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAt            >  1:2354615/1‑51 (MQ=255)
                                  aaaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAt            >  1:8948568/1‑51 (MQ=255)
                                   aaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAtt           >  1:6162286/1‑51 (MQ=255)
                                   aaaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAtt           <  1:2371441/51‑1 (MQ=255)
                                    aaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAttt          >  1:3824374/1‑51 (MQ=255)
                                    aaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAttt          <  1:8673162/51‑1 (MQ=255)
                                    aaaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAAttt          <  1:7018357/51‑1 (MQ=255)
                                     aaCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTa         >  1:2085181/1‑51 (MQ=255)
                                      aCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAatatat        <  1:5267719/51‑1 (MQ=255)
                                      aCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTAt        >  1:8471206/1‑51 (MQ=255)
                                      aCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTAt        <  1:2405661/51‑1 (MQ=255)
                                      aCGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTAt        >  1:6336874/1‑51 (MQ=255)
                                       cGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATc       <  1:6583681/51‑1 (MQ=255)
                                       cGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATc       >  1:7110818/1‑51 (MQ=255)
                                       cGACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATc       <  1:9587407/51‑1 (MQ=255)
                                        gACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCg      <  1:3396759/51‑1 (MQ=255)
                                        gACTGTCGCTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCg      >  1:2388572/1‑51 (MQ=255)
                                         actgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGc     >  1:10641317/1‑51 (MQ=255)
                                          ctgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGcc    <  1:7241975/51‑1 (MQ=255)
                                           tgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTAt        <  1:6330077/46‑1 (MQ=255)
                                           tgtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCa   >  1:6944490/1‑51 (MQ=255)
                                            gtcgctgtcCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCAt  >  1:1471898/1‑51 (MQ=255)
                                               |                                               
GATAAGCAGTTACTTAATTTAAGTGACGATCGCTAAAAACGACTGTCACTGTCCTAATCTTATACGACATCCGAATGAGATTAATTTATCGCCAT  >  REL606/1916062‑1916156

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