 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 3409 | 回復: 3 | ||
wangxueqin金蟲 (小有名氣)
|
[求助]
GATK生成的VCF轉(zhuǎn)換成FASTA格式后，不是alignment結(jié)果 已有1人參與
|
|
有哪位高手熟悉GATK的使用么？我之前用BOWTIE+SAMTOOLS來分析過數(shù)據(jù)，生成的vcf格式文件，轉(zhuǎn)換成FASTA之后直接就是alignment的結(jié)果，每個樣品的序列長度是一致的。但是不知道為什么用GATK生成的vcf格式文件，轉(zhuǎn)換成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 |
回復此樓» 猜你喜歡
不合理蛙科研實驗中的趣事：實驗器材的 “烏龍”
已經(jīng)有0人回復
-
不合理蛙科研實驗中的趣事：和實驗材料的 “斗智斗勇”
已經(jīng)有0人回復
-
化學工程及工業(yè)化學論文潤色/翻譯怎么收費?
已經(jīng)有214人回復
-
蛋白質(zhì)檢測：精準分析，解鎖生物分子的密碼
已經(jīng)有0人回復
-
不合理蛙科研實驗之小鼠實驗：嚴謹設(shè)計，解析生命機制的重要載體
已經(jīng)有0人回復
-
不合理蛙科研實驗之重金屬檢測：精準篩查，守護健康與環(huán)境的防線
已經(jīng)有0人回復
-
不合理蛙科研實驗之“蛙測重金屬我背鍋三千”
已經(jīng)有0人回復
-
不合理蛙科研實驗之“鼠逃三次我發(fā)三篇SCI”
已經(jīng)有0人回復
-
納米粒度分析
已經(jīng)有3人回復
-
林科院林化所招收材料與化工專業(yè)碩士，坐標南京
已經(jīng)有0人回復
zhongtong
鐵桿木蟲 (正式寫手)
- 應助: 5 (幼兒園)
- 金幣: 14855.5
- 散金: 15
- 紅花: 3
- 帖子: 907
- 在線: 95.7小時
- 蟲號: 267350
- 注冊: 2006-07-22
- 性別: MM
- 專業(yè): 腫瘤遺傳
2樓2015-09-16 12:30:23
zhongtong
鐵桿木蟲 (正式寫手)
- 應助: 5 (幼兒園)
- 金幣: 14855.5
- 散金: 15
- 紅花: 3
- 帖子: 907
- 在線: 95.7小時
- 蟲號: 267350
- 注冊: 2006-07-22
- 性別: MM
- 專業(yè): 腫瘤遺傳
3樓2015-09-16 12:35:19
4樓2020-07-07 09:17:06
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wangxueqin 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +5 | 陌の森林 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 10/500 |
|
|
|
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 9/450 |
|
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 080500，材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +3 | d如愿上岸 2026-03-13 | 3/150 |
|