---
title: "04-exploring-3.1-deeper"
output:
  md_document:
    variant: markdown_github
---

# We have a blast table and some Uniprot data... lets see what the R can do..

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

```{r}
library(stringr)
```



```{r}

blast <- read.csv("../data/cbai_transcriptome_v3.1.blastx.outfmt6.txt", header = FALSE, sep="\t") %>%
  separate(col = V2, into = c("sp", "Accession_ID", "Gene"), sep = "\\|")
  


head(blast)
```

```{r}
up <- read.csv("../data/uniprot-reviewed-20210901.tab", header = TRUE, sep="\t")

head(up)
```
```{r}
comb <- left_join(blast, up, by =c("Accession_ID" = "Entry") )
```


```{r}
comb %>% filter(str_detect(Gene.ontology..biological.process., "immune"))
```
