Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,567,261 G→A 5.9% intergenic (‑210/+20) mtlD ← / ← mtlA mannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding/fused mannitol‑specific PTS enzyme IIABC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,567,2610GA 5.9% 412.3 / 10.3 118intergenic (‑210/+20)mtlD/mtlAmannitol‑1‑phosphate dehydrogenase, NAD(P)‑binding/fused mannitol‑specific PTS enzyme IIABC components
Reads supporting (aligned to +/- strand):  ref base G (40/71);  new base A (4/3);  total (44/74)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 123 total aligned reads displayed.

AACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGGTGGT  >  minE/2567191‑2567328
                                                                      |                                                                   
aacaaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtg                                                                     <  1:781704/71‑1 (MQ=255)
 acaaATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgt                                                                    >  1:891699/1‑71 (MQ=255)
   aaaTTTTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTgg                                                                  <  1:2956612/71‑1 (MQ=255)
    aaTTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGga                                                                 <  1:3006938/71‑1 (MQ=255)
       tgtgATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                              >  1:4409376/1‑71 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                           <  1:3590455/70‑1 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                           <  1:3242477/70‑1 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                           <  1:1176472/70‑1 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTg                                                          <  1:298036/71‑1 (MQ=255)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGGAGAAGGt                                                           <  1:3091785/70‑1 (MQ=21)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGGAGAAGGt                                                           <  1:4439536/70‑1 (MQ=21)
           aTGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCCTCTCCACATGGAGAAGGt                                                           <  1:3588592/70‑1 (MQ=21)
            tGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGt                                                           <  1:3330848/69‑1 (MQ=255)
            tGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTgg                                                         <  1:2445853/71‑1 (MQ=255)
            tGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTgg                                                         <  1:1884574/71‑1 (MQ=255)
            tGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTg                                                          <  1:428253/70‑1 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                   <  1:1973435/59‑1 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                   <  1:2793174/59‑1 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCAtgtg                                                                   <  1:4003552/59‑1 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTggg                                                        <  1:282487/70‑1 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGgga                                                       <  1:3358786/71‑1 (MQ=255)
              cATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGgga                                                       <  1:556498/71‑1 (MQ=255)
                tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTgg                                                                  <  1:3697883/58‑1 (MQ=255)
                tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTgg                                                                  <  1:4256397/58‑1 (MQ=255)
                tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTgg                                                                  <  1:866629/58‑1 (MQ=255)
                tCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTgg                                                                  <  1:1227673/58‑1 (MQ=255)
                 ccTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                                    >  1:226292/1‑71 (MQ=255)
                  cTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattgg                                                   <  1:1041713/71‑1 (MQ=255)
                  cTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattgg                                                   <  1:1604559/71‑1 (MQ=255)
                  cTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                                    >  1:1883784/1‑70 (MQ=255)
                  cTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGggattg                                                    >  1:3163063/1‑70 (MQ=255)
                        gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                              <  1:4592/54‑1 (MQ=255)
                        gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGagag                                                              <  1:544862/54‑1 (MQ=255)
                        gCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                             <  1:1808628/71‑1 (MQ=255)
                          cTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGga                                                                 <  1:1507705/49‑1 (MQ=255)
                          cTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGga                                                                 <  1:1983927/49‑1 (MQ=255)
                           tGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttactt                                            <  1:1666934/69‑1 (MQ=255)
                            gATAAGCGTAGCGCATCAGGCAATTAATCCCACCCTCTCCACATGgagaa                                                              >  1:977454/1‑49 (MQ=18)
                            gATAAGCGTAGCGCATCAGGCAATTAATCCCACCCTCTCCACATGGAGAagg                                                            >  1:4367440/1‑52 (MQ=17)
                            gATAAGCGTAGCGCATCAGGCAATTAATCACACCCTCTCCACATGGAGAagg                                                            >  1:3493153/1‑52 (MQ=17)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTg                                                          <  1:4012893/53‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttacttac                                          <  1:2134512/69‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttacttac                                          <  1:2823364/69‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttact                                             <  1:633058/66‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACg                                         <  1:1377683/70‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACg                                         <  1:1199479/70‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACg                                         <  1:2495068/70‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACg                                         <  1:3929811/70‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        <  1:1307737/71‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        <  1:626989/71‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        <  1:2284931/71‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        <  1:2847608/71‑1 (MQ=255)
                             aTAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        >  1:2738317/1‑71 (MQ=255)
                              tAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGCGAGGGTGGGATTGGATTACTTACGAc                                       <  1:558813/71‑1 (MQ=255)
                              tAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGAc                                       <  1:2293646/71‑1 (MQ=255)
                                aGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                            >  1:607228/1‑48 (MQ=255)
                                aGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAggg                                                            >  1:1040382/1‑48 (MQ=255)
                                             aGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGAttac                                              >  1:3594054/1‑49 (MQ=255)
                                             aGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGAc                                       >  1:3279558/1‑56 (MQ=255)
                                             aGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTcc                         <  1:3560561/70‑1 (MQ=255)
                                              ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTcc                         <  1:262613/69‑1 (MQ=255)
                                              ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                       <  1:2390378/71‑1 (MQ=255)
                                              ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                       <  1:927702/71‑1 (MQ=255)
                                              ggCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                       <  1:2212452/71‑1 (MQ=255)
                                               gCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                       <  1:61595/70‑1 (MQ=255)
                                                  aTTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACt                   >  1:2306916/1‑71 (MQ=255)
                                                     aaCCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCa                >  1:4344839/1‑71 (MQ=255)
                                                     aaCCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCa                >  1:2728458/1‑71 (MQ=255)
                                                      aCCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCAt               <  1:918401/71‑1 (MQ=255)
                                                        cccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTc                 <  1:4164105/67‑1 (MQ=255)
                                                         ccACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATcc             >  1:2539323/1‑70 (MQ=255)
                                                          cACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCAc           <  1:1849186/71‑1 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        >  1:2198380/1‑41 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        >  1:288898/1‑41 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGa                                        >  1:1064229/1‑41 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                       <  1:1770434/58‑1 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                       <  1:1198366/58‑1 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAg                       <  1:1547238/58‑1 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACg          <  1:2705384/71‑1 (MQ=255)
                                                           aCCTTCTCCATGTGGAGAGGGTGGGATTGATTACTTACGACCTGCcagc                               <  1:2191409/49‑1 (MQ=255)
                                                            ccTTCTCCATGTGGAGAGGGTGGGATTGGAttacttac                                          >  1:830441/1‑38 (MQ=255)
                                                            ccTTCTCCATGTGGAGAGGGTGGGATTGGAttactta                                           >  1:1029954/1‑37 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:3036982/1‑58 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:1621516/1‑58 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:864015/1‑58 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:637194/1‑58 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:1082365/1‑58 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACt                   <  1:1898190/60‑1 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCAc           <  1:1863717/68‑1 (MQ=255)
                                                             cTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCt        <  1:1195300/71‑1 (MQ=255)
                                                              ttCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCt                                     >  1:1362051/1‑41 (MQ=255)
                                                              ttCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCtg       >  1:4141888/1‑71 (MQ=255)
                                                               tctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:2680011/1‑56 (MQ=255)
                                                               tctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:3746820/1‑56 (MQ=255)
                                                               tctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:3821952/1‑56 (MQ=255)
                                                               tctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:4154295/1‑56 (MQ=255)
                                                                ctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:3418365/1‑55 (MQ=255)
                                                                ctcCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCa                     >  1:4029507/1‑55 (MQ=255)
                                                                   cATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTggtggt  <  1:3642522/71‑1 (MQ=255)
                                                                   cATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTggtggt  <  1:3208455/71‑1 (MQ=255)
                                                                      |                                                                   
AACAAATTGTGATGCATCCTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTAACCCCACCTTCTCCATGTGGAGAGGGTGGGATTGGATTACTTACGACCTGCCAGCAGTTCCAGCACTTCATCCACGCTGGTGGT  >  minE/2567191‑2567328

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