Difference between revisions of "SK met"

From Crop Genomics Lab.
Jump to: navigation, search
Line 1: Line 1:
 +
== materials ==
 
  raw_data : /NGS/NGS/VignaRadiata/DNA/SK/
 
  raw_data : /NGS/NGS/VignaRadiata/DNA/SK/
TBI[7]=TN1509D0657--TCCGGAGA-AGGCTATA
+
/mount_test/sdd1/Mungbean_assemble/ver6/SNP/Vradi_ver6.reseqKJ.sorted.bam
sample[7]=SK010
+
{| cellpadding="20" cellspacing="0" border="1"
 
+
|SK010
TBI[57]=TN1509D0710--TAATGCGC-TAAGATTA
+
|TN1509D0657--TCCGGAGA-AGGCTATA
sample[57]=SK068
+
|-
 
+
|SK068
TBI[177]=TN1510D0932--TCTCGCGC-TCAGAGCC
+
|TN1509D0710--TAATGCGC-TAAGATTA
sample[177]=SK186
+
|-
 
+
|SK186
TBI[43]=TN1509D0695--GAATTCGT-ACGTCCTG
+
|TN1510D0932--TCTCGCGC-TCAGAGCC
sample[43]=SK049
+
|-
 
+
|SK049
TBI[83]=TN1509D0736--TCTCGCGC-GTCAGTAC
+
|TN1509D0695--GAATTCGT-ACGTCCTG
sample[83]=SK094
+
|-
 
+
|SK094
TBI[143]=TN1510D0898--CTGAAGCT-GCCTCTAT
+
|TN1509D0736--TCTCGCGC-GTCAGTAC
sample[143]=SK152
+
|-
 
+
|SK152
TBI[167]=TN1510D0922--TCCGCGAA-GCCTCTAT
+
|TN1510D0898--CTGAAGCT-GCCTCTAT  
sample[167]=SK176
+
|-
 
+
|SK176
TBI[147]=TN1510D0902--CTGAAGCT-TAAGATTA
+
|TN1510D0922--TCCGCGAA-GCCTCTAT
sample[147]=SK156
+
|-
 
+
|SK156
244 > 63 /data/haggui/raw_data/ 에 sample#.fq.gz 로 저장함
+
|TN1510D0902--CTGAAGCT-TAAGATTA
 +
|}
 +
*244 > 63 /data/haggui/raw_data/ 에 sample#.fq.gz 로 저장함
  
 
  ref : 244 /mount_test/sdd1/Mungbean_assemble/ver6/SNP/Vradi_ver6.fa  
 
  ref : 244 /mount_test/sdd1/Mungbean_assemble/ver6/SNP/Vradi_ver6.fa  
244 > 63 data/haggui/Vradi_ver6.fa
+
*244 > 63 data/haggui/Vradi_ver6.fa
1. bwa index <ref.fa>
+
# bwa index <ref.fa>
2. cat SK_met.list | parallel --gnu -j 3 "bwa mem Vradi_ver6.fa ./raw_data/{}_1.fq.gz ./raw_data/{}_2.fq.gz | samtools view -bS -q 30 -> {}.bam"
+
# cat SK_met.list | parallel --gnu -j 3 "bwa mem Vradi_ver6.fa ./raw_data/{}_1.fq.gz ./raw_data/{}_2.fq.gz | samtools view -bS -q 30 -> {}.bam"

Revision as of 05:44, 23 January 2019

materials

raw_data : /NGS/NGS/VignaRadiata/DNA/SK/
/mount_test/sdd1/Mungbean_assemble/ver6/SNP/Vradi_ver6.reseqKJ.sorted.bam
SK010 TN1509D0657--TCCGGAGA-AGGCTATA
SK068 TN1509D0710--TAATGCGC-TAAGATTA
SK186 TN1510D0932--TCTCGCGC-TCAGAGCC
SK049 TN1509D0695--GAATTCGT-ACGTCCTG
SK094 TN1509D0736--TCTCGCGC-GTCAGTAC
SK152 TN1510D0898--CTGAAGCT-GCCTCTAT
SK176 TN1510D0922--TCCGCGAA-GCCTCTAT
SK156 TN1510D0902--CTGAAGCT-TAAGATTA
  • 244 > 63 /data/haggui/raw_data/ 에 sample#.fq.gz 로 저장함
ref : 244 /mount_test/sdd1/Mungbean_assemble/ver6/SNP/Vradi_ver6.fa 
  • 244 > 63 data/haggui/Vradi_ver6.fa
  1. bwa index <ref.fa>
  2. cat SK_met.list | parallel --gnu -j 3 "bwa mem Vradi_ver6.fa ./raw_data/{}_1.fq.gz ./raw_data/{}_2.fq.gz | samtools view -bS -q 30 -> {}.bam"