 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 3411 | 回復(fù): 3 | |||
wangxueqin金蟲(chóng) (小有名氣)
|
[求助]
GATK生成的VCF轉(zhuǎn)換成FASTA格式后，不是alignment結(jié)果 已有1人參與
|
|
有哪位高手熟悉GATK的使用么？我之前用BOWTIE+SAMTOOLS來(lái)分析過(guò)數(shù)據(jù)，生成的vcf格式文件，轉(zhuǎn)換成FASTA之后直接就是alignment的結(jié)果，每個(gè)樣品的序列長(zhǎng)度是一致的。但是不知道為什么用GATK生成的vcf格式文件，轉(zhuǎn)換成FASTA之后,每個(gè)樣品的序列長(zhǎng)度不一樣？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 |
回復(fù)此樓» 猜你喜歡
不合理蛙科研實(shí)驗(yàn)中的趣事：實(shí)驗(yàn)器材的 “烏龍”
已經(jīng)有0人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)中的趣事：和實(shí)驗(yàn)材料的 “斗智斗勇”
已經(jīng)有0人回復(fù)
-
化學(xué)工程及工業(yè)化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有204人回復(fù)
-
蛋白質(zhì)檢測(cè)：精準(zhǔn)分析，解鎖生物分子的密碼
已經(jīng)有0人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)之小鼠實(shí)驗(yàn)：嚴(yán)謹(jǐn)設(shè)計(jì)，解析生命機(jī)制的重要載體
已經(jīng)有0人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)之重金屬檢測(cè)：精準(zhǔn)篩查，守護(hù)健康與環(huán)境的防線
已經(jīng)有0人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)之“蛙測(cè)重金屬我背鍋三千”
已經(jīng)有0人回復(fù)
-
不合理蛙科研實(shí)驗(yàn)之“鼠逃三次我發(fā)三篇SCI”
已經(jīng)有0人回復(fù)
-
納米粒度分析
已經(jīng)有3人回復(fù)
-
林科院林化所招收材料與化工專(zhuān)業(yè)碩士，坐標(biāo)南京
已經(jīng)有0人回復(fù)
zhongtong
鐵桿木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 14855.5
- 散金: 15
- 紅花: 3
- 帖子: 907
- 在線: 95.7小時(shí)
- 蟲(chóng)號(hào): 267350
- 注冊(cè): 2006-07-22
- 性別: MM
- 專(zhuān)業(yè): 腫瘤遺傳
2樓2015-09-16 12:30:23
zhongtong
鐵桿木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 14855.5
- 散金: 15
- 紅花: 3
- 帖子: 907
- 在線: 95.7小時(shí)
- 蟲(chóng)號(hào): 267350
- 注冊(cè): 2006-07-22
- 性別: MM
- 專(zhuān)業(yè): 腫瘤遺傳
3樓2015-09-16 12:35:19
radiomumm
鐵蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 8
- 帖子: 17
- 在線: 6.8小時(shí)
- 蟲(chóng)號(hào): 14090163
- 注冊(cè): 2019-02-15
- 專(zhuān)業(yè): 作物分子育種
4樓2020-07-07 09:17:06
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wangxueqin 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 271材料工程求調(diào)劑
+5
|
.6lL 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|
|
[考研] 085601材料工程專(zhuān)碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +8 | 大可digkids 2026-03-16 | 8/400 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|