Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 141,288 A→T 53.1% T112T (ACT→ACA yadN ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE141,2880AT53.1% 7.3 / 70.5 64T112T (ACT→ACAyadNpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base T (9/25);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 7.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GTGTAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGA  >  minE/141231‑141354
                                                         |                                                                  
gtgtAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGc                                                          <  1:948319/68‑1 (MQ=255)
gtgtAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGc                                                          <  1:556465/68‑1 (MQ=255)
gtgtAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGc                                                          <  1:1369399/68‑1 (MQ=255)
gtgtAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGc                                                          <  1:3382588/68‑1 (MQ=255)
gtgtAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCATTTGc                                                          <  1:2502055/68‑1 (MQ=255)
   tAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGt                                                       <  1:1372140/68‑1 (MQ=255)
   tAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGt                                                       <  1:88066/68‑1 (MQ=255)
   tAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGt                                                       <  1:465365/68‑1 (MQ=255)
        cGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTcc                                                                >  1:1428696/1‑54 (MQ=255)
        cGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTAc                                                    >  1:3108185/1‑66 (MQ=255)
        cGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTAc                                                    >  1:2496243/1‑66 (MQ=255)
        cGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAACTCCCTTTGCTGTTAc                                                    >  1:3241435/1‑66 (MQ=255)
            atTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGaa                                             >  1:2050118/1‑69 (MQ=255)
               aaCGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGaa                                             >  1:2135642/1‑66 (MQ=255)
               aaCGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGaa                                             >  1:1008238/1‑66 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:2157117/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:2473194/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:2542390/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:2713665/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:380009/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:2933179/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:3101334/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:3129938/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:215615/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:1089267/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:753468/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:1432358/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:1348289/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:664103/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:732243/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:92261/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:151001/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:753673/69‑1 (MQ=255)
                                 gttAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCAtttt                        <  1:13063/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCGGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:2416876/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:37950/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:3116142/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:909803/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:3066429/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:2979815/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:754957/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:2914283/69‑1 (MQ=255)
                                    aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:2021928/69‑1 (MQ=255)
                                     aTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:3060745/68‑1 (MQ=255)
                                     aTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:379633/68‑1 (MQ=255)
                                      tAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                     <  1:3409462/67‑1 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:747510/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:2886813/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:2931302/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:2866839/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:2139445/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:1545499/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:1460798/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:143881/1‑66 (MQ=255)
                                          cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt               >  1:1381223/1‑66 (MQ=255)
                                             gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:355765/1‑69 (MQ=255)
                                             gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:1805546/1‑69 (MQ=255)
                                             gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:1238580/1‑69 (MQ=255)
                                             gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:1353637/1‑69 (MQ=255)
                                              tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:3248453/1‑68 (MQ=255)
                                              tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:886182/1‑68 (MQ=255)
                                              tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:1686035/1‑68 (MQ=255)
                                              tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc            >  1:677764/1‑68 (MQ=255)
                                                      cAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCgg   <  1:1787470/69‑1 (MQ=255)
                                                        aaGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGa  <  1:1055308/68‑1 (MQ=255)
                                                         |                                                                  
GTGTAGAGCGATATTAACGCCGTCAGATGGGTTGTTAATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGA  >  minE/141231‑141354

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: