methylseq data analysis using bismark and low alignment rate

I have methylseq data and used bismark for the analysis. but looking at the bam files, the alignment rate is only 8 %. then I checked the alignment using fastq_screen and got 70% alignment rate. do you know what could be the problem with bismark?

next-gen

• 786 views • link

written 5.6 years ago by Sara

0 answers

No answers yet.

More posts like this

zero counts for all genes in RNAseq data of Ferret
written 2.8 years ago by Sara

I have bulk RNA-seq data from Ferret and trying to get counts per gene. to do so I used hisat2 and got the genome from …
MethylR on bam files Gives Error
written 3.4 years ago by anasjamshed

I used bismark to get bam files and now trying to implement processBismarkAln() on it but getting the following error: Using htslib. Error: Unexpected cigar: …
no read counts per allele while using gatk ASEReadCounter
written 4.2 years ago by Sara

I am using gatk ASEReadCounter to get the read counts per allele. To do so, I used the following command: gatk ASEReadCounter -R /path_to_genome/hg38_genome/GRCh38.p13.genome.fa -I …
featureCounts has low successfully assigned reads
written 4.8 years ago by tomas4482 •

After finishing STAR two-step alignment, I got 62% uniquely mapped reads. But featureCounts gives me only 17% successfully aligned rate. featureCounts -T 8 -F GTF …
low alignment rate of methyl-seq data
written 5.5 years ago by Sara

I have 8 samples and paired-end `methyl-seq` data for all of them. I have a serious problem with alignment of these files. when I aligned …
getting error when indexing bam files
written 7.1 years ago by Sara

I have aligned my fastq files to hg19 using this command: bwa mem hg19.fa /DATA/myfile.fastq.gz then I made bam files from sam files using: samtools …
Problem with methylKit and Bismark
written 8.7 years ago by ahmad mousavi

Hi I have used Bismark for aligning RRBS reads and now I should run methylKit on bismarks' output, like CpG and bedGraph and .Cov. How …
Does insert size short length affect mapping rate?
written 8.7 years ago by Tania

Hi All If the insert size is small, just few bases (5-8 in average) longer than a single read length, the reads overlap with 135 …
how to visualize the RNAseq data on IGV
written 8.8 years ago by Sara

I have aligned my `RNAseq` data and 3 mismatches were allowed and sorted the `bam` files (also have .bai files). I am trying to visualize …
Bowtie2 settings in Tophat
written 12.1 years ago by zioutisch •

<p>I tried to run Tophat with different bowtie2 options. I actually used the preset bowtie2 options (i.e. --b2-sensitive). Looking for differences between different runs, i …