Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,602,606 T→C 15.2% intergenic (+142/+98) yfeW → / ← yfeX predicted periplasmic esterase/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,602,6060TC15.2% 210.3 / 27.5 92intergenic (+142/+98)yfeW/yfeXpredicted periplasmic esterase/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/40);  new base C (8/6);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 102 total aligned reads displayed.

ACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGC  >  minE/1602536‑1602674
                                                                      |                                                                    
gcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGtt                                                                      <  1:2019165‑M1/22‑1 (MQ=255)
gcaaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGtt                                                                      <  1:3459302‑M1/22‑1 (MQ=255)
 cGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                      <  1:3969903/70‑1 (MQ=255)
  aaatgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCa                                                                    <  1:4205089‑M1/24‑1 (MQ=255)
    atgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCTTTCAcg                                                                  >  1:9523‑M1/46‑71 (MQ=255)
    atgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCAcg                                                                  >  1:2084738‑M1/46‑71 (MQ=255)
    atgcctgatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCAc                                                                   >  1:2198300‑M1/46‑70 (MQ=255)
    atgcctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCAcg                                                                  <  1:5272066‑M1/26‑1 (MQ=255)
     tgcctgatgcgctacgcttatcaggcctaaagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCAcgc                                                                 <  1:711982‑M1/27‑1 (MQ=255)
       tCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccg                                                               >  1:4105072/1‑71 (MQ=255)
       cctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCAcgccg                                                               >  1:4313221‑M1/43‑71 (MQ=255)
       cctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCAcgccg                                                               >  1:5410137‑M1/43‑71 (MQ=255)
       cctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCAcgccg                                                               >  1:77811‑M1/43‑71 (MQ=255)
        ctgatgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCAcgccgc                                                              >  1:2867373‑M1/42‑71 (MQ=255)
        cAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccg                                                               >  1:4164196/1‑70 (MQ=255)
        cAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccg                                                               >  1:3341650/1‑70 (MQ=255)
         ctgatgcgctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                             >  1:4208134‑M1/41‑71 (MQ=255)
         ctgatgcgctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                             >  1:4201390‑M1/41‑71 (MQ=255)
         aGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgccgc                                                              <  1:5149331/70‑1 (MQ=255)
          gatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                            <  1:3279022‑M1/32‑1 (MQ=255)
          gatgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                            <  1:3601042‑M1/32‑1 (MQ=255)
           atgcgctacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                           >  1:1533701‑M1/39‑71 (MQ=255)
           atgcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCa                                                             >  1:884991‑M1/39‑69 (MQ=255)
            ccTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                          >  1:5105273/1‑71 (MQ=255)
             gcgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCg                                                         <  1:534969‑M1/35‑1 (MQ=255)
              tACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                      <  1:2807604/57‑1 (MQ=255)
              tACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                        >  1:4429880/1‑71 (MQ=255)
              cgctacgcttatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCgg                                                        >  1:3046333‑M1/36‑71 (MQ=255)
               cgctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                        <  1:3979161‑M1/36‑1 (MQ=255)
                gctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGgca                                                      >  1:2570550‑M1/34‑70 (MQ=255)
                gctacgcttatcaggcctacaattctcgcgccaCTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                       >  1:5620472‑M1/34‑70 (MQ=255)
                cGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                                    >  1:2535100/1‑57 (MQ=255)
                cGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                        >  1:4481191/1‑69 (MQ=255)
                 tacgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCggcaa                                                     <  1:4788693‑M1/39‑3 (MQ=255)
                  tGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtg                                                    <  1:536118/71‑1 (MQ=255)
                  tGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtg                                                    <  1:5593013/71‑1 (MQ=255)
                   cgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGgca                                                      <  1:73329‑M1/38‑2 (MQ=255)
                   cgcttatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                       <  1:3955404‑M1/37‑1 (MQ=255)
                    gAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCgtgt                                                   >  1:2940791/1‑70 (MQ=255)
                      ttatcaggcctacaaattcccgcgccaTCCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                           <  1:2807260‑M1/33‑1 (MQ=255)
                      ttatcagacctacaaattctcacgccaTCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                        >  1:3392862‑M1/28‑63 (MQ=255)
                       tatcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAcat                                               <  1:768055‑M1/45‑2 (MQ=255)
                       tatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATcc                                                          >  1:1181507‑M1/27‑60 (MQ=255)
                       tatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCGgca                                                      <  1:1179752‑M1/38‑2 (MQ=255)
                       tatcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCGgca                                                      <  1:591948‑M1/38‑2 (MQ=255)
                       tatcaggcctacaaattcccgcgccaTCCGTAGGCCGGATAAGGCGTTCAcgccgc                                                              >  1:2548375‑M1/27‑56 (MQ=255)
                         tcaggcctacagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAacatg                                              >  1:769987‑M1/25‑68 (MQ=255)
                         tcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCg                                                         <  1:4204713‑M1/35‑1 (MQ=255)
                         tcaggcctacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATCCg                                                         <  1:2144493‑M1/35‑1 (MQ=255)
                         tcaggcctacaaattcccgcaccaTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                       >  1:4522534‑M1/25‑61 (MQ=255)
                         tcaggcctacaaattcccgcaccaTCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                       >  1:3267046‑M1/25‑61 (MQ=255)
                           caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                      >  1:1857770/1‑44 (MQ=38)
                           caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                      >  1:4214999/1‑44 (MQ=38)
                           caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                      >  1:841473/1‑44 (MQ=38)
                           caATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGtt                                                                      >  1:5003444/1‑44 (MQ=38)
                             atatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTc                                                                     >  1:2608166/1‑43 (MQ=255)
                              tatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTc                                                                     >  1:2734752/1‑42 (MQ=255)
                              tatatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCAc                                        >  1:5068593/1‑71 (MQ=255)
                               ctacaaattctcacgccaTCCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                           <  1:1625640‑M1/33‑1 (MQ=255)
                                tatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                                    >  1:1222205/1‑41 (MQ=255)
                                tatTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                             <  1:1007690/48‑1 (MQ=255)
                                tacagattctcgcgccaCTCGTAGGCCGGATAAGGCGTCCACGCCGCATc                                                           >  1:2597924‑M1/18‑50 (MQ=255)
                                tacaaattctcgcgccaCCCGTAGGCCGGATAAGGCGTTCAcgccgc                                                              >  1:5399880‑M1/18‑47 (MQ=255)
                                 atTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCAcgc                                                                 <  1:195599/43‑1 (MQ=255)
                                  ttGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                         >  1:801883/1‑50 (MQ=255)
                                  ttGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                         >  1:280150/1‑50 (MQ=255)
                                  cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCa                                                                    <  1:3717949‑M1/24‑1 (MQ=255)
                                  cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                             <  1:315120‑M1/31‑1 (MQ=255)
                                  cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                             <  1:4790849‑M1/31‑1 (MQ=255)
                                  cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                           <  1:4403124‑M1/33‑1 (MQ=255)
                                  cagattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCAACCggcaa                                                     <  1:2715978‑M1/39‑3 (MQ=255)
                                    gattcttgcgccaTTCGTAGGCCGGATAAGGCGTTCAcgccg                                                               <  1:870899‑M1/29‑1 (MQ=255)
                                    gAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTc                                                                     >  1:3506159/1‑36 (MQ=255)
                                    gAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCa                                                                    >  1:1994417/1‑37 (MQ=255)
                                      ttctcacgccaTCCGTAGGCCGGATAAGGCGTTCAcgcc                                                                <  1:2525628‑M1/28‑1 (MQ=255)
                                      ttctcacgccaTCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                        >  1:4501568‑M1/12‑47 (MQ=255)
                                         ttgtgccaTTCGTAGGCCGGATAAGGCGTTCAcgccgc                                                              <  1:1444150‑M1/30‑1 (MQ=255)
                                         ttgcgccaTTCGTAGGCCGGATAAGGCGTTCACGCCGCa                                                             <  1:2060699‑M1/31‑1 (MQ=255)
                                            aCGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCt                          >  1:2755046/1‑71 (MQ=255)
                                              gATTTCGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCt                          >  1:2185763/1‑69 (MQ=255)
                                                tttCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAgg                      <  1:2886075/71‑1 (MQ=255)
                                                  tCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                     <  1:5068809/70‑1 (MQ=255)
                                                  tCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGc                     >  1:1406879/1‑70 (MQ=255)
                                                   cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGcc                    <  1:2489035/70‑1 (MQ=255)
                                                       ggCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagt                 >  1:5048955/1‑69 (MQ=255)
                                                       ggCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCcagt                 >  1:1982972/1‑69 (MQ=255)
                                                          cGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCt                          <  1:4444136/57‑1 (MQ=255)
                                                          cGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTa            <  1:2578570/71‑1 (MQ=255)
                                                           ggATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTaa           <  1:526786/71‑1 (MQ=255)
                                                             aTAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATg         <  1:330856/71‑1 (MQ=255)
                                                              tAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCAc                                        <  1:4952622/39‑1 (MQ=255)
                                                              tAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGt        <  1:1501737/71‑1 (MQ=255)
                                                              tAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGt        >  1:3374618/1‑71 (MQ=255)
                                                                aGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGtt       <  1:874546/70‑1 (MQ=255)
                                                                 ggCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGAcagtcagt             <  1:950834/63‑1 (MQ=255)
                                                                  gCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTaa           >  1:1533119/1‑64 (MQ=255)
                                                                   cGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTg   <  1:1285640/71‑1 (MQ=255)
                                                                    gTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGc  <  1:1894343/71‑1 (MQ=255)
                                                                     ttCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTcaacaa                             <  1:3889928/43‑1 (MQ=255)
                                                                     ttCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTcaaca                              <  1:2062302/42‑1 (MQ=255)
                                                                      |                                                                    
ACGCTTATCAGGCCTACGTGGAGCCCACAATATATTGAATTAGCACGATTTCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCGTGTACAACGCGCACTTGGTCAACAATCTGAGGCCAGTCAGTAATGTTACTGC  >  minE/1602536‑1602674

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: