Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 2,965,801 | 0 | C | A | 32.7% | 97.0 / 38.2 | 52 | P428Q (CCG→CAG) | deoA | thymidine phosphorylase |
Reads supporting (aligned to +/- strand): ref base C (35/0); new base A (17/0); total (52/0) | |||||||||||
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.32e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAA > minE/2965757‑2965841 | ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACt > 1:100612/1‑61 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACt > 1:453989/1‑61 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:220577/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:633596/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:660975/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:474045/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:433155/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:679299/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:690979/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:766774/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:769878/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:197067/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:105307/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:119628/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:1171723/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:932720/1‑62 (MQ=255) ggcgAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACAGGAAAGCACACCAACTg > 1:1188750/1‑62 (MQ=255) aTTAACCTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:173085/1‑62 (MQ=255) aTTAAACTTGCGGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:184340/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGc > 1:481250/1‑57 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:587165/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:974280/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:638374/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:954581/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:67088/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:680792/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:93204/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:730385/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:877450/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:987919/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:799552/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:859015/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:245759/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:1007887/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:1016290/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:1055221/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:1193190/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:1193969/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:136474/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:195693/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:23055/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:565061/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:320220/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:333295/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:340736/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:343164/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:373353/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:490562/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:498802/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:53880/1‑62 (MQ=255) aTTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATa > 1:543447/1‑62 (MQ=255) aTTAAACTT‑CCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATaa > 1:200290/1‑62 (MQ=255) | GGCGAAAGCGGTGAAAGCGGCAATTAAACTTGCCGATAAAGCACCGGAAAGCACACCAACTGTCTATCGCCGTATCAGCGAATAA > minE/2965757‑2965841 |
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: ‑ |