Experiment: deep sequencing for mutants in 700nt fragment. the fragment of dna was preamplified by primers flanking the fragment followed by hiseq. per base coverage was calculated by coverageBed -d -abam in.bam -b ref.bed > out.cov

Observation: two distinct peaks in coverage at the ends as below plot.. coverage vs positions enter image description here the peaks are made from reads having part of primers..thus also show soft clipping at ends.. there is a huge difference in the calculations if i include such reads And if I exclude them.

Question: is there anyone who knows how to handle such a situation?


closed as unclear what you're asking by biogirl, jarlemag, Chris, Ilan, Daniel Standage Apr 15 '14 at 14:41

Please clarify your specific problem or add additional details to highlight exactly what you need. As it's currently written, it’s hard to tell exactly what you're asking. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.

  • 1
    $\begingroup$ It's absolutely not clear, how you prepared your library for sequencing, where the boundaries of your primers are and what exactly bothers you in this context. $\endgroup$ – har-wradim Apr 9 '14 at 11:08

Browse other questions tagged or ask your own question.