tmhmxg59
tmhmxg59
Did you read [35](https://github.com/humanlongevity/HLA/issues/35)? I think it might help.
I believe you are missing fna, as @cjieming noted. https://github.com/humanlongevity/HLA/issues/27 https://github.com/humanlongevity/HLA/issues/35
I think [bwt_restore_sa] SA-BWT inconsistency: primary is not the same. Abort! is your main issue. I think this is because you are missing .fna and/or proper bwa index.
I get those warnings too. I'm not sure if they are significant or not. I tried resorting and removing unpaired reads (not sure if I was successful), but the warnings...
Please forgive my naïveté, but could someone suggest steps to use a GRCh37 aligned BAM? Like is it best to start with fq and realign or is there an easier...
In my hands, the provided get-reads script worked well for conversion from GRCh37 to GRCh38. https://github.com/humanlongevity/HLA/issues/35
And/or check out https://github.com/humanlongevity/HLA/issues/35