How can I generate sequences from a bam/sam file starting at a specific position and remove what is before this position? Thank you!
2 answers
I'm deeply sorry if my question was obscure: what I want to do is to get, from an alignment, reads without nucleic bases before position 30 of the reference sequence. Is it possible to crop a bam file?
Thank you for your answer! However, I got reads matching within the interval given, not cropped reads inside this interval
Thank you for your answer, it means to me to have a return. My programming skills are not great, do you know a script that I can use as start basis to do my custom script?
Could you use an igv screenshot and a graphical program (e.g. MS paint) to clarify what you aim for?
curiousbiologist wants individual reads chopped so they start and end at a specific position i.e. nothing should extend to left or right of an interval a <--> b
Okay, makes me wonder "why" OP would want that, but fine. This is not a straightforward question, requires modification of CIGAR, sequence, qualities, start,...
Yes you got it genomax2. I want to have several (and switchable) windows of reads from different samples in order to compare them using different score calculations; stats, entropy (shannon entropy score). If I haven't same size of reads pieces, my results will be misrepresented
would it be possible to resolve this problem using an awk script? I was thinking of an alignment, conversion to fasta (with gap or x for non-aligned bases) and then trimming using awk or fastx_trimmer? Maybe there is something easier? how do I get gap or 'x' for non-aligned bases before and after each reads?
Log in to answer this question.
-
Extarct specific bases from BAM/SAM files
written by Ankit
-
Extract base information from aligned reads in a sam file
written by yliueagle
-
ref and alt allele count from BAM file
written by monaallouba •
-
Size different between two uBam file which generate from the same Fastq.
written by weichi •
-
How do we find the number of sequences with terminal mismatches?
written by MAPK
-
BAM/SAM to aligned FASTA file
written by ryan.englander •
-
alignment containing sequences from position a to b
written by curiousbiologist •
-
How to remove specifically non aligned regions
written by curiousbiologist •
-
How to efficiently remove a list of reads from BAM file?
written by Tao
-
How Can I Find Position Specific Weight Matrix From A Multi Fasta File
written by k.nirmalraman
To me, it's unclear what you are asking for. Is this the same as alignment containing sequences from position a to b ? What do you want to obtain? "Sequences"? Is that a read/fasta/fastq/reference/variant...?
I would like to obtain a bam "cropped" (all my reads aligned and starting at a defined position)
One of the options in this threads should do this: How to get the consensus sequence from a BAM alignment
I would like to play with all the selected reads after, I'm not sure consensus is a good approach