Closed nikitagambhir closed 6 years ago
The default minimum size is 32bp. This can be changed by specifying a CONFIGURATION_FILE and changing/including the following line in that file:
variantcalling.minSize = 32
I can see even 1 bp mutations in the bedpe file
GRIDSS calculates the size of an event based on to the difference between the two position of the two breakends plus the number of inserted bases. Note that this is not the same as the edit distance and an 16bp inversion will have a CIGAR of 16D16I and will be reported by the default 32bp threshold.
The 1bp events you see will either be a) insertions of least 32 bases, or b) fold-back inversions.
Thanks, I appreciate it! I just want to make sure that I follow you, so could you please comment on my explanation? In the following example, a piece of chromosome 1 got inserted into chromosome 4. So this can be classified as an inter-chromosomal SV (which might be a transposon movement).
But in the following example, one part of chromosome 1 broke and got inserted somewhere else in the same chromosome.
In this case, would it be fair to call it a simple deletion or insertion? Could it be a transposon movement? If so, can we go back and trace the sequence information?
I used GRIDSS and StructuralVariantAnnotation and I can see even 1 bp mutations in the bedpe file. If I want to change the minimum length, what would be the best route?