| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | Exported | 643,500 | G→C | 11.5% | T61S (ACC→AGC) | ompF ← | outer membrane porin 1a (Ia;b;F) |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | Exported | 643,500 | 0 | G | C | 11.5% | 227.1 / 15.0 | 78 | T61S (ACC→AGC) | ompF | outer membrane porin 1a (Ia;b;F) |
| Reads supporting (aligned to +/- strand): ref base G (43/26); new base C (3/6); total (46/32) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.49e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01 | |||||||||||
CCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTT > Exported/643432‑643569 | cccACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGt < 1:2949791/70‑1 (MQ=255) aCTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg > 1:4718141/1‑71 (MQ=255) aCTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATg > 1:4164276/1‑71 (MQ=255) tAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc > 1:3717188/1‑65 (MQ=255) tAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc > 1:844067/1‑65 (MQ=255) aaCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCa > 1:2067888/1‑68 (MQ=255) aaCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTg > 1:236655/1‑71 (MQ=255) aaCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTg > 1:5361736/1‑71 (MQ=255) aCCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGc < 1:5061059/71‑1 (MQ=255) cGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCa < 1:5320295/71‑1 (MQ=255) ggTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCAc > 1:4015265/1‑71 (MQ=255) ggTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCAc > 1:1943712/1‑71 (MQ=255) aGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGt < 1:3138920/70‑1 (MQ=255) aTCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCAt < 1:4611718/51‑1 (MQ=255) tCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCa > 1:2894089/1‑63 (MQ=255) cGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCa > 1:114233/1‑48 (MQ=255) cGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTc > 1:1219940/1‑52 (MQ=255) gAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt < 1:2414017/71‑1 (MQ=255) aaTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGttt < 1:1447785/71‑1 (MQ=255) aTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt < 1:4679742/71‑1 (MQ=255) aTTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt > 1:1169998/1‑71 (MQ=255) tGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa > 1:1289310/1‑62 (MQ=255) tGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTc > 1:4538312/1‑70 (MQ=255) tGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCa > 1:1560731/1‑71 (MQ=255) gATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCAc > 1:966891/1‑71 (MQ=255) tGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGt > 1:2483133/1‑70 (MQ=255) tGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGt < 1:2997788/70‑1 (MQ=255) aGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa < 1:5263873/70‑1 (MQ=255) aGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa < 1:39490/70‑1 (MQ=255) gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg < 1:335279/39‑1 (MQ=9) gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCg < 1:979285/39‑1 (MQ=9) gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct < 1:4658526/45‑3 (MQ=2) gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct < 1:4665175/45‑3 (MQ=2) gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct < 1:4325777/45‑3 (MQ=2) gTTTCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct < 1:670260/45‑3 (MQ=2) tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct > 1:4732126/1‑42 (MQ=2) tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct > 1:5123028/1‑42 (MQ=2) tttCGCCTTTAAAACCCAGACGCGCATAGCTCTGGTCGCCAtct > 1:1392188/1‑42 (MQ=2) tttCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCAt < 1:742732/42‑1 (MQ=255) tgccccTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGt < 1:4086236/48‑1 (MQ=255) cTTTGAATCCGATACGGGCATATGACATGTCGCCAt < 1:797672/36‑1 (MQ=255) cTTTGAAACCGATACGGGCATATGACATGTCGCCAt < 1:5139789/36‑1 (MQ=16) cTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTg < 1:4465853/50‑1 (MQ=255) cTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTAcc > 1:3924379/1‑64 (MQ=255) tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTg < 1:3289042/49‑1 (MQ=255) tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTc < 1:754256/54‑1 (MQ=255) tttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTc < 1:1441848/54‑1 (MQ=255) ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt > 1:3584994/1‑50 (MQ=255) ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt > 1:5041755/1‑50 (MQ=255) ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt > 1:4355507/1‑50 (MQ=255) ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtt > 1:4361800/1‑50 (MQ=255) ttAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaa > 1:1325758/1‑70 (MQ=255) tAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaaaa > 1:3530018/1‑71 (MQ=255) tAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaaa > 1:3906047/1‑70 (MQ=255) aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa > 1:2371831/1‑43 (MQ=255) aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa > 1:1926999/1‑43 (MQ=255) aaaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa > 1:2859589/1‑43 (MQ=255) aaaCCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAaat < 1:2786438/70‑1 (MQ=255) ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg < 1:1073922/71‑1 (MQ=255) ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg > 1:1483044/1‑71 (MQ=255) ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg < 1:3853084/71‑1 (MQ=255) ccAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg < 1:4367416/71‑1 (MQ=255) aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGaaaaa > 1:1557125/1‑64 (MQ=255) aaGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACACTTGGaaaaa > 1:1801515/1‑64 (MQ=255) aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa > 1:3435306/1‑36 (MQ=255) aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa > 1:5147641/1‑36 (MQ=255) aGACGGGCATAGGTCATGTCGCCATTGCCACCGTaa > 1:4179006/1‑36 (MQ=255) aGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCa < 1:3336347/70‑1 (MQ=255) gACGGGCATAGGTCATGTCACCATTGCCACCGTAACTGTTTTCACCGt < 1:3296958/48‑1 (MQ=255) aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt > 1:4928621/1‑41 (MQ=255) aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt > 1:1408658/1‑41 (MQ=255) aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt > 1:5432919/1‑41 (MQ=255) aCGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGtttt > 1:5410093/1‑41 (MQ=255) cGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGAc > 1:2221897/1‑70 (MQ=255) gggaTAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCaa < 1:4682230/68‑1 (MQ=255) gCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCaa > 1:1352152/1‑70 (MQ=255) cATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGa > 1:4410730/1‑65 (MQ=255) aGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCt > 1:505377/1‑71 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa > 1:4322509/1‑40 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa > 1:3942761/1‑40 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTa > 1:2099510/1‑40 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt > 1:1205247/1‑71 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt > 1:3284373/1‑71 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt > 1:5667624/1‑71 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt < 1:5544781/71‑1 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt > 1:5531783/1‑71 (MQ=255) ggTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCtt > 1:1396457/1‑71 (MQ=255) gTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg < 1:3724021/56‑1 (MQ=255) gTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATg < 1:2587379/56‑1 (MQ=255) | CCCACTGACCATAACCGGTCAGATCGGAATTGATTTGAGTTTCCCCTTTAAAACCAAGACGGGCATAGGTCATGTCGCCATTGCCACCGTAACTGTTTTCACCGTTACCCTTGGAAAAATAATGCAGACCAACAGCTT > Exported/643432‑643569 |
| 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: ‑ |