Question: Paired end reads merged to single end for mapping - should i allow multi-overlap similar to paired end reads?
gravatar for devikaparvathy
2.5 years ago by
devikaparvathy40 wrote:

I am using featurecounts after mapping of my reads. I have a paired-end sequencing data but Only a low percentage of reads were mapped when the paired-end reads were mapped separately and hence they were merged to a single-end read and mapped.

I know that when we use paired-end reads in featurecounts, we have to give paired-end option and also allow multi-overlap.

Does this stand true for merged single end reads too? Or should i give it as single-end read (since I have only a single bam file) and allow multi-overlaps?

The read counts are similar to ht-seq count when the multi-overlap option is set false.

ADD COMMENTlink modified 2.5 years ago • written 2.5 years ago by devikaparvathy40

multi-overlap == multi-mapped in this case?

You should not be counting reads that multi-map. It sounds like you have short inserts which are allowing your paired-end reads to largely overlap.

low percentage of reads were mapped when the paired-end reads were mapped separately

Why did you try to align the PE reads separately?

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by genomax83k

No, multimapping and multioverlapping are different. I checked this by running alternate featurecount analysis on the same dataset by changing each of the option to true and false and analysing the results. Multioverlapping when kept to false, produces many ambiguous, and the counts for lot of genes are given as zero (Same as HT-seq count). However, multimapping when kept to false produces a different result.

If multi-overlap is allowed (advantage of featurecounts over ht-seq count) then I get counts for some important genes (whose counts are zero otherwise) and these genes are actually known to play a role in the phenomenon that I am studying.

NOTE: I am working on a bacterial genome, so should not be multi-map reads counted?

I just tried and checked out initially the mapping percentage with the pared-end reads. Tried with different aligners but still the mapping percentage was low. then proceeded to merge them.

ADD REPLYlink written 2.5 years ago by devikaparvathy40
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1566 users visited in the last hour