#!/bin/bash
WORK_DIR="/sf/storage/gmo/hic/out/by_User/KhabarovaAnya/Pig/raw/Orlov_NGS_15_05_25/"
cd $WORK_DIR || { echo "Невозможно перейти в директорию $WORK_DIR"; exit 1; }

# Объявление образцов без R1 и R2, только общие части
samples=("20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_1_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_3_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_4_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_5_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_6_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_7_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_8_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_9_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_10_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_11_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_12_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_13_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_14_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_15_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_16_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_17_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_18_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_19_L00" "20250513130908_RA01230901018_5P241020135US2S2656DX_seqonc9924_subcl_1903_20_L00")

# Создание директории для выравнивания, если её ещё нет
mkdir -p ${WORK_DIR}/align_paired

for sample in "${samples[@]}"; do
    echo "Обрабатываем образец: ${sample}"
    echo "=========="

    # Выполняем выравнивание
    bwa mem /sf/storage/gmo/hic/out/by_User/KhabarovaAnya/pig_genome/GCF_000003025.6_Sscrofa11.1_genomic.fna ${WORK_DIR}/${sample}_R1.f*q.gz ${WORK_DIR}/${sample}_R2.f*q.gz > ${WORK_DIR}/align_paired/${sample}.sam

    # Преобразуем SAM в BAM
    samtools view -b -h ${WORK_DIR}/align_paired/${sample}.sam > ${WORK_DIR}/align_paired/${sample}.bam 

    # Сортируем BAM
    samtools sort -o ${WORK_DIR}/align_paired/${sample}.sorted.bam ${WORK_DIR}/align_paired/${sample}.bam 

    # Индексируем отсортированный BAM файл
    samtools index ${WORK_DIR}/align_paired/${sample}.sorted.bam 
done