Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 4,223,597 | 0 | C | T | 29.5% | 176.3 / 39.0 | 88 | R11K (AGA→AAA) | iclR | DNA‑binding transcriptional repressor IclR |
Reads supporting (aligned to +/- strand): ref base C (17/45); new base T (26/0); total (43/45) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.88e-11 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTG > NC_000913/4223549‑4223646 | aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc < 1:3489312/49‑1 (MQ=255) aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc < 1:2572243/49‑1 (MQ=255) aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc < 1:5709994/49‑1 (MQ=255) aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTc < 1:1668155/49‑1 (MQ=255) aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTg < 1:8097743/50‑1 (MQ=255) aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTg < 1:153304/50‑1 (MQ=255) aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTg > 1:5554294/1‑50 (MQ=38) aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTg > 1:5725824/1‑50 (MQ=38) aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGc < 1:760473/50‑1 (MQ=255) aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTgg > 1:3763480/1‑49 (MQ=38) gACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGcc < 1:6834958/50‑1 (MQ=255) aCTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcg > 1:4828880/1‑50 (MQ=38) cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCcgc < 1:1146999/50‑1 (MQ=255) cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGTcgc > 1:6895622/1‑50 (MQ=38) cTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCcgc > 1:7597003/1‑50 (MQ=38) tGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcgcg > 1:6658865/1‑50 (MQ=38) gAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCCGt > 1:7608458/1‑50 (MQ=38) aaCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCcgcgct > 1:4673084/1‑48 (MQ=37) aCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGttt > 1:5433212/1‑50 (MQ=38) cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcg < 1:6775865/50‑1 (MQ=255) cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcg < 1:2015457/50‑1 (MQ=255) cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGgtttg > 1:6596963/1‑48 (MQ=25) cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGgtttg > 1:8133371/1‑48 (MQ=37) cTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTcg > 1:3698930/1‑50 (MQ=39) tGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTtgc > 1:4767192/1‑50 (MQ=25) tCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTcgc < 1:4234470/48‑1 (MQ=255) cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCggg < 1:248691/49‑1 (MQ=255) cAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGa < 1:3229893/50‑1 (MQ=255) cAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTGCGCGGTTCGCggg > 1:1447648/1‑49 (MQ=25) aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCgggga > 1:7585229/1‑48 (MQ=255) aGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCgggg > 1:5446809/1‑48 (MQ=255) gTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTTGCGGGGAt > 1:2600125/1‑50 (MQ=38) gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGGTTCGCGGGGAt > 1:3385052/1‑50 (MQ=255) gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAAt < 1:3238814/50‑1 (MQ=255) gTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAAt < 1:1528631/50‑1 (MQ=255) tCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATg < 1:1369607/50‑1 (MQ=255) tCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATg < 1:2998199/50‑1 (MQ=255) cGCTGGTGCGGTGGCAACGGCGGGGTTTTTTTCCCGTTTTGCGGGGATgg > 1:1963089/1‑50 (MQ=16) cGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTCGCGGGGATgg > 1:4303269/1‑50 (MQ=255) gCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCGCGGTTCGCGGGGATtgg > 1:2659997/1‑50 (MQ=25) gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATggg < 1:7876877/50‑1 (MQ=255) gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATggg < 1:4403070/50‑1 (MQ=255) gCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATgg < 1:4153162/49‑1 (MQ=255) gCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGGTTTGCGGGGAAggg > 1:517991/1‑50 (MQ=255) gCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGCTTCGCGGGAATggg > 1:3792885/1‑50 (MQ=37) gCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGATggg > 1:5953791/1‑50 (MQ=255) gCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGGGAAggg > 1:1637315/1‑50 (MQ=255) cTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGt < 1:6477979/50‑1 (MQ=255) cTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTCCGCGGGAATGGGt < 1:7673466/50‑1 (MQ=255) cTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGt > 1:4886117/1‑50 (MQ=255) cTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGCGGGGATGGGt > 1:1798908/1‑50 (MQ=255) ggTGCGGTGGCAACGGCGGGTTTTTTTCCGCGTTTCCCGGGGATAGGGGc > 1:4222346/1‑50 (MQ=19) gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTg < 1:5890881/48‑1 (MQ=255) gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGc < 1:2814850/49‑1 (MQ=255) gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCg < 1:880248/50‑1 (MQ=255) gTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCg < 1:111783/50‑1 (MQ=255) tGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGa < 1:2870279/50‑1 (MQ=255) cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGAcc < 1:259592/50‑1 (MQ=255) cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGAcc < 1:3332734/50‑1 (MQ=255) ggTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCa < 1:1432568/50‑1 (MQ=255) gTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCAt < 1:7454228/50‑1 (MQ=255) gTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCAt < 1:5973117/50‑1 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGa < 1:4404083/49‑1 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGa < 1:2887654/49‑1 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAc < 1:1593132/50‑1 (MQ=255) gCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGAc < 1:5440991/50‑1 (MQ=255) aCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt < 1:4894489/50‑1 (MQ=255) aCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGt < 1:3568093/50‑1 (MQ=255) cggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtc < 1:1710235/50‑1 (MQ=255) ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtct > 1:1689691/1‑50 (MQ=255) gcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtctc < 1:136963/50‑1 (MQ=255) gcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGtctc < 1:837013/50‑1 (MQ=255) cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct > 1:3857229/1‑50 (MQ=255) cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct > 1:2298915/1‑50 (MQ=255) cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct < 1:7989724/50‑1 (MQ=255) cggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTctct > 1:8238211/1‑50 (MQ=255) gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt > 1:1660802/1‑50 (MQ=255) gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt > 1:130813/1‑50 (MQ=255) gggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtt > 1:844830/1‑50 (MQ=255) ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt < 1:765380/50‑1 (MQ=255) ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt > 1:3040590/1‑50 (MQ=255) ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt > 1:5055051/1‑50 (MQ=255) ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttt < 1:8078027/50‑1 (MQ=255) ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttttt < 1:6731052/50‑1 (MQ=255) ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttttt < 1:1700070/50‑1 (MQ=255) ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCttttt > 1:5559726/1‑50 (MQ=255) tttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCtttttt > 1:3145805/1‑50 (MQ=255) ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc > 1:7292687/1‑50 (MQ=255) ttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc > 1:2490530/1‑50 (MQ=255) tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc < 1:2789163/49‑1 (MQ=255) tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTc < 1:482721/49‑1 (MQ=255) tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt > 1:700868/1‑50 (MQ=255) cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt < 1:1776619/49‑1 (MQ=255) cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt < 1:5490392/49‑1 (MQ=255) cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCt < 1:2269911/49‑1 (MQ=255) cTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTg > 1:3470749/1‑50 (MQ=255) | AAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTG > NC_000913/4223549‑4223646 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |