Biostar Beta. Not for public use.
i want to convert .ht2 files into FASTA file
0
Entering edit mode
16 months ago

Hey,

I want to convert .ht2 files into FASTA file. as I have obtained .ht2 files from Hisat2 tool and now I need to convert these files into fasta format please help? these .hts2 files are obtained as I have indexed my reference genome. now, I am calling haplotypecaller from GATK tools for haplotypes and it needs INDEXED reference genome as I have .hts2 index files.GATK is not accepting those files as it only accepts fasta file. please guide me?

P.S- after running hisat2, there are 8 resulting files.

0
Entering edit mode

If you still have your original reference why do you need to do this conversion?

0
Entering edit mode

because GATK haplotypecaller needs fasta file of the indexed reference genome. if I am indexing my reference file again from samtools it says "not supported file format". so I think I need to convert the indexed files which I have obtained from hisat2 into FASTA file.

if you could suggest any other way to call haplotypes from vcf files. Please suggest?

0
Entering edit mode

if I am indexing my reference file again from samtools it says "not supported file format".

samtools index ref.fa is generating that error?

0
Entering edit mode

it generates a ".fai" file but when I use that file on GATK it says not supported file format.

0
Entering edit mode
0
Entering edit mode

I have created a ".dict" file from picard tools. but then I have the error

0
Entering edit mode

Just to be sure you used the same version of GATK to make the .dict file as one you plan to use?

0
Entering edit mode

okayyy...can you please guide me any other way to find haplotypes? is there any way to obtain haplotypes from a vcf file? please provide me with the command for making ".dict file"?

can you make me clear about the gatk command? I have to input ".dict file" "fai.file" while using the command or I have to just put these files into the same directory while running the command? P.S- I am working on ubuntu

0
Entering edit mode

Those files should just be in the same directory as your fasta file

0
Entering edit mode

in command which file will be used? the original reference genome or .dict file or .fai file?