Plot Sealer results of PG29 mitochondrion

Load libraries

library(ggplot2)
library(knitr)

Read data

stats <- read.delim('stats.tsv')
data <- read.delim('gapstats.tsv')

Table of contiguity stats

kable(stats)

n	n.500	n.N50	min	N80	N50	N20	E.size	max	sum	name
78	78	11	2229	80869	157089	399506	215523	519137	5942459	pg29mt-8.fa
78	78	11	2229	81206	157361	399766	216050	520537	5958723	run1_scaffold.fa
202	15	7	500	529	592	652	588	670	8744	run1_merged.fa
n	n:1	n:N50	min	N80	N50	N20	E-size	max	sum	name
202	202	69	106	204	285	464	329	670	56664	run1_merged.fa

Scatter plot of actual gap size vs. estimated gap size

ggplot(na.omit(data), aes(x = gap, y = fill)) +
    geom_point() +
    geom_smooth(method = 'loess') +
    xlab('Estimated gap size by ABySS (bp)') +
    ylab('Actual gap size by Sealer (bp)')