Hi!
First of all thanks to the community who often solve my problems before I even have time to ask the question.
I have a problem though. I’m working with datasets from 454 sequencing and I know one organism to be present and want to filter out those reads to find possible other organisms by analyzing the remaining reads.
From what I have managed to grasp running around in these forums and other forums it should be an easy task if I could extract the data from an alignment. MAQ as well as Bowtie have these functions built into them (to dump the unaligned sequences) but since I’m using 454 they won’t work.
So my question is, is this doable with BWA and samtools? I think it should be doable but I’m left with mostly incomplete clues in the forums and I’m a bit lost right now.
I know that samtools view with the flag –f should be able to dump the reads but I’m unsure of the flag and even more unsure on how to handle the output (ideally I will have an output that I can incorporate into MIRA for assembly (or velvet / Ray)).
First of all thanks to the community who often solve my problems before I even have time to ask the question.
I have a problem though. I’m working with datasets from 454 sequencing and I know one organism to be present and want to filter out those reads to find possible other organisms by analyzing the remaining reads.
From what I have managed to grasp running around in these forums and other forums it should be an easy task if I could extract the data from an alignment. MAQ as well as Bowtie have these functions built into them (to dump the unaligned sequences) but since I’m using 454 they won’t work.
So my question is, is this doable with BWA and samtools? I think it should be doable but I’m left with mostly incomplete clues in the forums and I’m a bit lost right now.
I know that samtools view with the flag –f should be able to dump the reads but I’m unsure of the flag and even more unsure on how to handle the output (ideally I will have an output that I can incorporate into MIRA for assembly (or velvet / Ray)).
Comment