Created
August 3, 2011 18:44
-
-
Save brantfaircloth/1123452 to your computer and use it in GitHub Desktop.
Getting singleton reads from a 454 assembly:
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # get singleton reads names from assembly stats: | |
| grep Singleton 454ReadStatus.txt > singles.txt | |
| # if you assembled using fasta+qual files, use faSomeRecords from kent source: | |
| faSomeRecords ../reads/schisto.454reads.combined.fsa singles.txt singles.fa | |
| # if you assembled using SFF files: | |
| sfffile -o singles.sff -i singles.txt sff/* | |
| sffinfo -s singles.sff > singles.fna |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment