| 5 | 1/1 | 返回列表 |
| 查看: 3410 | 回復: 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 |
鐵蟲 (初入文壇)

鐵桿木蟲 (正式寫手)
鐵桿木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 化學工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 0703化學調(diào)劑 ,六級已過,有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 0703化學調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學666 2026-03-16 | 6/300 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|