Entering edit mode
8.5 years ago
TheMaster
•
0
I'd like to understand how can I pick the amount of reads for a specific sequence.
E.g. I have a sequence of 200 bp, and need to know all reads from pair 6 to pair 70, or detect where the half of all reads (!=the middle of coordinates) for the sequence given.
I tried to read https://samtools.github.io/htsjdk/javadoc/htsjdk/ and picard documentation, but with zero result. Please, help me with an advice, I'm a total newcomer.
Thank you in advance!
Please consider editing your question to include more details. Could you provide an example of an input file and your desired output? Perhaps you can elaborate about your specific goal?