bp_scripts/remove_adaptor_454.sh

   1 #!/bin/sh
   2
   3 # Remove adaptor from sequences and produce a number of plots.
   4
   5 # Usage: read_fastq -i test.fq | remove_adaptor_454.sh | write_fastq -xo test_no_adaptor.fq
   6
   7 adaptor='TCGTATGCCGTCTTCTGCTTG' # 454 adaptor
   8 pid=$$
   9
  10 plot_lendist -t post -o 00_remove_adaptor_lendist_seq_before.$pid.ps -k SEQ_LEN |
  11 find_adaptor -a $adaptor -p |
  12 plot_histogram -t post -o 01_remove_adaptor_histogram_adaptor_pos.$pid.ps -k ADAPTOR_POS -s num |
  13 plot_lendist -t post -o 02_remove_adaptor_lendist_adaptor_len.$pid.ps -k ADAPTOR_LEN |
  14 plot_lendist -t post -o 03_remove_adaptor_lendist_seq_after.$pid.ps -k SEQ_LEN |
  15 analyze_vals -k ADAPTOR_POS,ADAPTOR_LEN -o 04_remove_adaptor_analyze_vals.$pid.txt |
  16 clip_adaptor