10XGenomics / scHLAcount

Count HLA alleles in single-cell RNA-seq data
MIT License
58 stars 12 forks source link

Panic occurred in file 'src/mapping.rs' at line 52 #13

Open danshu opened 3 years ago

danshu commented 3 years ago

Hello,

I'm trying to run scHLAcount but having this error.

Sharding sequences... Got 373889 sequence chunks Assembling 114 shards...

Done separate de Bruijn graph construction Starting merging disjoint graphs Merger of graphs complete Indexing de Bruijn graph Total 58019 kmers to process in dbg Making mphf of kmers Assigning offsets to kmers Well, this is embarrassing.

sc_hla_count had a problem and crashed. To help us diagnose the problem you can send us a crash report.

name = 'sc_hla_count' operating_system = 'unix:CentOS' crate_version = '0.1.0' explanation = ''' Panic occurred in file 'src/mapping.rs' at line 52 ''' method = 'Panic' backtrace = ''' stack backtrace: 0: backtrace::backtrace::libunwind::trace::h3b4dbc2ad001faa5 (0x555cc902bdb7) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/libunwind.rs:88 backtrace::backtrace::trace_unsynchronized::h77d6568b33e3258f at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/mod.rs:66 backtrace::backtrace::trace::h47c6435d0453bfe0 at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/mod.rs:53 1: backtrace::capture::Backtrace::create::ha7dcae36bc2f2bde (0x555cc902b315) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/capture.rs:163 backtrace::capture::Backtrace::new::h92bc83138b1926c0 at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/capture.rs:127 2: human_panic::report::Report::new::h1cbf997f6d47b7d3 (0x555cc901b749) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/report.rs:45 3: human_panic::handle_dump::h1add66df506345a8 (0x555cc901af13) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/lib.rs:197 4: sc_hla_count::main::{{closure}}::h4b47b3fa03a044b7 (0x555cc8d8b729) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/lib.rs:120 5: std::panicking::rust_panic_with_hook::he82f5d0644692441 (0x555cc9141668) at src/libstd/panicking.rs:530 6: rust_begin_unwind (0x555cc914123b) at src/libstd/panicking.rs:437 7: core::panicking::panic_fmt::h09c929f06bb87c98 (0x555cc9166891) at src/libcore/panicking.rs:85 8: core::panicking::panic::h7ece43057e5422d4 (0x555cc91667dd) at src/libcore/panicking.rs:50 9: core::option::Option::unwrap::h2b3d96d5cd337084 (0x555cc8d65298) at /rustc/04488afe34512aa4c33566eb16d8c912a3ae04f9/src/libcore/macros/mod.rs:10 sc_hla_count::mapping::BamSeqReader::fetch::h49bde93d68dcce94 at src/mapping.rs:52 sc_hla_count::mapping::mapping_wrapper::hebc162172d078d87 at src/mapping.rs:317 10: sc_hla_count::_main::h3253537751c2ae0e (0x555cc8d915e6) at src/main.rs:198 11: sc_hla_count::main::h6c74ee357aa88633 (0x555cc8d8b4cf) at src/main.rs:139 12: std::rt::lang_start::{{closure}}::h4a39b9fb10bdccb7 (0x555cc8d76bf3) at /rustc/04488afe34512aa4c33566eb16d8c912a3ae04f9/src/libstd/rt.rs:67 13: std::rt::lang_start_internal::{{closure}}::h5d3ea623498f5f43 (0x555cc9141ae3) at src/libstd/rt.rs:52 std::panicking::try::do_call::hac65e71be769a440 at src/libstd/panicking.rs:348 std::panicking::try::hd4706e264bcf6712 at src/libstd/panicking.rs:325 std::panic::catch_unwind::h948a0fb4a8b3ee82 at src/libstd/panic.rs:394 std::rt::lang_start_internal::h72cc068ed2d0ac53 at src/libstd/rt.rs:51 14: main (0x555cc8d93788) 15: __libc_start_main (0x7f4f89c2e505) 16: (0x555cc8d46dcf) 17: (0x0)'''

danshu commented 3 years ago

Here is the another report when running using the linux binary sc_hla_count_linux.

name = 'sc_hla_count' operating_system = 'unix:CentOS' crate_version = '0.1.0' explanation = ''' Panic occurred in file 'src/libcore/option.rs' at line 347 ''' method = 'Panic' backtrace = ''' stack backtrace: 0: backtrace::backtrace::libunwind::trace::hb73dfa724f78bc0d (0x5558f6e8fa07) at /mnt/home/patrick/.cargo/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/libunwind.rs:88 backtrace::backtrace::trace_unsynchronized::hd68d5720b8eeecb4 at /mnt/home/patrick/.cargo/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/mod.rs:66 backtrace::backtrace::trace::h07a8a1b16091cf6f at /mnt/home/patrick/.cargo/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/mod.rs:53 1: backtrace::capture::Backtrace::create::ha2012ed583b4dee0 (0x5558f6e8ef63) at /mnt/home/patrick/.cargo/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/capture.rs:163 backtrace::capture::Backtrace::new::h2402c178e424aa5f at /mnt/home/patrick/.cargo/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/capture.rs:127 2: human_panic::report::Report::new::h5e9fc17a20b5f651 (0x5558f6e7def7) at /mnt/home/patrick/.cargo/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/report.rs:45 3: human_panic::handle_dump::h6497bba18125ae8e (0x5558f6e7d5e7) at /mnt/home/patrick/.cargo/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/lib.rs:197 4: sc_hla_count::main::{{closure}}::h4d1b937fa8891909 (0x5558f6dcf129) at /mnt/home/pat/code/scHLAcount/<::human_panic::setup_panic macros>:18 5: std::panicking::rust_panic_with_hook::hffcefc09751839d1 (0x5558f6fb5da9) at src/libstd/panicking.rs:481 6: std::panicking::continue_panic_fmt::hc0f142c930c846fc (0x5558f6fb5842) at src/libstd/panicking.rs:384 7: rust_begin_unwind (0x5558f6fb5726) at src/libstd/panicking.rs:311 8: core::panicking::panic_fmt::h2daf88b2616ca2b2 (0x5558f6fd4a3d) at src/libcore/panicking.rs:85 9: core::panicking::panic::h2d0bc53a963fb996 (0x5558f6fd497c) at src/libcore/panicking.rs:49 10: core::option::Option::unwrap::hbcc241c67317b537 (0x5558f6d9d428) at /rustc/eae3437dfe991621e8afdc82734f4a172d7ddf9b/src/libcore/macros.rs:12 sc_hla_count::mapping::BamSeqReader::fetch::hf81cf3042ac7d942 at src/mapping.rs:52 sc_hla_count::mapping::mapping_wrapper::hd4dba6eb8a8c6fa1 at src/mapping.rs:320 11: sc_hla_count::_main::h60c6e0494337c80f (0x5558f6dd4dc3) at src/main.rs:198 12: sc_hla_count::main::he49cdf9c5df2b57b (0x5558f6dcef83) at src/main.rs:139 13: std::rt::lang_start::{{closure}}::h93883743bb8c8181 (0x5558f6daace3) at /rustc/eae3437dfe991621e8afdc82734f4a172d7ddf9b/src/libstd/rt.rs:64 14: std::rt::lang_start_internal::{{closure}}::h447d8812e3ee306d (0x5558f6fb56c3) at src/libstd/rt.rs:49 std::panicking::try::do_call::h4a61cb372364c745 at src/libstd/panicking.rs:296 15: __rust_maybe_catch_panic (0x5558f6fbcc4a) at src/libpanic_unwind/lib.rs:82 16: std::panicking::try::hdf71f938885bca42 (0x5558f6fb628d) at src/libstd/panicking.rs:275 std::panic::catch_unwind::h7e85dbf162b1611a at src/libstd/panic.rs:394 std::rt::lang_start_internal::h1e06cc26b9fc25ea at src/libstd/rt.rs:48 17: main (0x5558f6dd6788) 18: __libc_start_main (0x7f4f5c973505) /tmp/report-52281d61-da7b-4b9d-9d18-5dd0cc567aa4.toml

danshu commented 3 years ago

It seems that the error is related to "BamSeqReader". My bam file is actually produced using the most recent CR4.0, is this causing the error?

danshu commented 3 years ago

I may have figured out the cause of this error, which is the chromosomes names' update in the references (from 1 and MT to chr1 and chrM, https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/release-notes). I can successfully run scHLAcount on the test bam file, but when I have changed the chromosome names in the test bam file: samtools view -H test.bam | sed -e 's/SN:([0-9XY])/SN:chr\1/' -e 's/SN:MT/SN:chrM/' | samtools reheader - test.bam > test_chr.bam I can exactly reproduce the same error message: name = 'sc_hla_count' operating_system = 'unix:CentOS' crate_version = '0.1.0' explanation = ''' Panic occurred in file 'src/mapping.rs' at line 52 ''' method = 'Panic' backtrace = ''' stack backtrace: 0: backtrace::backtrace::libunwind::trace::h3b4dbc2ad001faa5 (0x55c46daf4db7) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/libunwind.rs:88 backtrace::backtrace::trace_unsynchronized::h77d6568b33e3258f at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/mod.rs:66 backtrace::backtrace::trace::h47c6435d0453bfe0 at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/backtrace/mod.rs:53 1: backtrace::capture::Backtrace::create::ha7dcae36bc2f2bde (0x55c46daf4315) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/capture.rs:163 backtrace::capture::Backtrace::new::h92bc83138b1926c0 at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/backtrace-0.3.35/src/capture.rs:127 2: human_panic::report::Report::new::h1cbf997f6d47b7d3 (0x55c46dae4749) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/report.rs:45 3: human_panic::handle_dump::h1add66df506345a8 (0x55c46dae3f13) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/lib.rs:197 4: sc_hla_count::main::{{closure}}::h4b47b3fa03a044b7 (0x55c46d854729) at /research/kevinyip12/jqcao/tools/HLA/scHLAcount/registry/src/github.com-1ecc6299db9ec823/human-panic-1.0.1/src/lib.rs:120 5: std::panicking::rust_panic_with_hook::he82f5d0644692441 (0x55c46dc0a668) at src/libstd/panicking.rs:530 6: rust_begin_unwind (0x55c46dc0a23b) at src/libstd/panicking.rs:437 7: core::panicking::panic_fmt::h09c929f06bb87c98 (0x55c46dc2f891) at src/libcore/panicking.rs:85 8: core::panicking::panic::h7ece43057e5422d4 (0x55c46dc2f7dd) at src/libcore/panicking.rs:50 9: core::option::Option::unwrap::h2b3d96d5cd337084 (0x55c46d8387ae) at /rustc/04488afe34512aa4c33566eb16d8c912a3ae04f9/src/libcore/macros/mod.rs:10 sc_hla_count::mapping::BamSeqReader::fetch::h49bde93d68dcce94 at src/mapping.rs:52 sc_hla_count::mapping::map_and_count_pseudo::h0307c0db42007237 at src/mapping.rs:658 10: sc_hla_count::_main::h3253537751c2ae0e (0x55c46d857bb2) at src/main.rs:227 11: sc_hla_count::main::h6c74ee357aa88633 (0x55c46d8544cf) at src/main.rs:139 12: std::rt::lang_start::{{closure}}::h4a39b9fb10bdccb7 (0x55c46d83fbf3) at /rustc/04488afe34512aa4c33566eb16d8c912a3ae04f9/src/libstd/rt.rs:67 13: std::rt::lang_start_internal::{{closure}}::h5d3ea623498f5f43 (0x55c46dc0aae3) at src/libstd/rt.rs:52 std::panicking::try::do_call::hac65e71be769a440 at src/libstd/panicking.rs:348 std::panicking::try::hd4706e264bcf6712 at src/libstd/panicking.rs:325 std::panic::catch_unwind::h948a0fb4a8b3ee82 at src/libstd/panic.rs:394 std::rt::lang_start_internal::h72cc068ed2d0ac53 at src/libstd/rt.rs:51 14: main (0x55c46d85c788) 15: __libc_start_main (0x7fa1e8b1f505) 16: (0x55c46d80fdcf) 17: (0x0)'''

Are you planning to verify and fix this bug recently? If not, I may need to change all my bam files and then fun scHLAcount.

Best, danshu