|
@@ -1 +1,26 @@
|
1
|
1
|
#!/bin/bash
|
|
2
|
+
|
|
3
|
+mkdir results
|
|
4
|
+
|
|
5
|
+# Question 1
|
|
6
|
+soft/bowtie2 --end-to-end --fast -x databases/all_genome.fasta -1 EchE_R1.fastq.gz -2 EchE_R2.fastq.gz -S results/alignment.sam
|
|
7
|
+
|
|
8
|
+# Question 2
|
|
9
|
+
|
|
10
|
+# a
|
|
11
|
+samtools view -b results/alignment.sam
|
|
12
|
+
|
|
13
|
+# b
|
|
14
|
+samtools sort -@ 4 results/alignment.bam -o results/alignment.sorted.bam
|
|
15
|
+
|
|
16
|
+# c
|
|
17
|
+samtools index results/alignment.sorted.bam
|
|
18
|
+
|
|
19
|
+# d
|
|
20
|
+samtools idxstats -@ 4 results/alignment.sorted.bam > results/alignment.idx
|
|
21
|
+
|
|
22
|
+# f
|
|
23
|
+grep ">" databases/all_genome.fasta|cut -f 2 -d ">" > results/association.tsv
|
|
24
|
+
|
|
25
|
+# Question 3
|
|
26
|
+
|