samtools mpileup include soft-clipped bases
0
0
Entering edit mode
4 months ago
noodle ▴ 590

Hi biostars,

Is someone aware of a convenient way to include soft-clipped bases with 'samtools mpileup'?

samtools NGS soft-clip mpileup • 592 views
ADD COMMENT
1
Entering edit mode

is there more context to this question or a reason you think they are not included by default?

ADD REPLY
0
Entering edit mode

Just want to compare how many bases I lose from clipping with different parameters

ADD REPLY
1
Entering edit mode

Hard to see how it's possible since soft-clipped bases technically don't have an aligned position... Maybe it would make sense to just input alignments from your different parameter runs to get the side-by-side coverage, e.g.

samtools pileup reads.param1.bam reads.param2.bam

and then you could do some custom summary if needed.

ADD REPLY
0
Entering edit mode

Ya, I just settled on doing 'endtoend' alignments instead of 'local' so the bases aren't clipped

ADD REPLY
1
Entering edit mode

don't you mean

is it possible to include the soft-clipped BASES in "samtools mpileup" ?

ADD REPLY

Login before adding your answer.

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