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.

ADD COMMENTlink
0
Entering edit mode

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

ADD REPLYlink
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?

ADD REPLYlink
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?

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink
0
0
Entering edit mode

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

ADD REPLYlink
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?

ADD REPLYlink
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

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

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

ADD REPLYlink
0
Entering edit mode

GATK uses the fasta file as input, and expects the other files to be present.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1