| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | W3110S.gb | 52,474 | 1 | . | A | 61.3% | 52.7 / 43.8 | 31 | coding (943/990 nt) | pdxA | 4‑hydroxy‑L‑threonine phosphate dehydrogenase, NAD‑dependent |
| Reads supporting (aligned to +/- strand): ref base . (7/5); new base A (19/0); total (26/5) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGG > W3110S.gb/52416‑52529 | gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat‑aa > 1:436375/1‑61 (MQ=255)gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat‑aa > 1:3026508/1‑61 (MQ=255)gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat‑aa > 1:2809299/1‑61 (MQ=255)gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat‑aa > 1:1503887/1‑61 (MQ=255)gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat‑aa > 1:2548737/1‑61 (MQ=255)gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat‑aa > 1:1886863/1‑61 (MQ=255)gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat‑aa > 1:71765/1‑61 (MQ=255)gACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGtaat > 1:1669611/1‑59 (MQ=255) cTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGc < 1:1998924/65‑1 (MQ=255) cAATCATTTTGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTcc < 1:1232097/64‑1 (MQ=255) cAATCATTTTGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTcc < 1:1976306/64‑1 (MQ=255) cAATCATTTTGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTcc < 1:2835243/64‑1 (MQ=255) cAATCATTTTGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTcc < 1:899612/64‑1 (MQ=255) ttGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGt > 1:1941517/1‑63 (MQ=255) tGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCa > 1:1982101/1‑65 (MQ=255) tGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCa > 1:509820/1‑65 (MQ=255) aGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGt > 1:2387198/1‑64 (MQ=255) aGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTg > 1:1990103/1‑65 (MQ=255) aGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTg > 1:1311371/1‑65 (MQ=255) aGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTg > 1:1954807/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:731458/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:898300/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:934175/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:956774/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:687041/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:1054184/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:999487/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:498246/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:3096109/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:3038231/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:282282/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:2702663/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:254872/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:2354247/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:2279572/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:2132599/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:1888714/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:1541767/1‑65 (MQ=255) cgcCGTAATAAAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTgg > 1:1462832/1‑65 (MQ=255) | GACTCGATTATTCATTGGGTGTTAACAATCATTTTGATGGCGAGATTAAGCGCCGTAAT‑AAAACTGCCGACATCGGCTTTGCCACGTCCCGCCAGTTCAAGCGCGGTGCCGTGG > W3110S.gb/52416‑52529 |
| 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: ‑ |