Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 859,054 | 0 | G | A | 32.3% | 35.5 / 18.0 | 31 | S3F (TCT→TTT) | ybiU | DUF1479 domain‑containing protein YbiU |
Reads supporting (aligned to +/- strand): ref base G (21/0); new base A (0/10); total (21/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.25e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATA > NC_000913/858992‑859115 | cgTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGaa > 1:3531861/1‑65 (MQ=255) cTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGGAAAAGGAGAAACCa > 1:2903856/1‑65 (MQ=255) tCTGACGGATAGCTGTTTTGTGATCGGCAGGAAATGTGTCGCTGGAAAAAGAAAAAGCCATAacc < 1:3042992/65‑1 (MQ=255) tCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATTTGTCGCTGGTAAAAGTAGAAGGCATaagc > 1:950342/1‑63 (MQ=255) tCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATaaac > 1:2791833/1‑63 (MQ=255) tCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacc > 1:346372/1‑65 (MQ=255) tCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacc > 1:3069507/1‑65 (MQ=255) tCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAacc > 1:907703/1‑65 (MQ=255) aCGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACct > 1:3759747/1‑65 (MQ=255) aCGGATAGCTGCTTTGTGATCGGCAGGCAATGGGTCGCTGGTAAAAGAAGAAGCCATAAGCACct > 1:2517846/1‑65 (MQ=255) cGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACctc > 1:1728194/1‑65 (MQ=255) ggttAGCTGTTTGGTGATGGGCAGGCAATGTGTCGCTGGAAAAAGAAAAAGCCATAACCACctct < 1:2069501/62‑1 (MQ=255) ggATGGCTGTTTTGTAAGGGGCAGGAAATGTGTCGCTGTAAAAAGAAAAAGCCATAACCACctct < 1:2536402/65‑1 (MQ=255) ggATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGGCATAACCACctct > 1:607944/1‑65 (MQ=255) aatgctcctTTGTGATGGGAGGGTAATGTGGCGTTGTAAAAATTAAAAGCCATAACCACCTCTCa < 1:2867006/62‑1 (MQ=255) aTACCTGTTTTGTGATCGGCAGGCAATGTGTCGCTGGAAAAAGAAAAAGCCATAACCACCTCTCa < 1:2823705/65‑1 (MQ=255) aGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCaat > 1:1245553/1‑63 (MQ=255) aGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCAtt > 1:654796/1‑65 (MQ=255) tttttGTGATCGGCAGGCAATGTGTCGCTGGAAAAAGTAAAAGCCATAACCACCTCTCATTAACt < 1:3902245/63‑1 (MQ=255) gctTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACt > 1:735555/1‑65 (MQ=255) gctTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACt > 1:345337/1‑65 (MQ=255) tttGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACttt > 1:398953/1‑65 (MQ=255) tgtgATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTcc > 1:2985985/1‑64 (MQ=255) tCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGa > 1:3403554/1‑65 (MQ=255) gcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGaacc > 1:1697929/1‑63 (MQ=255) gcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAAc > 1:3949759/1‑65 (MQ=255) gcaggcaATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAAc > 1:1628025/1‑65 (MQ=255) ggcaGTGTGTCGCTGAAAAAAGAAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCAt < 1:2799225/65‑1 (MQ=255) ggcaATGTTTCGCTGTAAAAAGAAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCAt < 1:1435879/65‑1 (MQ=255) aTGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTg > 1:2021745/1‑65 (MQ=255) tgtCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGAtt > 1:1692247/1‑65 (MQ=255) tgtCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGAtt > 1:3626795/1‑65 (MQ=255) gtCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGAttt > 1:4393763/1‑65 (MQ=255) gtCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGAttt > 1:2948501/1‑65 (MQ=255) ccttGGAAAAATTAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGc < 1:2537241/62‑1 (MQ=255) cGCTGAAAAAAATAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGc < 1:2428889/65‑1 (MQ=255) gCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCa > 1:3551350/1‑65 (MQ=255) gCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCa > 1:1635208/1‑65 (MQ=255) gtaaaagtaaaaGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg < 1:2312451/65‑1 (MQ=255) gaaaaaGAAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg < 1:3667242/63‑1 (MQ=255) gaaaaaGAAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg < 1:4063048/63‑1 (MQ=255) gaaaaaGAAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg < 1:2000139/63‑1 (MQ=255) gTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCAtgtg > 1:1087803/1‑65 (MQ=255) aaaGTAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACa < 1:2706823/65‑1 (MQ=255) aaaGGAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACa < 1:4306367/65‑1 (MQ=255) aaaGAAAAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACa < 1:1793655/65‑1 (MQ=255) aagaagaagCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACaa < 1:4593809/65‑1 (MQ=255) gTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAata < 1:4020542/65‑1 (MQ=255) gTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAata > 1:4395734/1‑65 (MQ=255) gTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAata > 1:4439334/1‑65 (MQ=255) gTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAata > 1:2931262/1‑65 (MQ=255) | CGTGCTTCATCTGACGGATAGCTGCTTTGTGATCGGCAGGCAATGTGTCGCTGGTAAAAGTAGAAGCCATAACCACCTCTCATTAACTTTCATTCGAAACCATTCTGATTTGCATGTGACAATA > NC_000913/858992‑859115 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |