-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathextract.sh
executable file
·39 lines (27 loc) · 908 Bytes
/
extract.sh
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
#!/bin/bash
gen=D_multi_fasta.fa
aa=$(cat $gen | grep "Fatima" | sed 's/[<>,]//g')
for f in $aa;
do
samtools faidx D_multi_fasta.fa $f >> D_RLG_Tura_Fatima_B_210J24-1.fa
done
aa=$(cat $gen | grep "Atau_Sabrina" | sed 's/[<>,]//g')
for f in $aa;
do
samtools faidx D_multi_fasta.fa $f >> D_RLG_Atau_Sabrina_AY368673-2.fa
done
aa=$(cat $gen | grep "Tura_Sabrina" | sed 's/[<>,]//g')
for f in $aa;
do
samtools faidx D_multi_fasta.fa $f >> D_RLG_Tura_Sabrina_210J24-1.fa
done
aa=$(cat $gen | grep "Tura_Angela" | sed 's/[<>,]//g')
for f in $aa;
do
samtools faidx D_multi_fasta.fa $f >> D_RLC_Tura_Angela_210J24-3.fa
done
aa=$(cat $gen | grep "Atau_Angela" | sed 's/[<>,]//g')
for f in $aa;
do
samtools faidx D_multi_fasta.fa $f >> D_RLC_Atau_Angela_AF497474-3.fa
done