Comments (1)
Dear Chris,
Thanks a lot for your questions and feedback.
In my last prep, Circle-Map assigned some very high circle scores to regions with no discordant pairs (see orange shading). Also coverage was highly variable across the region. So you might want to consider assigning more scoring weight to coverage uniformity of potential circles.
Regarding this question, I plan to improve substantially the scoring model during the next couple of months. The reason for the low priority is that we (in my lab) do not use the circle scores in our downstream analysis. However, If you (or anybody else reading this) consider that this feature is very important for the way you handle your data downstream of Circle-Map, inform me on this thread and I will try to push a bit the development of the model.
My confusion is that IGV says that the read is hard-clipped. I had the impression that Circle-Map would ignore hard-clipped reads. Is that correct? Or is it the case that Circle-Map will only looks for supplementary alignments, whether the read is hard or soft clipped?
This one is a very good question. The short answer is the following: Circle-Map uses the hard-clipped reads to construct the graph and guide the realignment of the soft-clipped reads. But it does not use hard-clipped reads for the detection.
The long answer: based on the SAM format, in some cases, BWA MEM is able to indentify split reads by itself. On those cases, BWA MEM represents the split reads (reads crossing the circular DNA breakpoint) as two separate alignments for the same read:
-
The primary alignment will be a soft-clipped read that will have an extra tag (the SA tag, you can find the information about the tag on SAM optional fields page 3) indicating the location of the alignment of the hard-clipped alignment.
-
The sencodary alignment(s): will be a hard-clipepd read with the tag indicating the location of the alignment of the soft-clipped alignment.
Hence, you can see that a single read can be represented by severall alignments on the SAM file. If Circle-Map encounters this cases, it will skip the hard-clipped alignment, but it will use the primary alignment of that read for detecting the circle.
I hope this is helpful
Best wishes,
Iñigo
from circle-map.
Related Issues (20)
- ignore unparsed_pd is empty
- Questions about Circle-Map bam2bam
- Why do the coordinates of the circle exceed the length of the contig? HOT 1
- glibc conflicts when installing
- Update Pypi Version to 1.1.5
- Could Circle-Map detect chimeric eccDNAs? HOT 1
- Circle_Map takes an unusually long time to read files. HOT 1
- ImportError: Bio.Alphabet has been removed from Biopython.
- Using bowtie2 with circle-map
- can I split alignment by chromosomes to speed up the last step
- How to plot the results in circular form?
- Circle-Map Repeats output file HOT 2
- Circle Map Realign
- some errors report about Circle-Map Repeats -i
- complicated circles
- input data HOT 1
- Circle-map Realign run too slow at the end.
- A problem of searching principle of Circle-Map HOT 3
- RuntimeWarning in Circle-Map Realign HOT 3
- Bio.Alphabet has been removed from Biopython
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 circle-map.