Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • arcolombo698
    replied
    Yes thank you... it is true that CNV analysis must have WGS or WES. thank you very much.

    Leave a comment:


  • pashadag
    replied
    CNVer was not designed to work with RNA seq data

    Leave a comment:


  • arcolombo698
    replied
    can CNVer be used on RNAseq (NGS short mated reads) ?

    Leave a comment:


  • pashadag
    replied
    You probably need to use an absolute filename inside of bam_files_to_use.txt

    Leave a comment:


  • dmacmillan
    replied
    Unfortunately not

    Leave a comment:


  • CNVboy
    replied
    Have you solved your problem?
    thx

    Originally posted by dmacmillan View Post
    I had a question regarding CNVer if anyone has some experience. I ran the cnver.pl with correct parameters, however, I get a ton of output on the screen and some of it doesn't looks good.

    Here is an example:

    Code:
    Exec at 12:44:30, Fri Jan 13, 2012 awk '{ if ($6 == 3) print $0 }'  /home/dmacmillan/copy_num_ana/CNVer/workdir/mapping_files/chr22.mmap  | sort -k4n,4 -k3n
    ,3 -k1n,1  | /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/cluster/cluster_matepairs --mean=398 --stdev=58 --colID=-1 --colDist=0 --colChr=1 --colLeft
    =2 --colRight=3 --colTemplate=4 --baseLenFactor=1 --mdJoinTolerance=6 --type=3 --concise=1 > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.t3...
    No matepairs to cluster.
    
    Exec at 12:44:30, Fri Jan 13, 2012 grep -h EDGE /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.t[0123] | /home/dmacmillan/copy_num_ana/CNVer/cnver
    -0.8.1/src/cluster/screen_contig_breaks --mean=398 --stdev=58 --tolerance=572 --breaksFile = /home/dmacmillan/copy_num_ana/CNVer/hg18comp/contig_breaks_fold
    er/chr22.cbs | awk '{ if (($5 >= 2) && ($8 < 10000000)) print $0 }' > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.links...
    
    Exec at 12:44:30, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/axt_to_binary /home/dmacmillan/copy_num_ana/CNVer/hg18comp/self_align
    ments_folder/chr22.axt 49691432 /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.edges...
    Axtbin
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/make_reference_graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr2
    2//chr22.edges /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.graph...
    Converged after 1 iterations
    num_edges: 0
    num_nodes: 0
    
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/fill_reference_graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr2
    2//chr22.graph 49691432...
    
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/make_donor_graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//c
    hr22.graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.links /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.scov /home/dmacmillan/cop
    y_num_ana/CNVer/workdir/chr22//chr22.masks /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.gc 1> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22/
    /chr22.problem 2> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.graphinfo...
    
    Exec at 12:44:31, Fri Jan 13, 2012 perl /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/flow_solve/make_cs2_graph.pl on coverage: 10 chr22.problem 2 0 2
    0 1 0 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/cs2-4.6/...
    Starting cs2 graph builder...
    0 read_arrival= nodecnt=1 totalarccnt=0
    Arc Count 0
    Starting cs2...
    
    line 1 of input - bad value of a parameter in the problem line
    Done with cs2.
    Done with cs2 graph builder...
    Running perl /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/flow_solve//affiner.pl 2  1 10 0 20 < outfile.txt > nf
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/report_cnvs /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.graphinfo /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.problem.out 100 0 2> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.raw 1> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.used_dgs...
    
    Exec at 12:44:31, Fri Jan 13, 2012 cat /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.raw | /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/intRemove2  /home/dmacmillan/copy_num_ana/CNVer/hg18comp/contig_breaks_folder/chr22.cbs > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.raw.screened...
    Read in 12 forbidden intervals.
    
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/smoother.sh 5 1000 chr22.cnvs.raw.screened > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.smoothed...
    
    Exec at 12:44:31, Fri Jan 13, 2012 cat /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.smoothed | awk -v OFS="\t" '{ if ($4 > 0) print $1, $2, $3, "gain"; else print $1, $2, $3, "loss" }' | /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/post_analysis/doc_walker /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.scov /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.gc | awk -v OFS="\t" '{ print $1, $2, $3, $4, $6 }' > /home/dmacmillan/copy_num_ana/CNVer/workdir/calls//chr22.cnvs...
    Here are the steps that I took to achieve this:
    Code:
    Steps:
    1. Untar the CNVer tar.gz source file
    2. cd into the directory
    3. type "./configure" (no quotations) and hit enter
    4. type "make" hit enter
    5. cd .. into folder containing cnver-0.8.1, hg18comp
    6. Created file "bam_files_to_use.txt" by doing
    	i. echo -e "file1.bam\nfile2.bam" > bam_files_to_use.txt
    7. type "export CNVER_FOLDER=/home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1" or wherever your cnver-0.8.1 folder is.
    8. ran the following command:
    ./cnver.pl --map_list /home/dmacmillan/copy_num_ana/CNVer/bams_to_use.txt --ref_folder /home/dmacmillan/copy_num_ana/CNVer/hg18comp --work_dir /home/dmacmillan/copy_num_ana/CNVer/workdir  --read_len 150 --mean_insert 398 --stdev_insert 58 --min_mps 2
    I do not receive a workdir/results folder in the end, and most, if not all the files in the chr[1-22]/x/y folders are 0 bytes. Any suggestions would be appreciated

    Leave a comment:


  • dmacmillan
    started a topic CNVer Help

    CNVer Help

    I had a question regarding CNVer if anyone has some experience. I ran the cnver.pl with correct parameters, however, I get a ton of output on the screen and some of it doesn't looks good.

    Here is an example:

    Code:
    Exec at 12:44:30, Fri Jan 13, 2012 awk '{ if ($6 == 3) print $0 }'  /home/dmacmillan/copy_num_ana/CNVer/workdir/mapping_files/chr22.mmap  | sort -k4n,4 -k3n
    ,3 -k1n,1  | /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/cluster/cluster_matepairs --mean=398 --stdev=58 --colID=-1 --colDist=0 --colChr=1 --colLeft
    =2 --colRight=3 --colTemplate=4 --baseLenFactor=1 --mdJoinTolerance=6 --type=3 --concise=1 > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.t3...
    No matepairs to cluster.
    
    Exec at 12:44:30, Fri Jan 13, 2012 grep -h EDGE /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.t[0123] | /home/dmacmillan/copy_num_ana/CNVer/cnver
    -0.8.1/src/cluster/screen_contig_breaks --mean=398 --stdev=58 --tolerance=572 --breaksFile = /home/dmacmillan/copy_num_ana/CNVer/hg18comp/contig_breaks_fold
    er/chr22.cbs | awk '{ if (($5 >= 2) && ($8 < 10000000)) print $0 }' > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.links...
    
    Exec at 12:44:30, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/axt_to_binary /home/dmacmillan/copy_num_ana/CNVer/hg18comp/self_align
    ments_folder/chr22.axt 49691432 /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.edges...
    Axtbin
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/make_reference_graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr2
    2//chr22.edges /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.graph...
    Converged after 1 iterations
    num_edges: 0
    num_nodes: 0
    
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/fill_reference_graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr2
    2//chr22.graph 49691432...
    
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/make_donor_graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//c
    hr22.graph /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.links /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.scov /home/dmacmillan/cop
    y_num_ana/CNVer/workdir/chr22//chr22.masks /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.gc 1> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22/
    /chr22.problem 2> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.graphinfo...
    
    Exec at 12:44:31, Fri Jan 13, 2012 perl /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/flow_solve/make_cs2_graph.pl on coverage: 10 chr22.problem 2 0 2
    0 1 0 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/cs2-4.6/...
    Starting cs2 graph builder...
    0 read_arrival= nodecnt=1 totalarccnt=0
    Arc Count 0
    Starting cs2...
    
    line 1 of input - bad value of a parameter in the problem line
    Done with cs2.
    Done with cs2 graph builder...
    Running perl /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/flow_solve//affiner.pl 2  1 10 0 20 < outfile.txt > nf
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/report_cnvs /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.graphinfo /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.problem.out 100 0 2> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.raw 1> /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.used_dgs...
    
    Exec at 12:44:31, Fri Jan 13, 2012 cat /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.raw | /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/intRemove2  /home/dmacmillan/copy_num_ana/CNVer/hg18comp/contig_breaks_folder/chr22.cbs > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.raw.screened...
    Read in 12 forbidden intervals.
    
    Exec at 12:44:31, Fri Jan 13, 2012 /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/smoother.sh 5 1000 chr22.cnvs.raw.screened > /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.smoothed...
    
    Exec at 12:44:31, Fri Jan 13, 2012 cat /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.cnvs.smoothed | awk -v OFS="\t" '{ if ($4 > 0) print $1, $2, $3, "gain"; else print $1, $2, $3, "loss" }' | /home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1/src/post_analysis/doc_walker /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.scov /home/dmacmillan/copy_num_ana/CNVer/workdir/chr22//chr22.gc | awk -v OFS="\t" '{ print $1, $2, $3, $4, $6 }' > /home/dmacmillan/copy_num_ana/CNVer/workdir/calls//chr22.cnvs...
    Here are the steps that I took to achieve this:
    Code:
    Steps:
    1. Untar the CNVer tar.gz source file
    2. cd into the directory
    3. type "./configure" (no quotations) and hit enter
    4. type "make" hit enter
    5. cd .. into folder containing cnver-0.8.1, hg18comp
    6. Created file "bam_files_to_use.txt" by doing
    	i. echo -e "file1.bam\nfile2.bam" > bam_files_to_use.txt
    7. type "export CNVER_FOLDER=/home/dmacmillan/copy_num_ana/CNVer/cnver-0.8.1" or wherever your cnver-0.8.1 folder is.
    8. ran the following command:
    ./cnver.pl --map_list /home/dmacmillan/copy_num_ana/CNVer/bams_to_use.txt --ref_folder /home/dmacmillan/copy_num_ana/CNVer/hg18comp --work_dir /home/dmacmillan/copy_num_ana/CNVer/workdir  --read_len 150 --mean_insert 398 --stdev_insert 58 --min_mps 2
    I do not receive a workdir/results folder in the end, and most, if not all the files in the chr[1-22]/x/y folders are 0 bytes. Any suggestions would be appreciated

Latest Articles

Collapse

  • seqadmin
    An Introduction to the Technologies Transforming Precision Medicine
    by seqadmin


    In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
    01-27-2025, 07:46 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 09:30 AM
0 responses
16 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-05-2025, 10:34 AM
0 responses
28 views
0 likes
Last Post seqadmin  
Started by seqadmin, 02-03-2025, 09:07 AM
0 responses
27 views
0 likes
Last Post seqadmin  
Started by seqadmin, 01-31-2025, 08:31 AM
0 responses
35 views
0 likes
Last Post seqadmin  
Working...
X