How to remove specifically non aligned regions
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!
• 1,985 views
•
link
updated
by
Brice Sarver
•
written
by
curiousbiologist •
•
written
by
curiousbiologist •
1 answer
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
• 0 views
•
link
written
by
Brice Sarver
Log in to answer this question.
More posts like this
-
Problem trimming SAM file to use with assembler
written by Emilio •
-
Extract base information from aligned reads in a sam file
written by yliueagle
-
BWA MEM output directory
written by sarastrafella.ss •
-
Skip orientation as there are not enough pairs by bwa mem
written by seta
-
Remove reads from fastq file
written by varsha619
-
alignment containing sequences from position a to b
written by curiousbiologist •
-
Mapping BAM/SAM file back to reference?!?!
written by ab.tsubaki •
-
BWA alignment issue in iPlant
written by aj123
-
removing sequencing gaps
written by cdwilliam524 •
-
Clarification For Running Bwa In Parallel
written by Mohamed •
