Skip to content

merge bams issue #170

@xuefenfei712

Description

@xuefenfei712

Hey there

Thanks for Graphtyper, I am call SNPs but merged bams from the same individuals, and got the error below, can I know if you have any solutions,
TMPDIR=/disk/regine/data2/xuexue/tmp/ ~/software/graphtyper.27 genotype ref.fasta --sams=bam --region=chr15 --threads=10 --sams_index=bai

[2026-01-12 16:20:42.956] caller.cpp:2534 We found file with multiple samples, sorry, this is currently not supported.

but if I rename the header, I still got the issue below

base=$(basename "$bam" .bam)
out=fixed_bams/${base}.fixed.bam
echo "[FIXING] $bam -> SM=${base}"
samtools addreplacerg -r ID:${base} -r SM:${base} -r PL:ILLUMINA -r LB:${base} -r PU:${base} -o "$out" "$bam"
samtools index "$out"

[2026-01-12 16:34:29.330] hts_parallel_reader.cpp:309 Reads with name=f# both have IS_FIRST_IN_PAIR=1

Thankyou

Xue

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions