Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 1,495,966 | 1 | . | T | 62.0% | 82.6 / 71.6 | 50 | coding (820/924 nt) | ydcI | predicted DNA‑binding transcriptional regulator |
Reads supporting (aligned to +/- strand): ref base . (3/16); new base T (17/14); total (20/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.97e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAG > W3110S.gb/1495896‑1496025 | atcaACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATg > 1:2199324/1‑71 (MQ=255) caACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑tt > 1:2672484/1‑70 (MQ=255) caACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑tt > 1:4100722/1‑70 (MQ=255) caACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑tt > 1:2266780/1‑70 (MQ=255) aaCCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATGGGTTCGCCTGCACCATGTTCCCGGaa > 1:3985329/1‑70 (MQ=255) aaCCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGaa > 1:2266169/1‑70 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg > 1:3323699/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg > 1:202032/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg > 1:2832614/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg > 1:1515927/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg > 1:3623910/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg > 1:2286819/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTg > 1:2659274/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTgg > 1:2814217/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTgg > 1:2990393/1‑71 (MQ=255) aTCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGACACCAGCGTgg > 1:3342588/1‑71 (MQ=255) tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc > 1:604469/1‑70 (MQ=255) tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa > 1:3215022/1‑71 (MQ=255) tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa > 1:3073073/1‑71 (MQ=255) tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa > 1:1363309/1‑71 (MQ=255) tACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCa > 1:636618/1‑71 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTgga < 1:2442608/68‑2 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc < 1:1536046/68‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc < 1:649081/68‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAAGAGGCAGGGCCACCAGCGTGGCAt < 1:851075/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:2682729/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:820210/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:1246941/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:1286454/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:3988610/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:1388366/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:3656852/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:1430965/69‑1 (MQ=255) cgcgGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCAt < 1:1431657/69‑1 (MQ=255) ggTCAGTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACg > 1:3914715/1‑70 (MQ=255) gTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtc < 1:2823350/71‑1 (MQ=255) gTATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtc < 1:1697375/71‑1 (MQ=255) tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcg < 1:1219474/71‑1 (MQ=255) tATTCCAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGgtcg < 1:2243241/71‑1 (MQ=255) tATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc < 1:1581956/64‑1 (MQ=255) tATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc < 1:1008290/64‑1 (MQ=255) tATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc < 1:537729/64‑1 (MQ=255) tATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc < 1:3639503/64‑1 (MQ=255) tATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc < 1:3005788/64‑1 (MQ=255) cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt < 1:3903870/71‑1 (MQ=255) cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt < 1:286187/71‑1 (MQ=255) cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt < 1:1841289/71‑1 (MQ=255) cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt < 1:3000188/71‑1 (MQ=255) cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt < 1:3381154/71‑1 (MQ=255) cAATCGGTTCGCCTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCttt < 1:1239624/71‑1 (MQ=255) gCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGc < 1:1610177/49‑1 (MQ=255) cTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCt > 1:1483486/1‑57 (MQ=255) cTGCACCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAg > 1:2453602/1‑70 (MQ=255) aCCATGTTCCCGGAACAGGCAGGGCCACCAGCGTGGc < 1:679084/37‑1 (MQ=37) | ATCAACTGGCAACCAGAAGAGAACGTCGCATCTACGCGGGTCAGTATTCCAATCGGTTCGCCTGCACCATG‑TCCCGGAACAGGCAGGGCCACCAGCGTGGCATGACGCAGGTCGTCTTTTACAGCGCCAG > W3110S.gb/1495896‑1496025 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |