Question: How to remove specifically non aligned regions
gravatar for curiousbiologist
4.0 years ago by
curiousbiologist40 wrote:

Hi all,

I would like to select a specific region from a fastq file (all read matching this region) and if they exceed the reference, trim only the non matching regions. I've done a bwa alignment and generated sam file but I'm not sure how remove the regions that are not matching to the reference without removing all the read

Thanks for your help!

sequencing alignment genome • 999 views
ADD COMMENTlink modified 4.0 years ago by Brice Sarver3.5k • written 4.0 years ago by curiousbiologist40
gravatar for Brice Sarver
4.0 years ago by
Brice Sarver3.5k
United States
Brice Sarver3.5k wrote:

You'll want to subset the BAM based on the region you want. This can be done easily with samtools.

Check out this Biostars post: Extract Reads From A Bam File That Fall Within A Given Region

ADD COMMENTlink written 4.0 years ago by Brice Sarver3.5k
Please log in to add an answer.


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