Hi everyone, I'm new in this forum. I'm a biology student from Italy.
I have to create a plot using a values obtained from an alignment from DNA sequences (from BLAST).
the results are in this format:
seq. id is the name of my sequence, coverage % is how much my DNA sequence covers on the reference, lenght is the lenght in nucleotides of my sequences paired on the reference, and start/stop are the numbers of the DNA bases of the reference where my sequence aligns.
My aim is to plot a thing like this:
To make this, I have to generate all the values from the start to the stop, and give these numbers the corrispective % of coverage (example with 90% of coverage the start is at 190 bp and the stop at 800 bp, so 190 90%,191 90%,192 90%,193 90%..... 798 90%, 799 %, 800 %, and write it in a new table).
How I can do it? I import my txt table on R using read.table, so I think that is a data.frame.
I'm going to use ggplot2 to plot it
thank you all in advance, greatings from Italy.