Questions tagged [structural-variation]
Refers to large-scale genomic variation in a population or between individuals, including inversions, translocations, and copy number variants
33
questions
3
votes
0
answers
70
views
How to use hgvs to project a variant list on a given protein sequence string?
I have the following protein structural variance list:
...
- 1,437
0
votes
0
answers
194
views
How to manually curate a genome assembly for sequence variation or error?
I have a PacBio HiFi assembly of 1.1 Gb from a heterozygous species. I have aligned this assembly against a reference genome which is around 0.9 Gb. I can see that there are quite a few INDELs, ...
- 75
4
votes
1
answer
105
views
Counting hexamers in fasta sequence and identify its structure (and interruptions)
I have a lot of fasta files, each one with thousand of reads containing the hexameric motif "CCCTCT". The hexameric motif is highly continuous in most cases but interruptions may occur. I ...
- 41
1
vote
2
answers
85
views
Solutions to reference bias issue
I've recently learned about the reference bias issue - inability to properly map NGS reads in certain genomic regions caused by the fact that our standard genomic references are linear and do not ...
- 111
3
votes
2
answers
1k
views
What is the exact definition of "breakend"?
I'm reading the Manta User Guide
→ it is a method to discover structural variants and indels from next generation sequencing data.
The word 'breakend' is used in several contexts
I have searched for a ...
- 280
4
votes
1
answer
804
views
Difference between genome assembly and genome sequence alignment to a reference to find structural variants
I'm trying to determine what the difference and benefits of genome assembly and genome sequence alignments are when trying to identify structural variants or transposons in populations.
I've been ...
- 41
3
votes
0
answers
85
views
Pindel has an extremely large output
I am wondering if I am doing something wrong here.
I ran Pindel (alongside two other structural variant [SV] callers) on four WGS samples. The two other callers ...
- 31
2
votes
1
answer
62
views
Analyzing proteins based on sequence similarity
Using NCBI Gene database, I retrieved list of 7015 protein which are implicated in autism pathogenesis. This list includes tax-ID, organism name, symbol, aliases, map location, chromosome, genomic ...
2
votes
1
answer
455
views
BED file from .bam alignment structure
I am using a pipeline from nanopore to detect structural variants (SVs) in a human sample with long-reads sequencing. The first steps of the pipeline are:
index the reference genome with ...
- 283
5
votes
1
answer
2k
views
Is there a standard tool used to convert a VCF to a BEDPE?
Many popular SV callers output a VCF. Unfortunately, there isn't a unified system at the present to label events with the same notation.
However, is there a standard method for converting these VCFs ...
- 1,413
1
vote
1
answer
472
views
Quality GL entries in vcf file
I have a vcf v4.2 file with structural variants, the 8 columns describing the variant are clear, but in the genotype information I got confused by values of ...
- 5,542
-4
votes
1
answer
109
views
Finding matched genes for each genomic range
I have called structural variants and now I have
...
- 2,148
0
votes
1
answer
32
views
Generate table for total number of SV events per sample
I have called and annotated structural variants and I have a table by name of samples and structural events like below
...
- 2,148
2
votes
1
answer
107
views
Delly finds identical deletions with different consensus sequences
I am using delly to call structural variants on a human WGS sample. I noticed that it finds certain variants that are essentially duplicates. For example:
...
- 10.1k
4
votes
1
answer
3k
views
How do I find split reads?
How can I detect a split read in a BAM file?
Is there any sign in the CIGAR string that describes split read?
2
votes
0
answers
38
views
Is there a risk in inference of structural variations using mate-pairs?
Inference of structural variants from short read data is always sneaky (for a question about pair-end reads, see this). Recently I bumped into a paper where they used mate-pairs for inference of ...
- 5,542
2
votes
0
answers
47
views
structural variants: why doesn't AS + AP = SU?
I have run smoove then duphold on 46 individuals, then pasted the duphold output back into a squared vcf with all 46 samples.
I have been using SU/AC as one of my filters, looking for variants that ...
- 305
2
votes
0
answers
23
views
Visual representation of detecting deletion by different read types
I'm having some trouble differentiating between CNV, split reads and read pairs and how each one elucidates the presence of a structural variation, say deletion.
Given a reference and donor sequence,...
- 21
3
votes
1
answer
84
views
RSVSim insertions from chr1 to chr2
edit explaining python tag
I would still rather have a solution based on the RSVSim package in R, but while waiting for someone who might have an answer I wanted to look at other solutions as well.
...
- 153
5
votes
1
answer
897
views
Simplest way to work out structural variant type?
In VCF 4.2, a structural variant (SV) can be described with the BND keyword in SVTYPE. For example, the following example is an ...
- 2,699
1
vote
0
answers
66
views
Complete Genomics Junction file for SV Detection
At present, I have vcf file for each patient generated from tumor data. I am using virtual normal correction method for somatic mutation detection because I don't have normal data. For SNP or small ...
- 530
2
votes
1
answer
703
views
Can I pilon-polish long reads with Illumina short reads to improve structural variant detection?
I have pacbio Sequel data at 50 x coverage for a strain of animal. I would like to find structural variants compared to the reference genome sequence. At the moment, I align my reads to the reference ...
- 2,459
1
vote
0
answers
179
views
Can Mauve export indels or is there an whole-genome aligner that can export indels?
I have aligned two whole genomes of the same species in Mauve using progressiveMauve. One of the genomes is 2 megabases longer than the other so I know it must contain some insertions. Looking at the ...
- 2,459
1
vote
2
answers
174
views
Can structural - equation modeling be used in Bioinformatics?
I am a bioinformatics student and I was learning the basics of structural equation modeling (SEM).
I was wondering, are there any applications of SEM on bioinformatics? Are there any strong sides ...
- 39
3
votes
1
answer
561
views
Why are inversions defined as the reverse complement and not just the reverse of the reference?
I can’t quite understand the way inversions are defined. In particular I expect an inversion to be only a reversed version of the reference and not its reverse complement. Most sources use diagrams ...
- 309
4
votes
2
answers
272
views
After artificially creating events in a FASTA file, how do I keep track of the old coordinates?
I'm beginning with the reference genome in FASTA format, hg19. I am reading the sequence into a Python dictionary with BioPython:
...
- 1,691
6
votes
1
answer
225
views
Detecting structural variants with MinION data
Working on various cancers I have an interest in detecting structural variation (SV) in human, we've successfully used various tools like Pindel, SVDetect, Manta, and LUMPY, to name a few for ...
- 1,069
6
votes
2
answers
162
views
Structural variant calling for low-coverage PacBio data
PacBio is selling ~10x PacBio SEQUEL long reads as an upgrade to Illumina data for SV discovery.
In a clinical setting, the main requirements are proper sensitivity and specificity but also the ...
- 588
13
votes
2
answers
403
views
Remapping genomic coordinates to account for indels
I'm interested in obtaining coding sequences of my favourite gene in all individuals from the 1000Genomes (and similar projects). I use GATK to get the right subset of variants, vcf-consensus to map ...
- 232
13
votes
2
answers
2k
views
How to read structural variant VCF?
The IGSR has a sample for encoding structural variants in the VCF 4.0 format.
An example from the site (the first record):
...
- 2,699
17
votes
5
answers
10k
views
Is there an easy way to create a summary of a VCF file (v4.1) with structural variations?
I got a bunch of vcf files (v4.1) with structural variations of bunch of non-model organisms (i.e. there are no known variants). I found there are quite a some tools to manipulate vcf files like ...
- 5,542
13
votes
3
answers
371
views
Given a VCF of a human genome, how do I assess the quality against known SNVs?
I'm looking for tools to check the quality of a VCF I have of a human genome. I would like to check the VCF against publicly known variants across other human genomes, e.g. how many SNPs are already ...
- 1,691
15
votes
2
answers
439
views
How can I call structural variants (SVs) from pair-end short read resequencing data?
I have a reference genome and now I would like to call structural variants from Illumina pair-end whole genome resequencing data (insert size 700bp).
There are many tools for SV calls (I made an ...
- 5,542