Skip to content

Commit

Permalink
Browse files Browse the repository at this point in the history
  • Loading branch information
ChrisMantegna committed Jun 2, 2023
2 parents dfc157d + 10862c3 commit 3e1ca51
Show file tree
Hide file tree
Showing 5 changed files with 2,608 additions and 82 deletions.
21 changes: 5 additions & 16 deletions code/05-kallisto.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,8 @@ knitr::opts_chunk$set(
```


**Data**\
# **Data**

Existing RNA-Seq data was retrieved from the following complete RNASeq data available in the [NCBI database](https://www.ncbi.nlm.nih.gov/):\
- SRR19782039 [Exposure to Valsartan & Carbamazepine](https://www.ncbi.nlm.nih.gov/sra/SRX15826291%5Baccn%5D)\
- SRR16771870 [Exposure to a synthetic hormone 17 a-Ethinylestradiol (EE2)](https://www.ncbi.nlm.nih.gov/sra/SRX12971792%5Baccn%5D)\
Expand All @@ -50,18 +51,13 @@ cp /home/shared/8TB_HDD_01/sr320/ncbi/ncbi_dataset/data/GCA_900618805.1/cds_from
```




### Pulling my reference genome from prefab code provided by NCBI
curled genome and copied just cds file

```{r, engine='bash', eval=TRUE}
head ../data/cds_from_genomic.fna
```

# Index Gene Set


# Create the index file to align my short read files to the genes from the MGAL_10
Create the index file to align my short read files to the genes from the MGAL_10
```{r, engine='bash'}
/home/shared/kallisto/kallisto \
index \
Expand All @@ -70,12 +66,7 @@ index \
```




# Where is data


#Runnning Kallisto
# Runnning Kallisto

```{r, engine='bash', eval=TRUE}
/home/shared/kallisto/kallisto quant
Expand All @@ -89,8 +80,6 @@ ls /home/shared/8TB_HDD_02/cnmntgna/GitHub/chris-musselcon/output/ncbi/
#mkdir ../output
#mkdir ../output/kallisto_01
#OLD CODE
find /home/shared/8TB_HDD_02/cnmntgna/GitHub/chris-musselcon/output/ncbi/*_1.fastq \
| xargs basename -s _1.fastq | xargs -I{} /home/shared/kallisto/kallisto \
quant -i ../data/MGAL_cds.index \
Expand Down
86 changes: 20 additions & 66 deletions code/05-kallisto.html

Large diffs are not rendered by default.

5 changes: 5 additions & 0 deletions code/07-blastingCDS.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -232,7 +232,12 @@ library(kableExtra)
library(knitr)
library(dplyr)
library(readr)
#new table process
full <-left_join(cdsDEG, uni, by = c("Column 2" = "Entry"))
write.table(full, file = "/home/shared/8TB_HDD_02/cnmntgna/GitHub/chris-musselcon/data/annoblast.tab", sep = "\t", row.names = FALSE)
```

```{r}
# Generate table
kbl(
head(
Expand Down
Original file line number Diff line number Diff line change
@@ -0,0 +1,11 @@
name: Document
title:
username:
account: rpubs
server: rpubs.com
hostUrl: rpubs.com
appId: https://api.rpubs.com/api/v1/document/1033647/813bc1a525b94585a7a2ddc1eb48233e
bundleId: https://api.rpubs.com/api/v1/document/1033647/813bc1a525b94585a7a2ddc1eb48233e
url: http://rpubs.com/publish/claim/1033647/6e7075d73dc34714a6321f17180b1d99
when: 1682452314.04458
lastSyncTime: 1682452314.04459
Loading

0 comments on commit 3e1ca51

Please sign in to comment.