Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,670 G→T 6.3% Y499* (TAC→TAA htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,6700GT 6.3% 371.5 / 5.2 107Y499* (TAC→TAAhtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base G (57/43);  new base T (4/3);  total (61/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 110 total aligned reads displayed.

CCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCAT  >  minE/138602‑138740
                                                                    |                                                                      
ccGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATTCGTGAGTa                                                                      <  1:3209585/71‑1 (MQ=255)
       cTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                                >  1:630053/1‑70 (MQ=255)
        ttttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCa                                                              >  1:3620740/1‑71 (MQ=255)
          ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                                <  1:479278/67‑1 (MQ=255)
          ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                                <  1:3159958/67‑1 (MQ=255)
          ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                            >  1:1435694/1‑71 (MQ=255)
           tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                                <  1:3165100/66‑1 (MQ=255)
           tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                            >  1:1757488/1‑70 (MQ=255)
           tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                            >  1:2420384/1‑70 (MQ=255)
            ttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATgg                                                           <  1:1177864/70‑1 (MQ=255)
            ttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATgg                                                           <  1:3049796/70‑1 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                     >  1:1004124/1‑71 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                     >  1:3719719/1‑71 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                     >  1:3094467/1‑71 (MQ=255)
                 aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                     >  1:3105421/1‑71 (MQ=255)
                    tttCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGtt                                                  >  1:3597348/1‑71 (MQ=255)
                    tttCAACGGTTGGTTAATACTGACCGTAACCTGATATTTCATGCGTGAGTAATAACGCATgg                                                           >  1:2151492/1‑62 (MQ=255)
                       cAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCa                                                              >  1:4112014/1‑56 (MQ=255)
                       cAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCa                                                <  1:1403787/70‑1 (MQ=255)
                        aaCGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                               <  1:2131423/70‑1 (MQ=255)
                          cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                     >  1:565893/1‑62 (MQ=255)
                          cGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                     >  1:2979259/1‑62 (MQ=255)
                           ggttggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                           <  1:1607005/71‑1 (MQ=255)
                            gttggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATgg                                                           <  1:2166734/54‑1 (MQ=255)
                              tggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                        >  1:2954923/1‑71 (MQ=255)
                              tggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                        >  1:1921383/1‑71 (MQ=255)
                               ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                        <  1:1868192/70‑1 (MQ=255)
                               ggttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTg                                        <  1:110519/70‑1 (MQ=255)
                                gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGtt                                      <  1:526344/71‑1 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAAt                                                                    >  1:3466460/1‑38 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                           >  1:155244/1‑63 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                           >  1:279796/1‑63 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                           >  1:383213/1‑63 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                   >  1:1623906/1‑71 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                   >  1:2334553/1‑71 (MQ=255)
                                   aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCg                                    <  1:713027/70‑1 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGAt                                            >  1:1416778/1‑61 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGAt                                            >  1:3366235/1‑61 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                   <  1:334817/70‑1 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                   <  1:2489645/70‑1 (MQ=255)
                                    aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                   <  1:1549132/70‑1 (MQ=255)
                                      aCTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                           <  1:4298124/60‑1 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                               >  1:228666/1‑55 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                               >  1:4057896/1‑55 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                           >  1:2110137/1‑59 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                           >  1:548630/1‑59 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                <  1:2195673/70‑1 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATg                               >  1:478779/1‑71 (MQ=255)
                                       cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATg                               >  1:2544496/1‑71 (MQ=255)
                                         gACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGtt                             >  1:1036854/1‑71 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                            <  1:3120584/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                            <  1:3365012/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                            <  1:758373/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                            <  1:3919660/71‑1 (MQ=255)
                                          aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                            <  1:1856634/71‑1 (MQ=255)
                                             gTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTc                                                 <  1:146049/47‑1 (MQ=255)
                                              tAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                        >  1:2555248/1‑71 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                >  1:1147256/1‑61 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                >  1:4023860/1‑61 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                           >  1:79171/1‑66 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                       >  1:747287/1‑70 (MQ=255)
                                                aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                       >  1:146286/1‑70 (MQ=255)
                                                    gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc                   >  1:4102601/1‑70 (MQ=255)
                                                     aTTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc                   <  1:4118998/69‑1 (MQ=255)
                                                      tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                               <  1:1033615/40‑1 (MQ=255)
                                                      tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                               <  1:1570073/40‑1 (MQ=255)
                                                        tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCtgg                                  <  1:3618505/51‑1 (MQ=255)
                                                        tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCggg                                  <  1:4000234/51‑1 (MQ=255)
                                                        tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAAtt               <  1:1475779/70‑1 (MQ=255)
                                                         ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa              <  1:1417112/70‑1 (MQ=255)
                                                         ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa              <  1:2923398/70‑1 (MQ=255)
                                                         ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa              <  1:2648552/70‑1 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                         >  1:2362988/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                         >  1:1660399/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                         >  1:1078763/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                         >  1:863713/1‑42 (MQ=255)
                                                          tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                         >  1:3215845/1‑42 (MQ=255)
                                                          ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                       <  1:406446/59‑1 (MQ=25)
                                                          ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                       <  1:4251833/59‑1 (MQ=25)
                                                          ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                       <  1:2730351/59‑1 (MQ=25)
                                                            aTGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt          <  1:1069039/71‑1 (MQ=255)
                                                            aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca                     >  1:2794561/1‑60 (MQ=37)
                                                            aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca                     >  1:1580655/1‑60 (MQ=37)
                                                            aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca                     >  1:3247666/1‑60 (MQ=37)
                                                             tGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt          >  1:65728/1‑70 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                           <  1:105678/49‑1 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                           <  1:783074/49‑1 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                           <  1:1743654/49‑1 (MQ=255)
                                                                 tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt          <  1:933718/66‑1 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGtt                             >  1:1269698/1‑45 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCAc                         >  1:3576342/1‑49 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCAc                         >  1:4424158/1‑49 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg       >  1:1790881/1‑67 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg       >  1:4249771/1‑67 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg       >  1:1865245/1‑67 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg       >  1:2711126/1‑67 (MQ=255)
                                                                   aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg       >  1:895935/1‑67 (MQ=255)
                                                                    gTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcat  >  1:4441266/1‑71 (MQ=255)
                                                                    gTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcat  >  1:2450160/1‑71 (MQ=255)
                                                                    gTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcat  >  1:3297332/1‑71 (MQ=255)
                                                                    |                                                                      
CCGTAATCTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCAT  >  minE/138602‑138740

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