You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Hi! I was trying to mask a multi-fasta file with a bed file, and discovered that bedtools provides a program that does exactly this, but at the same time not exactly. The program masks a fasta file with complete chromosomes. So if I have
test.fa
chr1
AAATGTTTAT
chr2
ATTTTGG
and test.bed
chr1 0 1
chr2 2 3
The result of bedtools maskfasta -fi test.fa -fo test_masked.fa -bed test.bed
Would be a file test_masked.fa
chr1 NAATGTTTAT
chr2
ATNTTGG
However, if I have a file like this
test.fa
chr1:1-5
AAATG
chr1:6-10
TTTAT
chr2:1-7
ATTTTGG
And I run the same line of code, I will obtain a file where the sequences are not changed, but at the same time no error is printed on the screen, saying that the program could not find a match between the fasta file and the bed file.
chr1:1-5
AAATG
chr1:6-10
TTTAT
chr2:1-7
ATTTTGG
I feel it would be very useful if bedtools would implement a masking also for non complete sequences, or at least print an error when no match is found
The text was updated successfully, but these errors were encountered:
Hi! I was trying to mask a multi-fasta file with a bed file, and discovered that bedtools provides a program that does exactly this, but at the same time not exactly. The program masks a fasta file with complete chromosomes. So if I have
test.fa
and test.bed
chr1 0 1
chr2 2 3
The result of
bedtools maskfasta -fi test.fa -fo test_masked.fa -bed test.bed
Would be a file test_masked.fa
However, if I have a file like this
test.fa
And I run the same line of code, I will obtain a file where the sequences are not changed, but at the same time no error is printed on the screen, saying that the program could not find a match between the fasta file and the bed file.
I feel it would be very useful if bedtools would implement a masking also for non complete sequences, or at least print an error when no match is found
The text was updated successfully, but these errors were encountered: