plotting roh from bcftools
0
0
Entering edit mode
3.3 years ago
gubrins ▴ 350

Heys,

I am following this small tutorial on how to calculate ROHs from a vcf file using bcftools (https://samtools.github.io/bcftools/howtos/roh-calling.html) and I am getting this txt file:

# This file was produced by: bcftools roh(1.10.2+htslib-1.10.2-3)
# The command line was: bcftools roh -G30 --AF-dflt 0.4 my_file.vcf
#
# RG    [2]Sample   [3]Chromosome   [4]Start    [5]End  [6]Length (bp)  [7]Number of markers    [8]Quality (average fwd-bwd phred score)
# ST    [2]Sample   [3]Chromosome   [4]Position [5]State (0:HW, 1:AZ)   [6]Quality (fwd-bwd phred score)
ST  Leopards_1  NC_018723.3 280 0   3.0
ST  Leopards_1  NC_018723.3 281 0   27.7
ST  Leopards_1  NC_018723.3 292 0   27.7
ST  Leopards_1  NC_018723.3 293 0   27.7
ST  Leopards_1  NC_018723.3 303 0   27.7

I am wondering how to interpret the results and how to plot them. From the plot they show I don't understand which are the ROHs areas and if there is a difference between ST and RG. Any type of help would be more than appreaciated!

bcftools roh plotting • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 1722 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6