Questions tagged [cnv]
Copy Number Variation. Questions about how to define, find, annotate, study CNVs. Algorithms.
40
questions
3
votes
0
answers
52
views
Issue creating CNV plot from WES data
I received Whole Exome Sequencing data from an NGS company (CARIS, specifically). I received R1 and R2 FASTQ files, a BAM file aligned to hg38, and a VCF file.
I used CNVPytor to create a CNV plot (...
0
votes
0
answers
14
views
Adding BAF and Log2 ratio to VCF file
I have access to pre-processed VCF files from a genetic dataset, and am trying to use the MoChA tool to calculate mosaic chromosomal alterations in this dataset. However, MoChA requires BAF and ...
2
votes
0
answers
15
views
consortium, biobank or browser with CNV eQTLs
I found that GTEx contains several SNPs eQTLs, but I'm looking for previously reported CNV-eQTLs, where can I find such data (freely available)?
- 307
2
votes
0
answers
15
views
datasets or biobanks with both (WGS or genotype intensity) and methylation data
My aim is to associate CNVs to DNAm, but I don't know where to find datasets which contain both (sequencing + methylation). Do you know any?
- 307
1
vote
2
answers
560
views
Plotting infercnv results
I'm working with matched single cell data, where we have treated and untreated samples for the same patient. I ran CNV analysis using the infercnv package.
I've ...
- 371
2
votes
1
answer
82
views
CNV Kit ` from . import commands ImportError: cannot import name 'commands' from '__main__'`
I am trying to run some code for my colleague in bash
...
- 380
-1
votes
1
answer
52
views
GATK g CNV code Aborts Before Running Main Code
I have been struggling to get the code I wrote for running GATK g CNV on my Linux machine (Steps Here: https://gatk.broadinstitute.org/hc/en-us/articles/360035531152 ) to work. Specifically it seems ...
- 380
1
vote
1
answer
27
views
what does ALT ID DUP refer to when ALT ID for copy number are already specified?
ALT ID CN2, CN3,CN4,CN5 are regions of elevated copy number relative to the reference, Then what does the ID =DUP separately cover? I am new to these analyses. please help
...
- 13
2
votes
1
answer
100
views
CNVKit does not output all the accessible regions in the targets bed file
This question was also asked on Biostars
I am using CNVkit on my data using hg38 as reference. The command that I am using is the following:
...
- 21
1
vote
0
answers
27
views
Creating msk file for Illumina Array Analysis CLI Copy number train
recently I'm researching on how to detect CNV from Illumina microarray data. I found that Illumina has Array Analysis CLI software for this task, and I need to train the CN model myself, since the ...
- 121
2
votes
0
answers
49
views
open source tools for CNV analysis from illumina SNP array
What are the best open source tools to analyse copy number variation using Illumina SNP array data...preferably R based tools, but open to other languages as well.
I'm working with leukemia sample and ...
- 21
2
votes
0
answers
173
views
Assess ploidy estimation through SNV and CNV results
I want to assess ploidy estimation (i.e., diploid vs tetraploid) in a set of samples for which I have both results from the variant caller (i.e., SNVs) and from the copy number caller.
What are the ...
- 1,783
4
votes
0
answers
46
views
How to summarize multiple exon copy numbers into copy number of the corresponding gene
I have a matrix, sample by exon, containing a copy number value for each pair (sample, exon). I would like to generate a second matrix, sample by gene, where the copy number of the exons is ...
- 1,783
1
vote
1
answer
41
views
Definition of "proband" in CNV analysis
I am reading this paper regarding de novo CNV analysis in Schizophrenia:
Accordingly, such CNVs must frequently occur as recurrent de novo mutations. In a sample of 662 schizophrenia proband–parent ...
- 1,783
1
vote
1
answer
30
views
A unified database for CNV, SNP, Indel and MSI
I am looking for a database or different databases where I can find information on different gene variants for a gene. As an example if I enter PPARG, I could be able to see SNPs, CNV, InDels and MSI.
...
- 227
0
votes
1
answer
82
views
polidy VCF file for Canvas somatic-wgs cnv calling
I am trying to run Illumina CANVAS cnv caller for Somatic-WGS. There is an option "--ploidy-vcf" which is mandatory to supply, but don't know what exactly that mean. I had supplied the CNV....
- 11
2
votes
1
answer
41
views
Cromwell: can't find file ID
I am using Cromwell to test a CWL (common workflow language) implementation of the GATK somatic CNV (copy number variation) pipeline. I am currently getting the following error (full log at the end):
...
-4
votes
1
answer
65
views
How to visualise the proportion of clonal genes to sub clonal genes [duplicate]
I have a data frame
...
- 1,981
-1
votes
2
answers
71
views
Calculating the number of probes for a given genomic range
I have a big data frame like this for copy number (exome seq)
...
- 1,981
2
votes
0
answers
52
views
Confusion about GISTIC threshold
Please, can you help me with a confusion?
In whole genome sequencing, if
Total copy number (major + minor allele copy number by SCAT R package) for SMAD4 in a ...
- 1,981
0
votes
1
answer
90
views
1
vote
1
answer
61
views
What does this statement mean?
I have a copy number segment file with this by SCAT R package
...
- 1,981
2
votes
2
answers
69
views
Detection of CNV(InDel) of intermediate size
My impression is that small InDel (a couple of bp) is identified through cigar string in BAM and typical CNV (at least thousands of bp) is detected through read depth.
What about InDel or CNV with ...
- 191
0
votes
1
answer
27
views
Arm-level SCNA does not change mRNA expression level?
It is almost common sense that local amplification of gene increase the its mRNA expression . I recently heard that large scale SCNA (such as arm-level amplification) does not increase the mRNA ...
- 191
2
votes
2
answers
47
views
Method to detect genome doubling
CNV is usually detected using either read depth or SNP VAF. However, neither of these methods can detect CNV in case of genome doubling (assuming genome doubling is a kind of genome scale CNV ...
- 191
3
votes
1
answer
61
views
Database for germline copy number variations in cancer
I am interested in looking at germline copy number variations in individuals that are at high risk of developing cancer. Are there any databases where I can look, if a CNV seen in our test case has ...
- 31
1
vote
1
answer
445
views
cn.MOPS fails with 'missing value where TRUE/FALSE needed'
I am trying to use cn.MOPS to call CNVs on a set of whole exome sequencing data bam files. My script (up to the point where it fails) is:
...
- 9,662
2
votes
1
answer
157
views
Save output to a specific folder and/or with a specific prefix in Cancer Genomics Cloud
I would like to run the CNVkit workflow on Cancer Genomics Cloud.
However, the two workflows that I'd like to use:
CNVkit Reference, for building a reference from normal samples
CNVkit CNV Calling, ...
- 1,783
1
vote
1
answer
101
views
Parallel version of cn.mops' bamDataRanges
I am trying to run cn.MOPS to detect CNVs in multiple bam files. Following the example code, I have run:
...
- 9,662
4
votes
1
answer
127
views
Simulate and test CNV workflow?
I'd like to evaluate a CNV project.
My aim is to evaluate if the scripts are sufficient for calling reasonable CNVs. I know they have a paper, but their scripts may be buggy... and all papers are ...
- 2,739
2
votes
1
answer
84
views
Comparing two files with overlapping regions and get their associated information
I have two file that have overlapping regions. One file has chromosomes, position, CNV data and samples.Other file has chromosomes, positions and ensemble ids of genes.
Here is some lines of two ...
- 105
1
vote
0
answers
101
views
Detect gene duplication events from a whole genome in FASTA format
Are there any tools for finding duplications from a whole genome in FASTA format?
Let's say the genome has gene A duplicated somewhere in the genome and gene B triplicated in tandem elsewhere in the ...
- 2,459
1
vote
0
answers
73
views
Normalize by counting CNV within a pool of VCF files
I’m (trying) using the GATK4 germline CNV calling pipeline. I successfully got 57 VCFs from my sample batch, called with segments (obtained by merging the contiguous intervals), like in a classic VCF :...
- 11
4
votes
1
answer
93
views
Compute copy number from cases and controls
I have some data on Copy Number Variation (SNP chip) for a population of samples.
In particular, I have a set of samples (considered as cases) which display a specific disease phenotype, and another ...
- 1,783
2
votes
0
answers
55
views
Methodologies for CNV analysis on SNP chip data
I'd like to analyze some Copy Number Variation data produced with SNP chip technologies.
Data is provided using contingency tables, i.e. for each gene it is provided:
...
- 1,783
2
votes
0
answers
113
views
how to plot LD between SNPs and CNVs
I have a bunch of CNVs which I have estimated LD between each allele of CNV with all SNPs across the entire human genome. I used my R-script to do a pairwise LD analysis.
The following is the first ...
- 516
2
votes
1
answer
152
views
how to count the number of each allele of CNVs?
I have a bunch of files (vcf format) each containing site and sample annotations of a copy number variant (following example is part of one of my files, I am showing 11 individuals but they are nearly ...
- 516
7
votes
2
answers
601
views
VCF merge containing CNV
How do I merge VCFs files containing CNVs?
I use vcf-merge, a VCFtools function, and after bgzip and tabix, SAMtools, to index and tab separate variants, but I don't know if it is the right way.
- 129
3
votes
1
answer
312
views
Database of copy number alterations in different cancer types
Is there a database of somatic copy number alterations (CNAs) in different cancer types? That is, annotations of which CNAs are more commonly seen in different cancer types.
- 2,324
6
votes
1
answer
59
views
Does GISTIC (v 2.0) estimate amplified/deleted probabilities on a single sample basis?
Does GISTIC 2.0 estimate the background model:
G = -log(Probability | Background)
by permuting within the sample or across all samples in the set?
The paper ...
- 183