Hello.
I am mapping 2x100 Illumina reads against hg19 using Bowtie using a tiered approach. I want separate output files for reads that map to 1 location and reads that map to 2 locations.
In Tier 1, I set m=1 and use --max to save the reads that mapped to more than 1 location. (eg: m 1 --max max.fastq)
In Tier 2, I set m=2 and map the reads from the max.fastq file. Usually, the number of unaligned reads in this step is 0 (as would be expected). Sometimes however there are unaligned reads. Does anyone know why that is?
Thanks.
I am mapping 2x100 Illumina reads against hg19 using Bowtie using a tiered approach. I want separate output files for reads that map to 1 location and reads that map to 2 locations.
In Tier 1, I set m=1 and use --max to save the reads that mapped to more than 1 location. (eg: m 1 --max max.fastq)
In Tier 2, I set m=2 and map the reads from the max.fastq file. Usually, the number of unaligned reads in this step is 0 (as would be expected). Sometimes however there are unaligned reads. Does anyone know why that is?
Thanks.
