Trimming a section of a scaffold from a genome
1
0
Entering edit mode
Karen • 0
@b10a14c9
Last seen 11 days ago
United States

I am having trouble using Biostrings. I need to locate a section of a specific scaffold by position in my genome assembly and remove it then save as two separate scaffolds. I am not sure how to do this using BioSstrings, if its even possible.

Thanks

```

Biostrings GenomeAssembly • 123 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 13 hours ago
United States

It's undoubtedly possible, but you don't provide enough information to go on. If I assume you are talking about an existing BSgenome package, you can do it like this.

> library(BSgenome.Hsapiens.UCSC.hg38)
>  getSeq(Hsapiens, GRanges(c("chr1:12345-54321","chr2:12345-55555")))
 DNAStringSet object of length 2:
     width seq
 [1] 41977 TCAGACCAGCCGGCTGGAGGGAGGGGCTCAGCAG...ACATACTTTTTTCATGTATTTGTTTCTACTGCT
 [2] 43211 AACAATAATAGTCTTACTTTGATTATGATCAAAA...TTGGTCAAAACAACTCTGGGTAAAGACTTTAGG

Login before adding your answer.

Traffic: 630 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