| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 287,335 | C→T | S933S (AGC→AGT) | eaeH → | attaching and effacing protein, pathogenesis factor |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 287,335 | 0 | C | T | 98.0% | 189.2 / ‑4.4 | 51 | S933S (AGC→AGT) | eaeH | attaching and effacing protein, pathogenesis factor |
| Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (23/27); total (23/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01 | |||||||||||
AACTCACCTTCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGCGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAATGGGGTCTACG > REL606/287262‑287405 | aaCTCACCTTCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTgg < 2:675427/76‑1 (MQ=255)aaCTCACCTTCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTgg < 2:513229/76‑1 (MQ=255)aaCTCACCTTCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTgg > 1:580421/1‑76 (MQ=255) aCCTTCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTgccg > 1:1254617/1‑76 (MQ=255) ttCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCa < 2:1326938/76‑1 (MQ=255) tCACCGTGAAGGATGCGTACGNTAACGCGGTCACCGGGGTGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCaa < 2:246632/76‑1 (MQ=255) tCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCaa < 2:1042598/76‑1 (MQ=255) aCCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAaca < 2:1145133/76‑1 (MQ=255) ccGTGAAGGATGCGTACGGGACCCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAacac < 2:736987/76‑1 (MQ=255) aaGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGa < 2:1148608/76‑1 (MQ=255) aGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAg > 2:1087762/1‑76 (MQ=255) aGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAg < 2:1149215/76‑1 (MQ=255) gATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTg < 2:1294375/76‑1 (MQ=255) aTGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGa < 1:804596/76‑1 (MQ=255) tGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAg < 1:172638/76‑1 (MQ=255) cGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGcn < 1:360533/76‑2 (MQ=255) cGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCg < 1:252951/76‑1 (MQ=255) cGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCg > 1:1338451/1‑76 (MQ=255) gTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGt > 1:980002/1‑76 (MQ=255) gTACGGGAACCCGGTCACCGGGATGAAGCCAGACGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGt > 1:187164/1‑76 (MQ=255) tACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTc > 1:550758/1‑76 (MQ=255) aCGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTcc > 2:1141197/1‑76 (MQ=255) aCGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTcc > 2:27079/1‑76 (MQ=255) cGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGAGCCGCCAACACGGGGAGTGAGCGTCCt > 2:927804/1‑76 (MQ=255) gggAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCtt > 2:118239/1‑76 (MQ=255) gggAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCtt < 1:31235/76‑1 (MQ=255) gggAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCtt < 2:1019596/76‑1 (MQ=255) ggAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTc < 1:768745/76‑1 (MQ=255) ggAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTc < 2:1313692/76‑1 (MQ=255) aaCCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTcag > 2:945635/1‑76 (MQ=255) ggTCCCCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGa < 1:259768/76‑1 (MQ=255) gTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGaa < 2:1202684/76‑1 (MQ=255) aCCGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACt > 1:1279503/1‑76 (MQ=255) ccGGGATGAAGCCAGATGCCCCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGGGCGTCCTTCAGCAGGAAACTg < 1:478756/76‑1 (MQ=255) cGGGATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTgg < 1:1028868/76‑1 (MQ=255) gggATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGa > 2:1491199/1‑76 (MQ=255) gggATGAAGCCAGATGCACCAGTGTTTAGCGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGa < 2:1333483/76‑1 (MQ=255) ggATGAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCACGCAACTGGAc > 1:754582/1‑76 (MQ=255) gAAGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACagag > 1:393501/1‑76 (MQ=255) gAAGCCAGATGCACCAGTGTTTAGTGGTACCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACagag > 1:1171532/1‑76 (MQ=255) aaGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAgaga < 2:1131436/76‑1 (MQ=255) aaGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAgaga > 2:938730/1‑76 (MQ=255) aaGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAgaga < 1:995917/76‑1 (MQ=255) aGCCAGATGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGaa < 1:1169005/76‑1 (MQ=255) tGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAAt > 2:1177600/1‑76 (MQ=255) tGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAAt > 1:1207352/1‑76 (MQ=255) tGCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGAAAt > 1:1177242/1‑76 (MQ=255) gCACCAGTGTTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAATg < 2:1338451/76‑1 (MQ=255) gtgtTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAATGGGGTCt > 1:968652/1‑76 (MQ=255) gtgtTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAATGGGGTCt < 1:1279016/76‑1 (MQ=255) gtTTAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAATGGGGTCTAc > 1:336466/1‑76 (MQ=255) tttAGTGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAATGGGGTCTACg < 2:1281531/76‑1 (MQ=255) | AACTCACCTTCACCGTGAAGGATGCGTACGGGAACCCGGTCACCGGGATGAAGCCAGATGCACCAGTGTTTAGCGGTGCCGCCAACACGGGGAGTGAGCGTCCTTCAGCAGGAAACTGGACAGAGAAAGGTAATGGGGTCTACG > REL606/287262‑287405 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |