Ns126:Calendar/NOTES/2016-1-5: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Shicheng
>Shicheng
Line 22: Line 22:
===Methylation haploinfo===
===Methylation haploinfo===


bismark_methylation_extractor \
--single-end \
--bedGraph \
--buffer_size 2G \
--remove_spaces \
--zero_based \
--merge_non_CpG \
--comprehensive \
--output ../methyfreq \
.bam


===Methylation Haplotype Load Matrix===
===Methylation Haplotype Load Matrix===

Revision as of 09:29, 7 January 2016

Collection RRBS Dataset From Richard Myers, HAIB (Encode Project)

Data Description

Fastq Download

xargs -n 1 curl -O -L < haib.download.files.txt

Check Phred Score

perl ~/bin/checkphred.pl *fastq

Fastq Quality Control

trim_galore --phred64 --fastqc --illumina --non_directional --rrbs *.fastq

Fastq Alignment

bismark --bowtie2 --phred64-quals --fastq -L 30 -N 1 /home/shg047/db/aligndb/hg19/bismark -1 ENCFF000LUN.fastq_qual_trimmed.fastq -o ../bam

Methylation haploinfo

bismark_methylation_extractor \
--single-end \
--bedGraph \
--buffer_size 2G \
--remove_spaces \
--zero_based \
--merge_non_CpG \
--comprehensive \
--output ../methyfreq \
.bam

Methylation Haplotype Load Matrix

Evaluation of the data quality

Merge with our RRBS data