gymrek-lab / LongTR

Tandem repeat genotyping with long reads
GNU General Public License v2.0
22 stars 0 forks source link

segfault when using --skip-assembly on pacbio HG002 bam #14

Closed bw2 closed 3 months ago

bw2 commented 3 months ago

I'm having trouble running

LongTR --skip-assembly --bams HG002.subset.bam  --fasta Homo_sapiens_assembly38.fasta --regions repeats.bed --log out.log --tr-vcf out.vcf.gz

on these files: https://www.dropbox.com/scl/fi/ik6umjhkjujwjlpnp0ia4/files.zip?rlkey=nfl6lg4a4g9jgzjryxqs4t897&st=v4cccuat&dl=0

but don't see an obvious reason:

image

It completed successfully on another pacbio bam from a different source.

heliziii commented 3 months ago

Hello,

We made the skip-assembly tag as default as it was causing so much confusion, please remove it from your tags. I have updated the document, I apologize for the inconvenience.

Best, Helia

bw2 commented 3 months ago

Will do, thanks.