In addition to Weibo, there is also WeChat
Please pay attention
WeChat public account
Shulou
2025-09-23 Update From: SLTechnology News&Howtos shulou NAV: SLTechnology News&Howtos > Internet Technology >
Share
Shulou(Shulou.com)06/01 Report--
This article is about how to use R language to use vcf format files to calculate nucleotide diversity. I think it is very practical, so I share it with you. I hope you can get something after reading this article. Let's take a look at it.
The first is to use bcftools software to manipulate the vcf file to split the vcf file according to the chromosome bcftools view snp.vcf.gz scaffold_1 > popgenome-vcf/scaffold_1
Bcftools view snp.vcf.gz scaffold_2 > popgenome-vcf/scaffold_2
If there are only vcf files in the current directory, you will encounter an error Failed to open .vcf.gz: could not load index. Please refer to https://www.cnblogs.com/chenwenyan/p/11945445.html.
Tabix-p vcf snp.vcf.gz
If there is no popgenome-vcf in the current directory, you also need to create a new directory.
Mkdir popgenome-vcf
Today's reference article says In theory, the r PopGenome can readVCF files directly, using the readVCF function. However, because our samples are haploid, we need to use a different function, r readData, which requires a folder with a separate VCF for each scaffold. Why is this?
The next step is to read in the data # install.packages ("PopGenome") in R language.
Library (PopGenome)
Getwd ()
Setwd ("VCF/")
Snp
Welcome to subscribe "Shulou Technology Information " to get latest news, interesting things and hot topics in the IT industry, and controls the hottest and latest Internet news, technology news and IT industry trends.
Views: 0
*The comments in the above article only represent the author's personal views and do not represent the views and positions of this website. If you have more insights, please feel free to contribute and share.
The market share of Chrome browser on the desktop has exceeded 70%, and users are complaining about
The world's first 2nm mobile chip: Samsung Exynos 2600 is ready for mass production.According to a r
A US federal judge has ruled that Google can keep its Chrome browser, but it will be prohibited from
Continue with the installation of the previous hadoop.First, install zookooper1. Decompress zookoope
About us Contact us Product review car news thenatureplanet
More Form oMedia: AutoTimes. Bestcoffee. SL News. Jarebook. Coffee Hunters. Sundaily. Modezone. NNB. Coffee. Game News. FrontStreet. GGAMEN
© 2024 shulou.com SLNews company. All rights reserved.