Skipping Reads Bowtie
1
1
Entering edit mode
8.2 years ago
roll ▴ 330

Hi,

I have a question about the alignment of chip-seq data.

I have got sample with 100 bp reads. After QC and trimming, these got shorter. I used bowtie to do the alignment. I got lots of warning message like

Warning: Skipping read (myData.110237305 VADER:39:C0NL4ACXX:1:2316:19783:100655 length=101) because it had length 0

What does this warning message mean? and how can i overcome it? is there an option in bowtie that will deal with this? (i check the options but could not see one!)

chip-seq bowtie • 2.9k views
ADD COMMENT
2
Entering edit mode
8.2 years ago
brentp 24k

I would guess that your trimmer trimmed some reads so that they have no bases left (length 0). It does not discard them so as to maintain pairing between paired-end reads. If yours are single-end, there is likely an option to discard those. Still, this looks like a warning so removing reads before aligning won't affect your results.

ADD COMMENT
0
Entering edit mode

Thanks. that make sense. I went through the bowtie parameter options but did not find which one would discard these. Do you know the name or which parameter i should set?

ADD REPLY

Login before adding your answer.

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