library(dplyr)
library(ggplot2)
library(knitr)
library(readr)
Read data
kmerstream_data <- read_tsv("HG004/sequence.index.AJtrio_Illumina_2x250bps_02192016.bfc.kmerstream.tsv")
abyss_data <- read_tsv("abyss-bfc-counts.tsv")
data <- rbind(abyss_data, kmerstream_data %>% transmute(Tool = "KmerStream 1.0", k, F0))
Plot F0 vs k
ggplot(data) + aes(x = k, y = F0, colour = Tool) + xlim(0, NA) + ylim(0, NA) + geom_point(shape = 1) +
geom_vline(xintercept = c(64, 128, 192), alpha = 0.2)
Tabulate data
kable(data)
| ABySS 1.9.0 |
32 |
7175140434 |
| ABySS 1.9.0 |
64 |
9097982371 |
| ABySS 1.9.0 |
96 |
9422137746 |
| ABySS 1.9.0 |
104 |
9380437212 |
| ABySS 1.9.0 |
112 |
9305498249 |
| ABySS 1.9.0 |
120 |
9199641105 |
| ABySS 1.9.0 |
128 |
9063966262 |
| ABySS 1.9.0 |
136 |
8899021013 |
| ABySS 1.9.0 |
140 |
8805573981 |
| ABySS 1.9.0 |
142 |
8756146067 |
| ABySS 1.9.0 |
143 |
8730751835 |
| ABySS 1.9.0 |
144 |
8704901295 |
| ABySS 1.9.0 |
145 |
8678609309 |
| ABySS 1.9.0 |
146 |
8651856109 |
| ABySS 1.9.0 |
148 |
8596966512 |
| ABySS 1.9.0 |
152 |
8481787855 |
| ABySS 1.9.0 |
160 |
8229897407 |
| ABySS 1.9.0 |
176 |
7639075163 |
| ABySS 1.9.0 |
192 |
6931146429 |
| ABySS 1.9.0 |
224 |
5097216308 |
| KmerStream 1.0 |
31 |
6583267269 |
| KmerStream 1.0 |
32 |
6833603685 |
| KmerStream 1.0 |
33 |
6933594507 |
| KmerStream 1.0 |
63 |
9108946311 |
| KmerStream 1.0 |
64 |
142381287401 |
| KmerStream 1.0 |
65 |
141238890517 |
| KmerStream 1.0 |
95 |
119036750375 |
| KmerStream 1.0 |
96 |
118010306557 |
| KmerStream 1.0 |
97 |
117448489335 |
| KmerStream 1.0 |
127 |
94776863611 |
| KmerStream 1.0 |
128 |
9200409063 |
| KmerStream 1.0 |
129 |
9138756101 |
| KmerStream 1.0 |
159 |
8394983190 |
| KmerStream 1.0 |
160 |
8313046852 |
| KmerStream 1.0 |
161 |
8317101431 |
| KmerStream 1.0 |
191 |
7102119285 |
| KmerStream 1.0 |
192 |
44943131884 |
| KmerStream 1.0 |
193 |
44192417397 |
| KmerStream 1.0 |
223 |
21204523750 |
| KmerStream 1.0 |
224 |
20400302262 |
| KmerStream 1.0 |
225 |
19638524129 |