commandlines and sequences
1
0
Entering edit mode
5.4 years ago

Which command can I use to grab the first 10 nucleotides in each line of my nucleotide sequence in a file?

sequence • 1.1k views
ADD COMMENT
2
Entering edit mode

Use google please! enter image description here

ADD REPLY
0
Entering edit mode

What have you tried? Otherwise this is a pure programming question.

ADD REPLY
0
Entering edit mode

I'd suggest starting with the biopython cookbook and tutorial. The solution to your problem is in there.

ADD REPLY
0
Entering edit mode

Refer to tips and tricks on trimming ngs data online.

ADD REPLY
1
Entering edit mode
5.4 years ago
GenoMax 141k

You can use reformat.sh from BBMap suite to do this.

reformat.sh in=infile.fq out=outfile.fq ftr=9

If you have paired-end data then you would want to use something like this

reformat.sh in1=infile_R1.fq in2=infile_R2.fq out1=outfile_R1.fq out2=outfile_R2.fq  ftr=9
ADD COMMENT

Login before adding your answer.

Traffic: 2800 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6