Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,028,431 C→T 8.7% A462A (GCC→GCT malX → fused maltose and glucose‑specific PTS enzyme IIBC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,028,4310CT 8.7% 473.4 / 25.2 149A462A (GCC→GCTmalXfused maltose and glucose‑specific PTS enzyme IIBC components
Reads supporting (aligned to +/- strand):  ref base C (64/72);  new base T (5/8);  total (69/80)
Fisher's exact test for biased strand distribution p-value = 7.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

AAGCCGTTGCCGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTA  >  minE/1028361‑1028501
                                                                      |                                                                      
aaGCCGTTGCCGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGcc                                                                        <  1:4427491/71‑1 (MQ=255)
   ccGTTGCCGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGa                                                                      >  1:4019434/1‑70 (MQ=255)
    cGTTGCCGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGTAGCATTAGGCGGTGCCGAc                                                                     >  1:2029915/1‑70 (MQ=255)
    cGTTGCCGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGAc                                                                     <  1:5987200/70‑1 (MQ=255)
     gTTGCCGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACaa                                                                   <  1:1335531/71‑1 (MQ=255)
          cGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGAc                                                                     <  1:5874601/64‑1 (MQ=255)
          cGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATAtt                                                               <  1:6510940/70‑1 (MQ=255)
          cGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTg                                                              <  1:3820353/71‑1 (MQ=255)
               cgcCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGc                                                         >  1:3685251/1‑71 (MQ=255)
                  cGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTc                                                      <  1:260361/71‑1 (MQ=255)
                     gTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGAt                                                   <  1:2560371/71‑1 (MQ=255)
                     gTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGAt                                                   <  1:5511302/71‑1 (MQ=255)
                      tAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATa                                                  <  1:5377023/71‑1 (MQ=255)
                          tCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCa                                                           <  1:540391/58‑1 (MQ=255)
                           cAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCa                                                           <  1:5461682/57‑1 (MQ=255)
                           cAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTg                                              >  1:4319881/1‑70 (MQ=255)
                            aGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGt                                                             >  1:6270877/1‑54 (MQ=255)
                            aGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCa                                            >  1:5732304/1‑71 (MQ=255)
                            aGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCa                                            >  1:4355116/1‑71 (MQ=255)
                            aGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCa                                            >  1:4140490/1‑71 (MQ=255)
                            aGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCa                                            >  1:6835701/1‑71 (MQ=255)
                             ggTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGAt                                                   <  1:2914532/63‑1 (MQ=255)
                              gTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCTGACAATATTGTcg                                                           <  1:5853152/54‑2 (MQ=255)
                              gTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCTGACAATATTGTcg                                                           <  1:6406630/54‑2 (MQ=255)
                              gTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCAt                                           >  1:4803458/1‑70 (MQ=255)
                              gTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCAt                                           >  1:3113305/1‑70 (MQ=255)
                              gTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCAt                                           >  1:6684698/1‑70 (MQ=255)
                               ttACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTa                                         >  1:433672/1‑71 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTTCTGACAATATTGTCGGc                                                         >  1:5165247/1‑54 (MQ=38)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCTGACAATATTGTCGGc                                                         >  1:2353696/1‑54 (MQ=38)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCTGACAATATTGTCGGc                                                         >  1:3216399/1‑54 (MQ=38)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTg                                                              <  1:6227677/49‑1 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTg                                                              <  1:749109/49‑1 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTg                                                              <  1:1215072/49‑1 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTg                                                              <  1:6060288/49‑1 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTg                                                              <  1:905259/49‑1 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTg                                                              <  1:927941/49‑1 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTAc                                        >  1:1855641/1‑71 (MQ=255)
                                tACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTAc                                        >  1:4617914/1‑71 (MQ=255)
                                 aCAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTAc                                        <  1:3008345/70‑1 (MQ=255)
                                  cAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTAccc                                      <  1:6429878/71‑1 (MQ=255)
                                  cAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTAccc                                      <  1:35826/71‑1 (MQ=255)
                                  cAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTAccc                                      <  1:1543135/71‑1 (MQ=255)
                                   aaCGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCg                                     >  1:4613200/1‑71 (MQ=255)
                                   aaCGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCg                                     >  1:947537/1‑71 (MQ=255)
                                   aaCGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCg                                     <  1:1686350/71‑1 (MQ=255)
                                   aaCGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCg                                     <  1:3707705/71‑1 (MQ=255)
                                     cGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTc                                   >  1:4230145/1‑71 (MQ=255)
                                     cGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTc                                   >  1:3325763/1‑71 (MQ=255)
                                       ttCCTGCAATCCTCGAAGCATTAGGCGGTGCTGACAATATTGTcg                                                           <  1:5878059/45‑2 (MQ=255)
                                       ttCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTa                                         <  1:4020323/63‑1 (MQ=255)
                                        tCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTg                                 >  1:6661526/1‑70 (MQ=255)
                                        tCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGc                                >  1:6411985/1‑71 (MQ=255)
                                        tCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGc                                >  1:1737459/1‑71 (MQ=255)
                                         ccTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCg                               >  1:5728817/1‑71 (MQ=255)
                                          cTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGt                              <  1:1387410/71‑1 (MQ=255)
                                           tGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGc                                                         <  1:6830783/43‑1 (MQ=255)
                                           tGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTc                                                      <  1:6770652/46‑1 (MQ=255)
                                           tGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTc                                                      <  1:3437834/46‑1 (MQ=255)
                                           tGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTc                                                      <  1:4287221/46‑1 (MQ=255)
                                           tGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTc                                                      <  1:2743247/46‑1 (MQ=255)
                                           tGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCAt                                           <  1:3448304/57‑1 (MQ=255)
                                             cAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGttt                            <  1:2664668/70‑1 (MQ=255)
                                             cAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGtt                             >  1:2187764/1‑69 (MQ=255)
                                             cAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTg                           >  1:1137175/1‑71 (MQ=255)
                                             cAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTg                           >  1:1986012/1‑71 (MQ=255)
                                                  cTCGAATCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCtgt                      >  1:6463831/1‑71 (MQ=255)
                                                  cTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTc                                   >  1:3922164/1‑58 (MQ=255)
                                                   tCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGAt                                                   <  1:6223749/41‑1 (MQ=255)
                                                    cGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCtgtg                     >  1:3036638/1‑70 (MQ=255)
                                                     gAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGc                                             >  1:3595056/1‑45 (MQ=255)
                                                     gAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGc                                             >  1:4170024/1‑45 (MQ=255)
                                                     gAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGc                                             >  1:1210797/1‑45 (MQ=255)
                                                     gAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGc                                             >  1:4358631/1‑45 (MQ=255)
                                                     gAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGc                                             >  1:1482563/1‑45 (MQ=255)
                                                     gAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCtgtg                     >  1:6958982/1‑69 (MQ=255)
                                                       aGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATaa                                                 >  1:211335/1‑39 (MQ=255)
                                                       aGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATaa                                                 >  1:4463063/1‑39 (MQ=255)
                                                        gCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGa                <  1:3117986/71‑1 (MQ=255)
                                                         cATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGat               <  1:5428986/71‑1 (MQ=255)
                                                           ttAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGaaa                  <  1:6960505/66‑1 (MQ=255)
                                                           ttAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAg                 <  1:3217399/67‑1 (MQ=255)
                                                            tAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGa                <  1:2485902/67‑1 (MQ=255)
                                                             aGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGa                <  1:38725/66‑1 (MQ=255)
                                                              ggCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTc          <  1:5207011/71‑1 (MQ=255)
                                                              ggCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTc          <  1:4997951/71‑1 (MQ=255)
                                                               gCGGTGCTGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCg         <  1:5757540/71‑1 (MQ=255)
                                                               gCGGTGCTGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCg         <  1:3318034/71‑1 (MQ=255)
                                                               gCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCg                               <  1:1872745/49‑1 (MQ=255)
                                                               gCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCg                               <  1:72436/49‑1 (MQ=255)
                                                               gCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCg                               <  1:793033/49‑1 (MQ=255)
                                                               gCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCg                               <  1:801154/49‑1 (MQ=255)
                                                                 ggTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGtt                             >  1:2227463/1‑49 (MQ=255)
                                                                 ggTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCt       <  1:2281271/71‑1 (MQ=255)
                                                                    gCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGt    >  1:6405182/1‑71 (MQ=255)
                                                                    gCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGt    >  1:6310455/1‑71 (MQ=255)
                                                                    gCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGt    >  1:6030500/1‑71 (MQ=255)
                                                                    gCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGt    >  1:1079371/1‑71 (MQ=255)
                                                                     ccGACAATATTGTCAGCATCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGt    >  1:1936904/1‑70 (MQ=255)
                                                                      cGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTa  <  1:4125647/71‑1 (MQ=255)
                                                                      |                                                                      
AAGCCGTTGCCGGTGCGCCGGGTAAATCAGGTTACAACGTTCCTGCAATCCTCGAAGCATTAGGCGGTGCCGACAATATTGTCAGCCTCGATAACTGCATTACCCGTCTGCGTTTGTCTGTGAAAGATATGTCGCTTGTTA  >  minE/1028361‑1028501

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: