Extracting Reads from paired end data in SAM file that map to a list of contig IDs
1 answer
ah, I see. Using awk, matching chromosomes RF01 or RF02
/samtools view -h S1.bam | awk -F '\t' 'function fun(C) { return C=="RF01" || C=="RF02";} /^@/ {print;next;} {if(fun($3) || fun($7)) print;}'
or using samjdk: http://lindenb.github.io/jvarkit/SamJdk.html
java -jar dist/samjdk.jar -e 'Predicate<String> f=C->C.equals("RF01") || C.equals("RF02"); return (!record.getReadUnmappedFlag() && f.test(record.getReferenceName())) || (!record.getMateUnmappedFlag() && f.test(record.getMateReferenceName()));' in.bam
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
Extract reads from BAM file with known variant
written by pablo
-
extract ALL mapped paired-end reads from BAM file
written by romain.derelle •
-
Contig extension using PRICE
written by A_heath •
-
Split BAM file from paired end read into forward and reverse strand w.r.t reference genome
written by Susmita Mandal
-
Extract reads matching variants from VCF file
written by sppearce •
-
Trying to get consensus sequence from paired end alignment
written by Sus •
-
Extract reads from specific target in .bam file
written by Picasa
-
Extracting reads that one pair is aligned and the other one is not
written by m.ghorbani1991 •
-
id in text file for sequence retrival from multifasta
written by tcf.hcdg
-
Extract Alignment By Read Id From A Sam File
written by Nicolas Rosewick
many duplicates. e.g: Extract ONLY chromosomes 1-22 from bam file - removing extraneous chr annotations ...
would the solution in that link also work for paired data? Thanks!
@max_19 why wouldn't it work for paired data ?
If one read in a pair matches the contig but the other read in that pair does not. Is there a way to include that pair even if only one read matches?