Description Usage Arguments Details Author(s)
Estimate Insert size library from pairs of sequence mapping on the same contig.
Estimate Insert size library from pairs of sequence mapping on two different contigs.
Estimate the assembly quality
Count the number of pairs of sequence that connect two contigs.
Plot method for the class ContigLinkViewer
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ## S4 method for signature 'ContigLinkViewer'
EstimateInsertSizeFromSame(object, ...)
## S4 method for signature 'ContigLinkViewer'
EstimateInsertSizeFromDiff(object, ...)
## S4 method for signature 'ContigLinkViewer'
EstimateAssQual(object, EstimateCutOff = TRUE,
Min = NA, Max = NA, SdNumber = NA)
## S4 method for signature 'ContigLinkViewer'
ContigLinkCount(object, MinCount)
## S4 method for signature 'ContigLinkViewer'
plot(object, x = NA, ...)
|
object |
An object of class
|
An |
object of class |
object |
An object of class
|
EstimateCutOff |
If TRUE, the interval for the insert size is [MeanIsize-(SdNumber*SdIsize), MeanIsize+(SdNumber*SdIsize)] |
Min |
If EstimateCutOff is FALSE The min value for the insert Size |
Max |
If EstimateCutOff is FALSE The max value for the insert Size |
SdNumber |
The number of sd to keep for computing the insert size interval. |
An |
object of class |
MinCount |
The minimum number of pairs connecting two contigs |
object |
An object of class
|
... |
Arguments to be passed to methods, such as
graphical parameters (see
|
x |
the coordinates of points in the plot.
Alternatively, a single plotting structure, function or
any R object with a |
Estimate Insert size library from pairs of sequence mapping on the same contig and having the good orientation. The well mapped paired-end sequence (-> <-) belong to this two categorism:
1 The strand of the first sequence in pair is + and its sam flag is 99. The strand of the second sequence in the pair is - and its flag is 147. The insert size is obtained by substracting the start position of the first sequence on the contig to the end position of the second sequence on the contig and by adding 1.
2 The strand of the first sequence in pair is - and its sam flag is 83. The strand of the second sequence in the pair is + and its flag is 163. The insert size is obtained by substracting the start position of the second sequence on the contig to the end position of the first sequence on the contig and by adding 1.
Estimate Insert size library from pairs of sequence mapping on two contigs. In this case the orientation of the sequence on the contig does not matter. To obtain the insert size, two cases:
The sequence is the first in the pair The insert size is the length of the first contig plus the end position of the second sequence on the second contig minus the start position of the first sequence on the first contig minus 1.
The sequence is the second in the pair The insert size is the length of the first contig plus the end position of the first sequence on the first contig minus the start position of the second sequence on the second contig minus 1.
This method creates the AssQual slot which contain 4 objects:
parameters A vector giving the insert size parameters which will be used to estimate the assembly quality. Min and Max or Mean and sd.
dfSameQualAssByContig A data.frame summarizing the repartition of the pairs according to the insert size range and the orientation of the pairs. This data.frame has 6 columns:
seqnames.1 A vector giving the name of the contig
strand.1 A vector giving the mapping strand of the first sequence in pair
strand.2 A vector giving the mapping strand of the second sequence in pair
outof A vector giving the number of pair that are outside the confidence interval for the insert size
nonav A vector giving the number of pair for which the insert size can not be computed
within A vector giving the number of pair with an insert size inside the confifence interval
dfSameQualAssTot A data.frame which sum the dfSameQualAssByContig by orientation only disregarding contig information.
dfDiffQualAssByContig A data.frame summarizing the repartition of the pairs according to the insert size range and the orientation of the contigs that are linked. This data.frame has 6 columns: for the insert size
seqnames.1 A vector giving the name of the first contig
seqnames.2 A vector giving the name of the second contig
Ori1 A vector giving the orientation of the first contig
Ori2 A vector giving the orientation of the second contig
outof A vector giving the number of pair that are outside the confidence interval for the insert size
within A vector giving the number of pair with an insert size inside the confifence interval
To get the quality, a confident interval has to be computed for the insert size.
This method creates the CLC slot which is a list of 2 objects:
MinCount A vector giving the minimum number of pairs binding two contigs
CLC A data.frame with 5 columns giving the number of pairs above the MinCount parameter that support the connection beween two contigs.
seqnames.1 A vector giving the name of the first contig
seqnames.2 A vector giving the name of the second contig
Ori1 A vector giving the orientation of the first contig
Ori2 A vector giving the orientation of the second contig
within A vector giving the number of pairs above the MinCount parameter that support the connection beween the two contigs.
This plot function allow to generate a graph to vizualize the contigs and their link using the Rgraphviz package. Contigs are the nodes and the number of pairs supporting the link between two are specified by the edge. The ContigLinkViewer object has to be complete.
Delphine Charif & Joseph Tran
Delphine Charif & Joseph Tran
Delphine Charif & Joseph Tran
Delphine Charif & Joseph Tran
Delphine Charif & Joseph Tran
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.