::: {.cell layout-align="center" hash='05-eval-assemblies_cache/html/unnamed-chunk-8_f9afc777b70680584d1d4daaac57f149'}

```{.r .cell-code}
# Read in FastA file
fasta_filegawn <- "../results/GAWN-gene.fa"
sequences <- readDNAStringSet(fasta_filegawn)

# Calculate lengths
lengths <- width(sequences)

# Create a data frame
df <- data.frame(length = lengths)

# Plot histogram with ggplot2
ggplot(df, aes(x = length)) +
  geom_histogram(binwidth = 1000, fill = "blue", color = "black") +
  xlab("Sequence Length") +
  ylab("Frequency") +
  ggtitle("Histogram of Sequence Lengths") +
  theme_minimal()
```
:::