More posts like this
-
Get bond angles for nucleotides in pdb/cif structure file
written by noodleHi Biostars, Is someone aware of a tool that can return the bond angles across all nucleotides in a chain from a PDB file? I'm …
-
PyMOL scripting with GPUs on LSF
written by noodleHi Biostars, I'm wondering if PyMOL should be compatible with GPUs when running from a command line as `pymol -c myscript.py`, or if there is …
-
How to set soft-clipped bases option in IGV batch script
written by JellyFrog •Hi, I would like to view soft clipped bases in BAMs. I know I can click View-Preference-Alignments-Show soft clipped bases for an individual BAM, but …
-
Mapping quality score of clipped reads
written by radwa.raed •Hi, I have a question please. For soft-clipped or hard-clipped reads, for example a read with a cigar string 100M50S, does the mapping quality score …
-
Calculate Depth with and without Soft-Clipping
written by ccnn •Hi guys, I'm trying to use GATK3 (since v4 doesn't have the DepthOfCoverage tool yet) to calculate the depth of coverage across some WGS BAM …
-
When Hard clipping is used instead of the Soft clipping ?
written by 2012secondseason •Hi all. The documentation for Samtools is minimal at best. I'm still confused on the concept of a clipped read. I know that there are …
-
Remove clipped reads
written by finswimmerHello, is there a tool or script which can remove reads from my bam files which contain soft clipped bases? It would be nice if …
-
Samtools mpileup of soft-clipped alignment behaves non-deterministically?
written by iovercast •I have a sorted, indexed bam of soft-clipped reads. I'm using samtools mpileup and the pileup file it's generating is malformed. For example I see …
-
Difference between Hard Clip and Soft Clip in Samtools CIGAR string
written by QVINTVS_FABIVS_MAXIMVSHi all, The documentation for Samtools is minimal at best. I'm still confused on the concept of a clipped read. 1. What is a clipped …
-
Why Is Bwa Soft-Clipping Exactly Matching Reads?
written by ktpetrie •I am using bwasw to align reads from a MiSeq run to a very short reference sequence. For some reason, [bwa](http://bio-bwa.sourceforge.net/) is soft-clipping (assigning S …
is there more context to this question or a reason you think they are not included by default?
Just want to compare how many bases I lose from clipping with different parameters
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.
and then you could do some custom summary if needed.
Ya, I just settled on doing 'endtoend' alignments instead of 'local' so the bases aren't clipped
don't you mean