 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3414 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
wangxueqin金蟲 (小有名氣)
|
[求助]
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)有269人回復(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ù)
-
林科院林化所招收材料與化工專業(yè)碩士，坐標(biāo)南京
已經(jīng)有0人回復(fù)
zhongtong
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 14855.5
- 散金: 15
- 紅花: 3
- 帖子: 907
- 在線: 95.7小時(shí)
- 蟲號(hào): 267350
- 注冊(cè): 2006-07-22
- 性別: MM
- 專業(yè): 腫瘤遺傳
3樓2015-09-16 12:35:19
zhongtong
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 14855.5
- 散金: 15
- 紅花: 3
- 帖子: 907
- 在線: 95.7小時(shí)
- 蟲號(hào): 267350
- 注冊(cè): 2006-07-22
- 性別: MM
- 專業(yè): 腫瘤遺傳
2樓2015-09-16 12:30:23
radiomumm
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 8
- 帖子: 17
- 在線: 6.8小時(shí)
- 蟲號(hào): 14090163
- 注冊(cè): 2019-02-15
- 專業(yè): 作物分子育種
4樓2020-07-07 09:17:06
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +6 | 墨墨漠 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
|
[考研] 281求調(diào)劑（0805） +9 | 煙汐憶海 2026-03-16 | 19/950 |
|
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|