Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 138,172 A→G 52.8% N665N (AAT→AAC htrE ← predicted outer membrane usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE138,1720AG52.8% 13.4 / 102.0 72N665N (AAT→AAChtrEpredicted outer membrane usher protein
Reads supporting (aligned to +/- strand):  ref base A (19/15);  new base G (22/16);  total (41/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

CAGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACATA  >  minE/138106‑138237
                                                                  |                                                                 
cagtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTg                                                               >  1:1292232/1‑71 (MQ=255)
 agtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGATAGAAACTTGACGGCTGTTATCGCTGTTTg                                                               <  1:3499757/70‑1 (MQ=255)
 agtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                               <  1:3218643/70‑1 (MQ=255)
 agtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTg                                                               <  1:4092560/70‑1 (MQ=255)
 agtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGc                                                              <  1:3928179/71‑1 (MQ=255)
  gtcCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCa                                                             <  1:4103344/71‑1 (MQ=255)
     cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGa                                                           <  1:4170510/70‑1 (MQ=255)
     cACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGa                                                           <  1:3821705/70‑1 (MQ=255)
      aCCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGaaa                                                         >  1:3498025/1‑71 (MQ=255)
        cGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGa                                                           <  1:326707/67‑1 (MQ=255)
        cGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGa                                                           <  1:3948405/67‑1 (MQ=255)
        cGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGa                                                           <  1:1681010/67‑1 (MQ=255)
         gCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTATGCAGAAAtt                                                       >  1:533748/1‑70 (MQ=255)
         gCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTg                                                      <  1:909614/71‑1 (MQ=255)
          cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGaaa                                                         >  1:3161689/1‑67 (MQ=255)
          cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGaaa                                                         >  1:2758349/1‑67 (MQ=255)
          cTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGaaa                                                         >  1:2428858/1‑67 (MQ=255)
             tGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAAcc                                                  >  1:156198/1‑71 (MQ=255)
               cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGaa                                                    >  1:3301215/1‑67 (MQ=255)
               cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTg                                                >  1:3589007/1‑71 (MQ=255)
               cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTg                                                >  1:4104627/1‑71 (MQ=255)
               cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTg                                                >  1:3238631/1‑71 (MQ=255)
               cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTg                                                >  1:1748220/1‑71 (MQ=255)
               cAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTg                                                >  1:1314292/1‑71 (MQ=255)
                 aTACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCa                                                             >  1:3048136/1‑56 (MQ=255)
                   aCAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAg                                            >  1:2345388/1‑71 (MQ=255)
                      aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                                     <  1:1341136/59‑1 (MQ=255)
                      aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                                     <  1:3498914/59‑1 (MQ=255)
                      aaaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGa                                                     <  1:1515305/59‑1 (MQ=255)
                       aaCCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGa                                                     <  1:2529611/58‑1 (MQ=255)
                          caccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTccc                                     >  1:402022/1‑71 (MQ=255)
                            ccGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCa                                   <  1:2064500/71‑1 (MQ=255)
                                cGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGt                                         >  1:4139003/1‑61 (MQ=255)
                                cGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGt                                         >  1:843812/1‑61 (MQ=255)
                                 ggTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGt                               <  1:2411293/70‑1 (MQ=255)
                                 ggTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGt                               <  1:2319205/70‑1 (MQ=255)
                                   tGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGc                                                              <  1:3923339/37‑1 (MQ=255)
                                   tGCTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAAtt                                                       <  1:1108859/44‑1 (MQ=255)
                                   tGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGAc                            >  1:460526/1‑71 (MQ=255)
                                     cTGAGAGAAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACt                           >  1:3000359/1‑70 (MQ=255)
                                     cTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                                  >  1:903193/1‑63 (MQ=255)
                                     cTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                                  >  1:4135773/1‑63 (MQ=255)
                                     cTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                                  >  1:1689442/1‑63 (MQ=255)
                                     cTGAGAGAAACTTGACGGATGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCAt                                  >  1:1093435/1‑63 (MQ=255)
                                      tGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGtt                                        <  1:1994358/56‑1 (MQ=255)
                                      tGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGtt                                        <  1:31212/56‑1 (MQ=255)
                                        agagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCa                         <  1:4005445/69‑1 (MQ=255)
                                        agagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa                       <  1:1009418/71‑1 (MQ=255)
                                        agagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa                       <  1:495778/71‑1 (MQ=255)
                                        agagaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATa                       <  1:2826117/71‑1 (MQ=255)
                                          agaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGc                     >  1:2358730/1‑71 (MQ=255)
                                          agaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGc                     >  1:3768169/1‑71 (MQ=255)
                                           gaAACTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGc                     <  1:2266926/70‑1 (MQ=255)
                                           gaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCa                                   >  1:4098461/1‑56 (MQ=255)
                                           gaAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCa                                   >  1:1658535/1‑56 (MQ=255)
                                               cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGc                 >  1:2083495/1‑70 (MQ=255)
                                               cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCa                >  1:770033/1‑71 (MQ=255)
                                               cTTGACGGCTGTTATCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCa                >  1:419195/1‑71 (MQ=255)
                                                  gACGTCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCAt               <  1:4072477/69‑1 (MQ=255)
                                                  gACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATaa             >  1:1193379/1‑71 (MQ=255)
                                                   aCGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCa                                   <  1:1685743/48‑1 (MQ=255)
                                                   aCGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCa                                   <  1:2995454/48‑1 (MQ=255)
                                                     ggCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATg                                 >  1:2943952/1‑48 (MQ=255)
                                                       cTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAAcccc         >  1:4105044/1‑70 (MQ=255)
                                                            aTCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCa                                   >  1:4127035/1‑39 (MQ=255)
                                                            aTCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATaa             >  1:101181/1‑61 (MQ=255)
                                                            aTCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACa    >  1:601597/1‑70 (MQ=255)
                                                            aTCGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACa    >  1:402950/1‑70 (MQ=255)
                                                            aTCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCa                         <  1:2356870/49‑1 (MQ=255)
                                                            aTCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACa    >  1:181454/1‑70 (MQ=255)
                                                            aTCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACa    >  1:4146025/1‑70 (MQ=255)
                                                              cGCTGTTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACATa  >  1:4008417/1‑70 (MQ=255)
                                                                  |                                                                 
CAGTCCACCGCTATGCAATACAAAACCACCGTCGGTGCTGAGAGAAACTTGACGGCTGTTATCGCTATTTGCAGAAATTGAACCTGCCAGCGTTCCCCATGGTGACTCATAGCTGGCATAACCCCCAACATA  >  minE/138106‑138237

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: