This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can I map certain values for gene locus against genome?

I have a data like below. I want to make a plot where I want to map each locus to the genome sequence (Chr 1..15 as X axis) and show the counts in Y axis. How can I create that type of plot? Is there any bioconductor package to do this?

locus   counts
SS1G_03009  40

SS1G_02499  10

SS1G_09248  10

SS1G_09768  22

SS1G_04468  22

SS1G_08838  11

SS1G_02828  23
r
• 1,611 views • link
written by MAPK MAPK

1 answer

Manhattan plot using R?

• 0 views • link
written by GenoMax GenoMax

Yes, that would do, but I need to get the chr, start and end information for the given loci. How can I extract that from NCBI or using genes and genome fasta files? Or using the package you have suggested, I would need BP variable, so I am not sure how I can obtain that.

• 0 views • link
written by MAPK MAPK

Recall the NCBI unix utils answers provided by @vvkodali to one of your past questions. You should be able to use that to get the info you are asking for above.

• 0 views • link
written by GenoMax GenoMax

Thank you so much, I just got it :)

• 0 views • link
written by MAPK MAPK

Log in to answer this question.

More posts like this