Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- /data/data/bowtie2-2.3.5.1-linux-x86_64/bowtie2-build /data/data/spike.fa spike
- #sumest i viena folderi spike paskui ir
- /data/data_fastq/reads_colapsed # export BOWTIE2_INDEXES=/data/data_fastq/reads_colapsed/spike/
- /data/data_fastq/reads_colapsed # for file in *.fa; do /data/data/bowtie2-2.3.5.1-linux-x86_64/bowtie2 -f -N 0 --no-unal -x spike -U $file -S ${file%fa}sam; done
- # for file in *.sam; do grep -v "@" $file | awk -F "\t|x" '{print $4, " ", $2}' | sort -k 1 | awk '{ seen[$1] += $2 } END { for (i in seen) print i, seen[i] }' | sort -k 1 > ${file}.out; done
- 983aaaec4dff /data/data_fastq/reads_colapsed # #for file in *.sam; do grep -v "@" $file | awk -F "\t|x" '{print $4, " ", $2}' | sort -k 1 > ${file}.out; done
- # tmp yra visur po 0 dar, kad susivienodintu eiluciu sk
- for file in *.out; do cat tmp >> $file; done
- for file in *.out; do cat $file | awk '{ seen[$1] += $2 } END { for (i in seen) print i, seen[i] }' | sort -k 1 > ${file}.out; done
- pirmas stulpelis tai tmp su nauje eilute
- for file in ../????; do echo ${file:3:8} > ${file:3:8}; done
- for file in *.out.out; do awk '{print $2}' $file >> ${file:0:4};done
- paste -d "," tmp ???? > composition_matrix_spike.csv
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement