How to calculate percent coverage of each read for nanopore

0

Hello, I am assigned to investigate the read the cover the reference more than 80%. I know that samtools can be used to calculate coverage and depth for the interested regions but what about each read?
For example, I have read 1-10 that cover region A, I want to know their percent coverage for all of each read.
Expect output: read 1 70% cover, leftmost position 1 (like sam format column)
Is there any tools or script that can be used for this?
Thanks in advance.


alignment


nanopore


sequencing

• 23 views



Source link