Bowtie mapping for single_end read
bowtie --threads 5 -X 1000 -m 1 -v 2 --best --strata --sam IndexedGenome ${DATA_DIR}${SRR_ID}.trim.fastq > ${SAM_DIR}${SRR_ID}.sam
Hi All,
I am using the above script to map single_end ChIP-Seq reads. The percentage of the aligned reads is around 20%. How can I modify the script to achieve a better alignment rate?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
CHIP-Seq data analysis
written by prs •
-
Filtering parasite reads from host reads
written by m.habib •
-
Mouse miRNA mapping
written by AellSidious •
-
Header error while converting sam file to bam
written by MS •
-
Finding the cause for corrupt ChIP-Seq data
written by Aspire
-
Modifying fastq base at specific reference location on different length reads
written by yryan •
-
bowtie (1) 0% alignment on paired-end RNAseq data
written by gulamaltab •
-
Different version of Tophat gives me different number of mapped reads
written by BioDH •
-
question about the result of MergeBamAlignment in picard
written by RC •
-
Convert Settings From Bowtie To Bowtie2
written by lsvijfhuizen •
Maybe you could use bowtie2 instead of bowtie
Thank you, Tao. We prefer bowtie to bowtie2 in this case.
I am sceptical any tweaking of parameters will make a notable difference. Parameters in bowtie are very good at default, and given it's ChIP-seq the most likely explanation is that the sample is mostly of poor quality. Make a bigwig file out of the mapped reads and inspect on a browser if you see any evidence for peaks. Make also sure that reference genome is correct.
Thank you, ATpoint!