Spørgsmål:
Flet 2 VCF'er fra forskellige opkaldsvarianter
Lot_to_learn
2018-04-27 19:58:12 UTC
view on stackexchange narkive permalink

Jeg arbejder med WES-data til påvisning af somatiske varianter, og jeg har brugt to variantopkald, fordi ingen variantopkaldere er komplette i sig selv. Jeg har brugt GATK Haplotypecaller til små varianter som snp, ins, del osv. Og til store varianter har jeg brugt pindel (specielt til SVs-detektion). Nu har jeg to vcf-filer til hver patient, en fra gatk og en anden fra pindel. Selvom jeg kender forskellige vcf-fusionsværktøjer, men er det rigtigt at flette to forskellige vcf-filer fra forskellige opkaldere. Giver de rigtige og korrekte fusionerede vcf-filer. Er disse vcf-fusionerende værktøjer SV-opmærksomme og betragter kun trivielt opkald til den samme variant repræsenteret forskelligt i forskellige opkaldere. Jeg ved, at dette er for mange spørgsmål til en enkelt tråd, men alle disse er relateret til hinanden. Jeg er i dilemma, om jeg skal fusionere dem eller bruge dem separat.

Tak.

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

Fletteværktøjet fra SURVIVOR ser ud til at være det, du leder efter. Jeg prøvede det ikke endnu, men det ser ud til at være specielt designet til at håndtere SV.

Den enkleste ting at gøre ville sandsynligvis først flette dine patientfiler sammen for hvert værktøj og have en multipatientfil pr. variant opkald. Derefter ville jeg fusionere multi-sample-filer fra de forskellige værktøjer.

Her er et detaljeret blogindlæg, der beskriver de forskellige udfordringer forbundet med at flette strukturelle varianter fra forskellige opkaldere og forklare trinnene de tog for at adressere det. Al koden til at gengive rørledningen er tilgængelig på en github repo.

Tak for dit svar og undskyld for sent svar. Survivor er til fletning af SV'er fra forskellige opkaldsvarianter, men kan vi flette små og store varianter af vcf-fil ved hjælp af denne?
Jeg er ikke sikker, men det ser ud til, at det håndterer både små og store varianter ved at se på [kildekoden] (https://github.com/fritzsedlazeck/SURVIVOR/blob/master/src/vcfs/Merge_VCF.cpp). De har også en [wiki] (https://github.com/fritzsedlazeck/SURVIVOR/wiki/Methods-and-Parameter), der beskriver metoderne. Jeg tror, ​​du måske vil prøve en testfil og se, om den gør, hvad du forventer.
Tak. Jeg prøver at fortælle dig, om det fungerer godt med begge varianter eller ej.
Hej!! Jeg ved, det er meget sent at svare, men jeg havde ret travlt med andre ting. Jeg prøvede dette værktøj og fandt ud af, at dette værktøj udpakkede SV'er fra begge VCF-filer og flettede dem sammen i en enkelt fil. Men jeg har brug for fletningsværktøj, der kombinerer alle varianter fra begge VCF-filer inklusive SV'er og flettes sammen end til en enkelt VCF-fil, ikke kun SV'er. Desværre fungerer det ikke. Enhver anden idé ... Tak


Denne spørgsmål og svar blev automatisk oversat fra det engelske sprog.Det originale indhold er tilgængeligt på stackexchange, som vi takker for den cc by-sa 3.0-licens, den distribueres under.
Loading...