1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/extensions-VCF.R
\name{breakpointRanges}
\alias{breakpointRanges}
\alias{breakpointRanges,VCF-method}
\alias{.breakpointRanges}
\title{Extracting the structural variants as a GRanges.}
\usage{
breakpointRanges(x, ...)
\S4method{breakpointRanges}{VCF}(x, ...)
.breakpointRanges(
vcf,
nominalPosition = FALSE,
placeholderName = "svrecord",
suffix = "_bp",
info_columns = NULL,
unpartneredBreakends = FALSE,
inferMissingBreakends = FALSE,
ignoreUnknownSymbolicAlleles = FALSE
)
}
\arguments{
\item{x}{A VCF object}
\item{...}{Parameters of \code{.breakpointRanges()}. See below.}
\item{vcf}{A VCF object.}
\item{nominalPosition}{Determines whether to call the variant at the
nominal VCF position, or to call the confidence interval (incorporating
any homology present). Default value is set to FALSE, where the interval is
called based on the CIPOS tag. When set to TRUE, the ranges field contains
the nominal variant position only.}
\item{placeholderName}{Variant name prefix to assign to unnamed variants.}
\item{suffix}{The suffix to append to variant names.}
\item{info_columns}{VCF INFO columns to include in the GRanges object.}
\item{unpartneredBreakends}{Determining whether to report unpartnered
breakends. Default is set to FALSE.}
\item{inferMissingBreakends}{Infer missing breakend records from ALT field of records without matching partners}
\item{ignoreUnknownSymbolicAlleles}{Ignore unknown symbolic alleles.
StructuralVariantAnnotation currently handles INS, INV, DEL, DUP as well as
the VCF specifications non-compliant RPL, TRA symbolic alleles.}
}
\value{
A GRanges object of SVs.
}
\description{
Extracting the structural variants as a GRanges.
.breakpointRanges() is an internal function for extracting structural
variants as GRanges.
}
\details{
Structural variants are converted to breakend notation.
Due to ambiguities in the VCF specifications, structural variants
with multiple alt alleles are not supported.
The CIPOS tag describes the uncertainty interval around the position
of the breakend. See Section 5.4.8 of
\url{https://samtools.github.io/hts-specs/VCFv4.3.pdf} for details of CIPOS.
If HOMLEN or HOMSEQ is defined without CIPOS, it is assumed that
the variant position is left aligned.
A breakend on the '+' strand indicates a break immediately after the given
position, to the left of which is the DNA segment involved in the breakpoint.
The '-' strand indicates a break immediately before the given position,
rightwards of which is the DNA segment involved in the breakpoint.
Unpaired variants are removed at this stage.
}
\section{Methods (by class)}{
\itemize{
\item \code{VCF}: Extracting structural variants as GRanges.
}}
\examples{
vcf.file <- system.file("extdata", "vcf4.2.example.sv.vcf",
package = "StructuralVariantAnnotation")
vcf <- VariantAnnotation::readVcf(vcf.file, "hg19")
breakpointRanges(vcf)
breakpointRanges(vcf, nominalPosition=TRUE)
}
|
