Skip to content
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

Error when path_to_bam have bam in middle #88

Open
sssimonyang opened this issue Oct 18, 2021 · 0 comments
Open

Error when path_to_bam have bam in middle #88

sssimonyang opened this issue Oct 18, 2021 · 0 comments

Comments

@sssimonyang
Copy link

Version:
bam-readcount 0.8 py37pl526h94a8ba4_4 bioconda

When path contains bam in middle, it will give index not found error.

$ bam-readcount -q 10 -b 10 -w 0 -f ${hg19} 1-bam/1-normal/Case1_normal.bam 17:7579300-7579300
Minimum mapping quality is set to 10
Expect library: Case1_NC in BAM
[bam_index_load] fail to load BAM index.
BAM indexing file is not available.

And cd 1-bam and rerun:

$ cd 1-bam/
$ bam-readcount -q 10 -b 10 -w 0 -f ${hg19} 1-normal/Case1_normal.bam 17:7579300-7579300
Minimum mapping quality is set to 10
Expect library: Case1_NC in BAM
17      7579300 A       239     =:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  A:239:60.00:32.53:0.00:134:105:0.55:0.00:2.84:134:0.50:145.08:0.43      C:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00  G:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00 T:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00   N:0:0.00:0.00:0.00:0:0:0.00:0.00:0.00:0:0.00:0.00:0.00

So the mistake comes from bam in middle in path_to_bam.

I didn't test 1.0 cause I can't find it in bioconda.

I can bypass the mistake but I think it should be fixed.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant