ANNOVAR 根据CHR_ID CHR_POS SNPS 怎么注释、
时间: 2024-06-04 22:13:51 浏览: 181
ANNOVAR可以使用以下命令根据CHR_ID、CHR_POS和SNPS对染色体位点进行注释:
annotate_variation.pl -buildver hg19 -downdb -webfrom annovar refGene humandb/
这个命令将会下载refGene数据库(如果没有下载的话),并将其与ANNOVAR脚本结合使用,以注释给定的染色体位点。其中,-buildver选项指定了参考基因组版本(这里是hg19),-downdb选项告诉ANNOVAR下载缺失的数据库,-webfrom annovar选项告诉ANNOVAR从ANNOVAR服务器上下载数据库。注释结果将存储在一个输出文件中。
相关问题
ANNOVAR 根据CHR_ID CHR_POS 怎么注释、
ANNOVAR可以根据CHR_ID和CHR_POS来进行注释。具体步骤如下:
1. 准备输入文件:将需要注释的变异信息按照以下格式保存在一个文本文件中:CHR_ID CHR_POS REF ALT。其中,CHR_ID表示染色体编号,CHR_POS表示变异位置,REF表示参考序列,ALT表示替换的序列。
2. 下载参考基因组文件和注释数据库文件:ANNOVAR需要参考基因组文件和注释数据库文件来进行注释。可以从ANNOVAR官网下载最新的参考基因组文件和注释数据库文件。
3. 运行ANNOVAR:在命令行中输入以下命令来运行ANNOVAR:
```
./annotate_variation.pl -downdb -buildver hg19 -webfrom annovar refGene humandb/
```
这个命令会下载并安装hg19版本的参考基因组文件和refGene注释数据库文件。
4. 注释变异信息:在命令行中输入以下命令来注释变异信息:
```
./annotate_variation.pl -out output_file -buildver hg19 input_file humandb/
```
其中,output_file是输出文件名,input_file是输入文件名,humandb/是注释数据库文件所在目录。
5. 查看注释结果:可以在输出文件中查看注释结果,包括变异类型、基因名称、转录本信息等。
这个代码有什么问题 #include<stdio.h> #include<string.h> int main() { int n; int i=0; char chr[200]; char chr_2[200]; char chr_3[200]; scanf("%d", &n); scanf("%[^\n]", chr); scanf("%[^\n]", chr_2); scanf("%[^\n]", chr_3); printf("%s\n", chr); printf("%s\n", chr_2); for(i=0; i<strlen(chr_3); i++){ if(chr_3[i] !=" "){ printf("%c\n", chr_3[i]); } printf("%c", chr_3[i]); } return 0; }
这段代码存在以下问题:
1. 在读取字符数组 chr_3 时,使用了错误的比较方式。应该使用单引号来表示字符,而不是双引号。因此,将 `if(chr_3[i] !=" "){` 改为 `if(chr_3[i] != ' '){`。
2. 在输出字符数组 chr_3 时,无论是否是空格都会输出一次。因此,应该将 `printf("%c", chr_3[i]);` 放到 if 语句的外面。
3. 在输入多个字符串时,应该在每个 scanf 函数前加上 getchar(),以消耗输入缓冲区中的换行符。因此,将代码修改为:
```c
#include <stdio.h>
#include <string.h>
int main() {
int n;
int i = 0;
char chr[200];
char chr_2[200];
char chr_3[200];
scanf("%d", &n);
getchar(); // 消耗第一行输入缓冲区中的换行符
scanf("%[^\n]", chr);
getchar(); // 消耗第二行输入缓冲区中的换行符
scanf("%[^\n]", chr_2);
getchar(); // 消耗第三行输入缓冲区中的换行符
scanf("%[^\n]", chr_3);
printf("%s\n", chr);
printf("%s\n", chr_2);
for(i = 0; i < strlen(chr_3); i++){
if(chr_3[i] != ' '){
printf("%c\n", chr_3[i]);
}
}
printf("%s", chr_3);
return 0;
}
```
阅读全文