Different result by bowtie2
0
0
Entering edit mode
7 hours ago
daffodil ▴ 10

I ran the following Bowtie2 command, but my results show fewer uniquely mapped reads compared to my supervisor's results, who obtained 22,000,000 uniquely mapped reads. The command I used is:

bowtie2 -x /sw/data/iGenomes/Mus_musculus/UCSC/mm10/Sequence/Bowtie2Index/genome \
  -U cut_trimmed.fastq.gz \
  --no-unal \
  -X 700 \
  -S cut.sam

I have also used cutadapt for trimming of single read. What could be the reason for this discrepancy in the number of uniquely mapped reads?

Bowtie2 • 76 views
ADD COMMENT
0
Entering edit mode

There are many reasons you might get different number of reads with a tool like Bowtie2. The tool uses heuristics to be fast and efficient, so even two runs with most of the same parameters will differ.

How large is the difference between your two datasets? Are all the other parameters identical in your runs?

ADD REPLY

Login before adding your answer.

Traffic: 1881 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6