Difference between revisions of "2017 Haneul Lab note"
From Crop Genomics Lab.
(→1/18) |
(→Mungbean pacbio assembly) |
||
Line 33: | Line 33: | ||
coping sorted.bam file from 244 server to 63 server | coping sorted.bam file from 244 server to 63 server | ||
− | variant calling (244, 63) | + | variant calling (244, 63) (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 -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 == | == 1/18 == | ||
=== Mungbean pacbio assembly === | === Mungbean pacbio assembly === | ||
variant calling Kyoungki Jarae #5 with pacbio falcon scaffold | variant calling Kyoungki Jarae #5 with pacbio falcon scaffold |
Revision as of 02:21, 18 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, 63) (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