How to extract ~2M short sequences based on coordinates from a 3G fasta file?
(Sorry for the Bioinformatics 101 question!)
I have a file with ~2 million sets of coordinates (chromosome<tab>begin-position<tab>end-position), corresponding to short (~50nt) human genomic sequences (hg19). I want to extract the actual sequences from a human genome assembly 19 fasta file (~3.0G).
I imagine this is a relatively common task, and therefore, that there must be standard tools to carry it out efficiently.
Sadly, my Google fu has not been up to the task of finding them.
I would appreciate not only the name of a tool to use, but also the command line one would use, especially if there are important flags and options I should be aware of when using such a tool.
1 answer
Log in to answer this question.
More posts like this
-
Computing pairwise AA global sequence alignment between all pairs in a vector of sequences
written by rubic
-
Annotate a list of rsIDs with minor allele freq for specific populations
written by petrstep •
-
On the conventions for annotating the position of insertion and deletion variants
written by kynnjo •
-
How to use plink to consolidate multiple single-sample vcf files?
written by kynnjo •
-
Filter out all known genes & regulatory elements for a given genome in a local blast search
written by Proteus00 •
-
Consensus by aligning and extending the sequences
written by BDK_compbio
-
Inserting millions of records into MySQL database using Python
written by Prasad •
-
randomly sequence cleavage - average length constraint
written by Nifaste •
-
Extract reads with variant at given position
written by philip.jonsson •
-
Human Whole Genome Sequence Data With Certain Read Length And Coverage
written by Steffi
