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

GCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACAT  >  minE/1394792‑1394927
                                                                     |                                                                   
gcgcCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCg                                                                    >  1:3430849/1‑70 (MQ=255)
gcgcCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGa                                                                   <  1:4009994/71‑1 (MQ=255)
 cgcCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGaa                                                                  <  1:2321964/71‑1 (MQ=255)
    cAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCg                                                               >  1:379873/1‑71 (MQ=255)
     aGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCg                                                               <  1:3690740/71‑1 (MQ=255)
      gCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGcc                                                             <  1:1702335/71‑1 (MQ=255)
      gCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGcc                                                             <  1:2116708/71‑1 (MQ=255)
       cATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCct                                                            >  1:906502/1‑71 (MQ=255)
          gCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCctct                                                          >  1:2040592/1‑71 (MQ=255)
           ccATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCtcttc                                                        <  1:2150024/71‑1 (MQ=255)
           ccATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCtcttc                                                        <  1:2145236/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCtcttct                                                       <  1:1023612/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCtcttct                                                       <  1:2867687/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCtcttct                                                       <  1:1921657/71‑1 (MQ=255)
            cATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCtcttct                                                       <  1:1627056/71‑1 (MQ=255)
              ttGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCCTcttctt                                                      >  1:389881/1‑71 (MQ=255)
              ttGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAAt                                                                 <  1:3924367/59‑1 (MQ=255)
                gACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGt                                                   >  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)
                     aGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAAt                                                                 >  1:135369/1‑52 (MQ=255)
                     aGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCg                                              >  1:3674063/1‑71 (MQ=255)
                       ccTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGcc                                                             <  1:3402644/54‑1 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGCCTCTTCTTTGTCTCCGGcc                                           >  1:872940/1‑71 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCgga                                                                  <  1:2782013/47‑3 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCgga                                                                  <  1:203149/47‑3 (MQ=255)
                         ttCTGCTACATTACCGGCCTGCACCAGAGCCGCAGCGACCTGGCgga                                                                  <  1:1131341/47‑3 (MQ=39)
                          tCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGc                                                              >  1:2388243/1‑51 (MQ=39)
                          tCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGGAATCGc                                                              >  1:577833/1‑51 (MQ=39)
                           cTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCttt                                        >  1:992700/1‑71 (MQ=255)
                           cTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCtt                                         <  1:1992016/70‑1 (MQ=255)
                              cTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCttt                                        <  1:1972627/68‑1 (MQ=255)
                              cTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCttt                                        <  1:158506/68‑1 (MQ=255)
                                aCATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCttt                                                     >  1:3555977/1‑53 (MQ=255)
                                 cATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGccc                                  >  1:2350104/1‑71 (MQ=255)
                                 cATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGc                                    <  1:2788936/69‑1 (MQ=255)
                                 cATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGc                                    <  1:2302581/69‑1 (MQ=255)
                                  aTTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCttt                                        >  1:2301547/1‑64 (MQ=255)
                                  aTTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCg                                 <  1:3578222/71‑1 (MQ=255)
                                    tACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCgga                               >  1:2410680/1‑71 (MQ=255)
                                     aCCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCggat                              >  1:3936789/1‑71 (MQ=255)
                                      ccGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCggatg                             <  1:1413954/71‑1 (MQ=255)
                                       cGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCggatgg                            <  1:3677689/71‑1 (MQ=255)
                                          ccTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGgatggatg                         <  1:3912591/71‑1 (MQ=255)
                                          ccTGCACCAGCGCCGCAGCGACATGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGgatggatg                         <  1:3307176/71‑1 (MQ=255)
                                           cTGCACCAGCGCCGCAGCGACCTGGCGAATCGCct                                                            >  1:2973884/1‑35 (MQ=255)
                                               aCCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTg                    <  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)
                                                   gcgcCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACg                >  1:3732378/1‑71 (MQ=255)
                                                   gcgcCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACg                >  1:433586/1‑71 (MQ=255)
                                                     gccgcAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGa              >  1:2111573/1‑71 (MQ=255)
                                                      ccgcAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGAt             <  1:3028805/71‑1 (MQ=255)
                                                            cGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:3973457/71‑1 (MQ=255)
                                                            cGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:3974690/71‑1 (MQ=255)
                                                            cGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:297623/71‑1 (MQ=255)
                                                            cGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:270939/71‑1 (MQ=255)
                                                            cGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTgg       <  1:1018714/71‑1 (MQ=255)
                                                              aCCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGaa     <  1:2746775/71‑1 (MQ=255)
                                                                cTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACa   <  1:3312796/71‑1 (MQ=255)
                                                                cTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACa   <  1:424021/71‑1 (MQ=255)
                                                                  ggCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGaa     >  1:371874/1‑67 (MQ=255)
                                                                  ggCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGaa     >  1:3286524/1‑67 (MQ=255)
                                                                   gCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACAt  <  1:2075435/69‑1 (MQ=255)
                                                                     |                                                                   
GCGCCAGCATGCCATTGACGTAGCCTTCTGCTACATTACCGGCCTGCACCAGCGCCGCAGCGACCTGGCGAATCGCCTCTTCTTTGTCTCCGGCCTTTTCGCCCGGATGGATGTCCTGTACGGATAACTGGAACAT  >  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: