| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,180,341 | 0 | G | T | 24.6% | 140.6 / 19.4 | 61 | G200C (GGT→TGT) | ygiB | DUF1190 domain‑containing protein YgiB |
| Reads supporting (aligned to +/- strand): ref base G (30/16); new base T (15/0); total (45/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.36e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTAC > NC_000913/3180292‑3180390 | aaTGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTg > 1:7581569/1‑50 (MQ=255) tGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTTGt > 1:7483087/1‑50 (MQ=255) cACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTGaa > 1:3555032/1‑50 (MQ=255) cACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTgga > 1:3488659/1‑48 (MQ=255) cACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGaa > 1:1884622/1‑50 (MQ=255) cACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGaa < 1:3852698/50‑1 (MQ=255) ccAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTGAATc > 1:5959322/1‑50 (MQ=255) ccAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑ATTTGGTGAATc > 1:6110550/1‑50 (MQ=255) cAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCt < 1:602891/50‑1 (MQ=255) cAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCt < 1:6747453/50‑1 (MQ=255) cAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCt > 1:3051696/1‑50 (MQ=255) aaaaCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTTTTGAAACta > 1:172908/1‑49 (MQ=255) aaaaCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTGGATCtt > 1:3983620/1‑49 (MQ=255) aaaaCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTGAATCtt > 1:3196385/1‑49 (MQ=255) aaaaCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGGGAATCTg > 1:2284564/1‑50 (MQ=255) aaaCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTGAATCTGt > 1:6965826/1‑50 (MQ=255) aaaCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGt > 1:3808265/1‑50 (MQ=255) aaCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGtt < 1:5754429/50‑1 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGGATCTGttt > 1:177694/1‑49 (MQ=255) ggCGACCACCACTACCGTTACCCGTGGCGGCTTTTGGGGAA‑CTTTTGccc > 1:5514676/1‑49 (MQ=255) cGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGGATCTGTTGCcaa > 1:5646341/1‑49 (MQ=255) cGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCcaaa < 1:5667494/50‑1 (MQ=255) gACCACCACTACCGTTACCCGTGGCGGTTTTTGGGGAAT‑TGTTGCccaa > 1:6841902/1‑49 (MQ=255) gACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTTAACCTGTTTCcaaa > 1:1927424/1‑49 (MQ=255) gACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTGAATCTGTTGCccaa > 1:3069759/1‑49 (MQ=39) gACCACCACTACCGTTACCCGTGGCGG‑TTTTTGTGAATCTGTTGCcaaa > 1:3079276/1‑49 (MQ=255) gACCACCACTACCGTTACCCGTGGCGG‑TTTTTGGGAATCTGTTGCcaaa > 1:1286247/1‑49 (MQ=255) gACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCcaaa > 1:939680/1‑49 (MQ=255) accaccACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCcaaaca < 1:2964163/50‑1 (MQ=255) accaccACTACCGTTACCCGTGGCGG‑CTTTTGTGAATCCGGTGCccaac > 1:6882032/1‑49 (MQ=255) ccaccaCTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCcaaacaa < 1:8266617/50‑1 (MQ=255) caccacTACCGTTACCCGTGGCGG‑TTTTTGTGGAACTGTTTCccaacaa > 1:72128/1‑49 (MQ=25) caccacTACCGTTACCCGTGGCGG‑TTTTTGTGAATCTGTTGCccaacaa > 1:7213406/1‑49 (MQ=255) caccacTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCcaaacaaa < 1:4985093/50‑1 (MQ=255) caccacTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCcaaacaaa < 1:2409368/50‑1 (MQ=255) ccacTACCGTTACCCGTGGCGG‑TTTTGGTGGATTTGTTGCcaaacaaaa > 1:4465324/1‑48 (MQ=38) ccacTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTTCcaaacaaaa > 1:3718539/1‑48 (MQ=255) ccacTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCcaaaaaaaa > 1:1878957/1‑48 (MQ=255) ccacTACCGTTACCCGTGGCGG‑TTTTGGGGGAACTGTTTCcaaacaaaa > 1:5138207/1‑48 (MQ=255) acTACCGTTACCCGTGGCGG‑TTTTGGTGGAACTTTTGTCAAACAAAacac > 1:7995709/1‑50 (MQ=25) aCCGTTACCCGTGGCGG‑TTTTTGTGAATCTGTTGCCCAACAAAACACTAt > 1:7930674/1‑50 (MQ=255) gTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCa < 1:6634927/50‑1 (MQ=255) cccGTGGCGG‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCg < 1:1173828/49‑1 (MQ=255) tGGCGG‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTg > 1:6222996/1‑50 (MQ=255) ggCGG‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTg > 1:8270564/1‑49 (MQ=255) cGG‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCaa < 1:1206098/50‑1 (MQ=255) gg‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAAc < 1:6213298/50‑1 (MQ=255) gg‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAAc > 1:6692213/1‑50 (MQ=255) gg‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAAc > 1:1765634/1‑50 (MQ=255) gg‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAAc > 1:1910187/1‑50 (MQ=255) gg‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAAc < 1:2524944/50‑1 (MQ=255) g‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAAcc < 1:2757665/50‑1 (MQ=255) ttttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCg > 1:7901795/1‑50 (MQ=255) ttttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCg < 1:8571719/50‑1 (MQ=255) tttGGTGAATCTGTTGCCAAACAACGCACTATGCAGCGTAGTGCAACCgg > 1:5758000/1‑50 (MQ=255) tttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:7323909/1‑50 (MQ=255) tttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:6082622/1‑50 (MQ=255) tttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:4823381/1‑50 (MQ=255) tttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:5027440/1‑50 (MQ=255) ttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:7523657/1‑49 (MQ=255) ttGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:3952555/1‑49 (MQ=255) tGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTa > 1:324444/1‑50 (MQ=255) ggTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTAc > 1:6476952/1‑50 (MQ=255) | AATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGG‑TTTTGGTGAATCTGTTGCCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTAC > NC_000913/3180292‑3180390 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |