diff --git a/bcftools-man.html b/bcftools-man.html index 87533badf..eb8b1cbe9 100644 --- a/bcftools-man.html +++ b/bcftools-man.html @@ -50,7 +50,7 @@

DESCRIPTION

VERSION

-

This manual page was last updated 2024-12-16 09:31 GMT and refers to bcftools git version 1.21-58-g6559a12a+.

+

This manual page was last updated 2024-12-28 22:19 GMT and refers to bcftools git version 1.21-61-g68f13f2e+.

@@ -178,10 +178,13 @@

LIST OF SCRIPTS

@@ -4380,7 +4383,8 @@

bcftools reheader [OPTIONS] file.vcf.gz

bcftools roh [OPTIONS] file.vcf.gz

A program for detecting runs of homo/autozygosity. Only bi-allelic sites -are considered.

+are considered. The results can be visualized in HTML/JavaScript using +the roh-viz script.

The HMM model:

@@ -5114,6 +5118,65 @@

plot-vcfstats [OPTIONS] file.vchk [R

+
+

roh-viz [OPTIONS]

+
+

Script for visualizing the output of bcftools roh. It creates +a HTML/JavaScript document which can be interactively viewed in your favorite +web browser.

+
+
+
+
-i, --RoH-file FILE
+
+

Output of bcftools/roh

+
+
-l, --min-length NUM
+
+

Mimimum length of ROH [1e6]

+
+
-o, --output FILE
+
+

HTML output file

+
+
-r, --regions LIST
+
+

List of chromosomes/regions

+
+
-s, --samples LIST
+
+

List of samples to show

+
+
-v, --VCF-file FILE
+
+

VCF file to determine homozygosity rate

+
+
+
+
+

-h, -?, --help + This help message

+
+
+

Example:

+
+
+
+
+
+
# Find RoH regions
+bcftools roh --AF-dflt 0.5 -G 30 -Or -o roh.txt file.bcf
+
+
+
+
+
# Create HTML/JavaScript visualization
+roh-viz -r roh.txt -v file.bcf -o output.html
+
+
+
+
+
@@ -5654,7 +5717,7 @@

COPYING

diff --git a/bcftools.html b/bcftools.html index 87533badf..eb8b1cbe9 100644 --- a/bcftools.html +++ b/bcftools.html @@ -50,7 +50,7 @@

DESCRIPTION

VERSION

-

This manual page was last updated 2024-12-16 09:31 GMT and refers to bcftools git version 1.21-58-g6559a12a+.

+

This manual page was last updated 2024-12-28 22:19 GMT and refers to bcftools git version 1.21-61-g68f13f2e+.

@@ -178,10 +178,13 @@

LIST OF SCRIPTS

@@ -4380,7 +4383,8 @@

bcftools reheader [OPTIONS] file.vcf.gz

bcftools roh [OPTIONS] file.vcf.gz

A program for detecting runs of homo/autozygosity. Only bi-allelic sites -are considered.

+are considered. The results can be visualized in HTML/JavaScript using +the roh-viz script.

The HMM model:

@@ -5114,6 +5118,65 @@

plot-vcfstats [OPTIONS] file.vchk [R

+
+

roh-viz [OPTIONS]

+
+

Script for visualizing the output of bcftools roh. It creates +a HTML/JavaScript document which can be interactively viewed in your favorite +web browser.

+
+
+
+
-i, --RoH-file FILE
+
+

Output of bcftools/roh

+
+
-l, --min-length NUM
+
+

Mimimum length of ROH [1e6]

+
+
-o, --output FILE
+
+

HTML output file

+
+
-r, --regions LIST
+
+

List of chromosomes/regions

+
+
-s, --samples LIST
+
+

List of samples to show

+
+
-v, --VCF-file FILE
+
+

VCF file to determine homozygosity rate

+
+
+
+
+

-h, -?, --help + This help message

+
+
+

Example:

+
+
+
+
+
+
# Find RoH regions
+bcftools roh --AF-dflt 0.5 -G 30 -Or -o roh.txt file.bcf
+
+
+
+
+
# Create HTML/JavaScript visualization
+roh-viz -r roh.txt -v file.bcf -o output.html
+
+
+
+
+
@@ -5654,7 +5717,7 @@

COPYING

diff --git a/howtos/plugin.split-vep.html b/howtos/plugin.split-vep.html index 37573f653..40def73e0 100644 --- a/howtos/plugin.split-vep.html +++ b/howtos/plugin.split-vep.html @@ -4,7 +4,7 @@ - + Plugin split-vep diff --git a/howtos/roh-calling.html b/howtos/roh-calling.html index f52f98880..c1846862e 100644 --- a/howtos/roh-calling.html +++ b/howtos/roh-calling.html @@ -131,7 +131,8 @@

Example 1

-
bcftools roh -G30 --AF-dflt 0.4 file.vcf
+
bcftools roh -G30 --AF-dflt 0.4 file.vcf -o out.txt
+roh-viz -i out.txt -v test.vcf.gz -o rmme.html
diff --git a/howtos/roh-calling.txt b/howtos/roh-calling.txt index 344cb8878..f8a346a8a 100644 --- a/howtos/roh-calling.txt +++ b/howtos/roh-calling.txt @@ -45,7 +45,8 @@ and there is just one sample so it cannot be estimated on the fly. Hence we use a default AF value 0.4. ---- -bcftools roh -G30 --AF-dflt 0.4 file.vcf +bcftools roh -G30 --AF-dflt 0.4 file.vcf -o out.txt +roh-viz -i out.txt -v test.vcf.gz -o rmme.html ---- Please check this link:roh-calling/usage-example.tgz[usage example] for details @@ -148,7 +149,7 @@ program. For example in this run many sites were filtered: Number of lines: total/processed: 599218/37730 ---- -If the number of the processed sites is too low, check what was the reason for excluding +If the number of the processed sites is too low, check what was the reason for excluding them. This command should give the number of sites that were processed: ---- diff --git a/howtos/roh-calling/roh-example.png b/howtos/roh-calling/roh-example.png index 6fe6440b7..0e6869ebf 100644 Binary files a/howtos/roh-calling/roh-example.png and b/howtos/roh-calling/roh-example.png differ diff --git a/howtos/variant-calling.html b/howtos/variant-calling.html index 93e4ecf92..ef822a947 100644 --- a/howtos/variant-calling.html +++ b/howtos/variant-calling.html @@ -4,7 +4,7 @@ - + Variant calling