Skip to content

Commit 21519c9

Browse files
authored
Update SARS-CoV-2_get_ngs.sh
1 parent 60e5c14 commit 21519c9

File tree

1 file changed

+2
-11
lines changed

1 file changed

+2
-11
lines changed

SARS-CoV-2_get_ngs.sh

+2-11
Original file line numberDiff line numberDiff line change
@@ -104,20 +104,11 @@ echo ""
104104

105105
b= ls -1 *.fastq.gz.fastp
106106
for b in *.fastq.gz.fastp; do minimap2 -ax sr ${2} ${b} > ${b}.sam -t ${3}
107+
samtools sort ${b}.sam > ${b}.sam.sorted.bam -@ ${3}
108+
rm ${b}.sam
107109
rm ${b}
108110
done
109111

110-
###################################################
111-
# Sorting SAM files, using n threads (-@ option) #
112-
###################################################
113-
114-
echo "Sorting SAM files, using n threads."
115-
echo ""
116-
117-
c= ls -1 *.sam
118-
for c in *.sam; do samtools sort ${c} > ${c}.sorted.bam -@ ${3}
119-
done
120-
121112
######################
122113
# Renaming BAM files #
123114
######################

0 commit comments

Comments
 (0)