Comments (1)
Here is the reasoning behind the extra clips you see:
Paired-end sequences from each end of a DNA fragment. Based on that, the start position of the forward strand corresponds with one end of the DNA fragment. If the reverse strand extends beyond the forward strand's start position, it is actually going past the end of the DNA fragment being sequenced, so ClipOverlap clips that extra sequence. Similar logic applies to clipping the part of a forward strand that goes past the end of the reverse strand.
Reading past the end of the fragment results in adaptors being read, which you don't want in your SAM/BAM file, since they are not part of the DNA sequence, so we clip anything that goes beyond the "length" of the fragment as indicated by the start of the forward & reverse reads.
Does this help answer your question?
from bamutil.
Related Issues (20)
- Building bamUtil returns error: call of overloaded 'to_string(int&)' is ambiguous HOT 4
- Compile error HOT 4
- clipoverlap on raw output straight from bwa HOT 3
- Building bamUtil returns error: undefined reference HOT 4
- Cut middle of reads HOT 2
- `bam diff` should exit with a return code
- Dedup does not work from std/in HOT 2
- Cannot Compile NonPrimaryDedup Branch HOT 5
- cannot consume -.ubam on a pipe HOT 2
- Makefile should not ignore system LDFLAGS
- `bam diff` seems to return the reads that look identical between A and B HOT 1
- the clang compiler does not support -pg option on versions of OS X 10.9 and later HOT 4
- Error building from source: libStatGen HOT 3
- BamUtil: recab error: failed to open reference genome HOT 6
- bam2FastQ generates empty fastq files
- bash -x testBam2FastQ.sh and bash -x testMergeBam.sh failed on centos8_aarch64
- can not install with conda HOT 2
- bamutil diff - Get record for every read being compared HOT 2
- build is broken due to deprecated git protocol HOT 1
- Odd stats reported for clipOverlap
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from bamutil.