DiltheyLab / HLA-LA

Fast HLA type inference from whole-genome data
GNU General Public License v3.0
120 stars 40 forks source link

if I have fastq paired read data, what should i do? #33

Closed PavitaKae closed 4 years ago

PavitaKae commented 4 years ago

I noticed that your code extract bam to fastq first, and mapping with your reference. if I have fastq already, Should i skip extract bam file and use my fastq as input to your program? How?

AlexanderDilthey commented 4 years ago

There is an un-tested code branch in the C++ component of HLA*LA that accepts FASTQ files as input, but we currently don't support this option. In any case there are not going to be any great performance advantages, as, even if FASTQ is supplied, the first step will always consist of a complete mapping of all reads.