Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | W3110S.gb | 1,467,504 | G→C | 8.3% | pseudogene (399/2559 nt) pseudogene (399/8622 nt) |
ydbA → ydbA → |
ECK1398:JW5802:b1401; predicted outer membrane protein, N‑ter fragment ECK1398:JW5802+JW1402:b4492; predicted outer membrane protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | W3110S.gb | 1,467,504 | 0 | G | C | 8.3% | 159.8 / 4.2 | 60 | pseudogene (399/2559 nt) pseudogene (399/8622 nt) | ydbA ydbA | ECK1398:JW5802:b1401; predicted outer membrane protein, N‑ter fragment ECK1398:JW5802+JW1402:b4492; predicted outer membrane protein |
Reads supporting (aligned to +/- strand): ref base G (36/19); new base C (2/3); total (38/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.46e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01 |
AATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTT > W3110S.gb/1467434‑1467569 | aaTCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACg > 1:4208939/1‑71 (MQ=255) tCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGcc < 1:766111/71‑1 (MQ=255) cAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGAAAATACCGGACGA‑‑‑TCCGGATGATACg > 1:2199225/1‑68 (MQ=255) aGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTa > 1:496090/1‑71 (MQ=255) gacgacACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTac < 1:3885763/71‑1 (MQ=255) cgacACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTa > 1:1201336/1‑68 (MQ=255) ccAGTTCCTCCAACTCCTGGCGGAGATGAGATAATACCAGACGATGGTCCGGATGATACCCCGacgc > 1:352620/1‑65 (MQ=255) ccAGTTCCTCCAACTCCTGGCGGAGATGAGATAATACCAGACGATGGTCCGGATGATACCCCGacgc > 1:712010/1‑65 (MQ=255) ccAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCa > 1:603488/1‑70 (MQ=255) cAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCaaa > 1:858956/1‑71 (MQ=255) cAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCaaa > 1:2523998/1‑71 (MQ=255) cAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCaaa > 1:52347/1‑71 (MQ=255) cAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCaa < 1:689211/70‑1 (MQ=255) aGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAAc < 1:243892/71‑1 (MQ=255) cAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTcc > 1:3174064/1‑58 (MQ=255) cAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑CCCGGATGATACGCCTACACCTCCCAAACCGGTTTCg > 1:559033/1‑71 (MQ=255) tCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGc > 1:3146289/1‑44 (MQ=255) tCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAAc > 1:2628887/1‑59 (MQ=255) ccTGGGGGAGATGAGATAATACCGGACGT‑‑‑TCCGGATGATACGCCTACACCTCCCaaa > 1:300360/1‑57 (MQ=255) ccTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCaaa > 1:4179810/1‑57 (MQ=255) cTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaat < 1:3288888/71‑1 (MQ=255) ggCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaataa > 1:1846704/1‑71 (MQ=255) ggCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaataa > 1:2472051/1‑71 (MQ=255) ggCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTAAACCTCCCAAACCGGTTTCGTTtaataa > 1:2141491/1‑71 (MQ=255) gCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGtt > 1:2053777/1‑64 (MQ=255) gCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaata < 1:1685429/69‑1 (MQ=255) gCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTaataat > 1:1605352/1‑71 (MQ=255) agatgagatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACg > 1:3076105/1‑71 (MQ=255) agatgagatAATACCAGACGATGGTCCGGATGATACCCCGACGCCTCCCAAACCGATTTCGTTTaataat < 1:1148578/70‑1 (MQ=255) gatgagatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGt > 1:2920928/1‑71 (MQ=255) atgagatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGttt < 1:3605089/59‑1 (MQ=255) atgagatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGt > 1:3110227/1‑70 (MQ=255) gagatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTAt > 1:475171/1‑71 (MQ=255) agatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTcc > 1:799477/1‑40 (MQ=255) agatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTcc > 1:2549320/1‑40 (MQ=255) agatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCg < 1:2993940/47‑1 (MQ=255) agatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTAtt > 1:2814206/1‑71 (MQ=255) agatAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTAt > 1:4369820/1‑70 (MQ=255) ataataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCg < 1:2860811/51‑1 (MQ=255) taataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATtctc < 1:3438765/71‑1 (MQ=255) aataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATtctc < 1:1098451/70‑1 (MQ=255) aataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATtctc < 1:1610063/70‑1 (MQ=255) aataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATtctc > 1:191183/1‑70 (MQ=255) aataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATtctc > 1:3707966/1‑70 (MQ=255) ataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAAc > 1:2310083/1‑40 (MQ=255) ataCCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTaataat > 1:3501356/1‑57 (MQ=255) ccGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGa > 1:2720103/1‑68 (MQ=255) ccGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATa < 1:2601986/70‑1 (MQ=255) cGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTta < 1:1255248/48‑1 (MQ=255) ggACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATaaaa < 1:821544/71‑1 (MQ=255) cGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaataa > 1:3329677/1‑48 (MQ=255) cGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATaaaa > 1:3493544/1‑68 (MQ=255) a‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTa < 1:2567326/54‑1 (MQ=255) tCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTtaa > 1:454652/1‑42 (MQ=255) tCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTa < 1:3486774/53‑1 (MQ=255) tCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATtc < 1:2096253/56‑1 (MQ=255) tCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATaa < 1:4358052/63‑1 (MQ=255) ccGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGaaa > 1:215066/1‑70 (MQ=255) atgatACCCCGACGCCTCCCAAACCGATTTCGTTTAATAATGACGTTAtcc < 1:4056020/51‑3 (MQ=255) atgatACCCCGACGCCTCCCAAACCGATTTCGTTTAATAATGACGTTAtcc < 1:816305/51‑3 (MQ=255) tgatACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGt > 1:2535955/1‑71 (MQ=255) atACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGtt > 1:2677586/1‑70 (MQ=255) | AATCAGACGACACTCCAGTTCCCCCAACTCCTGGCGGAGATGAGATAATACCGGACGA‑‑‑TCCGGATGATACGCCTACACCTCCCAAACCGGTTTCGTTTAATAATGACGTTATTCTCGATAAAACAGAAAAAACGTT > W3110S.gb/1467434‑1467569 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |