Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,906,419 C→T 33.7% G382G (GGG→GGA treB ← fused trehalose(maltose)‑specific enzyme IIBC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,906,4190CT33.7% 112.1 / 74.5 136G382G (GGG→GGAtreBfused trehalose(maltose)‑specific enzyme IIBC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (48/42);  new base T (23/23);  total (71/65)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.12e-01

Only 100 of 152 total aligned reads displayed.

ATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAT  >  minE/2906357‑2906487
                                                              |                                                                    
atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                                  >  1:1838635/1‑67 (MQ=255)
atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                                  >  1:3048368/1‑67 (MQ=255)
atcGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCaa                                                                  >  1:666157/1‑67 (MQ=255)
  cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGcc                                                               >  1:1303382/1‑68 (MQ=255)
  cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                                              >  1:1736987/1‑69 (MQ=255)
  cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                                              >  1:328803/1‑69 (MQ=255)
  cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                                              >  1:1568588/1‑69 (MQ=255)
  cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                                              >  1:759515/1‑69 (MQ=255)
  cgcgCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                                              >  1:2317671/1‑69 (MQ=255)
  cgcgCACAGCATCGGGAAGAGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCa                                                              >  1:1019005/1‑69 (MQ=255)
      cacaGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCg                                                          <  1:1531533/69‑1 (MQ=255)
        caGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                                        <  1:1287043/69‑1 (MQ=255)
        caGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                                        <  1:3215992/69‑1 (MQ=255)
         aGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                                        <  1:874449/68‑1 (MQ=255)
          gCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                                        >  1:2778165/1‑67 (MQ=255)
          gCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                                        >  1:1210329/1‑67 (MQ=255)
           cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGAt                                                     >  1:1437456/1‑69 (MQ=255)
           cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGAt                                                     <  1:10319/69‑1 (MQ=255)
           cATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGAt                                                     >  1:1772158/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:2466910/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:2911731/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:3143102/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:2266500/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:1717678/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:402579/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:404390/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:849153/1‑69 (MQ=255)
            aTCGGGAAGCGATATTTCAGGTTGATGCCGTACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATa                                                    >  1:555006/1‑69 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                >  1:2794604/1‑54 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                <  1:107721/54‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                <  1:2357954/54‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                >  1:2176941/1‑54 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                <  1:2657429/54‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                >  1:2717927/1‑54 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                <  1:1372707/54‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                >  1:3277111/1‑54 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                <  1:534425/54‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                <  1:3418075/54‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGc                                                                <  1:3264554/54‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                                                 <  1:2569471/69‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                                                 <  1:3108588/69‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                                                 <  1:966132/69‑1 (MQ=255)
               gggAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCg                                                 <  1:1426276/69‑1 (MQ=255)
                        tatTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                                        <  1:1832175/53‑1 (MQ=255)
                        tatTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGga                                                        <  1:1245686/53‑1 (MQ=255)
                            tCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGgag                                     <  1:2854057/68‑1 (MQ=255)
                            tCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGgag                                     <  1:2936930/68‑1 (MQ=255)
                                gTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGAtct                                 <  1:1053174/68‑1 (MQ=255)
                                 ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg                               >  1:1750821/1‑69 (MQ=255)
                                 ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg                               >  1:2726201/1‑69 (MQ=255)
                                 ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg                               >  1:321772/1‑69 (MQ=255)
                                 ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg                               >  1:2237150/1‑69 (MQ=255)
                                 ttGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcg                               >  1:1452155/1‑69 (MQ=255)
                                  tGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTcgc                              >  1:2981947/1‑69 (MQ=255)
                                  tGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACCGAGATCTcgc                              >  1:1395337/1‑69 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:1860471/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:1786529/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:1876738/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:3068068/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:3347156/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:3350979/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:3371411/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:1615740/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:3423736/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:3431678/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:1520685/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:1449906/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:1990295/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:867190/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:2071896/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:2686109/69‑1 (MQ=255)
                                          tACATAGCAGGCTCAGTGACTCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATtg                      <  1:2507899/69‑1 (MQ=255)
                                            cATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTTATtgtg                    <  1:440542/69‑1 (MQ=255)
                                                gCAGGCTCTGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTg                >  1:554153/1‑69 (MQ=255)
                                                gCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCgg                                                >  1:497543/1‑37 (MQ=255)
                                                gCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTg                >  1:2262724/1‑69 (MQ=255)
                                                gCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTg                >  1:2075633/1‑69 (MQ=255)
                                                gCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTg                >  1:2168630/1‑69 (MQ=255)
                                                gCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTg                >  1:2190398/1‑69 (MQ=255)
                                                gCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTg                >  1:433398/1‑69 (MQ=255)
                                                gCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTg                >  1:2245403/1‑69 (MQ=255)
                                                   ggCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCg              <  1:2189114/68‑1 (MQ=255)
                                                       cAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTgg         <  1:249703/69‑1 (MQ=255)
                                                       cAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTgg         <  1:2354199/69‑1 (MQ=255)
                                                       cAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTgg         <  1:832290/69‑1 (MQ=255)
                                                           gACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAAt     >  1:491460/1‑69 (MQ=255)
                                                           gACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAAt     >  1:2015228/1‑69 (MQ=255)
                                                           gACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAAt     >  1:2993861/1‑69 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:1021776/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:141981/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:785360/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:652558/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:460464/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:1677560/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:1802614/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:2848181/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:2362394/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:2306194/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:2287183/69‑1 (MQ=255)
                                                              cccAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAt  <  1:2280088/69‑1 (MQ=255)
                                                              |                                                                    
ATCGCGCACAGCATCGGGAAGCGATATTTCAGGTTGATGCCGTACATTGCAGGCTCAGTGACCCCAAGCCAGGCGGAGATAGCGGCAGGCACGGAGATCTCGCGTTCATTGTGCTTGCGGCTGGAAATGAT  >  minE/2906357‑2906487

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: