Question: How to remove specifically non aligned regions
0
gravatar for curiousbiologist
2.4 years ago by
France
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 • 665 views
ADD COMMENTlink modified 2.4 years ago by Brice Sarver2.5k • written 2.4 years ago by curiousbiologist40
1
gravatar for Brice Sarver
2.4 years ago by
Brice Sarver2.5k
United States
Brice Sarver2.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 2.4 years ago by Brice Sarver2.5k
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: 1183 users visited in the last hour