Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,673 A→C 6.4% N498K (AAT→AAG) ‡ htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,6730AC 6.4% 298.2 / 16.3 110N498K (AAT→AAG) ‡htrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (60/43);  new base C (4/3);  total (64/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.24e-01

Only 100 of 119 total aligned reads displayed.

CTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCAT  >  minE/138609‑138743
                                                                |                                                                      
cTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                                   >  1:630053/1‑70 (MQ=255)
   ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                                   <  1:3159958/67‑1 (MQ=255)
   ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACg                                                                   <  1:479278/67‑1 (MQ=255)
   ttttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                               >  1:1435694/1‑71 (MQ=255)
    tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                               >  1:2420384/1‑70 (MQ=255)
    tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                               >  1:1757488/1‑70 (MQ=255)
    tttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATg                                                               >  1:4126721/1‑70 (MQ=255)
     ttCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATgg                                                              <  1:1177864/70‑1 (MQ=255)
          aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                        >  1:1004124/1‑71 (MQ=255)
          aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                        >  1:3094467/1‑71 (MQ=255)
          aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                        >  1:3105421/1‑71 (MQ=255)
          aaaTTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCg                                                        >  1:3719719/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:2979259/1‑62 (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)
                         gttAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGtt                                         <  1:526344/71‑1 (MQ=255)
                            aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAAt                                                                       >  1:3466460/1‑38 (MQ=255)
                            aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                              >  1:383213/1‑63 (MQ=255)
                            aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                              >  1:155244/1‑63 (MQ=255)
                            aaTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                      >  1:2924489/1‑71 (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)
                             aTACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCggg                                     <  1:1579855/71‑1 (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:4057896/1‑55 (MQ=255)
                                cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATc                                              >  1:548630/1‑59 (MQ=255)
                                cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                   <  1:2195673/70‑1 (MQ=255)
                                cTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATg                                  >  1:478779/1‑71 (MQ=255)
                                  gACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGtt                                >  1:1036854/1‑71 (MQ=255)
                                   aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                               <  1:3365012/71‑1 (MQ=255)
                                   aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                               <  1:3120584/71‑1 (MQ=255)
                                   aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                               <  1:321674/71‑1 (MQ=255)
                                   aCCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGttt                               <  1:3919660/71‑1 (MQ=255)
                                      gTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTc                                                    <  1:146049/47‑1 (MQ=255)
                                       tAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                           >  1:2555248/1‑71 (MQ=255)
                                       tAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACt                           >  1:3627984/1‑71 (MQ=255)
                                         aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                   >  1:1147256/1‑61 (MQ=255)
                                         aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                   >  1:4023860/1‑61 (MQ=255)
                                         aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGAt                                   >  1:1919919/1‑61 (MQ=255)
                                         aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                              >  1:79171/1‑66 (MQ=255)
                                         aCCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                          >  1:747287/1‑70 (MQ=255)
                                             gATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc                      >  1:4102601/1‑70 (MQ=255)
                                               tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                                  <  1:1033615/40‑1 (MQ=255)
                                               tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCaa                                                  <  1:1570073/40‑1 (MQ=255)
                                               tttttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtcatc                      <  1:244486/68‑1 (MQ=255)
                                                 tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCtgg                                     <  1:3618505/51‑1 (MQ=255)
                                                 tttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAAtt                  <  1:1475779/70‑1 (MQ=255)
                                                 atccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAAt                   >  1:616378/4‑69 (MQ=25)
                                                  ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa                 <  1:2923398/70‑1 (MQ=255)
                                                  ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa                 <  1:1417112/70‑1 (MQ=255)
                                                  ttCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTa                 <  1:2648552/70‑1 (MQ=255)
                                                   tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                            >  1:1660399/1‑42 (MQ=255)
                                                   tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                            >  1:2362988/1‑42 (MQ=255)
                                                   tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                            >  1:1078763/1‑42 (MQ=255)
                                                   tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                            >  1:3215845/1‑42 (MQ=255)
                                                   tCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCtt                                            >  1:863713/1‑42 (MQ=255)
                                                   ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                          <  1:2730351/59‑1 (MQ=25)
                                                   ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                          <  1:4251833/59‑1 (MQ=25)
                                                   ccATGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACtt                          <  1:406446/59‑1 (MQ=25)
                                                     aTGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt             <  1:1069039/71‑1 (MQ=255)
                                                     aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca                        >  1:3247666/1‑60 (MQ=37)
                                                     aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca                        >  1:2794561/1‑60 (MQ=37)
                                                     aTGCGTAATTACTCACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTtca                        >  1:1580655/1‑60 (MQ=37)
                                                      tGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGt             >  1:65728/1‑70 (MQ=255)
                                                          tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                              <  1:783074/49‑1 (MQ=255)
                                                          tGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTc                              <  1:105678/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:4424158/1‑49 (MQ=255)
                                                            aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCAc                            >  1:3576342/1‑49 (MQ=255)
                                                            aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg          >  1:4249771/1‑67 (MQ=255)
                                                            aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg          >  1:895935/1‑67 (MQ=255)
                                                            aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg          >  1:2711126/1‑67 (MQ=255)
                                                            aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg          >  1:1865245/1‑67 (MQ=255)
                                                            aGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAg          >  1:1790881/1‑67 (MQ=255)
                                                             gTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcat     >  1:2450160/1‑71 (MQ=255)
                                                             gTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcat     >  1:4441266/1‑71 (MQ=255)
                                                             gTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcat     >  1:3297332/1‑71 (MQ=255)
                                                              tAATTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                      <  1:2177230/37‑1 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                      >  1:481588/1‑35 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                      >  1:3706993/1‑35 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                      >  1:1118508/1‑35 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCgg                                      >  1:3524552/1‑35 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcatcat  >  1:2038196/1‑71 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcatcat  >  1:1255225/1‑71 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcatcat  >  1:2846594/1‑71 (MQ=255)
                                                                aTTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGcatca   <  1:3056988/70‑1 (MQ=255)
                                                                |                                                                      
CTTTTTTCTCAAATTTCAACGGTTGGTTAATACTGACCGTAACCTGATTTTTCATGCGTGAGTAATTACGCATGGATTTCGGTTCAAGATCTTGTTCGGGATGTTTCACTTCATCAATTAGAGTTAGTGCATCAT  >  minE/138609‑138743

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: