echo "" && echo "Locating and counting MID tags ... "
read_fastq -i $fq_file |
progress_meter |
+ clip_seq |
find_mids |
write_tab -o $table_mid -c -k MID_NUM,MID_SEQ,MID_COUNT -x
echo "" && echo "Locating and counting MID tags for sequences longer than 250 ... "
read_fastq -i $fq_file |
progress_meter |
+ clip_seq |
grab -e 'SEQ_LEN >= 250' |
find_mids |
write_tab -o $table_mid_len -c -k MID_NUM,MID_SEQ,MID_COUNT -x
echo "" && echo "Locating and counting MID tags for sequences longer than 250 and mean score above 20 ... "
read_fastq -i $fq_file |
progress_meter |
+ clip_seq |
grab -e 'SEQ_LEN >= 250' |
mean_scores |
grab -e 'SCORES_MEAN >= 20' |