bash - trim sequences and quality in fastq file -

-

i have bunch of fastq files in directory , want trim sequence 2 nucleotides , quality(if read has 51 base pairs , ends-with ctg or ttg).

here wrote shell script getting errors,need new shell scripting

input:

@hwi-st1072:187:c35yuacxx:7:1101:1609:1983 1:n:0:acagtg nggagaaagagagtgtgtttttagggggagatttttaaaatggttgttttg + #0<bfffffffff<bfffiifffffiiibfffffiifiiiiiffbffffff @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttattcgggaggttgagctg + #0<bffffffffffiibffiiiiiifiiiffiifiiifiifiiffffiiff @hwi-st1072:187:c35yuacxx:7:1101:9351:2210 1:n:0:acagtg cggttttgttttattttgtatgattaggagggttttggaggtttagttacc + bbbffffffffffiiiiiffiifiiiiiiiiiffiififiiffiiifiiii @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttat + #0<bffffffffffiibffiiiiiifiiiffiifi

output:

@hwi-st1072:187:c35yuacxx:7:1101:1609:1983 1:n:0:acagtg nggagaaagagagtgtgtttttagggggagatttttaaaatggttgttt + #0<bfffffffff<bfffiifffffiiibfffffiifiiiiiffbffff @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttattcgggaggttgagc + #0<bffffffffffiibffiiiiiifiiiffiifiiifiifiiffffii @hwi-st1072:187:c35yuacxx:7:1101:9351:2210 1:n:0:acagtg cggttttgttttattttgtatgattaggagggttttggaggtttagttacc + bbbffffffffffiiiiiffiifiiiiiiiiiffiififiiffiiifiiii @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttat + #0<bffffffffffiibffiiiiiifiiiffiifi

script:

for sample in *.fastq;do     name=$(echo ${sample} | sed 's/.fastq//')     while read line;do         if [ ${line:0:1} == "@" ] ;                 head="${line}"                 $echo $head         elif [ "${head}" ] && [ "${line}" ] ;                 length=${#line}                 if [ "${length}" = 51 -a "${line}" =~ *ctg|*ttg ] ;                         sequence= substr($line,0,49)                         #echo $sequence                 fi         elif [ ${line:0:1} == "+" ] ;                 plus="${line}"                 #echo $plus         elif [ "${plus}" ] && [ "${line}" ] ;                 quality= substr($line,0,49)                 #echo $quality         fi         print "${head}\n${sequence}\n${plus}\n${quality}" > ${name}_new.fq    done < $sample done

don't 100% understand you're doing, fixed few things. try below

#!/bin/bash sample in *.fastq;   name="${sample/.fastq/}"   while read -r line;     if [[ $line == '@'* ]];       head="$line" && echo "$head" >> "${name}_new.fq"     elif [[ -n $head && ${#line} == 51 && $line =~ (ctg|ttg)$ ]];       sequence="${line:0:49}" && echo "$sequence" >> "${name}_new.fq"     elif [[ $line == '+'* ]];       plus="$line" && echo "$line" >> "${name}_new.fq"     else       quality="$line" && echo "$quality" >> "${name}_new.fq"     fi   done < "$sample" done

example output

> cat sample_new.fq  > cat sample.fastq @hwi-st1072:187:c35yuacxx:7:1101:1609:1983 1:n:0:acagtg nggagaaagagagtgtgtttttagggggagatttttaaaatggttgttttg + #0<bfffffffff<bfffiifffffiiibfffffiifiiiiiffbffffff @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttattcgggaggttgagctg + #0<bffffffffffiibffiiiiiifiiiffiifiiifiifiiffffiiff @hwi-st1072:187:c35yuacxx:7:1101:9351:2210 1:n:0:acagtg cggttttgttttattttgtatgattaggagggttttggaggtttagttacc + bbbffffffffffiiiiiffiifiiiiiiiiiffiififiiffiiifiiii @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttat + #0<bffffffffffiibffiiiiiifiiiffiifi  > ./abovescript  > cat sample_new.fq @hwi-st1072:187:c35yuacxx:7:1101:1609:1983 1:n:0:acagtg nggagaaagagagtgtgtttttagggggagatttttaaaatggttgttt + @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttattcgggaggttgagc + @hwi-st1072:187:c35yuacxx:7:1101:9351:2210 1:n:0:acagtg cggttttgttttattttgtatgattaggagggttttggaggtttagttacc + bbbffffffffffiiiiiffiifiiiiiiiiiffiififiiffiiifiiii @hwi-st1072:187:c35yuacxx:7:1101:1747:1995 1:n:0:acagtg nggttgtggtggtgggtatttgtagttttatttat +

Comments

Post a Comment

Popular posts from this blog

python - Subclassed QStyledItemDelegate ignores Stylesheet -

-
i'm trying subclass qstyleditemdelegate remove focus rectangle in qcombobox . even though i'm calling base implementation of paint function , nothing else, result different. looks if parts of style-sheet influence bounding box of item taken consideration. class pstyleditemdelegate(qstyleditemdelegate): def __init__(self, *args, **kwds): super(pstyleditemdelegate, self).__init__(*args, **kwds) def paint(self, *args, **kwargs): qstyleditemdelegate.paint(*args, **kwargs) what have paint unmodified qstyleditemdelegate ? as suggested tried replacing pyside pyqt4 , works, appears bug. update pyside 1.1.2 1.2.1 result same. unfortunately switch breaks other parts of code, if there no other suggestions i'm going accept answer. edit bug tracked here
Read more

java - HttpClient 3.1 Connection pooling vs HttpClient 4.3.2 -

-
i have used httpclient3.1(java) connection pooling below : import org.apache.commons.httpclient.httpclient; import org.apache.commons.httpclient.multithreadedhttpconnectionmanager; import org.apache.commons.httpclient.params.httpclientparams; import org.apache.commons.httpclient.params.httpmethodparams; import org.apache.log4j.logger; public class httpclientmanager { private static logger logger = logger.getlogger(httpclientmanager.class); private static multithreadedhttpconnectionmanager connectionmanager = new multithreadedhttpconnectionmanager(); private static httpclientparams clientparams = new httpclientparams(); static { try { clientparams.makelenient(); clientparams.setauthenticationpreemptive(false); clientparams.setparameter(httpmethodparams.so_timeout, new integer(30000)); // set so_timeout clientparams.setparameter(httpmethodparams.head_body_check_timeout, new integer(3000)); // head responses ...
Read more

node.js - StackOverflow API not returning JSON -

-
i using node.js script gather data stackoverflow. know responses gzipped, put code in should take care of that. script follows: var request = require('request'); var zlib = require('zlib'); function getstackoverflowresponse(url, callback){ request(url, {encoding: null}, function(err, response, body){ if(response.headers['content-encoding'] == 'gzip'){ zlib.gunzip(body, function(err, dezipped) { callback(dezipped); }); } else { callback(body); } }); } var url = "https://api.stackexchange.com/docs/questions#pagesize=2&order=desc&min=2014-01-04&max=2014-02-02&sort=activity&tagged=apigee&filter=default&site=stackoverflow&run=true"; getstackoverflowresponse(url, function(questions) { console.log(questions); }); instead of getting json output, i'm getting following response: buffer 0d 0a 0d 0a 0d 0a 0d 0a 3c 21 44 4f 43 54 59 50 45 20 48...
Read more