-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathsort_samfile.sh
More file actions
22 lines (18 loc) · 765 Bytes
/
sort_samfile.sh
File metadata and controls
22 lines (18 loc) · 765 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
#!/bin/sh
#############################
# File Name : sort_samfile.sh
#
# Purpose : [???]
#
# Creation Date : 05-04-2024
#
# Last Modified : Mon 08 Apr 2024 02:38:23 PM PDT
#
# Created By : Roozbeh Dehghannasiri
#
##############################
cut -f2 sig_ancors_in_normal_or_immune_mapped_to_peter_cdna_compactors.tsv | grep -v "comp" | makefasta > compactors.fasta
/oak/stanford/groups/horence/Roozbeh/software/minimap2/minimap2 -N 0 -a /oak/stanford/groups/horence/Roozbeh/Sponge_project/Granny_Data/Jacob_assembly/Peters_cDNAs_11-20.fasta compactors.fasta > compactors_minimap.sam
samtools view -bS compactors_minimap.sam > compactors_minimap.bam
samtools sort compactors_minimap.bam compactors_minimap_sorted
samtools index compactors_minimap_sorted.bam