Generate genome location and read ID where read has a variant (different from the reference)?
Hi, I'm trying to write some code that generates both the genome location and read ID when there is a variant read(insertion or deletion) to the reference genome. My intution tells me to use samtools view but since i'm pretty new to the pysam library the documentation doesn't seem to indicate this. Any insight is appreciated!
• 583 views
•
link
0 answers
No answers yet.
Log in to answer this question.