Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,700,888 G→A 29.5% A166T (GCG→ACG)  hofQ → predicted fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,700,8880GA29.5% 95.3 / 41.8 95A166T (GCG→ACG) hofQpredicted fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base G (31/35);  new base A (6/22);  total (37/58)
Fisher's exact test for biased strand distribution p-value = 2.27e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTATGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCATAT  >  minE/2700824‑2700952
                                                                  |                                                                
gAGTATGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAgcgcgc                                                               >  1:1006401/1‑68 (MQ=255)
gAGTATGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAgcgcgc                                                               >  1:699625/1‑68 (MQ=255)
    aTGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTg                                                            >  1:2698819/1‑67 (MQ=255)
    aTGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTg                                                            >  1:1674724/1‑67 (MQ=255)
    aTGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTg                                                            >  1:2057856/1‑67 (MQ=255)
    aTGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTg                                                            >  1:3392705/1‑67 (MQ=255)
    aTGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTg                                                            >  1:2438660/1‑67 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:245621/69‑1 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:2000599/69‑1 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:2629128/69‑1 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:156352/69‑1 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:354339/69‑1 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:2226092/69‑1 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:240072/69‑1 (MQ=255)
     tGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:2800914/69‑1 (MQ=255)
       aCCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAAc                                                         <  1:456585/67‑1 (MQ=255)
       aCCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTAGAACAg                                                       >  1:91811/1‑69 (MQ=255)
       aCCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTAGAACAg                                                       >  1:2564120/1‑69 (MQ=255)
       aCCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTAGAACAg                                                       >  1:1662412/1‑69 (MQ=255)
                 aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:2644669/1‑68 (MQ=255)
                 aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:1769290/1‑68 (MQ=255)
                 aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:1264964/1‑68 (MQ=255)
                 aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:2829894/1‑68 (MQ=255)
                   cGCACCAATCGCCTTTTGCTACGAGATAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:786411/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:1393252/64‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCTCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:465953/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:958835/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:222361/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:908678/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2471275/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:1333908/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:3003443/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:3125659/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2888634/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:1871284/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:63953/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2287720/69‑1 (MQ=255)
                        cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2307828/69‑1 (MQ=255)
                            cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTTGGTAGCGCa                                            <  1:2268296/59‑1 (MQ=255)
                            cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:3431320/59‑1 (MQ=255)
                            cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:865499/59‑1 (MQ=255)
                            cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:1824147/59‑1 (MQ=255)
                            cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:1669134/59‑1 (MQ=255)
                             gCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGc                                 >  1:2150913/1‑69 (MQ=255)
                             gCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGc                                 >  1:2147037/1‑69 (MQ=255)
                             gCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGc                                 >  1:2128239/1‑69 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGTGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2447236/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGTGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1661521/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:810776/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:3262287/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1081550/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:3005422/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1372492/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2978163/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1378146/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2898164/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1581047/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2835528/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2238419/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2769379/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:269075/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1807701/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1939390/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2627743/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:257306/65‑1 (MQ=255)
                                     ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:19860/65‑1 (MQ=255)
                                          agaTAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:1743896/67‑1 (MQ=255)
                                          agaTAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:2888497/67‑1 (MQ=255)
                                          agaTAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:2892537/67‑1 (MQ=255)
                                          agaTAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:2617337/67‑1 (MQ=255)
                                          agaTAACAAACACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:1275075/67‑1 (MQ=255)
                                          agaTAACAAAACGGCGTTAAGGACGCTTGAACAGTGGGTATCGCa                                            >  1:3046740/1‑45 (MQ=255)
                                          agaTAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAg                      <  1:196613/67‑1 (MQ=255)
                                          agaTAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCa                                            >  1:1037354/1‑45 (MQ=255)
                                          agaTAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCa                                            >  1:1249524/1‑45 (MQ=255)
                                          agaTAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCa                                            >  1:1161263/1‑45 (MQ=255)
                                          agaTAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCa                                            >  1:1800256/1‑45 (MQ=255)
                                          agaTAACAAAACGGCGTTAAGCACGCTTGAACAGAGGGTat                                                >  1:1992495/1‑40 (MQ=38)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTg                  >  1:3370788/1‑67 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:552349/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:1725884/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:731688/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:1826250/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:2030420/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:2656020/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:2948365/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:3392754/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:307552/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:3331683/1‑69 (MQ=255)
                                              aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGCTCGGGCAGGTTGAg                >  1:1496537/1‑69 (MQ=255)
                                                       cGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:518800/52‑1 (MQ=255)
                                                         tcGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:1297810/37‑1 (MQ=255)
                                                             tAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1683602/46‑1 (MQ=255)
                                                              aaGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCatat  >  1:993389/1‑69 (MQ=255)
                                                              aaGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCatat  >  1:3244023/1‑69 (MQ=255)
                                                              aaGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCatat  >  1:76879/1‑69 (MQ=255)
                                                              aaGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCatat  >  1:580015/1‑69 (MQ=255)
                                                                  |                                                                
GAGTATGACCGTCGATAAACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCATAT  >  minE/2700824‑2700952

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: