---
title: "01-genome-scaff-size"
output: html_document
---

```{bash}
pwd
```

```{bash}
awk '$0 ~ ">" {print c; c=0;printf substr($0,2,100) "\t"; } $0 !~ ">" {c+=length($0);} END { print c; }' O_lurida/data/Olurida_v081.fa \
> O_lurida/analyses/scaff_size.tab


```

```{r}
scaff <- read.csv("../analyses/scaff_size.tab", sep = "\t", header = FALSE)
```

```{r}
library(tidyverse)
```


```{r}
shortscaff <- scaff %>% filter(V2 < 50000)
```

