User Tools

Site Tools


scaffolding

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
scaffolding [2016/02/16 14:24] – [개요] hyjeongscaffolding [2021/03/17 13:09] (current) – external edit 127.0.0.1
Line 14: Line 14:
   $ SolexaQA++ dynamictrim SB_iontor.fastq –-torrent (default: P = 0.05)   $ SolexaQA++ dynamictrim SB_iontor.fastq –-torrent (default: P = 0.05)
   $ sed -n '/^@BYI6V.*\/1$/{N;N;N;p;}' SB_iontor.fastq.trimmed > SB-trimmed_1.fastq   $ sed -n '/^@BYI6V.*\/1$/{N;N;N;p;}' SB_iontor.fastq.trimmed > SB-trimmed_1.fastq
-  $ sed -n '/^@BYI6V.*\/1$/{N;N;N;p;}' SB_iontor.fastq.trimmed > SB-trimmed_1.fastq+  $ sed -n '/^@BYI6V.*\/2$/{N;N;N;p;}' SB_iontor.fastq.trimmed > SB-trimmed_2.fastq
   [optional] SolexaQA++ lengthsort SB-trimmed_1.fastq SB-trimmed_2.fastq -l 40   [optional] SolexaQA++ lengthsort SB-trimmed_1.fastq SB-trimmed_2.fastq -l 40
  
Line 30: Line 30:
  
 ===== SSPACE를 이용한 scaffolding ===== ===== SSPACE를 이용한 scaffolding =====
-A5-miseq에 포함된 구 버전의 SSPAE(v1-1)는 paired end(-> <-, 0)와 mate-pair(<- ->, 1)만을 지원하므로, Ion Torrent에서 만들어진 mate-pair library read(RR, <- <-)을 사용하기 위해서는 read file 중 하나를 reverse complementary로 전환해야 +A5-miseq에 포함된 구 버전의 SSPAE(v1-1)는 paired end(-> <-, 0)와 mate-pair(<- ->, 1)만을 지원하므로, Ion Torrent에서 만들어진 mate-pair library read(RR, <- <-)를 사용하기 위해서는 read file 중 하나를 reverse complementary로 전환해야 한다. 또한 라이브러리 정보 파일의 형식도 다르다.
  
   $ seqtk seq -r SB-trimmed_1.fastq > 1rc.fastq   $ seqtk seq -r SB-trimmed_1.fastq > 1rc.fastq
Line 43: Line 43:
   $cat library2.txt   $cat library2.txt
   PGM_mate-pair-library 1rc.fastq SB-trimmed_2.fastq 3000 0.5 0   PGM_mate-pair-library 1rc.fastq SB-trimmed_2.fastq 3000 0.5 0
 +
 +==== SSPACE v3.0의 library file 사례====
 +  Lib1 bwa file1.1.fasta file1.2.fasta 400 0.25 FR
 +  Lib1 bowtie file2.1.fasta file2.2.fasta 400 0.25 FR
 +  Lib2 bwasw file3.1.fastq file3.2.fastq 4000 0.5 RF
 +  Lib2 TAB file4.tab 4000 0.5 RF
 +  Lib3 TAB file5.tab 10000 0.5 RF
 +  unpaired bowtie unpaired_reads1.fasta
 +  unpaired bwasw unpaired_longreads1.gz
  
  
scaffolding.1455600259.txt.gz · Last modified: 2021/03/17 13:09 (external edit)