| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 825,519 | C→T | 50.8% | intergenic (‑241/‑85) | ychN ← / → ychP | conserved hypothetical protein/predicted invasin |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 825,519 | 0 | C | T | 50.8% | ‑4.8 / 87.4 | 59 | intergenic (‑241/‑85) | ychN/ychP | conserved hypothetical protein/predicted invasin |
| Reads supporting (aligned to +/- strand): ref base C (13/16); new base T (13/17); total (26/33) | |||||||||||
| 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.72e-01 | |||||||||||
AGTTACCGACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAACAATGATG > minE/825450‑825586 | aGTTACCGACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCg > 1:2089944/1‑71 (MQ=255) gTTACCGACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTgg < 1:2439016/71‑1 (MQ=255) gTTACCGACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCgg < 1:979064/71‑1 (MQ=255) gTTACCGACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCgg < 1:2590660/71‑1 (MQ=255) aCCGACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTa > 1:414337/1‑70 (MQ=255) gACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTa < 1:3842343/67‑1 (MQ=255) aCGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACa > 1:2349534/1‑68 (MQ=255) aCGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCt < 1:1344545/71‑1 (MQ=255) gTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTaa > 1:3459495/1‑71 (MQ=255) gTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTa > 1:2688275/1‑70 (MQ=255) ttttGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTa < 1:2359413/69‑1 (MQ=255) tGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACg < 1:1236157/69‑1 (MQ=255) aGCCGTTTCGATCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGcaca < 1:1375532/71‑1 (MQ=255) ttCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTAc > 1:2085836/1‑71 (MQ=255) tCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAAt > 1:1568782/1‑66 (MQ=255) tCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTAcc > 1:2983218/1‑71 (MQ=255) tCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTAc < 1:2354268/70‑1 (MQ=255) tCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCAAcc > 1:4026530/1‑71 (MQ=255) ttCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCgg < 1:2954715/46‑1 (MQ=255) tCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCtt > 1:4113124/1‑69 (MQ=255) tCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCt > 1:2621481/1‑68 (MQ=255) cTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAg < 1:3042949/48‑1 (MQ=255) cTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAg < 1:1223329/48‑1 (MQ=255) gCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCa > 1:3947183/1‑70 (MQ=255) cATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCaaa < 1:1733232/71‑1 (MQ=255) ttATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAAt < 1:2477762/54‑1 (MQ=255) ttATCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAAt < 1:1151337/54‑1 (MQ=255) tATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTaa > 1:3152442/1‑45 (MQ=255) tATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGcaca > 1:2602439/1‑51 (MQ=255) tATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGcaca > 1:967026/1‑51 (MQ=255) aTCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAg < 1:1757753/70‑1 (MQ=255) aTCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGc > 1:3115193/1‑71 (MQ=255) tCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCt < 1:3085292/41‑1 (MQ=255) tCCCGTTTTATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGc < 1:3589817/70‑1 (MQ=255) ttttATTTACTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCt < 1:220887/48‑1 (MQ=255) ttatttaCTCTTGCTTGTTGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATc > 1:3970380/1‑71 (MQ=255) ttatttaCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCaaa > 1:551937/1‑59 (MQ=255) ttatttaCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAAt < 1:2523927/70‑1 (MQ=255) atttaCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCg > 1:3267298/1‑62 (MQ=255) tttaCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCt < 1:169791/47‑1 (MQ=255) tttaCTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCt < 1:3602897/71‑1 (MQ=255) taCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCg > 1:2467429/1‑48 (MQ=255) taCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCaaa > 1:1359887/1‑64 (MQ=255) taCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCttt > 1:1768372/1‑71 (MQ=255) taCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCt > 1:1776430/1‑69 (MQ=255) aCTCTTGCTTGTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAAt > 1:952680/1‑64 (MQ=255) ttgcttgTAGCAGGCGGTACAGCTAACGCACAATCt < 1:806198/36‑1 (MQ=255) gcttgTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGataat < 1:3005445/70‑1 (MQ=255) gcttgTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGataa < 1:1325370/69‑1 (MQ=255) cttgTAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGataaaa < 1:2723383/70‑3 (MQ=255) cttgTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGataata < 1:3619437/70‑1 (MQ=255) tgTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAAc > 1:338566/1‑70 (MQ=255) gTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAAc < 1:3904290/69‑1 (MQ=255) gTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAAc > 1:3303368/1‑69 (MQ=255) tAGCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAACAat < 1:3935650/71‑1 (MQ=255) gCAGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAACAatga < 1:1632793/71‑1 (MQ=255) gCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAACa < 1:3031391/67‑1 (MQ=255) aGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAACAatgatg < 1:3541476/71‑1 (MQ=255) aGGTGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAACAatgatg < 1:3266005/71‑1 (MQ=255) | AGTTACCGACGTTTTGAGCCGTTTCGTTCCTCGCATTATCCCGTTTTATTTACTCTTGCTTGTAGCAGGCGGTACAGCTAACGCACAATCTACCTTCGAGCAAAAAGCGGCAAATCCCTTTGATAATAACAATGATG > minE/825450‑825586 |
| 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: ‑ |