Chris:LabNotes/Microbiome/2011/2011-5-2: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Cjwei
(Created page with '=Slide Shortening Assembly Test= ==Background== *Sam tried to increase SAG contig length by using recruited sequences from metagenomic contigs (http://genome-tech.ucsd.edu/LabNot…')
 
>Cjwei
Line 10: Line 10:
**Test found in: </media/disk-2/cjwei/slide-shortening_Assembly_Test>
**Test found in: </media/disk-2/cjwei/slide-shortening_Assembly_Test>
**Used various sliding cut lengths:
**Used various sliding cut lengths:
**Test 1: >1000bp contigs cut; step size = 200bp
***Test 1: >1000bp contigs cut; step size = 200bp
**Test 2: >2000bp contigs cut; step size = 100bp (same as "Whole Genome Amplification and Denovo Assembly of Single Bacterial Cells")
***Test 2: >2000bp contigs cut; step size = 100bp (same as "Whole Genome Amplification and Denovo Assembly of Single Bacterial Cells")
***Test 3:


==Procedure==
==Procedure==

Revision as of 21:11, 2 May 2011

Slide Shortening Assembly Test

Background

  • Sam tried to increase SAG contig length by using recruited sequences from metagenomic contigs (http://genome-tech.ucsd.edu/LabNotes/index.php/Sam:LabNotes/Microbiome-new/2011-4-19)
    • But results showed some decreased contig lengths after recruitment (ID2, ID9, ID11)
    • Dr. Zhang suggested that the decreased contig lengths may be fixed by slide shortening of contigs
  • On 4/29, I wrote perl script that could help perform the slide shortening of contigs
  • Here, I test the script/method of slide shortening to increase the contig length after recruitment
  • SAG sequence: ID9 because showed the lower contig length but showed relatively good contig length
    • Sequences found in: </media/disk-1/samchiang/Sub_tests_folder_imp/2_Extending_Mouse_Bac_SAG_contigs/BLAT_Meta2SAG_040311/Step4_pooled_together/SAG_A_ID9>
    • Test found in: </media/disk-2/cjwei/slide-shortening_Assembly_Test>
    • Used various sliding cut lengths:
      • Test 1: >1000bp contigs cut; step size = 200bp
      • Test 2: >2000bp contigs cut; step size = 100bp (same as "Whole Genome Amplification and Denovo Assembly of Single Bacterial Cells")
      • Test 3:

Procedure

  • Use Sliding Cut for combined SAG and recruited metagenome contigs:
perl contig_sliding_cut.pl SAG_A_ID9_and_recruted_Meta_all >SAG_A_ID9_and_recruted_Meta_all_slide_cut
  • Run Velveth
    • k-mer 27, -short
/media/disk-2/cjwei/Velvet/velveth SAG_A_ID9_and_recruted_Meta_all_slide_cut_Assembly/ 27 -fasta -short SAG_A_ID9_and_recruted_Meta_all_slide_cut
  • Run Velvetg
    • fixed setting as cov_cutoff 2, min-contig lgth 100
/media/disk-2/cjwei/Velvet/velvetg SAG_A_ID9_and_recruted_Meta_all_slide_cut_Assembly/ -cov_cutoff 2 -min_contig_lgth 100

Results

n       n:100   n:N50  min    median mean   N50    max    sum
7871    7871    737    100    196    465    1000   20981  3662057 contigs_1000_200.fa
819     7819    710    100    197    473    1057   21081  3703868 contigs_2000_100.fa
7861    7861    726    100    196    469    1039   21081  3689257 contigs_1000_100.fa
7869    7869    718    100    196    468    1041   21091  3690448 contigs_100_10.fa
14828   14828   1140   100    170    403    902    38317  5988190 ../SAG_A_ID9_and_recruted_Meta_all
  • Based on results, notice that