Skip to content

Issue with subsampling normal sample #1

Description

@KJHeub

Hello, i have issues with using this scirpt "downsample_normal_gatk.sh"
The issues i found :

the Singularity image checking threw errors so i changed it to this :
" # Local image names with correct .sif filenames
singularity_local_images=()
for image in "${singularity_images[@]}"; do
image_name=$(basename "$image" | tr ':' '_')
singularity_local_images+=("$image_name.sif")
done

for i in "${!singularity_images[@]}"; do
image_uri="${singularity_images[$i]}"
image_file="${singularity_local_images[$i]}"

if [ ! -f "$image_file" ]; then
echo "Pulling Singularity image: $image_file"
singularity pull docker://"$image_uri"
else
echo "Singularity image: $image_file already exists. Skipping pull."
fi
done
"

"-f "not duplicated" seems to be an invalid syntax. https://github.com/biod/sambamba/wiki/%5Bsambamba-view%5D-Filter-expression-syntax

this command produced empty compressed files :
" $bedtools bamtofastq -i $OUTPUT_FILE -fq $OUTPUT_FQ | gzip > "$OUTPUT_FQ.gz" "
so i changed it to this :
" $bedtools bamtofastq -i "$OUTPUT_FILE" -fq "$OUTPUT_FQ"
gzip -c "$OUTPUT_FQ" > "$OUTPUT_FQ.gz" "

the problem is that the resulting subsamples give this error during pon build " [E::aux_parse] Incomplete aux field samtools sort: truncated file. Aborting" and " [E::sam_parse1] SEQ and QUAL are of different length samtools sort: truncated file. Aborting". I hope you can help me figure out why. Thank you in advance.

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