| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,180,341 | 0 | G | T | 26.6% | 168.3 / 21.7 | 65 | G200C (GGT→TGT) | ygiB | DUF1190 domain‑containing protein YgiB |
| Reads supporting (aligned to +/- strand): ref base G (13/34); new base T (17/0); total (31/34) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 8.68e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
ATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTAC > NC_000913/3180293‑3180390 | aTGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTgg > 1:5263775/1‑50 (MQ=255)aTGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGCTTTgg > 1:2740644/1‑50 (MQ=255) tGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTggg > 1:6436617/1‑49 (MQ=255) tGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTTGt > 1:6138900/1‑50 (MQ=255) tGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGt < 1:3025898/50‑1 (MQ=255) ggCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTTGt > 1:296082/1‑49 (MQ=255) gCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTtgtg > 1:7035389/1‑49 (MQ=255) aCCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTTGGGAAt > 1:2951791/1‑50 (MQ=255) ccAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTTGGGaaac > 1:5036824/1‑48 (MQ=255) cAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTTGGGAATCt > 1:520021/1‑50 (MQ=255) aaaaCCGGCGACCACCACTACCGTTACCCGTGGCGGCTTTTGTGAATCTg > 1:7235913/1‑50 (MQ=255) aaaaCCGGCGACCACCACTACCGTTACCCGTGGCGGCTTTTGTGAAACTg > 1:3601439/1‑50 (MQ=255) aaCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGtt < 1:7109920/50‑1 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTTGTGAATCTGTTg > 1:1776683/1‑50 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTTGTGAAACTGttt > 1:2842334/1‑49 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTg < 1:6347049/50‑1 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTg < 1:7252452/50‑1 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTg < 1:5953745/50‑1 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTg < 1:5164080/50‑1 (MQ=255) aCCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTg < 1:518918/50‑1 (MQ=255) ccGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑c < 1:766483/50‑1 (MQ=255) cGGCGACCACCACTACCGTTACCCGTGGCGGGTTTTGTGAAACTGGTggc > 1:5147545/1‑48 (MQ=255) ggCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑Cca < 1:3933745/50‑1 (MQ=255) ggCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑Cca < 1:7374362/50‑1 (MQ=255) ggCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑Cca < 1:2706088/50‑1 (MQ=255) gCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTT‑Cca > 1:2648433/1‑49 (MQ=255) gCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAAACTTTTT‑Cca > 1:6641660/1‑49 (MQ=38) cGACCACCACTACCGTTACCCGTGGCGGTTTTGGGGAATCTGTTG‑Ccaaa > 1:2356514/1‑50 (MQ=255) cGACCACCACTACCGTTACCCGTGGCGGTGTTCGGGCATCTGGTG‑Cccaa > 1:825209/1‑50 (MQ=255) gACCACCACTACCGTTACCCGTGGCGGTTTTTGGGAATCTGTTG‑Ccaaa > 1:5367887/1‑49 (MQ=255) gACCACCACTACCGTTACCCGTGGCGGTTTTGGGGAAACTGTTG‑Cccaa > 1:3354290/1‑49 (MQ=255) accaccACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑Ccaaaca < 1:6334657/50‑1 (MQ=255) ccaccaCTACCGTTACCCGTGGCGGTTTTTGTGAATTTGTTG‑Ccaaaca > 1:459593/1‑49 (MQ=39) ccaccaCTACCGTTACCCGTGGCGGTTTTGGTGAAACTGTTGTCcaaaa > 1:5204838/1‑48 (MQ=38) accacTACCGTTACCCGTGGCGGTTTTTGTGAAACTGTTG‑CCCAAAAAAg > 1:5305805/1‑50 (MQ=255) ccacTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGc < 1:5869609/50‑1 (MQ=255) cacTACCGTTACCCGTGGCGGTTTTTGTGGAACTCGTG‑CCAAATAAAgaa > 1:4678161/1‑48 (MQ=255) cacTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCa < 1:4239907/50‑1 (MQ=255) acTACCGTTACCCGTGGCGGTTTTTGTGAATCTGTTG‑CCCAACAAAACaa > 1:134556/1‑49 (MQ=255) acTACCGTTACCCGTGGCGGTTTTTGTGAATCTGTTG‑CCAAACAAAacac > 1:3257583/1‑50 (MQ=255) acTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCAc < 1:1902727/50‑1 (MQ=255) ccGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTa < 1:2225191/48‑1 (MQ=255) ccGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTAt < 1:5576480/49‑1 (MQ=255) gTTACCCGTGGCGGTTTTGGTGAATTTGTTG‑CCAAAAAAAACACTATg > 1:1103268/1‑48 (MQ=38) gTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCa < 1:6429492/50‑1 (MQ=255) gTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCa < 1:2793454/50‑1 (MQ=255) gTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCa < 1:3193396/50‑1 (MQ=255) ttACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAg < 1:5184015/50‑1 (MQ=255) aCCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGc < 1:4257945/49‑1 (MQ=255) cccGTGGCGGTTTTTGTGAATCTGTTG‑CCAAACAAAACACTATGTAGCg > 1:2151533/1‑49 (MQ=255) cccGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCg < 1:4848019/49‑1 (MQ=255) cccGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCg < 1:4274913/49‑1 (MQ=255) cccGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGt < 1:1788957/50‑1 (MQ=255) cccGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGt < 1:3778862/50‑1 (MQ=255) ccGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTa < 1:1545612/50‑1 (MQ=255) gTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGt < 1:5444021/50‑1 (MQ=255) tGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTg > 1:5316009/1‑50 (MQ=255) tGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTg < 1:1507769/50‑1 (MQ=255) ggCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTg > 1:7486909/1‑49 (MQ=255) ggCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTg > 1:190551/1‑49 (MQ=255) gCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGc > 1:809952/1‑49 (MQ=255) gCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCa < 1:910308/50‑1 (MQ=255) gCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCa < 1:4135085/50‑1 (MQ=255) gCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCa < 1:5166489/50‑1 (MQ=255) ggTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAAc < 1:779972/50‑1 (MQ=255) gTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAAcc > 1:172894/1‑50 (MQ=255) gTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAAcc < 1:4475024/50‑1 (MQ=255) tttGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:3172152/1‑50 (MQ=255) tttGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAACCgg > 1:5363772/1‑50 (MQ=255) ggTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTa < 1:4063269/49‑1 (MQ=255) ggTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTAc > 1:3854790/1‑50 (MQ=255) ggTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTAc > 1:4420749/1‑50 (MQ=255) | ATGGCACCAAAACCGGCGACCACCACTACCGTTACCCGTGGCGGTTTTGGTGAATCTGTTG‑CCAAACAAAGCACTATGCAGCGTAGTGCAACCGGTAC > NC_000913/3180293‑3180390 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |