Hello everyone,
We have troubles with mapping color-space SOLiD reads using bowtie. So I'm looking for your help to solve this weird issue.
We have yeast mRNA reads from SOLiD 4. We mapped these reads with bowtie using index from the mRNA sequences. About 30% of the reads mapped to these mRNAs. We checked with ABI rep, they said that 30% is normal for these libraries. So we went on with our analysis.
Now that we have some spare time. We want to know what the unmapped reads are, so we first mapped them to the yeast genome.
Surprisingly, ~50% of the unmapped reads mapped to the genome (presumably the non-mRNA regions). When we looked at the location of these genome-mapped reads, 10% are located in mRNA region in the yeast genome.
We went back and forth to check and it seems these genome-mapped are real and actually mapping to mRNAs.
Our question is how bowtie seems to miss these genome-mapped reads when mapping to the mRNAs. We also wonder how many of the remaining unmapped reads are actually unmapped, since it's not a viable option to convert from color-space to base-space and then mapping them.
We appreciate any advice on why this happens.
Thanks,
Hoa
We have troubles with mapping color-space SOLiD reads using bowtie. So I'm looking for your help to solve this weird issue.
We have yeast mRNA reads from SOLiD 4. We mapped these reads with bowtie using index from the mRNA sequences. About 30% of the reads mapped to these mRNAs. We checked with ABI rep, they said that 30% is normal for these libraries. So we went on with our analysis.
Now that we have some spare time. We want to know what the unmapped reads are, so we first mapped them to the yeast genome.
Surprisingly, ~50% of the unmapped reads mapped to the genome (presumably the non-mRNA regions). When we looked at the location of these genome-mapped reads, 10% are located in mRNA region in the yeast genome.
We went back and forth to check and it seems these genome-mapped are real and actually mapping to mRNAs.
Our question is how bowtie seems to miss these genome-mapped reads when mapping to the mRNAs. We also wonder how many of the remaining unmapped reads are actually unmapped, since it's not a viable option to convert from color-space to base-space and then mapping them.
We appreciate any advice on why this happens.
Thanks,
Hoa
Comment