Question: Calculating read count with paired end sequencing
gravatar for sb2108
19 months ago by
sb21080 wrote:

Hey all

I have paired-end sequencing data of RNAi target fragments from a genome wide RNAi library screen. Aligning the reads to the genome has proved relatively simple, but calculating the read depth with paired ends has been more difficult as my approach has not recognized the area between 2 ends on a single read as being covered. Instead only the region corresponding to the 70 bases read directly by the sequencing are considered covered on the genome.

Is there any way to calculate read depth in which the area between 2 ends of a fragment that is not directly sequenced are added to dead depth?

Alignment was done with Bowtie2 and the read depth calculated in artemis after processing with samtools.

I can provide more details if necessary.


ADD COMMENTlink modified 19 months ago • written 19 months ago by sb21080

Try featureCounts, it is part of the subread aligner software.

ADD REPLYlink written 19 months ago by Benn8.0k

Thanks, I'll see if that works

ADD REPLYlink written 19 months ago by sb21080
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: 820 users visited in the last hour