gravatar for mailard

2 hours ago by

Hi! I mapped methylome of HepG2 to reference human genome with bismark.
Is report of bismark okey? Or i have some troubles in mapping? And I have the question about coverage of WGBS. Is it okey that mean coverage of only one chromosome is 10X? I use this code for coverage: samtools coverage -h int.bam

Thanks for your answer!

Sequence pairs analysed in total: 695666796
Number of paired-end alignments with a unique best hit: 133717953
Mapping efficiency: 19.2%
Sequence pairs with no alignments under any condition: 548099782
Sequence pairs did not map uniquely: 13849061
Sequence pairs which were discarded because genomic sequence could not be extracted: 219

Number of sequence pairs with unique best (first) alignment came from the bowtie output:
CT/GA/CT: 66665199 ((converted) top strand)
GA/CT/CT: 0 (complementary to (converted) top strand)
GA/CT/GA: 0 (complementary to (converted) bottom strand)
CT/GA/GA: 67052535 ((converted) bottom strand)

Number of alignments to (merely theoretical) complementary strands being rejected in total: 0

Final Cytosine Methylation Report

Total number of C's analysed: 5669261795

Total methylated C's in CpG context: 153773530
Total methylated C's in CHG context: 5634249
Total methylated C's in CHH context: 54245515
Total methylated C's in Unknown context: 62135

Total unmethylated C's in CpG context: 205329168
Total unmethylated C's in CHG context: 1449704634
Total unmethylated C's in CHH context: 3800574699
Total unmethylated C's in Unknown context: 2788852

C methylated in CpG context: 42.8%
C methylated in CHG context: 0.4%
C methylated in CHH context: 1.4%
C methylated in unknown context (CN or CHN): 2.2%


modified 1 hour ago

2 hours ago

mailard 0

Source link