Comments (8)
do you use window=0 or window=500?
readCountThreshold is applied to the number of read starts in a window or exon in the control dataset. If you see a strange behavior of FREEC, can you send me an example via email please?
from freec.
I use window=0. I think that means the regions should match the BED file. With 500, regions should be 500 bp windows.
I think the point about the read starts may explain my problem. I am using the capture probes as the regions BED file. Most DNA fragments should be larger then the capture probes. Thus, most of the read starts would be outside my given regions. Should I be padding the regions to avoid this? Then I would have a lot of overlapping regions due to adjacent probes. Is that acceptable or should they be merged?
from freec.
In the previous versions of FREEC, with window=0, regions were extended by something like 150 so that read starts fall within regions surrounding exons. However, I could not find this piece of code in the latest version anymore. Apparently, it was removed by one of my collaborators when he was rewriting some functions.. :-/ I am really sorry for this. I must be back to the lab on Monday, I will put it back. Thanks a lot for pointing this out.
from freec.
So I read the code of v11.0, here is how it works now:
- in case of singe end BAM: a read is assigned to an exon if it overlap this exon (the whole read, not just its start is considered)
- in case of paired-end BAM: a read is assigned to the exon if the corresponding fragment overlaps the exon.
- in case of pileUp input file (instead of BAM), the read is considered to be 150bp long, and it is assigned to an exon, if it overlaps it.
If you run FREEC with mateOrientation=0
, all PE reads are considered as SE.
If you think that this is not what happens on your data, please share the corresponding BAM file and .cnp file and I will have a look.
from freec.
Is there a good way to share the BAMs since they are so large?
Here is an example screenshot (maybe you can see a problem from that alone):
The first two regions are present with 25x cutoff, but not with 50x cutoff. All of them have coverage above 50x or even 100x.
from freec.
Igor, and I believe that in the corresponding _control.cnp you have values below 50 for these regions, right?
Do you use mateOrientation=0
?
from freec.
The first region has coverage below 25x at the beginning, but it still passes the 25x threshold. The second region has coverage above 50x everywhere, but only passes the 25x threshold.
These are paired-end reads and mateOrientation
is set to FR
.
Is it possible that overlapping mate pairs are counted as 2 for the IGV coverage track, but counted as 1 by Control-FREEC?
from freec.
with mateOrientation=0 the reads in a pair are counted as 2. And with mateOrientation=FR each pair is counted as 1. I think I should mention this behavior in the README.
from freec.
Related Issues (20)
- no colors in makeGraph2.0.R script HOT 1
- Questions about the output of CNVs and the use of assess_significance HOT 3
- run with multiple control files HOT 2
- failed to run segmentation on chr*
- Error: to calculate BAF values HOT 2
- New release including makeGraph 2.0? HOT 1
- Make graph error HOT 1
- Problem making plots HOT 2
- Somatic cnv issue
- All zero reads with segfault HOT 1
- freec can't extract reads from bam file HOT 1
- questions about the output HOT 7
- missing p values in assess_significance.R
- Some settings for WGS explaind in manual are not included in provided config file for WGS.
- error in running gccount HOT 1
- No normalization although it is requested HOT 2
- questions about ratio.txt HOT 3
- Understanding _ratio.txt and _BAF.txt for counting informative markers HOT 1
- What is the meaning of suclonal information in the _ratio.txt file HOT 1
- get the float copy number values HOT 2
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 freec.