Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 141,333 A→T 35.8% A97A (GCT→GCA yadN ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE141,3330AT35.8% 56.5 / 52.0 67A97A (GCT→GCAyadNpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base A (31/12);  new base T (9/15);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 9.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTA  >  minE/141267‑141399
                                                                  |                                                                  
aaTAGACATGTCGTTGTTCGGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:2416876/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:2979815/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:909803/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:754957/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:3066429/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:3116142/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:2021928/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:2914283/69‑1 (MQ=255)
aaTAGACATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:37950/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:1460798/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:2866839/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:747510/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:2139445/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:1381223/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:2931302/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:1545499/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:2886813/1‑66 (MQ=255)
      cATGTCGTTGTTCAGTGTTCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGCAGtagt                                                            >  1:143881/1‑66 (MQ=255)
         gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:1805546/1‑69 (MQ=255)
         gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:1353637/1‑69 (MQ=255)
         gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:355765/1‑69 (MQ=255)
         gTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:1238580/1‑69 (MQ=255)
          tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:1686035/1‑68 (MQ=255)
          tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:677764/1‑68 (MQ=255)
          tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:3248453/1‑68 (MQ=255)
          tCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTAcc                                                         >  1:886182/1‑68 (MQ=255)
                  cAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCgg                                                <  1:1787470/69‑1 (MQ=255)
                    aaGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGa                                               <  1:1055308/68‑1 (MQ=255)
                       ttccTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:2645817/44‑1 (MQ=255)
                             tGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTTGc                                                                  <  1:2855215/40‑1 (MQ=255)
                             tGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGAtttgctttgc                                     >  1:3128209/1‑69 (MQ=255)
                             tGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGAtttgctttgc                                     >  1:3265060/1‑69 (MQ=255)
                             tGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGAtttgctttgc                                     >  1:3300595/1‑69 (MQ=255)
                             tGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGAtttgctttgc                                     >  1:973231/1‑69 (MQ=255)
                             tGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGAtttgctttgc                                     >  1:1864451/1‑69 (MQ=255)
                                    aCCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTc                              >  1:1453195/1‑69 (MQ=255)
                                    aCCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTc                              >  1:2620124/1‑69 (MQ=255)
                                    aCCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTc                              >  1:1110365/1‑69 (MQ=255)
                                       aaaGAAAACAGAACCGAAGGTCATTTTTGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAAc                           <  1:1475025/69‑1 (MQ=255)
                                       aaaGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAAc                           <  1:1012957/69‑1 (MQ=255)
                                       aaaGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAAc                           <  1:2676475/69‑1 (MQ=255)
                                       aaaGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAAc                           <  1:2501006/69‑1 (MQ=255)
                                       aaaGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAAc                           <  1:2830622/69‑1 (MQ=255)
                                       aaaGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAAc                           <  1:1504164/69‑1 (MQ=255)
                                        aaGAAAACAGAACCAAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACg                          <  1:2355820/69‑1 (MQ=255)
                                          gAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGt                        <  1:676489/69‑1 (MQ=255)
                                          gAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGt                        <  1:1823208/69‑1 (MQ=255)
                                          gAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGt                        <  1:2681918/69‑1 (MQ=255)
                                                      gAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACg                          <  1:1188512/55‑1 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:726692/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:1347489/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:965616/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:1431730/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:150042/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:654914/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:1436203/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:1826776/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:2284005/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:2376603/1‑69 (MQ=255)
                                                            cATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagc      >  1:2465029/1‑69 (MQ=255)
                                                                ttAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTc                              >  1:511547/1‑41 (MQ=255)
                                                                ttAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagcttta  >  1:1110248/1‑69 (MQ=255)
                                                                ttAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagcttta  >  1:2511866/1‑69 (MQ=255)
                                                                ttAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGtttagcttt   >  1:2846256/1‑68 (MQ=255)
                                                                  |                                                                  
AATAGACATGTCGTTGTTCAAAGTCCCTTTGCTGTTACCAAAGAAAACAGAACCGAAGGTCATTTTAGCCGTGCTACCTGGATTCGGATTTGCTTTGCTGCAGTCAACGGTGATGCTGAAAGGTTTAGCTTTA  >  minE/141267‑141399

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: