Split reference genome into callable regions by splitting on NNNN stretches
I am looking for a tool or piece of code that can take in a reference genome fasta file and output a bed file with X (equally sized) callable regions. The output can also be anything similar to the content of such a bed file.
These callable regions should be bordered upstream and downstream by any of these:
- The start of the chromosome
- The end of the chromosome
- a stretch of 250bp unknown (N) nucleotides (or e.g. 1000bp N)
This to get regions that can be variant called in parallel, without the risk of sequencing reads and variants going over the borders of the callable regions.
• 1,352 views
•
link
updated
by
Pierre Lindenbaum
•
written
by
William
•
written
by
William
1 answer
- use picard ScatterIntervalsByNs to extract the non-poly N regions.
- convert this interval list to bed with awk.
- run bedtools makewindows to generate your regions.
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Find reference genome regions spanned by only mapping quality 0 reads in multiple WGS samples
written by William
-
Reduce poly-N regions longer than 200bp back to 200bp in a reference genome FASTA
written by William
-
How to add transcript and cpg site information to a bed file
written by K.patel5 •
-
Count all (non N) nucleotides within regions defined in BED file
written by William
-
Best way to find intervals for parallelization of joint variant calling
written by William
-
How to use a bash variable in awk calculation
written by nadarajan_v •
-
Confused about TSS coordinates
written by srhic •
-
How to replace a stretch of a sequence with N's in a fasta file
written by Fluorine
-
bedop: output overlap for each region in a file
written by User 7754
-
Extract Coverage Track From Bam To Bed File
written by William
