ChenDepp

Results 11 comments of ChenDepp

@brentp thanks for your reply, I had solved the problem , good luck. i hava an another question about the five file *.disc.bam, *.disc.bam.orig.bam, *.histo, *.split.bam, *.split.bam.orig.bam. the first question...

@christian-lanius Thanks you, I will dissect the mm_idx_load and mm_idx_dump functions, But i am not good at C language, If you can provide the corresponding python code,will save me a...

hi when i running partition step, find it can't get result. The report message is as follows: ![aaf0fad14f461581b86f28bd72691ac](https://user-images.githubusercontent.com/48479509/189023225-cc614520-9f2a-4ef6-9d87-bebc46b2e417.png)

> Hi @ChenDepp > This is probably an error with HTSJDK's CSI decoder for BAM files. Can you reindex your bam file with BAI index and try again. Looks like...

> Hmn, we might have an off by one error somewhere... @lbergelson hi,Are you saying it's gatk's problem, or is it my problem when dealing with the file?

hi @gokalpcelik I used GenomicsDBImport to replace CombinedGVCFs, but it has new problems, GenotypeGVCFs for GenomicsDB is so slow, can only get 900K interval vcf in 9 hours. how can...

when i running vg surject convert gam to bam, it report warnning information as below ![Image](https://github.com/user-attachments/assets/8a6e061d-0ce7-4bff-923f-08e0a64731a2) Is this the reason why gatk cannot run successfully?

@adamnovak As you said, the chromosome names of the pan-genome are chr1、chr2. so the problem may be that gatk can't accept the bam converted by vg-surject. Thank you for your...

@adamnovak I have another question. A sample has 35016715 paired-end reads, but using vg stat to analyze the gam file gives the following information: ![Image](https://github.com/user-attachments/assets/9b0d25a0-2a4e-43f3-b997-afa94a28264f) there hava 37157187 reads,how to...