-
Notifications
You must be signed in to change notification settings - Fork 52
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
run get-reads-alt-unmap.sh in container #35
Comments
Hi, xHLA only works with hg38 reference. The get-reads-alt-unmap.sh script needs to be used when you have a BAM aligned on hg38 reference with alt contigs. Best, |
Conversion from 37 to 38 worked for me for one bam. I used: I did have to remove the memory detection part of the script as it gave an error at line 17. On a second data set, though, I'm loosing about half of the reads somewhere in the script. My first go I added a few reads as expected. The only thing I can think of is: when I ran the first time I was using a 37 bam PrintRead to chr6. The first step in the script is to trim, though, so I'm not sure what the issue is. |
Second data set had chromosome names as number only, so modified get-reads-alt-unmap.sh to 6 instead of chr6. |
I'm struggling to get typing information from hg19 aligned bams. I think I'd like to run get-reads-alt-unmap.sh in the docker container - I don't have all the dependencies on my machine.
I think I need to change the entrypoint, but I'm not sure the syntax. The coordinates in the .sh script appear to cover all genes in both references.
I tried several variations of:
$docker run -itv
pwd
:pwd
pwd
-w humanlongevity/hla --entrypoint /bin/bash lsdocker: invalid reference format: repository name must be lowercase.
See 'docker run --help'.
But doesn't seem like I've successfully changed the entrypoint.
The text was updated successfully, but these errors were encountered: