I want to do some HLA typing, most of the tools require the bam file is aligned to primary genome without alternative read handling. From the International Genome Sample Resource project, I can get the high coverage file, but they are all aligned to primary genome with hla, decoy and alternative reads. I want to make the alternative reads into unmapped reads.

I think about unmapping the whole file to fastq, and then realigning according to different assembly, but it's time-consuming and also costs a lot of storage, so I don't want to follow that approach.

I also tried this script https://github.com/humanlongevity/HLA/blob/master/bin/get-reads-alt-unmap.sh, but so far it's not working (I have all requirement packages installed).

Any help, tool or insight for this problem is warmly welcome.

Thank you

I want to do some HLA typing, most of the tools require the bam file is aligned to primary genome without alternative read handling. From the International Genome Sample Resource project, I can get the high coverage file, but they are all aligned to primary genome with hla, decoy and alternative reads. I want to make the alternative reads into unmapped reads.

I think about unmapping the whole file to fastq, and then realigning according to different assembly, but it's time-consuming and also costs a lot of storage, so I don't want to follow that approach.

I also tried this script, but so far it's not working (I have all requirement packages installed).

Any help, tool or insight for this problem is warmly welcome.

Thank you

I want to do some HLA typing, most of the tools require the bam file is aligned to primary genome without alternative read handling. From the International Genome Sample Resource project, I can get the high coverage file, but they are all aligned to primary genome with hla, decoy and alternative reads.

I think about unmapping the whole file to fastq, and then realigning according to different assembly, but it's time-consuming and also costs a lot of storage, so I don't want to follow that approach.

I also tried this script https://github.com/humanlongevity/HLA/blob/master/bin/get-reads-alt-unmap.sh, but so far it's not working (I have all requirement packages installed).

Any help, tool or insight for this problem is warmly welcome.

Thank you

I want to do some HLA typing, most of the tools require the bam file is aligned to primary genome without alternative read handling. From the International Genome Sample Resource project, I can get the high coverage file, but they are all aligned to primary genome with hla, decoy and alternative reads. I want to make the alternative reads into unmapped reads.

I think about unmapping the whole file to fastq, and then realigning according to different assembly, but it's time-consuming and also costs a lot of storage, so I don't want to follow that approach.

I also tried this script https://github.com/humanlongevity/HLA/blob/master/bin/get-reads-alt-unmap.sh, but so far it's not working (I have all requirement packages installed).

Any help, tool or insight for this problem is warmly welcome.

Thank you