Bowtie mapping for single_end read
0
0
Entering edit mode
3 months ago
Rita • 0
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?

bowtie • 405 views
ADD COMMENT
1
Entering edit mode

Maybe you could use bowtie2 instead of bowtie

ADD REPLY
0
Entering edit mode

Thank you, Tao. We prefer bowtie to bowtie2 in this case.

ADD REPLY
0
Entering edit mode

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.

ADD REPLY
0
Entering edit mode

Thank you, ATpoint!

ADD REPLY

Login before adding your answer.

Traffic: 1976 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6