troubleshooting bowtie2 low mapping rates
1
0
Entering edit mode
8.8 years ago
zizigolu ★ 4.3k

Hi,

I am mapping my illumina reads on the transcriptom by bowtie2 but the rate of alignment is too LOW(0.9). When I tried my reads on genome the alignment increased 3 times. anyway I have to use transcriptome by bowtie2 (I should not use tophat), then with which option I can truncate my reads to 50 bases?

I tried this option but

bowtie2 -x CDS -N 2 -L 28 -X 50 -U SRR1688549-ribo.fastq -S SRR1688549.sam

Thank you for your help

RNA-Seq • 3.6k views
ADD COMMENT
0
Entering edit mode

Firstly why are you mapping to the transcriptome?

ADD REPLY
0
Entering edit mode

we are working on ribosome drop-off then the R script we are using to calculating the slop needs three input file, bed file from ribo-seq, and two txt files from ribo-seq and rna-seq contains the number of reads mapped on each gene. Anyway I tried tophat but no way to get such an input there then I should used bowtie2 but because I am working with yeast I should consider interons then I should use transcriptome. Anyway the mapping rate is too vary and sometimes too low. then I used --local option in bowtie2 even the mapping rate increased three times I didn't try if the input file are allowed for my script

ADD REPLY
0
Entering edit mode

Ribo-seq. Good to know. Did you clip the adapter sequences?

ADD REPLY
0
Entering edit mode

I trimmed the adapter by universal illumina adapter sequence but the rate of trimmed reads was zero.

ADD REPLY
1
Entering edit mode

Is your library tagged with universal illumina or something else? Small RNA, CAGE and Ribo sequencing protocols often use custom adapters. You can diagnose adapter sequences with Minion (Kraken suite)

ADD REPLY
0
Entering edit mode

thank you for your comment

ADD REPLY
2
Entering edit mode
8.8 years ago
zizigolu ★ 4.3k

use --local option then the rate of mapping will increased 3 times!!!

ADD COMMENT

Login before adding your answer.

Traffic: 1241 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