| 查看: 3577 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
wangxueqin金虫 (小有名气)
|
[求助]
GATK生成的VCF转换成FASTA格式后，不是alignment结果 已有1人参与
|
|
|
有哪位高手熟悉GATK的使用么？我之前用BOWTIE+SAMTOOLS来分析过数据，生成的vcf格式文件，转换成FASTA之后直接就是alignment的结果，每个样品的序列长度是一致的。但是不知道为什么用GATK生成的vcf格式文件，转换成FASTA之后,每个样品的序列长度不一样？gatk使用的命令如下： java -jar GenomeAnalysisTK.jar -T HaplotypeCaller -R ./pe.scaffolds.db.fasta -I ./AEM.dedupped.bam -o ./AEM.raw_variants.g.vcf . --genotyping_mode DISCOVERY -stand_emit_conf 10 -stand_call_conf 30 --allow_potentially_misencoded_quality_scores --emitRefConfidence GVCF -variant_index_type LINEAR -variant_index_parameter 128000 java -jar GenomeAnalysisTK.jar \ -T GenotypeGVCFs \ -R pe.scaffolds.db.fasta \ --variant 123480.raw_variants.g.vcf \ --variant 123651.raw_variants.g.vcf \ -o 56species_raw.vcf \ java -jar GenomeAnalysisTK.jar \ -T SelectVariants \ -R pe.scaffolds.db.fasta \ -V 56species_raw.vcf \ -L PH01000000 \ -ef \ -o 56species_test.vcf |
回复此楼» 猜你喜欢
PRMT5·MTA 抑制剂 MRTX1719 的科研应用
已经有0人回复
-
SW033291 科研使用解析：分子特征、储存与实验获取要点
已经有0人回复
-
化学工程及工业化学论文润色/翻译怎么收费?
已经有83人回复
-
泛 RAS 抑制新思路：大环肽拟似物如何同时牵住突变型与野生型 RAS？
已经有0人回复
-
TYK2/JAK1 双靶点工具化合物 Brepocitinib：激酶选择性、研究进展与实验应用全解析
已经有0人回复
-
矢车菊素-3-O-葡萄糖苷质控实录：从一张HPLC图谱读懂纯度与稳定性
已经有0人回复
-
恩诺沙星（Enrofloxacin）科研解析：分子特征、LC-MS/MS 检测与实验室使用要点
已经有0人回复
-
Calderasib理化性质与实验操作技术笔记
已经有0人回复
-
RMC-6291怎么用？RAS(ON) 三重复合物工具化合物的实验室实操要点
已经有0人回复
-
苯并呋咱全解析：NBD 荧光衍生化试剂母核的化学特性与实验室应用
已经有0人回复
-
生物质炭与活化剂氢氧化钾混合注意事项
已经有0人回复
zhongtong
铁杆木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 14856.5
- 散金: 15
- 红花: 3
- 帖子: 907
- 在线: 96.3小时
- 虫号: 267350
- 注册: 2006-07-22
- 性别: MM
- 专业: 肿瘤遗传
3楼2015-09-16 12:35:19
zhongtong
铁杆木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 14856.5
- 散金: 15
- 红花: 3
- 帖子: 907
- 在线: 96.3小时
- 虫号: 267350
- 注册: 2006-07-22
- 性别: MM
- 专业: 肿瘤遗传
2楼2015-09-16 12:30:23
4楼2020-07-07 09:17:06