Dear all,
A strange idea has came to my mind recently.
There is no direct way to convert vcf into fasta. But
I've read that through Galaxy I can convert vcf to maf (multiple align),
and then it may be possible to convert maf to fasta.
How much information from the original vcf-file will I loose using this way,
and is it possible to avoid the loss?
THANK YOU!
Natasha
But this option killed me.
You definitely know. PLEASE, help me!
Thank you very much indeed!
Natasha
is -L required ?
Dear Pierre,
this is optional parameter. But when I omit it, I have just the full fasta-fail for corresponding chromosome, nothing else. I thought that GATK will allow me to cut a fragment corresponding.to vcf-file. But I have to give it "one or more genomic intervals over which to operate". It seems much more complicated. How to find these genomic intervals?
I have vcf-files, refs and chromosome sequences. I made fai-files, dict-files, but it still doesn't give me any hints to the fasta alignment. What else should be done?
Many thanks for your help!
Natasha
"one or more genomic intervals over which to operate": I don't understand your problem. You can provide a BED file or even a VCF file: http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_CommandLineGATK.html#--intervals
Great. And this will work??? " In order to perform the analysis
at specific positions based on the records present in the file
(e.g. -L file.vcf)" - this is exactly what I need. I will try, that's a miracle.
I was very poor in reading the manual... THANK YOU, Pierre!
I guess I missed the goal of OPs question, I just was not sure what they wanted in the fasta file. I presumed that they didn't have a preexisting fasta file with their reference sequences.