extract first and last n bp from fasta file from multiple fasta file in R
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Covert multiple short sequence into list of IUPAC motif
written by praasu •
-
comparing the intron/ intron location with the reference transcript in the orthologous transcripts
written by praasu •
-
How can I identity Branch point sequence from Intron?
written by praasu •
-
Validating splice site using protein sequence
written by praasu •
-
Software for multiple sequence alignment?
written by Ankit
-
Where I can find splice-leader (SL) sequence set ?
written by praasu •
-
How to extract the first and last 2 nt from multiple sequences in a FASTA file?
written by praasu •
-
how to extract intron coordinates (in bed) from bam
written by praasu •
-
How to extract the first and last N bases from a read in a fastq file?
written by lovegenetics •
-
Extracting Intron-Exon Reads from bam files
written by praasu •
Its not an answer but suggestion. Seqkit's subseq command does this. e.g. https://bioinf.shenwei.me/seqkit/tutorial/#play-with-mirna-hairpins
If you really want to do it in R you would want to loop over as you stream the file or else it would hog your memory really bad. Or you can call seqkit with system :).
https://stackoverflow.com/questions/42492351/stream-processing-large-csv-file-in-r
https://stackoverflow.com/questions/12626637/read-a-text-file-in-r-line-by-line
OP asked the exact same question 2 years ago and got the seqkit subseq answer from shenwei himself. Really odd behavior.