This is a test version of Biostars. For the public version, visit https://www.biostars.org.
creating BAM? alignment files for IGV

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?

r genome alignment

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:

chromName     gene1StartPos      gene1EndPos
chromName     gene2StartPos      gene2EndPos
...

And yes, to display on IGV you do need the reference genome.

0 answers

No answers yet.

Log in to answer this question.