Question: BLASR: How to subsample alignments from a region
0
gravatar for shwethacm
5.1 years ago by
shwethacm220
Seattle, WA
shwethacm220 wrote:

Hello! 

I have an alignment of ~20x PacBio long reads to my reference genome and I am trying to extract alignments within only a certain region, say 2mb-10mb (sort of like samtools view ) but I am not sure how to do this.

I have alignments in the m4 format.

I think I should be specifying a region table but I am not sure.

Can anybody help?  

ADD COMMENTlink modified 5.1 years ago • written 5.1 years ago by shwethacm220

Thanks, rhall.

I did eventually parse it myself but I was looking to check if I could be lazy about it and do it quickly using a preexisting tool :)

Well, now I know!

Thanks a ton anyway.

Shwetha

ADD REPLYlink modified 13 months ago by _r_am32k • written 5.1 years ago by shwethacm220
0
gravatar for rhall
5.1 years ago by
rhall160
United States
rhall160 wrote:

I'm not aware of any software to process m4 files, but they are straightforward to parse with simple scripts, https://github.com/PacificBiosciences/blasr/wiki/Blasr-Output-Format. If you are more comfortable with sam format, blasr can output sam instead of m4, the -sam parameter.

ADD COMMENTlink modified 13 months ago by _r_am32k • written 5.1 years ago by rhall160
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1663 users visited in the last hour
_