Entering edit mode
21 months ago
m166108
•
0
Hi I have a RNA-seq data BAM file. I want to check the Defective viral genome, and I can check the read junction in IGV. But I need a file that record the Reads' breakpoint.
EX: Readname | breakpoint position(5') | breakpoint position(3')
I used to the bedtool bamtobed but it just record the reads position in column 2 & 3