How to extract variable nucleotide regions from a list of contigs
I have a list of about 6000 NCBI contig accession numbers and I'd like to extract a specific 30kb region from each contig. I have, in a separate file, a list of "begin" and "end" indices that represent the region of interest for each contig. Is there a way to retrieve a fasta file of these trimmed contigs using an Entrez query? Given that I'm not a very skilled programmer, I would have just put the sequences in an excel spreadsheet and trimmed accordingly, but the sequence strings are too long to fit in an excel cell. Is there a simple way to do this via E-Utilities?
• 2,393 views
•
link
updated
by
lakhujanivijay
•
written
by
kayrouz.1 •
•
written
by
kayrouz.1 •
1 answer
To the point explanation of bedtools getfasta is here
• 0 views
•
link
written
by
lakhujanivijay
Log in to answer this question.
More posts like this
-
How to retrieve the corresponding genome assembly accession ID from a nucleotide accession ID using…
written by Mani •
-
I need to retrieve a set of protein and mRNA sequences
written by george •
-
retrieve all nt sequences for a taxid list
written by pe_se •
-
Efficiently retrieve the molecular type for a list of NCBI accessions
written by john
-
Batch Retrieval of Conserved Domains from a List of IDs?
written by kayrouz.1 •
-
Batch Annotation of Bacterial Gene Clusters
written by kayrouz.1 •
-
Extracting a Genomic Region From Nucleotide Accession
written by 13dsc •
-
Annotating gene symbols with BioPython
written by bojingjia •
-
Calculating N50 in Excel?
written by fantasticfaces2 •
-
Converting Hgnc Gene Symbol To Gi Number
written by hallister01 •
Do you have a reference genome? If you do you could use
bedtools getfasta.