rm HEK293_snp_our_data.txt
while read line;do
    samtools mpileup -C 0 -r $line HEK293_aligned.sorted.bam >> HEK293_snp_our_data.txt
done < HEK_snp_paper_hg38_simple.txt
