Ns126:Encode Methylation
Jump to navigation
Jump to search
RRBS Data Analysis to Encode Project
Aim
- Aim 1: methylation block with RRBS dataset
Background
- 1, Encode: http://genome.ucsc.edu/ENCODE/downloads.html
- 2, Encode|RRBS1: http://hgdownload.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeHaibMethylRrbs/
- 3, Encode|RRBS2: http://genome.ucsc.edu/cgi-bin/hgTrackUi?hgsid=437674359_aUhx08DjchWwtBjyCYv61EB7Yy8S&c=chr1&g=wgEncodeHaibMethylRrbs
- big problem for the fastq files from the above link. 1) not fastq but rather RAS file 2) reads is not true reads.
- 4, Encode|Methy450K: http://hgdownload.cse.ucsc.edu/goldenPath/hg19/encodeDCC/wgEncodeHaibMethyl450/
- 5, Encode|software: https://www.encodeproject.org/software/
Method
Data Download
Fastq and bed files can be downloaded from Encode Project. 101 RRBS data (bed files) were downloaded. 2,646,999 CpG loci were covered by 101 RRBS data while 866,979 CpG loci (32.8%) were detected in at least 80% samples.
Result
Alignment
- Maybe walltime time is too short , not all the samples were aligned completely, so I extend the walltime to 72 hours.
#PBS -q glean #PBS -l nodes=1:ppn=8 #PBS -l walltime=72:00:00 bismark --bowtie2 --phred64-quals --fastq -L 30 -N 1 --multicore 2 /home/shg047/db/hg19/meth/bismark ../fastq_trim/HOT197_trimmed.fq.gz -o ../bam2
Sample | N(reads) | N(mapped) | P(mapping) | N(C) | N(MCPG) | N(MCHG) | N(MCHH) | N(UCPG) | N(UCHG) | N(UCHH) | P(MCPG) | P(MCHG) | P(MCHH) |
ENCFF000LUP | 12830405 | 2122731 | 16.50% | 23121452 | 1755272 | 39237 | 63026 | 4964667 | 5412746 | 10886504 | 26.10% | 0.70% | 0.60% |
ENCFF000LUQ | 5561765 | 2658766 | 47.80% | 29121454 | 1939714 | 57696 | 100525 | 7299273 | 6782758 | 12941488 | 21.00% | 0.80% | 0.80% |
ENCFF000LVA | 8018858 | 1683010 | 21.00% | 15937339 | 1377038 | 59505 | 190243 | 2691809 | 3785119 | 7833625 | 33.80% | 1.50% | 2.40% |
ENCFF000LVB | 11054693 | 2404800 | 21.80% | 22228173 | 1455733 | 40643 | 119125 | 3161543 | 4987120 | 12464009 | 31.50% | 0.80% | 0.90% |
ENCFF000LVE | 10736789 | 3048286 | 28.40% | 32772970 | 1631989 | 49822 | 88996 | 8206874 | 7496120 | 15299169 | 16.60% | 0.70% | 0.60% |
ENCFF000LVK | 8542224 | 2399071 | 28.10% | 27152898 | 1893547 | 51170 | 70824 | 6965949 | 6095807 | 12075601 | 21.40% | 0.80% | 0.60% |
ENCFF000LVN | 13332915 | 4317747 | 32.40% | 31477310 | 1393770 | 79434 | 352105 | 2252186 | 6961310 | 20438505 | 38.20% | 1.10% | 1.70% |
ENCFF000LVO | 16368796 | 2638210 | 16.10% | 20542040 | 994068 | 43080 | 182778 | 2631559 | 4373546 | 12317009 | 27.40% | 1.00% | 1.50% |
ENCFF000LVR | 16091088 | 2615346 | 16.30% | 22490847 | 1312651 | 47566 | 188989 | 3319980 | 4807814 | 12813847 | 28.30% | 1.00% | 1.50% |
ENCFF000LVU | 13225098 | 1999398 | 15.10% | 18717042 | 1158340 | 37399 | 98519 | 2836054 | 4402056 | 10184674 | 29.00% | 0.80% | 1.00% |
ENCFF000LVV | 11857985 | 3571402 | 30.10% | 38309965 | 2844978 | 73249 | 122725 | 8543749 | 9108762 | 17616502 | 25.00% | 0.80% | 0.70% |
ENCFF000LWE | 21908830 | 2413683 | 11.00% | 24666773 | 1746092 | 71933 | 124011 | 4886835 | 5721680 | 12116222 | 26.30% | 1.20% | 1.00% |
ENCFF000LWL | 16598050 | 3530912 | 21.30% | 32715178 | 1675159 | 59309 | 192601 | 5947741 | 7094553 | 17745815 | 22.00% | 0.80% | 1.10% |
ENCFF000LWP | 16429464 | 2721119 | 16.60% | 28744997 | 1773533 | 43976 | 80430 | 6212197 | 6572632 | 14062229 | 22.20% | 0.70% | 0.60% |
ENCFF000LWS | 6325597 | 1224397 | 19.40% | 13284708 | 1257780 | 29276 | 46149 | 2987604 | 3195673 | 5768226 | 29.60% | 0.90% | 0.80% |
ENCFF000LWW | 16821394 | 2133899 | 12.70% | 22603245 | 1722157 | 39949 | 68530 | 4343561 | 5191728 | 11237320 | 28.40% | 0.80% | 0.60% |
ENCFF000LWX | 21208447 | 13218301 | 62.30% | 145215571 | 9668358 | 223882 | 356827 | 34200783 | 33840207 | 66925514 | 22.00% | 0.70% | 0.50% |
ENCFF000LXB | 18551059 | 3078942 | 16.60% | 33082876 | 1775486 | 66300 | 106053 | 8382345 | 7598213 | 15154479 | 17.50% | 0.90% | 0.70% |
Sample | N(reads) | N(mapped) | P(mapping) | N(C) | N(MCPG) | N(MCHG) | N(MCHH) | N(UCPG) | N(UCHG) | N(UCHH) | P(MCPG) | P(MCHG) | P(MCHH) |
ENCFF000LUP | 12830405 | 2122731 | 16.50% | 23121452 | 1755272 | 39237 | 63026 | 4964667 | 5412746 | 10886504 | 26.10% | 0.70% | 0.60% |
ENCFF000LUQ | 5561765 | 2658766 | 47.80% | 29121454 | 1939714 | 57696 | 100525 | 7299273 | 6782758 | 12941488 | 21.00% | 0.80% | 0.80% |
ENCFF000LVA | 8018858 | 1683010 | 21.00% | 15937339 | 1377038 | 59505 | 190243 | 2691809 | 3785119 | 7833625 | 33.80% | 1.50% | 2.40% |
ENCFF000LVB | 11054693 | 2404800 | 21.80% | 22228173 | 1455733 | 40643 | 119125 | 3161543 | 4987120 | 12464009 | 31.50% | 0.80% | 0.90% |
ENCFF000LVE | 10736789 | 3048286 | 28.40% | 32772970 | 1631989 | 49822 | 88996 | 8206874 | 7496120 | 15299169 | 16.60% | 0.70% | 0.60% |
ENCFF000LVK | 8542224 | 2399071 | 28.10% | 27152898 | 1893547 | 51170 | 70824 | 6965949 | 6095807 | 12075601 | 21.40% | 0.80% | 0.60% |
ENCFF000LVN | 13332915 | 4317747 | 32.40% | 31477310 | 1393770 | 79434 | 352105 | 2252186 | 6961310 | 20438505 | 38.20% | 1.10% | 1.70% |
ENCFF000LVO | 16368796 | 2638210 | 16.10% | 20542040 | 994068 | 43080 | 182778 | 2631559 | 4373546 | 12317009 | 27.40% | 1.00% | 1.50% |
ENCFF000LVR | 16091088 | 2615346 | 16.30% | 22490847 | 1312651 | 47566 | 188989 | 3319980 | 4807814 | 12813847 | 28.30% | 1.00% | 1.50% |
ENCFF000LVU | 13225098 | 1999398 | 15.10% | 18717042 | 1158340 | 37399 | 98519 | 2836054 | 4402056 | 10184674 | 29.00% | 0.80% | 1.00% |
ENCFF000LVV | 11857985 | 3571402 | 30.10% | 38309965 | 2844978 | 73249 | 122725 | 8543749 | 9108762 | 17616502 | 25.00% | 0.80% | 0.70% |
ENCFF000LWE | 21908830 | 2413683 | 11.00% | 24666773 | 1746092 | 71933 | 124011 | 4886835 | 5721680 | 12116222 | 26.30% | 1.20% | 1.00% |
ENCFF000LWL | 16598050 | 3530912 | 21.30% | 32715178 | 1675159 | 59309 | 192601 | 5947741 | 7094553 | 17745815 | 22.00% | 0.80% | 1.10% |
ENCFF000LWP | 16429464 | 2721119 | 16.60% | 28744997 | 1773533 | 43976 | 80430 | 6212197 | 6572632 | 14062229 | 22.20% | 0.70% | 0.60% |
ENCFF000LWS | 6325597 | 1224397 | 19.40% | 13284708 | 1257780 | 29276 | 46149 | 2987604 | 3195673 | 5768226 | 29.60% | 0.90% | 0.80% |
ENCFF000LWW | 16821394 | 2133899 | 12.70% | 22603245 | 1722157 | 39949 | 68530 | 4343561 | 5191728 | 11237320 | 28.40% | 0.80% | 0.60% |
ENCFF000LWX | 21208447 | 13218301 | 62.30% | 145215571 | 9668358 | 223882 | 356827 | 34200783 | 33840207 | 66925514 | 22.00% | 0.70% | 0.50% |
ENCFF000LXB | 18551059 | 3078942 | 16.60% | 33082876 | 1775486 | 66300 | 106053 | 8382345 | 7598213 | 15154479 | 17.50% | 0.90% | 0.70% |