Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,700,906 G→T 36.9% A172S (GCG→TCG)  hofQ → predicted fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,700,9060GT36.9% 69.6 / 51.5 73A172S (GCG→TCG) hofQpredicted fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base G (19/27);  new base T (5/22);  total (24/49)
Fisher's exact test for biased strand distribution p-value = 7.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCATAT  >  minE/2700841‑2700952
                                                                   |                                              
aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:2644669/1‑68 (MQ=255)
aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:2829894/1‑68 (MQ=255)
aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:1264964/1‑68 (MQ=255)
aaCGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAgcg                                              >  1:1769290/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:1871284/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2888634/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2287720/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2307828/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:3003443/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:3125659/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:2471275/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:63953/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:1333908/69‑1 (MQ=255)
       cAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGa                                      <  1:908678/69‑1 (MQ=255)
           cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTTGGTAGCGCa                                            <  1:2268296/59‑1 (MQ=255)
           cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:1669134/59‑1 (MQ=255)
           cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:3431320/59‑1 (MQ=255)
           cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:1824147/59‑1 (MQ=255)
           cGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCa                                            <  1:865499/59‑1 (MQ=255)
            gCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGc                                 >  1:2128239/1‑69 (MQ=255)
            gCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGc                                 >  1:2147037/1‑69 (MQ=255)
            gCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGc                                 >  1:2150913/1‑69 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGTGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2447236/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGTGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1661521/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1081550/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2627743/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:269075/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2769379/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2835528/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2898164/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2978163/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:3005422/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:3262287/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:810776/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:257306/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1581047/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1378146/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1807701/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1372492/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:1939390/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:19860/65‑1 (MQ=255)
                    ctgcGAGATAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCAAATGGATCTGCCGGTCgg                          <  1:2238419/65‑1 (MQ=255)
                         agaTAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:2888497/67‑1 (MQ=255)
                         agaTAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:2892537/67‑1 (MQ=255)
                         agaTAACAAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCa                       <  1:1743896/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:1161263/1‑45 (MQ=255)
                         agaTAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCa                                            >  1:1037354/1‑45 (MQ=255)
                         agaTAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCa                                            >  1:1249524/1‑45 (MQ=255)
                         agaTAACAAAACGGCGTTAAGCACGCTTGAACAGTGGGTATCGCa                                            >  1:1800256/1‑45 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTg                  >  1:3370788/1‑67 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:307552/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:1725884/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:3331683/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:3392754/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:2656020/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:552349/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:2948365/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:1826250/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:731688/1‑69 (MQ=255)
                             aaCAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAg                >  1:2030420/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:76879/1‑69 (MQ=255)
                                             aaGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCatat  >  1:580015/1‑69 (MQ=255)
                                             aaGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCatat  >  1:3244023/1‑69 (MQ=255)
                                                                   |                                              
AACGCACCAATCGCCTTTTGCTACGAGATAACAAAACGGCGTTAAGCGCGCTTGAACAGTGGGTAGCGCAAATGGATCTGCCGGTCGGGCAGGTTGAGCTGTCGGCGCATAT  >  minE/2700841‑2700952

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