-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathOrdering FASTA files to match the geneID order in a .txt file
37 lines (32 loc) · 2.58 KB
/
Ordering FASTA files to match the geneID order in a .txt file
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
# This is when I sorted the sig_down_isoforms
cd Trialling_Ayos_code/
head -10 down_sig.fasta
cat TRINITY_sig_down_isoform_seqs_copy.fasta | seqkit seq -n -i -s
cat TRINITY_sig_down_isoform_seqs_copy.fasta | sed -e 's/(^>TRINITY_*_*_*)_.*/$1/g' | head -5
cat TRINITY_sig_down_isoform_seqs_copy.fasta | seqkit seq -n -i -s
cat TRINITY_sig_down_isoform_seqs_copy.fasta | seqkit seq | head -10
cat TRINITY_sig_down_isoform_seqs_copy.fasta | paste - - | cut -f 1
cat TRINITY_sig_down_isoform_seqs_copy.fasta | seqkit fx2tab
cat TRINITY_sig_down_isoform_seqs_copy.fasta | seqkit seq -n -i -s
list=$(cat TRINITY_sig_down_isoforms_copy.txt)
cat $list
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do echo ${i}; done
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do seqkit grep -p ${i} TRINITY_sig_down_isoform_seqs_copy.fasta; done
cat TRINITY_sig_down_isoform_seqs_copy.fasta | seqkit seq -n -i -s
cat TRINITY_sig_down_isoform_seqs_copy.fasta | seqkit seq -n -i -s > TRINITY_sig_down_isoform_seqs_copy2.fasta
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do seqkit grep -p ${i} TRINITY_sig_down_isoform_seqs_copy2.fasta; done
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do seqkit grep -p ${i} TRINITY_sig_down_isoform_seqs_copy2.fasta; head -1; done
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do seqkit grep -p ${i} TRINITY_sig_down_isoform_seqs_copy2.fasta; grep -i '^>'; done
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do base=$i; seqkit grep -p '$base' TRINITY_sig_down_isoform_seqs_copy2.fasta; done
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do base=$i; seqkit grep -p '${base}' TRINITY_sig_down_isoform_seqs_copy2.fasta; done
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do base=$i; seqkit grep -p '${i}' TRINITY_sig_down_isoform_seqs_copy2.fasta; done
for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do base=$i; seqkit grep -p ${i} TRINITY_sig_down_isoform_seqs_copy2.fasta; done
# for i in $(cat TRINITY_sig_down_isoforms_copy.txt); do base=$i; seqkit grep -p ${i} TRINITY_sig_down_isoform_seqs_copy2.fasta; wc -l ; done
# parallel -j 1 "seqkit grep -f {} TRINITY_sig_down_isoform_seqs_copy2.fasta" :::: TRINITY_sig_down_isoforms_copy.txt
# seqkit grep -f {} TRINITY_sig_down_isoform_seqs_copy2.fasta | head -1
cat TRINITY_sig_down_isoform_seqs_copy2.fasta | head -1
parallel -j 1 "cat TRINITY_sig_down_isoform_seqs_copy2.fasta | seqkit grep -p '{}' >> down_sig_all.fasta" :::: TRINITY_sig_down_isoforms_copy.txt
tail -10 down_sig_all.fasta
tail down_sig_all.fasta | grep -i '^>'
cat down_sig_all.fasta | grep -i '^>' | tail
tail TRINITY_sig_down_isoforms_copy.txt