bedtools2
bedtools2 copied to clipboard
maskfasta/maskFastaFromBed: No option for minus (-) strand
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.