"sequenza" - chromosomes are removed when I run 'sequenza.extract'.
Hi.
I'm trying to run a "sequenza", but something is wrong. Before I made a pileup file, the bam file has chromosome 1 - 22 and X, Y but when I run a sequenza.extract, there is only chr1.
test <- sequenza.extract("out_PAAD_small.seqz.gz") Processing chr1: 42 variant calls; 2970 heterozygous positions; 357535 homozygous positions. names(test) [1] "BAF" "ratio" "mutations" "segments" "chromosomes" [6] "gc" "avg.depth" test$chromosomes [1] "chr1"
I don't know what happend. Could anyone let me know how can I get all chrs as an output?
• 2,030 views
•
link
written
by
a01039012387 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
output same result using different imputed data and same reference panel
written by sh •
-
Problem while working with sequenza - Chromosomes out of order
written by Camilo Andres •
-
Linkage disequilibrium from VCF file of variants in bacterial metagenomes
written by matthias.hotzinger •
-
PureCN issue: having problems running a small example
written by bboy •
-
Tutorial: Yes .. BBMap can do that! - Part II randomreads (read simulation), demuxbyname/filterbyna…
written by GenoMax
-
Filtering for individual allelic balance in multisample vcf using GATK4
written by mernster •
-
Encountering errors with PAGIT icorn
written by Ace
-
Perl script for reformatting a plink/seq output file
written by plink_9857
-
A simple question regarding homozygous and heterozygous variation.
written by mangfu100
-
F(St) Test? Coorelation Between Genetic Distance And Geographical Distance Between Isolates?
written by Adrian Pelin
What is in the file: out_PAAD_small.seqz.gz? Try and see if all the chromosomes are in the file. I tried sequenza on the test data set and it worked fine, so it is possible that there is a problem with the input.