SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Circular contig in multicontig reference remapping (SMRTPipe, BLASR) Gestahl Pacific Biosciences 1 08-30-2016 07:10 AM
Remapping MurielGB Genomic Resequencing 4 10-18-2013 07:33 AM
Remapping 1000 Genomes Data with BWA nexgengirl Bioinformatics 4 12-13-2012 03:09 PM
duplicate reads in Illumina short, single end reads of RNAseq data inbarpl Bioinformatics 4 05-22-2012 08:36 AM
Extending contigs by remapping AdrianP Bioinformatics 8 03-07-2012 01:23 PM

Reply
 
Thread Tools
Old 01-22-2018, 06:03 AM   #1
ekageyama
Junior Member
 
Location: Tuebingen

Join Date: Sep 2012
Posts: 3
Default Remapping reads

hello!
Currently I am working with a large set of STAR mapped reads, the output is in sam format. I would like to remap, but was wondering if its possible to do this using the SAM files without converting the sam files to fastq and then mapping. Any help would be appreciated.
ekageyama is offline   Reply With Quote
Old 01-22-2018, 06:12 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,780
Default

You could use unix pipe (e.g. reformat.sh in=your.sam out=stdout.fq | some_other_aligner commands sequence input from stdin ). reformat.sh comes from BBMap suite. If you have paired-end reads then things may become more complicated.

Any particular reason you want to do this? Are you not happy with STAR's output?
GenoMax is offline   Reply With Quote
Old 01-22-2018, 06:38 AM   #3
ekageyama
Junior Member
 
Location: Tuebingen

Join Date: Sep 2012
Posts: 3
Default

Actually, I would prefer to stay with STAR, the reason is I need to map against a different version of the genome, and the fastq files were deleted for storage reasons, so I would like to avoid converting each sam to fastq all together, since im working with a fairly large dataset.
ekageyama is offline   Reply With Quote
Old 01-22-2018, 06:41 AM   #4
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,780
Default

Then you could use STAR once again in place of "some_other_aligner" above.
GenoMax is offline   Reply With Quote
Old 02-07-2018, 09:29 PM   #5
Dario1984
Senior Member
 
Location: Sydney, Australia

Join Date: Jun 2011
Posts: 163
Default

If the new version of the genome is of higher quality and more complete, you would have lost all of those reads which did not map to the previously used genome but could map to the new genome, because the default value of --outSAMunmapped is None. I don't think that mapping again is worthwhile if you can't access all of the reads which were stored in the FASTQ files.
Dario1984 is offline   Reply With Quote
Reply

Tags
alignment, star, star aligner

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:35 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO