Difference between revisions of "2017 Haneul Lab note"

From Crop Genomics Lab.
Jump to: navigation, search
(1 / 10)
(1 / 9)
Line 4: Line 4:
 
  using lgcombine.py(63:/data/skyts0401/Mungbean/MY_UV/)
 
  using lgcombine.py(63:/data/skyts0401/Mungbean/MY_UV/)
 
  find linkage group 2 map is wrong, construct map newly
 
  find linkage group 2 map is wrong, construct map newly
 +
  
 
=== Mungbean synchronous QTL ===
 
=== Mungbean synchronous QTL ===
Line 9: Line 10:
 
  missing > 10%, hetero > 10, depth < 3 marker is filtered
 
  missing > 10%, hetero > 10, depth < 3 marker is filtered
 
  while grouping them, find vr03, vr04 is combined in a group and vr05 is splited 2 groups, check it.
 
  while grouping them, find vr03, vr04 is combined in a group and vr05 is splited 2 groups, check it.
 +
  
 
=== Mungbean pacbio assembly ===
 
=== Mungbean pacbio assembly ===

Revision as of 06:00, 23 January 2017

Contents

1 / 9

Minyoung_UV_QTL

parsing genotype data(Joinmap) and phenotype data to ICImapping format(bip format)
using lgcombine.py(63:/data/skyts0401/Mungbean/MY_UV/)
find linkage group 2 map is wrong, construct map newly


Mungbean synchronous QTL

make loc file(244:/home/skyts0401/reseq/chr/Mungbean_chr_coseq_parse_seg_dist.loc)
missing > 10%, hetero > 10, depth < 3 marker is filtered
while grouping them, find vr03, vr04 is combined in a group and vr05 is splited 2 groups, check it.


Mungbean pacbio assembly

moving SAM data(align reseq data on pacbio-scaffold) from NICEM server to 244 server (244:/kev8305/SK3/)

1 / 10

Minyoung_UV_QTL

QTL analysis by using IciMapping


Mungbean synchronous QTL

construct genetic map (JoinMap 4.1), just using chr 3, 4 combined and chr 5 splited linkage group.

ML method, Haldane algorithm


Mungbean pacbio assembly

convert SAM format to BAM format (244:/kev8305/SK3/)

./convertbam.sh

1/ 11

Mungbean synchronous QTL

QTL analysis by using RQTL(desktop:/Users/sky/desktop/Mungbean_syn_RQTL.csv)
just for checking locus

1/16

Mungbean pacbio assembly

coping sorted.bam file from 244 server to 63 server

variant calling (244:/kev8305/SK3/, 63:/data/skyts0401/Mungbean/mapping/resequencing/)

samtools mpileup -f falcon_500_sspace.final.scaffolds.fasta -v -t DP,AD,ADF,ADR,SP,INFO/AD,INFO/ADF,INFO/ADR -b bam_list | bcftools call -v -m -O v > variants.vcf
samtools mpileup -f falcon_500_sspace.final.scaffolds.fasta -I -v -t DP,AD,ADF,ADR,SP,INFO/AD,INFO/ADF,INFO/ADR -b bam_list | bcftools call -v -m -O v > variants_snp.vcf

1/18

Mungbean pacbio assembly

variant calling Kyoungki Jarae #5 with pacbio falcon scaffold (63:/data/skyts0401/Mungbean/mapping/resequencing/)

bwa index falcon_500_sspace.final.scaffolds.fasta
bwa mem -t 10 falcon_500_sspace.final.scaffolds.fasta KJ-C_1.fastq.gz KJ-C_2.fastq.gz > KJ-pe_falcon_scaffold.sam

1/19

Mungbean pacbio assembly

variant calling Kyoungki Jarae #5 with pacbio falcon scaffold (63:/data/skyts0401/Mungbean/mapping/resequencing/)

samtools view -Sb KJ-pe_falcon_scaffold.sam > KJ-pe_falcon_scaffold.bam
samtools sort KJ-pe_falcon_scaffold.bam -o KJ-pe_falcon_scaffold.sorted.bam
samtools index KJ-pe_falcon_scaffold.sorted.bam
samtools mpileup -f falcon_500_sspace.final.scaffolds.fasta -I -v -t DP,AD,ADF,ADR,SP,INFO/AD,INFO/ADF,INFO/ADR -u KJ-pe_falcon_scaffold.sorted.bam | bcftools call -v -m -O v > KJ_falcon_scaffold_variants_snp.vcf