Hi Everyone,
I tried Tophat v 2.0.6 and GSNAP v 2013-02-05 to align Illumina paired-end reads. However, the numbers of proper-paired reads were very different. I also used cufflinks for assembly, but got similar number of genes/isoforms. Is there anything wrong with my commands? Why did Tophat result in such low proper-paired read rates? Thanks!
To get proper-paired reads:
samtools view -f 0x2 accepted_hits.bam | cut -f1 | sort | uniq | wc -l
Tophat comands:
>> tophat -G genes.gtf -o tophat --no-novel-juncs genome read_1 read_2
Results: 1,431,500 (6.34%) proper-paired reads; 3746 genes (cufflinks).
GSNAP command
>> gsnap -A sam -N 0 -D ~/Software/gmap-2013-02-05/ -d mm10 -s mm10.splicesites.iit read_1 read_2
Results: 20,543,759 (90.9%) proper-paired reads;3959 genes (cufflinks).
-
Tophat/GSNAP: proper-paired reads
Last edited by riziai; 04-26-2013, 12:40 PM.
-
Spatial biology is an exciting field that encompasses a wide range of techniques and technologies aimed at mapping the organization and interactions of various biomolecules in their native environments. As this area of research progresses, new tools and methodologies are being introduced, accompanied by efforts to establish benchmarking standards and drive technological innovation.
3D Genomics
While spatial biology often involves studying proteins and RNAs in their...01-01-2025, 07:30 PM