Masking repetitive sequences ...can I then map using BWA?

0

Entering edit mode

10.5 years ago

ab.tsubaki ▴ 50

Hi there

In need of help please!

I'm masking the repetitive regions in my wheat sequencing data using RepeatMasker. This generates a .masked file. Next I want to map or compare this masked file to other genome sequences but BWA uses fastq if I remember correctly. How can I convert the .masked file to a useable format? Or can it be used as-is???

Appreciate any comments, suggestions or previous experience :-)

Anandi

next-gen blast • 3.5k views

ADD COMMENT • link updated 3.1 years ago by Ram 44k • written 10.5 years ago by ab.tsubaki ▴ 50

0

Entering edit mode

Do you have to use BWA? It's intended for short reads, not whole genome comparisons. Why don't you just use a tool designed to compare genomes?

ADD REPLY • link 10.5 years ago by Devon Ryan 104k

0

Entering edit mode

Well, I'm not QUITE comparing a whole genome. My contigs are quite short and that's what I'll be mapping first.

Any advice on other mapping programs will be appreciated? I know BWA and Bowtie and I don't think either of them makes mention of .masked files?

ADD REPLY • link updated 4.9 years ago by Ram 44k • written 10.5 years ago by ab.tsubaki ▴ 50

0

Entering edit mode

The .masked file is just a fasta file, at least that's my recollection.

ADD REPLY • link 10.5 years ago by Devon Ryan 104k

0

Entering edit mode

Thanks! I looked at the output from RepeatMasker and saw that it's a fasta file just now! Clever me!

ADD REPLY • link 10.5 years ago by ab.tsubaki ▴ 50

Login before adding your answer.