Why Bowtie2 and Lastz give so different results?
0
0
Entering edit mode
8.4 years ago
gbdias ▴ 150

Hey guys,

I'm mapping 50 bp reads from several ChIP-seq datasets to a short reference (170bp). I want to measure the relative enrichment of this specific consensus in the different IP conditions. It's repetitive DNA so I guess it's ok to use a consensus.

To do so, I used Bowtie2 (with the --end-to-end and --local modes) and Lastz (with the 12of19 seed and no filtering).

The number of reads mapping to my reference in all datasets (about 30 ngs files) was always bigger when I used Lastz. I can understand why bt2 end-to-end VS local give different results. But shouldn't Lastz and bt2 --local give comparable results?

Should I apply filtering to the Lastz mapping?

bowtie2 lastz alignment ChIP-Seq • 1.7k views
ADD COMMENT

Login before adding your answer.

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