| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | minE | 983,519 | G→A | 19.7% | intergenic (+431/‑137) | ydfC → / → dicB | hypothetical protein/cell division inhibition protein |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | minE | 983,519 | 0 | G | A | 19.7% | 127.2 / 23.7 | 61 | intergenic (+431/‑137) | ydfC/dicB | hypothetical protein/cell division inhibition protein |
| Reads supporting (aligned to +/- strand): ref base G (26/23); new base A (7/5); total (33/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGC > minE/983451‑983579 | tCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAg > 1:3054534/1‑71 (MQ=255) cAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAgc > 1:3741604/1‑71 (MQ=255) cAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAgc > 1:1048103/1‑71 (MQ=255) aCTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcgc < 1:3643948/70‑1 (MQ=255) cGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa > 1:3093459/1‑70 (MQ=255) tGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAg > 1:3683709/1‑71 (MQ=255) aCTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGt > 1:3041768/1‑70 (MQ=255) ctctGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTaaaa > 1:44310/1‑71 (MQ=255) ctGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAAc > 1:3472261/1‑70 (MQ=255) ccgccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTg < 1:3661696/44‑1 (MQ=255) ccgccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaa < 1:701072/71‑1 (MQ=255) gccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa > 1:3810154/1‑71 (MQ=255) gccCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa > 1:2859085/1‑71 (MQ=255) cccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa > 1:3801619/1‑70 (MQ=255) cccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaa < 1:2318810/67‑1 (MQ=255) cccTTTTTAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCAcg > 1:4553311/1‑49 (MQ=38) cccTTTTTAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCAcg > 1:5018031/1‑49 (MQ=38) cccTTTTTAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCAcg > 1:209816/1‑49 (MQ=38) ccTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAg > 1:2177498/1‑56 (MQ=255) tttttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCgg < 1:1431139/49‑1 (MQ=255) tttttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCgg < 1:2351362/49‑1 (MQ=255) tttAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAg > 1:5355539/1‑38 (MQ=37) tttAAAGTGAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTggca > 1:2288969/1‑47 (MQ=37) ttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaa < 1:1153227/63‑1 (MQ=255) ttAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAgtgt < 1:3432727/70‑1 (MQ=255) ttAAAGAGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAgtgt > 1:60501/1‑70 (MQ=255) tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa > 1:4790084/1‑47 (MQ=255) tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa > 1:3522728/1‑47 (MQ=255) tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGaa > 1:2777883/1‑47 (MQ=255) tAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGtt < 1:4752168/70‑1 (MQ=255) aaaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAAc > 1:4377844/1‑63 (MQ=255) aaaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAAc > 1:5115125/1‑63 (MQ=255) aaGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTAt < 1:1607623/70‑1 (MQ=255) tGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg < 1:3577234/37‑1 (MQ=255) tGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg < 1:4286088/37‑1 (MQ=255) tGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCAcg < 1:4858309/37‑1 (MQ=255) gAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaaaa > 1:4340373/1‑50 (MQ=255) gAATTTTGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaaaa > 1:1109611/1‑50 (MQ=255) aaTTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTg > 1:1791464/1‑70 (MQ=255) ttttGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa < 1:5535740/45‑1 (MQ=25) ttttGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa < 1:3726905/45‑1 (MQ=25) ttttGTGATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa < 1:1661474/45‑1 (MQ=25) ttGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGAtt > 1:5099483/1‑70 (MQ=255) gtgATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATtct < 1:3279245/70‑1 (MQ=255) gATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa < 1:4210935/39‑1 (MQ=25) gATGCGGTGAATGCGGCTAAGCGCACGTGGCACAGTTaa < 1:902769/39‑1 (MQ=25) tggtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGtt < 1:4610845/51‑1 (MQ=255) tggtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGtt < 1:1989506/51‑1 (MQ=255) gtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCg < 1:4843624/71‑1 (MQ=255) gtgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCg > 1:4111513/1‑71 (MQ=255) tgAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCgg < 1:1709916/71‑1 (MQ=255) gAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGc > 1:2829307/1‑71 (MQ=255) aaTGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa < 1:454490/39‑1 (MQ=255) aaTGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa < 1:2236473/39‑1 (MQ=255) aaTGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaaaa < 1:5371144/39‑1 (MQ=255) aaTGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGc < 1:2943790/70‑1 (MQ=255) tGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaa < 1:3934297/35‑1 (MQ=255) tGCGGCTGAGCGCACGCGGAACAGTTAAAACCaaa < 1:5357137/35‑1 (MQ=255) tGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAAc > 1:4410330/1‑38 (MQ=255) tGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGa > 1:4063963/1‑53 (MQ=255) tGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGa > 1:2614304/1‑53 (MQ=255) tGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGa > 1:4336838/1‑53 (MQ=255) | TCAGTTTTACTCCGTGACTGCTCTGCCGCCCTTTTTAAAGTGAATTTTGTGATGTGGTGAATGCGGCTGAGCGCACGCGGAACAGTTAAAACCAAAAACAGTGTTATGGGTGGATTCTCTGTATCCGGC > minE/983451‑983579 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |