Many thanks for these links. Your solution with flex work well.
• 0 views
•
link
written
by
marcofraca •
How Do I create a BED file from an hard masked fasta?
Hi to everyone,
I would like to create a BED file with the position hard masked in a fasta file.
This is the hard masked fasta:
>seq1
ACGTNNNNAGTCAGNNNNACGTC
>seq2
NNNNACGCATGNNNACGT
And this is the output bed file that I want:
seq1 5 9
seq1 16 19
seq2 1 4
seq2 12 14
• 3,769 views
•
link
updated
by
hewm2008 •
•
written
by
marcofraca •
•
written
by
marcofraca •
2 answers
duplicate : Identifying long stretches of Ns ; Annotating/Finding "NNNNNN" scaffold gaps in genome assembly
• 0 views
•
link
written
by
Pierre Lindenbaum
This tools : https://github.com/BGI-shenzhen/Reseqtools by this : iTools Fatools findN -InPut in.fa -OutPut N_Region.bed
• 0 views
•
link
written
by
hewm2008 •
Log in to answer this question.
More posts like this
-
NcbiblastpCommandline alignment results are different from blast webpage
written by laila.jafari •
-
Creating column that groups data from another columns
written by Bioinfo •
-
Blastx: FASTA-Reader: Ignoring invalid residues at position(s): On line 25975584
written by nikhilshinde0909 •
-
Sort Sam file by chromosome and position in python
written by hectorbaigessanchez •
-
MGRA returns error
written by little_more •
-
blast using python
written by santhoshhegde278 •
-
coverageBed problem: bam contains chromosome while query file does not
written by Chico Fernandez •
-
Extract only sequences that have N nucleotides
written by waqasnayab
-
how to compare sets using python (dealing with PDB file)
written by Jason Lin •
-
Split Seven Rows In Per Column
written by 2011101101 •
the
twoBitInfocommand from kentUtils (available here http://hgdownload.soe.ucsc.edu/admin/exe/linux.x86_64/) has an option-nBedto get a bed file of N stretches. You will need to convert the fasta to twobit by another commandtwoBitToFaalso available at the same site.