Dear all,
I have been looking in vain for documentation on how the "fraction of indel reads" is calculated that appears in the INFO field IS in the vcf output after samtools/bcftools (version 0.1.19).
It seems that this number is not based on the "maximum number of reads supporting the indel" divided by any of the depth values (DP, DP4, sample_DP).
Does anyone know how this fraction value is calculated?
Also, do you know, how the "maximum" number of reads is determined? Is it the highest number of indel reads found in any one sample in the analysis?
This is all it says in the vcf header about the IS field:
##INFO=<ID=IS,Number=2,Type=Float,Description="Maximum number of reads supporting an indel and fraction of indel reads">
Any comment or explication is highly appreciated.
Many thanks,
cheers,
Sophia
-
Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?
Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating...10-14-2026, 09:55 AM