Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,750,816 | 0 | A→C | 4.4% | 0.7 ‑30.1 0.0 | 232 | Y540S (TAC→TCC) | fhlA | DNA‑binding transcriptional activator |
Rejected: E-value exceeds prediction threshold. | ||||||||||
Rejected: Polymorphism prediction not supported by minimum number of reads on both strands. | ||||||||||
Rejected: Prediction has biased read strand distribution supporting polymorphism. | ||||||||||
Reads supporting (aligned to +/- strand): new base (0/16): ref base (190/26): total (190/42) | ||||||||||
Fisher's exact test strand distribution p-value = 8.40e-14 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
Only 100 of 262 total aligned reads displayed.
GTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCA > REL606/2750782‑2750850 | gTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTa < 1:8784092/35‑1 (MQ=255) gTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTa > 1:19860446/1‑35 (MQ=255) gTCGCCGACCGTGAGTTCCGTAGAGATCTCTATTa < 1:2195168/35‑1 (MQ=255) tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc > 1:3655817/1‑35 (MQ=255) tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc < 1:4930739/35‑1 (MQ=255) tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc < 1:6563476/35‑1 (MQ=255) tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc > 1:19770100/1‑35 (MQ=255) tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc < 1:18921557/35‑1 (MQ=255) tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc < 1:1008013/35‑1 (MQ=255) tCGCCGACCGTGAGTTCCGTAGCGATCTCTATTAc < 1:17772135/35‑1 (MQ=255) cGCCGACCGTTAGTTCCGTAGCGATCTCTATTAcc > 1:3146076/1‑35 (MQ=255) cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc > 1:3381714/1‑35 (MQ=255) cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc > 1:5150833/1‑35 (MQ=255) cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc > 1:13937704/1‑35 (MQ=255) cGCCGACCGTGAGTTCCGTAGCGATCTCTATTAcc < 1:13829340/35‑1 (MQ=255) gCCGACCGTGAGTTCCGTAGCGATCTCTATTACCg < 1:7468148/35‑1 (MQ=255) gCCGACCGTGAGTTCCGTAGCGATCTCTATTACCg > 1:754058/1‑35 (MQ=255) gCCGACCGTGAGTTCCGTAGCGATCTCTATTACCg < 1:6194970/35‑1 (MQ=255) ccGCCCGTGAGTTCCGTAGCGATCTCTATTACCGc < 1:19105114/35‑1 (MQ=255) ccGACCGTGAGTTCCGTAGCGATCTCTATTACCGc > 1:4795532/1‑35 (MQ=255) ccGACCGTGAGTTCCGTAGCGATCTCTATTACCGc > 1:6017864/1‑35 (MQ=255) cGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc > 1:3093451/1‑35 (MQ=255) cGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc > 1:11777490/1‑35 (MQ=255) cGACCGTGAGTTCCGTAGCGATCTCTATTACCGcc > 1:7848754/1‑35 (MQ=255) gACCGTGAGTTCCGTAGCGATCTCTATTACCGCCt > 1:4711865/1‑35 (MQ=255) gACCGTGAGTTCCGTAGCGATCTCTATTACCGCCt > 1:12212220/1‑35 (MQ=255) gACCGTGAGTTCCGTAGCGATCTCTATTACCGCCt > 1:8003024/1‑35 (MQ=255) aCCGTGAGTTCCGTAGCGATCTCTATTACCGCCTg > 1:390874/1‑35 (MQ=255) aCCGTGAGTTCCGTAGCGATCTCTATTACCGCCTg > 1:5433256/1‑35 (MQ=255) aCCGTGAGTTCCGTAGCGATCTCTATTACCGCCTg > 1:20421212/1‑35 (MQ=255) ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa > 1:6654030/1‑35 (MQ=255) ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa > 1:9554751/1‑35 (MQ=255) ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa > 1:10902437/1‑35 (MQ=255) ccGTGAGTTCCGTAGCGATCTCTATTACCGCCTGa > 1:16911773/1‑35 (MQ=255) cGTGAGTTCCGTATCGATCTCTATTACCGCCTGaa > 1:14652934/1‑35 (MQ=255) cGTGAGTTCCGTAGCGATCTCTATTACCGCCTGaa > 1:18403051/1‑35 (MQ=255) gTGAGTTCCGTAGCGATCTCTATTACCGCCTGAAc > 1:558058/1‑35 (MQ=255) gTGAGTTCCGTAGCGATCTCTATTACCGCCTGAAc > 1:3312353/1‑35 (MQ=255) gTGAGTTCCGTAGCGATCTCTATTACCGCCTGAAc > 1:8529125/1‑35 (MQ=255) tGAGTTCCGTAGCGATCTCTATTACCGCCTGAACg > 1:11944487/1‑35 (MQ=255) tGAGTTCCGTAGCGATCTCTATTACCGCCTGAACg > 1:8166572/1‑35 (MQ=255) gAGTTCCGTAGCGATCTCTATTACCGCCTGAACGt > 1:17854657/1‑35 (MQ=255) gAGTTCCGTAGCGATCTCTATTACCGCATGAACGt < 1:225565/35‑1 (MQ=255) aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa > 1:18529035/1‑35 (MQ=255) aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa > 1:12417063/1‑35 (MQ=255) aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa > 1:1603349/1‑35 (MQ=255) aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa > 1:7390132/1‑35 (MQ=255) aGTTCCGTAGCGATCTCTATTACCGCCTGAACGTa > 1:2984608/1‑35 (MQ=255) ttCCGTAGCGATCTCTATTACCGCCTGAACGTAtt > 1:2763191/1‑35 (MQ=255) ccGTAGCGATCTCTATTACCGCCTGAACGTATTcc > 1:9320687/1‑35 (MQ=255) ccGTAGCGATCTCTATTACCGCCTGAACGTATTcc > 1:17577953/1‑35 (MQ=255) cGTAGCGATCTCTATTACCGCCTGAACGTATTccc > 1:19496014/1‑35 (MQ=255) cGTAGCGATCTCTATTACCGCCTGAACGTATTccc > 1:16186466/1‑35 (MQ=255) gTAGCGATCTCTATTACCGCCTGAACGTATTCCCg > 1:11058281/1‑35 (MQ=255) gTAGCGATCTCTATTACCGCCTGAACGTATTCCCg > 1:5853831/1‑35 (MQ=255) tAGCGATCTCTATTACCGCCTGAACGTATTCCCGa > 1:10693319/1‑35 (MQ=255) tAGCGATCTCTATTACCGCCTGAACGTATTCCCGa > 1:1305870/1‑35 (MQ=255) tAGCGATCTCTATTACCGCCTGAACGTATTCCCGa > 1:6144740/1‑35 (MQ=255) aGCGTTCTTTTTTCCCGCCTGAACGTATTCCCGAt < 1:4255953/35‑1 (MQ=16) aGCGATCTCTATTACCGCCTGAACGTATTCCCGAt > 1:17228024/1‑35 (MQ=255) aGCGATCTCTATTACCGCCTGAACGTATTCCCGAt > 1:12037561/1‑35 (MQ=255) aGCGATCTCTATTACCGCCTGAACGTATTCCCGAt > 1:4320537/1‑35 (MQ=255) gCGATCTCTATTCCAGCCTGAACGTATTCCCGAtt < 1:18266408/35‑1 (MQ=37) gCGATCTCTATTACCGCCTGAACGTATTCCCGAtt > 1:6839034/1‑35 (MQ=255) gCGATCTCTATTACCGCCTGAACGTATTCCCGAtt > 1:14350798/1‑35 (MQ=255) cGATCTCTATTCCCGCCTGAACGTATTCCCGATTc < 1:7099606/35‑1 (MQ=255) cGATCTCTATTACCGCCTGAACGTATTCCCGATTc > 1:11343185/1‑35 (MQ=255) cGATCTCTATTACCGCCTGAACGTATTCCCGATTc > 1:2368706/1‑35 (MQ=255) gATTTCTTTTCCCGCCTGAACGTATTCCCGATTCa < 1:19974452/35‑1 (MQ=21) gATCTCTATTACCGCCTGAACGTATTCCCGATTCa > 1:6397641/1‑35 (MQ=255) gATCTCTATTACCGCCTGAACGTATTCCCGATTCa > 1:4423550/1‑35 (MQ=255) gATCTCTATTACCGCCTGAACGTATTCCCGATTCa > 1:18072988/1‑35 (MQ=255) tcttctTTTCCCGCCTGAACGTATTCCCGATTCAc < 1:12967297/32‑1 (MQ=21) aTCTCTATTACCGCCTGAACGTATTCCCGATTCAc > 1:6920583/1‑35 (MQ=255) aTCTCTATTACCGCCTGAACGTATTCCCGATTCAc > 1:10002627/1‑35 (MQ=255) tctcTATTACCGCCTGAACGTATTCCCGATTCAcc > 1:10595106/1‑35 (MQ=255) ctttATTCCCGCCTGAACGTATTCCCGATTCACCt < 1:9115727/32‑1 (MQ=25) ctttATTCCCGCCTGAACGTATTCCCGATTCACCt < 1:11151828/32‑1 (MQ=25) ctctATTCCCGCCTGAACGTATTCCCGATTCACCt < 1:19413268/35‑1 (MQ=255) ctctATTACCGCCTGAACGTATTCCCGATTCACCt > 1:8268196/1‑35 (MQ=255) ctctATTACCGCCTGAACGTATTCCCGATTCACCt > 1:8677001/1‑35 (MQ=255) ctctATTACCGCCTGAACGTATTCCCGATTCACCt > 1:2127141/1‑35 (MQ=255) tctATTACCGCCTGAACGTATTCCCGATTCACCTg > 1:17975079/1‑35 (MQ=255) cctATTACCGCCTGAACGTATTCCCGATTCACCTg > 1:16465849/2‑35 (MQ=255) ctATTACCGCCTGAACGTATTCCCGATTCACCTgc > 1:5942607/1‑35 (MQ=255) ctATTACCGCCTGAACGTATTCCCGATTCACCTgc > 1:12485529/1‑35 (MQ=255) ctATTACCGCCTGAACGTATTCCCGATTCACCTgc > 1:20283211/1‑35 (MQ=255) ttttACCGCCTGAACGTATTCCCGATTCACGTgcc < 1:14963496/33‑1 (MQ=37) tATTACCGCCTGAACGTATTCCCGATTCACCTgcc > 1:13444336/1‑35 (MQ=255) tATTACCGCCTGAACGTATTCCCGACTCATCTgcc > 1:16109855/1‑35 (MQ=25) aTTACCGCCTGAACGTATTCCCGATTCACCTgccg > 1:3752167/1‑35 (MQ=255) ttACCGCCTGAACGTATTCCCGATTCACCTgccgc > 1:15857354/1‑35 (MQ=255) tACCGCCTGAACGTATTCCCGATTCACCTgccgcc > 1:17397111/1‑35 (MQ=255) tACCGCCTGAACGTATTCCCGATTCACCTgccgcc > 1:13703072/1‑35 (MQ=255) aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa > 1:17379054/1‑35 (MQ=255) aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa > 1:7274107/1‑35 (MQ=255) aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa > 1:2706620/1‑35 (MQ=255) aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa > 1:15370106/1‑35 (MQ=255) aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa > 1:11206063/1‑35 (MQ=255) aCCGCCTGAACGTATTCCCGATTCACCTGCCGCCa > 1:5317774/1‑35 (MQ=255) | GTCGCCGACCGTGAGTTCCGTAGCGATCTCTATTACCGCCTGAACGTATTCCCGATTCACCTGCCGCCA > REL606/2750782‑2750850 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |
Reads not counted as support for junction |
read_name Not counted due to insufficient overlap past the breakpoint. |
read_name Not counted due to not crossing MOB target site duplication. |