Fråga:
Slå samman 2 VCF: er från olika varianter
Lot_to_learn
2018-04-27 19:58:12 UTC
view on stackexchange narkive permalink

Jag arbetar med WES-data för detektering av somatiska varianter och jag har använt två variantuppringare eftersom ingen variantuppringare är komplett i sig. Jag har använt GATK Haplotypecaller för små varianter som snp, ins, del etc. Och för stora varianter har jag använt pindel (speciellt för SV-detektering). Nu har jag två vcf-filer för varje patient, en från gatk och en annan från pindel. Även om jag känner till olika vcf-sammanslagningsverktyg men är det rätt att slå samman två olika vcf-filer från olika uppringare. Ger de rätt och korrekt sammanslagen vcf-fil. Är dessa vcf-sammanslagningsverktyg SV-medvetna och anser endast trivialt samtal för samma variant representerat annorlunda i olika uppringare. Jag vet att det här är för många frågor för en enda tråd men alla dessa är relaterade till varandra. Jag är i dilemma om jag ska slå samman dem eller använda dem separat.

Tack.

Ett svar:
cmdoret
2018-04-27 20:18:52 UTC
view on stackexchange narkive permalink

Sammanfogningsverktyget från SURVIVOR verkar vara det du letar efter. Jag försökte inte det ännu, men det verkar vara speciellt utformat för att hantera SV.

Det enklaste att göra skulle förmodligen först slå samman dina patientfiler för varje verktyg och ha en multipatientfil per variantuppringare. Då skulle jag slå samman flera exempelfiler för de olika verktygen.

Här är ett detaljerat blogginlägg som beskriver de olika utmaningarna som är förknippade med att slå samman strukturvarianter från olika uppringare och förklara stegen de tog för att ta itu med det. All kod för att reproducera rörledningen finns på en github repo.

Tack för ditt svar och ledsen för sent svar. Survivor är för att slå samman SV: er från olika varianter, men kan vi slå samman små och stora varianter av vcf-fil med den här?
Jag är inte säker, men det verkar som att det hanterar både små och stora varianter genom att titta på [källkoden] (https://github.com/fritzsedlazeck/SURVIVOR/blob/master/src/vcfs/Merge_VCF.cpp). De har också en [wiki] (https://github.com/fritzsedlazeck/SURVIVOR/wiki/Methods-and-Parameter) som beskriver metoderna. Jag tror att du kanske vill prova en testfil och se om den gör vad du förväntar dig.
Tack. Jag ska försöka låta dig veta om det fungerar bra med båda varianterna eller inte.
Hallå!! Jag vet att det är väldigt sent att svara men jag var ganska upptagen med andra saker. Jag försökte det här verktyget och upptäckte att det här verktyget extraherar SV: er från båda VCF-filerna och slår ihop dem i en enda fil. Men jag behöver ett sammanslagningsverktyg som kombinerar alla varianter från båda VCF-filerna inklusive SV-filer och slås samman än i en enda VCF-fil, inte bara SV-filer. Tyvärr fungerar det inte. Någon annan idé ... Tack


Denna fråga och svar översattes automatiskt från det engelska språket.Det ursprungliga innehållet finns tillgängligt på stackexchange, vilket vi tackar för cc by-sa 3.0-licensen som det distribueras under.
Loading...