Read alignment evidence... | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
seq id | position | change | freq | score | reads | annotation | genes | product | ||
* | REL606 | 2,767,883 | 0 | G→T | 6.5% | 20.7 ‑30.3 0.0 | 242 | L106M (CTG→ATG) | ygbE | conserved inner membrane protein |
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 (199/27): total (199/43) | ||||||||||
Fisher's exact test strand distribution p-value = 6.66e-14 | ||||||||||
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00 |
Only 100 of 298 total aligned reads displayed.
GAACAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCT > REL606/2767849‑2767917 | gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg > 1:6685169/1‑35 (MQ=255) gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg > 1:3850510/1‑35 (MQ=255) gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg > 1:8885956/1‑35 (MQ=255) gAACAGAGGAGACAAGAAAGTACAAATTAGCCCAg > 1:11554805/1‑35 (MQ=255) aaCAGAGGAGACAAGAAAGTACAAATTAGCCCAGt > 1:10113505/1‑35 (MQ=255) aaCAGAGGAGACAAGAAAGTACAAATTAGCCCAGt > 1:17357875/1‑35 (MQ=255) aaCAGAGGAGACAAGAAAGTACAAATTAGCCCAGt > 1:16316467/1‑35 (MQ=255) aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa > 1:7962821/1‑35 (MQ=255) aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa < 1:13685673/35‑1 (MQ=255) aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa > 1:15514304/1‑35 (MQ=255) aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa > 1:8602726/1‑35 (MQ=255) aCAGAGGAGACAAGAAAGTACAAATTAGCCCAGTa > 1:7705485/1‑35 (MQ=255) cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg > 1:4599382/1‑35 (MQ=255) cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg > 1:552971/1‑35 (MQ=255) cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg > 1:4192942/1‑35 (MQ=255) cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg > 1:1340492/1‑35 (MQ=255) cAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAg < 1:2893026/35‑1 (MQ=255) agagGAGACAAGAAAGTACAAATTAGCCCAGTAGc > 1:580501/1‑35 (MQ=255) agagGAGACAAGAAAGTACAAATTAGCCCAGTAGc < 1:16500379/35‑1 (MQ=255) agagGAGACAAGAAAGTACAAATTAGCCCAGTAGc > 1:18999800/1‑35 (MQ=255) gaggagACAAGAAAGTACAAATTAGCCCAGTAGcc > 1:4469134/1‑35 (MQ=255) gaggagACAAGAAAGTACAAATTAGCCCAGTAGcc < 1:10956169/35‑1 (MQ=255) agtagaCAAGAAAGTACAAATTAGCCCAGTAGCca > 1:9254041/4‑35 (MQ=255) aggagACAAGAAAGTACAAATTAGCCCAGTAGCca > 1:563373/1‑35 (MQ=255) aggagACAAGAAAGTACAAATTAGCCCAGTAGCca < 1:1951284/35‑1 (MQ=255) aggagACAAGAAAGTACAAATTAGCCCAGTAGCca < 1:14549755/35‑1 (MQ=255) aggagACAAGAAAGTACAAATTAGACCAGTAGCca < 1:11281854/35‑1 (MQ=255) gagaCAAGAAAGTACAAATTAGCCCAGTAGCcaca > 1:19932222/1‑35 (MQ=255) gagaCAAGAAAGTACAAATTAGCCCAGTAGCcaca > 1:18018284/1‑35 (MQ=255) agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt > 1:7488390/1‑35 (MQ=255) agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt > 1:983755/1‑35 (MQ=255) agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt < 1:18152709/35‑1 (MQ=255) agaCAAGAAAGTACAAATTAGCCCAGTAGCCACAt > 1:1058544/1‑35 (MQ=255) gaCAAGAAAGTACAAATTAGCCCAGTATCCACATa > 1:18698736/1‑35 (MQ=255) gaCAAGAAAGTACAAATTAGCCCAGTAGCCACATa > 1:4512041/1‑35 (MQ=255) aCAAGAAAGTACAAATTAGCCCAGTAGCCACATaa > 1:16682748/1‑35 (MQ=255) cAAGAAAGTACAAATTAGCCCAGTAGCCACATaaa > 1:19188399/1‑35 (MQ=255) aaGAAAGTACAAATTAGCCCAGTAGCCACATAAAc > 1:12687806/1‑35 (MQ=255) aGAAAGTACAAATTAGCCCAGTAGCCACATAAACa > 1:3300224/1‑35 (MQ=255) aGAAAGTACAAATTAGCCCAGTAGCCACATAAACa > 1:17172202/1‑35 (MQ=255) aGAAAGTACAAATTAGCCCAGTAGCCACATAAACa > 1:11765550/1‑35 (MQ=255) gAAAGTACAAATTAGCCCAGTAGCCACATAAACAg > 1:15807045/1‑35 (MQ=255) aaaGTACAAATTAGCCCAGTAGCCACATAAACAGt > 1:17877530/1‑35 (MQ=255) aaaGTACAAATTAGCCCAGTAGCCACATAAACAGt > 1:12004815/1‑35 (MQ=255) aaaGTACAAATTAGCCCAGTAGCCACATAAACAGt > 1:6569415/1‑35 (MQ=255) aaGTACAAATTAGCCCAGTAGCCACATAAACATTg > 1:8378580/1‑35 (MQ=255) aaGTACAAATTAGCCCAGTAGCCACATAAACATTg > 1:7819488/1‑35 (MQ=255) aaGTACAAATTAGCCCAGTAGCCACATAAACAGTg > 1:3493029/1‑35 (MQ=255) aGTACAAATTAGCCCAGTAGCCACATAAACAGTgc > 1:16053606/1‑35 (MQ=255) aGTACAAATTAGCCCAGTAGCCACATAAACAGTgc > 1:10473074/1‑35 (MQ=255) aGTACAAATTAGCCCAGTAGCCACATAAACAGTgc > 1:5427321/1‑35 (MQ=255) gTACAAATTAGCCCAGTAGCCACATAAACAGTgcg > 1:1147387/1‑35 (MQ=255) taaaaaTTACCCCTTTAGCCACATAAACAGTgcgc < 1:18475089/32‑1 (MQ=16) tACAAATTACCCCATTAGCCACATAAACAGTgcgc < 1:11202975/35‑1 (MQ=37) aCAAATTAGCCCAGTAGCCACATAAACAGTGCGcc > 1:9875537/1‑35 (MQ=255) aCAAATTAGCCCAGTAGCCACATAAACAGTGCGcc > 1:13519170/1‑35 (MQ=255) aaaTTAGCCTGGTAGCCACATAAACAGTGCGCCaa < 1:5943052/35‑1 (MQ=37) aaaTTAGCCCAGTAGCCACATAAACAGTGCGCCaa > 1:5094951/1‑35 (MQ=255) aaaTTAGCCCAGTAGCCACATAAACAGTGCGCCaa > 1:17716308/1‑35 (MQ=255) aaaTTACCCCATTAGCCACATAAACAGTGCGCCaa < 1:10040698/35‑1 (MQ=37) aaaTTACCCCATTAGCCACATAAACAGTGCGCCaa < 1:1032273/35‑1 (MQ=37) aaTTAGCCCATTAGCCACATAAACAGTGCGCCaaa < 1:16788812/35‑1 (MQ=255) aaTTAGCCCAGTAGCCACATAAACAGTGCGCCaaa > 1:16078584/1‑35 (MQ=255) aTTAGCCCATTAGCCACATAAACAGTGCGCCAAAc < 1:4103995/35‑1 (MQ=255) aTTAGCCCAGTAGCCACATAAACAGTGCGCCAAAc > 1:1313836/1‑35 (MQ=255) aTTAGCCAAGTAGCCACATAAACAGTGCGCCAAAc < 1:13369509/35‑1 (MQ=255) aTTACCCCATTAGCCACATAAACAGTGCGCCAAAc < 1:17490940/35‑1 (MQ=37) aTTACCCCAGTAGCCACATAAACAGTGCGCCAAAc < 1:18210587/35‑1 (MQ=255) ttAGCCCATTAGCCACATAAACAGTGCGCCAAACa < 1:12772930/35‑1 (MQ=255) ttAGCCCAGTAGCCACATAAACAGTGCGCCAAACa > 1:9540316/1‑35 (MQ=255) ttAGCCCAGTAGCCACATAAACAGTGCGCCAAACa < 1:10689521/35‑1 (MQ=255) tAGCCCATTAGCCACATAAACAGTGCGCCAAACAt < 1:2389524/35‑1 (MQ=255) tAGCCCAGTAGCCACATAAACAGTGCGCCAAACAt > 1:14255330/1‑35 (MQ=255) aaGCCCATTAGCCACATAAACAGTGAGCCAAACAt < 1:1472984/34‑1 (MQ=37) aGCCCAGTAGCCACATAAACAGTGCGCCAAACATa > 1:6902303/1‑35 (MQ=255) aGCCCAGTAGCCACATAAACAGTGCGCCAAACATa > 1:4916155/1‑35 (MQ=255) aGCCCAGTAGCCACATAAACAGTGCGCCAAACATa > 1:20501943/1‑35 (MQ=255) gCCCAGTAGCCACATAAACAGTGCGCCAAACATaa > 1:3374377/1‑35 (MQ=255) ccccATTAGCCACATAAACAGTGCGCCAAACATaa < 1:13090568/34‑1 (MQ=255) cccATTAGCCACATAAACAGTGCGCCAAACATAAt < 1:4399644/35‑1 (MQ=255) cccAGTAGCCACATAAACAGTGCGCCAAACATAAt > 1:6976202/1‑35 (MQ=255) cccAGTAGCCACATAAACAGTGCGCCAAACATAAt > 1:12256571/1‑35 (MQ=255) cccAGTAGCCACATAAACAGTGCGCCAAACATAAt > 1:16571716/1‑35 (MQ=255) cccAGTAGCCAAATAAACAGTGTGCCAAACATAAt < 1:1850144/35‑1 (MQ=37) cccAGTAGCAACATAAACAGTGCGCCAAACATAAt < 1:7360432/35‑1 (MQ=255) ccATTAGCCACATAAACAGTGCGCCAAACATAATg < 1:15244096/35‑1 (MQ=255) ccAGTAGCCACATAAACAGTGCGCCAAACATAATg > 1:8096664/1‑35 (MQ=255) ccAGTAGCCACATAAACAGTGCGCCAAACATAATg > 1:2088564/1‑35 (MQ=255) ccAGTAGCCACATAAACAGTGCGCCAAACATAATg > 1:18832731/1‑35 (MQ=255) ccAGTAGCCACATAAACAGTGCGCCAAACATAATg > 1:309486/1‑35 (MQ=255) ccAGTAGCCACATAAACAGTGCG‑CAAACATAATGc < 1:10025851/35‑1 (MQ=37) cAGTAGCCACATAAACAGTGCGCCAAACATAATGc > 1:14623231/1‑35 (MQ=255) cAGTAGCCACATAAACAGTGCGCCAAACATAATGc > 1:2666476/1‑35 (MQ=255) cAGTAGCCACATAAACAGTGCGCCAAACATAATGc > 1:19599034/1‑35 (MQ=255) aaGTAGCCACATAAACAGTGCGCCAAACATAATGc < 1:14463687/34‑1 (MQ=255) aGTAGCCACATAAACAGTGCGCCAAACATAATGcc > 1:3944413/1‑35 (MQ=255) gTAGCCACATAAACAGTGCGCCAAACATAATGCCt > 1:625763/1‑35 (MQ=255) gTAGCCACATAAACAGTGCGCCAAACATAATGCCt > 1:15083888/1‑35 (MQ=255) gTAGCCACATAAACAGTGCGCCAAACATAATGCCt > 1:15311097/1‑35 (MQ=255) gTACCCACATAAACAGTGCGCCAAACATAATGCCt > 1:20102210/1‑35 (MQ=255) | GAACAGAGGAGACAAGAAAGTACAAATTAGCCCAGTAGCCACATAAACAGTGCGCCAAACATAATGCCT > REL606/2767849‑2767917 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ 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. |