Does anyone know how I can pull just one read from a BAM file, when I know the read-id (and sequence) of that read? I'm doing this because I've just got an error saying that one of the reads has no quality information. Since the original BAM is huge, I don't want to convert the whole thing to SAM or fastq and find that read-id. Any suggestions (also on BAM validation) welcome!
exact duplicate of
Extracting subsets of reads from a BAM file
quickest way to pull a list of reads out of a BAM file?
Sorry, should have spotted that