Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,394,860:1 +G 19.4% coding (68/1131 nt) fruB ← fused fructose‑specific PTS enzyme IIA component and HPr component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,394,8601.G19.4% 134.5 / 64.6 62coding (68/1131 nt)fruBfused fructose‑specific PTS enzyme IIA component and HPr component
Reads supporting (aligned to +/- strand):  ref base . (17/33);  new base G (8/4);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 5.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

GCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACAT  >  minE/1394792‑1394927
                                                                     |                                                                   
gcgcCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑g                                                                    >  1:3430849/1‑70 (MQ=255)
gcgcCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑Ga                                                                   <  1:4009994/71‑1 (MQ=255)
 cgcCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑Gaa                                                                  <  1:2321964/71‑1 (MQ=255)
    cAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCg                                                               >  1:379873/1‑71 (MQ=255)
     aGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCg                                                               <  1:3690740/71‑1 (MQ=255)
      gCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGcc                                                             <  1:1702335/71‑1 (MQ=255)
      gCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGcc                                                             <  1:2116708/71‑1 (MQ=255)
       cATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCct                                                            >  1:906502/1‑71 (MQ=255)
          gCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCctct                                                          >  1:2040592/1‑71 (MQ=255)
           ccATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCtcttc                                                        <  1:2150024/71‑1 (MQ=255)
           ccATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCtcttc                                                        <  1:2145236/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCtcttct                                                       <  1:1023612/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCtcttct                                                       <  1:2867687/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCtcttct                                                       <  1:1921657/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCtcttct                                                       <  1:1627056/71‑1 (MQ=255)
              ttGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCCTcttctt                                                      >  1:389881/1‑71 (MQ=255)
              ttGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAAt                                                                 <  1:3924367/59‑1 (MQ=255)
                gACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGt                                                   >  1:1372039/1‑71 (MQ=255)
                     aGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTcc                                               >  1:3871869/1‑71 (MQ=255)
                     aGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTcc                                               >  1:2218867/1‑71 (MQ=255)
                     aGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTcc                                               >  1:3460295/1‑71 (MQ=255)
                     aGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAAt                                                                 >  1:135369/1‑52 (MQ=255)
                     aGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCg                                              >  1:3674063/1‑71 (MQ=255)
                       ccTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGcc                                                             <  1:3402644/54‑1 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTCCGGcc                                           >  1:872940/1‑71 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑gga                                                                  <  1:2782013/47‑3 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑gga                                                                  <  1:203149/47‑3 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGAGCCGCAGCGACCTGGC‑gga                                                                  <  1:1131341/47‑3 (MQ=39)
                          tCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGc                                                              >  1:2388243/1‑51 (MQ=39)
                          tCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGc                                                              >  1:577833/1‑51 (MQ=39)
                           cTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCttt                                        >  1:992700/1‑71 (MQ=255)
                           cTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCtt                                         <  1:1992016/70‑1 (MQ=255)
                              cTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCttt                                        <  1:1972627/68‑1 (MQ=255)
                              cTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCttt                                        <  1:158506/68‑1 (MQ=255)
                                aCATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCttt                                                     >  1:3555977/1‑53 (MQ=255)
                                 cATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGccc                                  >  1:2350104/1‑71 (MQ=255)
                                 cATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGc                                    <  1:2788936/69‑1 (MQ=255)
                                 cATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGc                                    <  1:2302581/69‑1 (MQ=255)
                                  aTTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCttt                                        >  1:2301547/1‑64 (MQ=255)
                                  aTTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCg                                 <  1:3578222/71‑1 (MQ=255)
                                    tACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCgga                               >  1:2410680/1‑71 (MQ=255)
                                     aCCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCggat                              >  1:3936789/1‑71 (MQ=255)
                                      ccGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCggatg                             <  1:1413954/71‑1 (MQ=255)
                                       cGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCggatgg                            <  1:3677689/71‑1 (MQ=255)
                                          ccTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGgatggatg                         <  1:3912591/71‑1 (MQ=255)
                                          ccTGCACCAGCGCCGCAGCGACATGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGgatggatg                         <  1:3307176/71‑1 (MQ=255)
                                           cTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCct                                                            >  1:2973884/1‑35 (MQ=255)
                                               aCCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTg                    <  1:2340820/71‑1 (MQ=255)
                                                ccAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTCCGGc                                            <  1:2038121/47‑1 (MQ=38)
                                                ccAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTCCGGc                                            <  1:529985/47‑1 (MQ=38)
                                                ccAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTCCGGc                                            <  1:1251377/47‑1 (MQ=38)
                                                   gcgcCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACg                >  1:3732378/1‑71 (MQ=255)
                                                   gcgcCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACg                >  1:433586/1‑71 (MQ=255)
                                                     gccgcAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGa              >  1:2111573/1‑71 (MQ=255)
                                                      ccgcAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGAt             <  1:3028805/71‑1 (MQ=255)
                                                            cGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:3973457/71‑1 (MQ=255)
                                                            cGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:3974690/71‑1 (MQ=255)
                                                            cGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:297623/71‑1 (MQ=255)
                                                            cGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:270939/71‑1 (MQ=255)
                                                            cGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:1018714/71‑1 (MQ=255)
                                                              aCCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGaa     <  1:2746775/71‑1 (MQ=255)
                                                                cTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACa   <  1:3312796/71‑1 (MQ=255)
                                                                cTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACa   <  1:424021/71‑1 (MQ=255)
                                                                  ggC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGaa     >  1:371874/1‑67 (MQ=255)
                                                                  ggC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGaa     >  1:3286524/1‑67 (MQ=255)
                                                                   gC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACAt  <  1:2075435/69‑1 (MQ=255)
                                                                     |                                                                   
GCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGC‑GAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACAT  >  minE/1394792‑1394927

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑