Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,418,849 A→G 38.6% F388F (TTT→TTC wzyE ← predicted Wzy protein involved in ECA polysaccharide chain elongation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,418,8490AG38.6% 55.9 / 77.8 114F388F (TTT→TTCwzyEpredicted Wzy protein involved in ECA polysaccharide chain elongation
Reads supporting (aligned to +/- strand):  ref base A (22/48);  new base G (18/26);  total (40/74)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 127 total aligned reads displayed.

AGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCA  >  minE/2418782‑2418909
                                                                   |                                                            
aggccACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:2455915/4‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:1004192/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:10279/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:858640/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:1094149/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:2201054/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:2292263/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:2952833/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:2912690/1‑69 (MQ=255)
aGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCaa                                                             >  1:14508/1‑69 (MQ=255)
       gcgTGAGACAAACGAATCCATCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2611573/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2083405/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:211468/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2367072/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2498887/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2700946/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2758712/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:298572/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:300952/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:331740/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:3339499/68‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGaa                                                      <  1:252243/69‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGaa                                                      <  1:2891688/69‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGaa                                                      <  1:2484761/69‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGaa                                                      <  1:2911481/69‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGaa                                                      <  1:215203/69‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGaa                                                      <  1:970376/69‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCACAAAGCAGaa                                                      <  1:816521/69‑1 (MQ=255)
       gcgTGAGACAAACGAATCCAGCCCTTCACGCGC‑‑‑CACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2587393/65‑1 (MQ=255)
         gTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:1154592/66‑1 (MQ=255)
                           gCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCtt                                  <  1:2147073/69‑1 (MQ=255)
                             ccTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCtt                                  >  1:1910767/1‑67 (MQ=255)
                             ccTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCtt                                  >  1:705346/1‑67 (MQ=255)
                             ccTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCtt                                  >  1:213412/1‑67 (MQ=255)
                              cTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:1598488/45‑1 (MQ=255)
                                 cACGCGCCAGCACGATCATATTGAAGATCGCCCCGAAGCAGa                                                       <  1:2380343/42‑1 (MQ=255)
                                 cACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGa                                                       <  1:1649780/42‑1 (MQ=255)
                                 cACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:578266/1‑69 (MQ=255)
                                    gcgcCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGAt                         <  1:1823713/69‑1 (MQ=255)
                                    gcgcCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGAt                         <  1:1713981/69‑1 (MQ=255)
                                    gcgcCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGAt                         <  1:2482366/69‑1 (MQ=255)
                                    gcgcCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGAt                         <  1:1030816/69‑1 (MQ=255)
                                    gcgcCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGAt                         <  1:2369906/69‑1 (MQ=255)
                                    gcgcCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGAt                         <  1:1560283/69‑1 (MQ=255)
                                        cAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAgg                     <  1:2717515/69‑1 (MQ=255)
                                        cAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAgg                     <  1:3324369/69‑1 (MQ=255)
                                         aGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAgg                     <  1:953096/68‑1 (MQ=255)
                                         aGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGc                    >  1:1217242/1‑69 (MQ=255)
                                         aGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGc                    >  1:2861628/1‑69 (MQ=255)
                                                 cATATTGAAGATCGCCGCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:31780/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:3139590/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:70118/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:1869662/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:3062788/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:1122923/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:2927222/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:3246833/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:1401118/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:1898364/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:2178073/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:2103309/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:1786531/49‑1 (MQ=255)
                                                 cATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:2458322/49‑1 (MQ=255)
                                                    atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc          <  1:1370078/68‑1 (MQ=255)
                                                    atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc          <  1:1317808/68‑1 (MQ=255)
                                                    atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc          <  1:694908/68‑1 (MQ=255)
                                                    atTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGc          <  1:1186567/68‑1 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAg                             >  1:1499886/1‑48 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:1058103/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:1678668/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:556194/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:2198559/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:2328238/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:728818/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:832739/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:1914798/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:913718/1‑49 (MQ=255)
                                                     ttGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGc                            >  1:1431290/1‑49 (MQ=255)
                                                      tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTcg                 >  1:1462487/1‑59 (MQ=255)
                                                      tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc        >  1:202159/1‑68 (MQ=255)
                                                      tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc        >  1:1223014/1‑68 (MQ=255)
                                                      tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc        >  1:2671785/1‑68 (MQ=255)
                                                      tGAAGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGccc        >  1:541787/1‑68 (MQ=255)
                                                         aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc     <  1:840325/68‑1 (MQ=255)
                                                         aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc     <  1:101314/68‑1 (MQ=255)
                                                         aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc     <  1:3149250/68‑1 (MQ=255)
                                                         aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc     <  1:2378216/68‑1 (MQ=255)
                                                         aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc     <  1:2498539/68‑1 (MQ=255)
                                                         aGATCGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGc     <  1:2672245/68‑1 (MQ=255)
                                                           aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTta                                <  1:3023527/39‑1 (MQ=255)
                                                           aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa  <  1:3255101/69‑1 (MQ=255)
                                                           aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa  <  1:244504/69‑1 (MQ=255)
                                                           aTCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCa  <  1:2587716/69‑1 (MQ=255)
                                                             cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa  >  1:352404/1‑67 (MQ=255)
                                                             cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa  >  1:1103389/1‑67 (MQ=255)
                                                             cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa  >  1:3243062/1‑67 (MQ=255)
                                                             cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa  >  1:320916/1‑67 (MQ=255)
                                                             cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa  >  1:1921868/1‑67 (MQ=255)
                                                             cGCCCCGAAGCAGAAACTGTGCAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa  >  1:1258592/1‑67 (MQ=255)
                                                             cGCCCCGAAGCAGAAACTGTGAAATATCGCAGCCTTATAGCGATTAGTCTCGCGGTTGCCCAGCTCa  >  1:1473443/1‑67 (MQ=255)
                                                                   |                                                            
AGACCACGCGTGAGACAAACGAATCCAGCCCTTCACGCGCCAGCACGATCATATTGAAGATCGCCCCAAAGCAGAAACTGTGCAATATCGCAGCTTTATAGCGATTAGGCTCGCGGTTGCCCAGCTCA  >  minE/2418782‑2418909

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: