bedtools2 icon indicating copy to clipboard operation
bedtools2 copied to clipboard

maskfasta/maskFastaFromBed: No option for minus (-) strand

Open githubgig opened this issue 5 years ago • 0 comments

Hello,

I'd like to use maskFastaFromBed to softmask the genome, but it doesn't consider strand information. Example run and output is below.

Input fasta file: genome.fasta

Chr1 ATTGACAGAAATATCATCACATCTATTCTTTCTCTCCCCTAGTTTAGCAAAT Chr2 GACATATAAATAATAGTGGGAAAGAGACCGGATGAAACCTCAACTGTGGCTTTCATTAACAGATCA

Input repeat file: genome.bed

Chr1 0 17 for 1 + Chr2 0 17 rev 1 -

Command: maskFastaFromBed -soft -fi genome.fasta -bed genome.bed -fo genome_softmasked.fasta

Output

Chr1 attgacagaaatatcatCACATCTATTCTTTCTCTCCCCTAGTTTAGCAAAT Chr2 gacatataaataatagtGGGAAAGAGACCGGATGAAACCTCAACTGTGGCTTTCATTAACAGATCA

Is there a way to include strand information in the command?

Thank you.

githubgig avatar Feb 05 '20 05:02 githubgig