This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Sequence extract from a big fast file

Dear All,

After sequencing reaction and extensive filtering, we converted ".fastq" files into ".fasta" files and each fasta file has approximately ~67 million reads. My question is that is there any script for extracting first 30 million reads, then remaining 37 million reads with sequential manner.

Thank you for all your help!

next-gen-sequencing rna-seq
• 3,164 views • link
updated by Ram Ram • written by vahapel vahapel

2 answers

Assuming no linebreaks in sequences, this is as simple as:

head -n x file > output

Where x is number of seqs times 2 (one line for header, one for sequence). Similarly, you can get the last x sequences utilizing tail.

• 0 views • link
written by 5heikki 5heikki

Thak you so much, 5heikki. It seems very practical way for such a purpose !

• 0 views • link
written by vahapel vahapel

split fasta file

• 0 views • link
written by GouthamAtla GouthamAtla

Hi, Geek_y, "split fasta file" scripts will be useful for my project, thank you for your help!

• 0 views • link
written by vahapel vahapel

Hi,

I am a novice using R to split a FASTA file with 300 000 contigs into 6 file of less than 50 000 contigs. I have seen many options but would anyone advise anything that I could be used in R? Thank you A

• 0 views • link
written by alexandra.lanot • alexandra.lanot

Log in to answer this question.

More posts like this