No I currently do not, to save time. I will try it with the GEM file though. I'll let you know how it goes!
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
If it looks like:
Code:hs1 1631000 1731000 3.953692 hs1 1731000 1831000 3.953692 hs1 1831000 1931000 3.953692 hs1 1931000 2031000 3.953692
Code:<plot> show = yes type = scatter file = yourFile.CNP.circos.txt glyph = rectangle glyph_size = 4 fill_color = dgrey stroke_color = dgrey stroke_thickness = 1 min = 0 max = 6 #r0 = 0.7r #r1 = 0.975r r0 = 0.76r r1 = 0.975r background = no background_color = vvlgrey background_stroke_color = black background_stroke_thickness = 1 axis = yes axis_color = lgrey axis_thickness = 1 axis_spacing = 1 <rules> <rule> importance = 100 condition = var(value) >= 6 value=6 glyph = rectangle glyph_size = 4 fill_color = red stroke_color = red stroke_thickness = 1 </rule> <rule> importance = 100 condition = var(value) > 2.5 && var(value) < 6 glyph = rectangle glyph_size = 4 fill_color = orange stroke_color = orange stroke_thickness = 1 </rule> <rule> importance = 85 condition = var(value) < 1.5 glyph_size = 4 fill_color = blue stroke_color = blue stroke_thickness = 1 </rule> <rule> show=0 importance = 85 condition = var(value) < 0 color = blue </rule> </rules> </plot>
Last edited by valeu; 07-30-2014, 07:00 AM.
Comment
-
Calling CNA and BAF for exome seq data without a control
Hi, I'm trying to use FREEC to call CNA and BAF in some exome seq samples, but haven't been able to make it work. The tool ran fine on the chr19 test data.
The problems I'm seeing are
- the tool runs but doesn't generate any of the output files except mpileup.txt_sample.cpn
It exits with the error message:
..There is no control sample!!!
You have to use a matching control sample to get adequite results since GC-bias is not the only bias in targeted sequencing
- even the .cpn file doesn't seem right because it tries to call copy number outside the target regions I provided via the TruSeq_exome_targeted_regions.hg19.bed
and not surprisingly calls them all as 0 copy number.
I'm running the tool on a pileup generated using samtools v.0.1.12. I'm using the following config file:
[general]
chrLenFile=hg19.len
window = 3000
step = 1000
ploidy = 2
#breakPointThreshold = -.001
#GCcontentProfile = GC_profile.cnp
intercept=1
minMappabilityPerWindow = 0.7
outputDir = .
sex=XY
breakPointType=4
#degree=3
#coefficientOfVariation = 0.05
#gemMappabilityFile = /hg19/out76_hg19.gem
chrFiles = FREEC_Linux64/chromosomes/
[sample]
mateFile=mpileup.txt
#mateCopyNumberFile=mpileup.txt_sample.cpn
inputFormat = pileup
mateOrientation = FR
[control]
[BAF]
SNPfile = hg19_snp131.SingleDiNucl.1based.txt
minimalCoveragePerPosition = 5
[target]
captureRegions=TruSeq_exome_targeted_regions.hg19.bed
Comment
-
Hi!
I think I should right a script to create a .WIG file from _ratio.txt (http://www.broadinstitute.org/igv/WIG)
If you code in perl, you could do it yourself and share the script. Otherwise, I will do it when I have time
Comment
-
Version 5.7 of Control-FREEC can now create BedGraph tracks on demand. Just set
Code:[general] BedGraphOutput=TRUE
Comment
-
Hi,
I have not foreseen an option to create a GC-content profile outside of FREEC. But you could make a trick: something like providing a read file with 10 reads and setting a window size.
I will think about adding a separate function to calculate this profile.
Comment
-
Hi,
so now you can generate a GC-content profile independently from FREEC:
http://bioinfo-out.curie.fr/projects...orial.html#GCP
Comment
-
a unusual problem in FREEC
I meet an unusual problem when I run FREEC. If I provide a mateCopyNumberFile, and then qsub, it will run sucessfully, but in fact I normally don't have mateCopyNumberFile, but if I don't provide it, it can't run, the error message is:
sh:samtools:command not found
Error: FREEC was not able to extract reads from /database/chenxi/task/cancerProgram/pipline/bwa/BPMICS4/sortedbam/PD1W_XXL_BPMICS4_NoIndex_modified.bam
Check your parameters: inputFormat and matesOrientation
Use "matesOrientation=0" if you have single end reads
Check the list of possible input formats at http://bioinfo-out.curie.fr/projects...al.html#CONFIG
I can't figure out what happens here. My parameters are absolutely correct, and I have samtools in my working directory and environmental PATH.
Can anyone give me some helps? I am not able to run FREEC because of this.
Comment
-
I think I know what happens. When you provide mateCopyNumberFile, FREEC does not try to read your .BAM file. When there is no mateCopyNumberFile available FREEC tries to read your file by calling "samtools" in the command line. And it seems that "samtools" does not point to anything.
I would suggest two solutions:
- in the bash script that you qsub, add something like
Code:
export PATH=$PATH:/pathToSamtools #check whether samtools exists samtools 2>samtools.usage.txt #run FREEC freec -conf myCofig.txt
- transform you .BAM file to .SAM. Then, FREEC will not need samtools.
If you prefer, you can write to me directly to [email protected]
Comment
- in the bash script that you qsub, add something like
Latest Articles
Collapse
-
by seqadmin
While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...-
Channel: Articles
06-06-2024, 07:15 AM -
-
by seqadmin
Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.
Somatic Genomics
“We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...-
Channel: Articles
05-24-2024, 01:16 PM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, Today, 02:20 PM
|
0 responses
9 views
0 likes
|
Last Post
by seqadmin
Today, 02:20 PM
|
||
Started by seqadmin, 06-07-2024, 06:58 AM
|
0 responses
181 views
0 likes
|
Last Post
by seqadmin
06-07-2024, 06:58 AM
|
||
Started by seqadmin, 06-06-2024, 08:18 AM
|
0 responses
228 views
0 likes
|
Last Post
by seqadmin
06-06-2024, 08:18 AM
|
||
Started by seqadmin, 06-06-2024, 08:04 AM
|
0 responses
184 views
0 likes
|
Last Post
by seqadmin
06-06-2024, 08:04 AM
|
Comment