-
Notifications
You must be signed in to change notification settings - Fork 4
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
unmatched chromosome name get error #5
Comments
abother question is we know wgs always give big fragments of cnv. so why here the configure file, the window size is 500, people seems to use 1M instead of 500bp |
We recommend remap your samples to the ref genome provided by us to avoid some unexpected behaviour. Or, you can generate your own ref data according to https://www.yfish.org/display/PUB/Accucopy#Accucopy-3.7Makeyourownreferencegenomepackage |
You probably need to watch some videos or read some reviews/tutorials to understand how DNA is extracted from a cell, fragmented, and PCRed before it is put on a DNA sequencing machine. 500bp is NOT the CNA length. It is the average length of DNA fragments to be sequenced by a high-throughput DNA sequencer, i.e. Illumina HiSeq or NovaSeq. These so-called next-gen sequencers can only sequence 100-150bp for one fragment, not from start to end of a chromosome. Anyhow, you need to get familiar with what a next-gen sequencer can and cannot do.
|
thanks a lot, so can this tool accept fastq file instead of bam? |
thanks a lot. so how to understand here is 500 for segmentation |
No, Accucopy accepts bam file only.
500bp for segmentation is just a proper parameter base on our testing and you can set other value. |
Dear professor,
thansk for such a accurate software.
when I am using it, it raise errors like follows. MAYbe caused by chrM and chrMT.
what is more, since your genome.dict has many patch chromosome names, like
SN:GL000207.1
SN:GL000226.1
SN:GL000229.1
SN:GL000231.1
, but when users input bam, they may use a different version of genome, even for hg19, the patch chromosome seems to be different, so why not the input is a fastq, but a bam? ,can you give me some suggestions
ERRORS screenshot like following
The text was updated successfully, but these errors were encountered: