/data/hasindu/realmuc1/install/bedpull/target/release/bedpull --bam /home/hasindu/scratch/realmuc1/KDG0246PR/KDG0246PR.muc1_analysis_out/KDG0246PR.muc1_flvntr-hap1_chm13.bam --bed /data/hasindu/realmuc1/data/ref_files/muc1_vntr_region_slop.chm13.bed -o a.fa
Opts {
bam: "/home/hasindu/scratch/realmuc1/KDG0246PR/KDG0246PR.muc1_analysis_out/KDG0246PR.muc1_flvntr-hap1_chm13.bam",
reference: "None",
bed: "/data/hasindu/realmuc1/data/ref_files/muc1_vntr_region_slop.chm13.bed",
paf: "None",
query_ref: "None",
output: "a.fa",
use_paf_index: true,
fastq: false,
}
Reading bed file
BedRecord { chrom: "chr1", start: 154327978, end: 154330908, name: None }
region: Region { name: [99, 104, 114, 49], interval: Interval { start: Some(Position(154327978)), end: Some(Position(154330908)) } }
Bam mode
Extracting sequences
===============================
Analysing region: chr1:154327978-154330908, chr1:154327978-154330908
===============================
thread 'main' (895344) panicked at src/main.rs:86:65:
Couldn't read header: Custom { kind: InvalidData, error: InvalidRecord(InvalidValue(InvalidProgram(InvalidOther(Other("VN"), Missing)))) }
stack backtrace:
0: __rustc::rust_begin_unwind
1: core::panicking::panic_fmt
2: core::result::unwrap_failed
3: bedpull::main
note: Some details are omitted, run with `RUST_BACKTRACE=full` for a verbose backtrace.
GTA100
/data/hasindu/realmuc1/install/bedpull/target/release/bedpull --bam /home/hasindu/scratch/realmuc1/KDG0246PR/KDG0246PR.muc1_analysis_out/KDG0246PR.muc1_flvntr-hap1_chm13.bam --bed /data/hasindu/realmuc1/data/ref_files/muc1_vntr_region_slop.chm13.bed -o a.fa