Entering edit mode
9.1 years ago
ad
▴
30
Hi, I have coordinate data for several sequences I want to align against genomes in IGV. The preferred format for alignment is said to be BAM. Unfortunately I haven't found a straightforward way of creating BAMs straight from coordinates. I assume I'll need the reference genome as well. Is there any simple way to do it from R/Bioconductor?
Do you have sequences you want to align, or do you have sequences (whose coordinates you know already) which you to display? If you have the genomic coordinates, all you need is a tab-delimited .BED file with:
And yes, to display on IGV you do need the reference genome.