[BioC] deqseq_count and BWA-based SAM files
Steve Lianoglou
mailinglist.honeypot at gmail.com
Tue Dec 13 22:02:01 CET 2011
I'd rather use `samtools sort ...`
You probably want to (1) convert your sam to bam, (2) then sort, (3)
then index, (4) then use.
>From the command line:
## convert to bam
samtools view -bS your.sam > your-tmp.bam
## sort bam
samtools sort your-tmp.bam your
## index bam
samtools index your.bam
Now you can junk your sam file and save mucho HD space to boot.
-steve
On Tue, Dec 13, 2011 at 10:41 AM, wang peter <wng.peter at gmail.com> wrote:
> you can sort them by linux command
>
> sort -T . -S 2G -k 1,1 -k 3,3 $target_fa.psl.sam > $target_fa.psl.nameSorted.sam
> --
> shan gao
> Room 231(Dr.Fei lab)
> Boyce Thompson Institute
> Cornell University
> Tower Road, Ithaca, NY 14853-1801
> Office phone: 1-607-254-1267(day)
> Official email:sg839 at cornell.edu
> Facebook:http://www.facebook.com/profile.php?id=100001986532253
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
--
Steve Lianoglou
Graduate Student: Computational Systems Biology
| Memorial Sloan-Kettering Cancer Center
| Weill Medical College of Cornell University
Contact Info: http://cbio.mskcc.org/~lianos/contact
More information about the Bioconductor
mailing list