Last updated: 2018-10-08

workflowr checks: (Click a bullet for more information)

• ✔ R Markdown file: up-to-date

  Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.

• ✔ Repository version: 50c8b76

  Great! You are using Git for version control. Tracking code development and connecting the code version to the results is critical for reproducibility. The version displayed above was the version of the Git repository at the time these results were generated.
  
  Note that you need to be careful to ensure that all relevant files for the analysis have been committed to Git prior to generating the results (you can use wflow_publish or wflow_git_commit). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:

  
  Ignored files:
      Ignored:    .DS_Store
      Ignored:    .Rhistory
      Ignored:    .Rproj.user/
      Ignored:    output/.DS_Store
  
  Untracked files:
      Untracked:  KalistoAbundance18486.txt
      Untracked:  analysis/genometrack_figs.Rmd
      Untracked:  analysis/ncbiRefSeq_sm.sort.mRNA.bed
      Untracked:  analysis/snake.config.notes.Rmd
      Untracked:  analysis/verifyBAM.Rmd
      Untracked:  data/18486.genecov.txt
      Untracked:  data/APApeaksYL.total.inbrain.bed
      Untracked:  data/NuclearApaQTLs.txt
      Untracked:  data/RNAkalisto/
      Untracked:  data/TotalApaQTLs.txt
      Untracked:  data/Totalpeaks_filtered_clean.bed
      Untracked:  data/YL-SP-18486-T-combined-genecov.txt
      Untracked:  data/YL-SP-18486-T_S9_R1_001-genecov.txt
      Untracked:  data/apaExamp/
      Untracked:  data/bedgraph_peaks/
      Untracked:  data/bin200.5.T.nuccov.bed
      Untracked:  data/bin200.Anuccov.bed
      Untracked:  data/bin200.nuccov.bed
      Untracked:  data/clean_peaks/
      Untracked:  data/comb_map_stats.csv
      Untracked:  data/comb_map_stats.xlsx
      Untracked:  data/comb_map_stats_39ind.csv
      Untracked:  data/combined_reads_mapped_three_prime_seq.csv
      Untracked:  data/ensemble_to_genename.txt
      Untracked:  data/filtered_APApeaks_merged_allchrom_refseqTrans.closest2End.bed
      Untracked:  data/filtered_APApeaks_merged_allchrom_refseqTrans.closest2End.noties.bed
      Untracked:  data/first50lines_closest.txt
      Untracked:  data/gencov.test.csv
      Untracked:  data/gencov.test.txt
      Untracked:  data/gencov_zero.test.csv
      Untracked:  data/gencov_zero.test.txt
      Untracked:  data/gene_cov/
      Untracked:  data/joined
      Untracked:  data/leafcutter/
      Untracked:  data/merged_combined_YL-SP-threeprimeseq.bg
      Untracked:  data/mol_overlap/
      Untracked:  data/nom_QTL/
      Untracked:  data/nom_QTL_opp/
      Untracked:  data/nom_QTL_trans/
      Untracked:  data/nuc6up/
      Untracked:  data/other_qtls/
      Untracked:  data/peakPerRefSeqGene/
      Untracked:  data/perm_QTL/
      Untracked:  data/perm_QTL_opp/
      Untracked:  data/perm_QTL_trans/
      Untracked:  data/reads_mapped_three_prime_seq.csv
      Untracked:  data/smash.cov.results.bed
      Untracked:  data/smash.cov.results.csv
      Untracked:  data/smash.cov.results.txt
      Untracked:  data/smash_testregion/
      Untracked:  data/ssFC200.cov.bed
      Untracked:  data/temp.file1
      Untracked:  data/temp.file2
      Untracked:  data/temp.gencov.test.txt
      Untracked:  data/temp.gencov_zero.test.txt
      Untracked:  output/picard/
      Untracked:  output/plots/
      Untracked:  output/qual.fig2.pdf
  
  Unstaged changes:
      Modified:   analysis/28ind.peak.explore.Rmd
      Modified:   analysis/39indQC.Rmd
      Modified:   analysis/PeakToGeneAssignment.Rmd
      Modified:   analysis/cleanupdtseq.internalpriming.Rmd
      Modified:   analysis/dif.iso.usage.leafcutter.Rmd
      Modified:   analysis/diff_iso_pipeline.Rmd
      Modified:   analysis/explore.filters.Rmd
      Modified:   analysis/overlapMolQTL.Rmd
      Modified:   analysis/overlap_qtls.Rmd
      Modified:   analysis/peakOverlap_oppstrand.Rmd
      Modified:   analysis/pheno.leaf.comb.Rmd
      Modified:   analysis/test.max2.Rmd
      Modified:   code/Snakefile
  
  

  Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.

Expand here to see past versions:

File	Version	Author	Date	Message
html	7369f8e	Briana Mittleman	2018-10-02	Build site.
Rmd	f66679b	Briana Mittleman	2018-10-02	overlap plots at peak level
html	2e0d959	Briana Mittleman	2018-10-01	Build site.
Rmd	9d6ee03	Briana Mittleman	2018-10-01	add 4su plots
html	da0f70e	Briana Mittleman	2018-10-01	Build site.
Rmd	fecb046	Briana Mittleman	2018-10-01	fix overlapfile
html	ac983db	Briana Mittleman	2018-10-01	Build site.
Rmd	35142fb	Briana Mittleman	2018-10-01	overlap QTL plots
html	149d033	Briana Mittleman	2018-09-26	Build site.
html	aaed5fd	Briana Mittleman	2018-09-26	Build site.
Rmd	eda266e	Briana Mittleman	2018-09-26	test peak to gene transcript dist
html	47cd4bf	Briana Mittleman	2018-09-26	Build site.
html	b1bcf99	Briana Mittleman	2018-09-25	Build site.
html	d3bb287	Briana Mittleman	2018-09-24	Build site.
Rmd	f7934ce	Briana Mittleman	2018-09-24	wflow_publish(c(“index.Rmd”, “39indQC.Rmd”))
html	f92a58f	Briana Mittleman	2018-09-06	Build site.
Rmd	46b7343	Briana Mittleman	2018-09-06	add overlap analysis with code to subset
html	b4bda1b	Briana Mittleman	2018-08-30	Build site.
html	e000583	Briana Mittleman	2018-08-30	Build site.
Rmd	a5f5276	Briana Mittleman	2018-08-30	initialize diff iso pipeline
html	c6dc97b	brimittleman	2018-08-28	Build site.
html	5564e25	brimittleman	2018-08-20	Build site.
Rmd	6b1b51c	brimittleman	2018-08-20	start qtl analsis, add to index
html	c67380e	brimittleman	2018-08-14	Build site.
html	38bfbaf	brimittleman	2018-08-09	Build site.
Rmd	03761cb	brimittleman	2018-08-09	start peak explore analysis for 54 libraries
html	e5a8da6	Briana Mittleman	2018-07-30	Build site.
Rmd	422a428	Briana Mittleman	2018-07-30	add peak cove pipeline and combined lane qc
html	0398fdb	Briana Mittleman	2018-07-25	Build site.
Rmd	2d41f11	Briana Mittleman	2018-07-25	add expand smash to index
html	3a5a8fe	Briana Mittleman	2018-07-25	Build site.
html	a34f9c1	Briana Mittleman	2018-07-17	Build site.
Rmd	ca9cee9	Briana Mittleman	2018-07-17	add smash test to index
html	c3ca148	Briana Mittleman	2018-07-17	Build site.
html	dc4a51a	Briana Mittleman	2018-07-16	Build site.
Rmd	4fede84	Briana Mittleman	2018-07-16	add eval brain analysis
html	cfc8f43	Briana Mittleman	2018-07-13	Build site.
Rmd	4a84769	Briana Mittleman	2018-07-13	add analysis for coparing coverage of RNAseq and 3’ seq
html	210cd2b	Briana Mittleman	2018-07-05	Build site.
Rmd	8b622c3	Briana Mittleman	2018-07-05	add denovo to index
html	ee777df	Briana Mittleman	2018-06-26	Build site.
Rmd	789d8ef	Briana Mittleman	2018-06-26	start test macs analysis. download package
html	30b415e	Briana Mittleman	2018-06-19	Build site.
html	08b5934	Briana Mittleman	2018-06-13	Build site.
Rmd	642faf0	Briana Mittleman	2018-06-13	start PAS enrichment analysis
html	4ef7d85	Briana Mittleman	2018-06-11	Build site.
Rmd	03c7b14	Briana Mittleman	2018-06-11	start A analysis
html	062808d	Briana Mittleman	2018-06-07	Build site.
html	86cfd9f	Briana Mittleman	2018-06-06	Build site.
Rmd	a19683c	Briana Mittleman	2018-06-06	start dif isoform analysis
html	710cf6a	Briana Mittleman	2018-05-29	Build site.
Rmd	d58bc13	Briana Mittleman	2018-05-29	start 200 bp analysis
html	96ecbb1	Briana Mittleman	2018-05-26	Build site.
Rmd	2076ce9	Briana Mittleman	2018-05-26	initial commit, gene level analysis
html	b5cdd59	Briana Mittleman	2018-05-26	Build site.
Rmd	9a49459	Briana Mittleman	2018-05-26	start map qc analysis
html	924952b	Briana Mittleman	2018-05-09	Build site.
Rmd	f8ce0fc	Briana Mittleman	2018-05-09	Update to 1.0
html	f8ce0fc	Briana Mittleman	2018-05-09	Update to 1.0
html	3dd80ef	Briana Mittleman	2018-05-08	Build site.
Rmd	519447d	Briana Mittleman	2018-05-08	Add total.nuc.rnaseq analysis
html	889c9c6	Briana Mittleman	2018-05-02	Build site.
html	9d196cf	Briana Mittleman	2018-05-02	Build site.
html	972ef3e	Briana Mittleman	2018-05-02	Build site.
Rmd	6cc5321	Briana Mittleman	2018-05-02	add home page information
html	e12dd6c	Briana Mittleman	2018-05-01	Build site.
Rmd	0e95115	Briana Mittleman	2018-05-01	add home page information
html	4c26d4c	Briana Mittleman	2018-05-01	Build site.
Rmd	539892a	Briana Mittleman	2018-05-01	Start workflowr project.

---

Briana Mittleman

This is the site I will use for my human 3’ seq project. I will have my workflow and analysis here. The three prime seq protocol I am using is Lexogen, QuantSeq 3’ mRNA-Seq library kit. The first analysis will be on the total and nucleic fractions of the LCL used in the net-seq analysis found at https://brimittleman.github.io/Net-seq/

Initial analysis on RNA seq of nuclear and total fractions:

• Initial data QC

Initial analysis on first 32 3’ seq libraries:

• Initial mapping QC
• Initial gene leve QC
  
• Compare coverage in Three prime seq to RNAseq
  
• Coverage in 200bp bins
  
• Differential isoform usage between fractions
  
• Explore filter cuttoffs
  
• Enrichment at PAS
  
• Second Lane data

Explore and filter mis-primming:

• Look at number of As in bins
  
• Create script to filter mapped A’s
  
• Test cleanUpdTSeq for calling internal priming
  
• Pipeline to create peaks and clean them

Peak calling:

• Testing macs2 on 3’seq data
  
• Call peak regions de-novo
• Compare peaks in Derti Brain set
  
• Test SMASH smoothing
  
• Expand SMASH analysis for use on midway
  
• Data Processing Figure design

Peak coverage/Nominal QTLs/Diff iso:

• Start evaluating peak quantificaiton
  
• Prepare the data for leafcutter QTL calling
• Call nominal QTLs
  
• Differential Isoform pipeline

Rerun analysis opp strand:

• Rerun QTL analysis with off strand

Evaluate and learn from QTLs:

• Overlap QTLs with other molecular QTLs

Analysis for all 39 indivduals:

• Read and peak QC plots

• APAqtl for full ind set

• Rethink Peak to Gene assignment

• Condition Total and Nuclear QTLs on each other

Molecular QTL analysis and Overlap:

• Rerun molecular QTLs

• Make overlap QTL plots

• Plot top apaQTLs in other phenotypes

---

---

---

This reproducible R Markdown analysis was created with workflowr 1.1.1

---