Hi!
I'm new to both RNA-seq and bioinformatics itself.
In the project that we're working on we sequenced bacterial RNA from different samples and one of the things we're trying to find from it are deletions and insertions resulting from mistakes in ligation of the RNA.
I have chosen to use bowtie2 for alignment as the argument about which aligner is better seems to be never-ending and also dependant on the situation. But I have seen that people have noted that bowtie2 might prefer a gapped alignment to another non-gapped alignment in some situations.
We had a crazy idea of first using bowtie to map all reads and then filter out all the reads that bowtie did not map(as bowtie doesn't map reads with indels) and maps them with bowtie2. It seems like this way I could calm myself about getting false positives due to alignment.
Is this a bad idea? I haven't seen anyone propose it, so I suppose there's something wrong with it.
Thank you for your answer.
I forgot to write it in the OP, but I planned on allowing none to one mismatch while mapping with bowtie. In this case it should only map the reads that are perfectly mappable somewhere on the genome and leave the rest unmapped for bowtie2 to find gapped alignments. I would only filter out the reads that have a perfect alignment somewhere on the genome just in case bowtie2 doesn't find it.
Ah, allowing a very low edit distance alleviates some of my concerns. BTW, here's a bwa mem-based pipeline that would presumably work and can perform realignment with soft-clipped alignments. I've never used it, but that might be a nicer alternative.