Eagle 2 - Phasing
Hello All,
I am working on a .vcf file with 200 samples and it looks like as follows
#CHROM POS ID REF ALT QUAL FILTER INFO FORMAT SRR2098177 SRR2098178
1 10146 . AC A 16738.1 VQSRTrancheINDEL99.00to99.90 . GT:AD:DP:GQ:PL 0/1:13,53:66:99:613,0,135 0/0:44,0:44:0:0,0,635 0/0:43,0:43:0:0,0,650 0/0:29
The file hold 1 to 23 + y chromosome information. My problem is, i wanted to conduct Phasing for this vcf file which is biallelic, using Eagle2, i have used forllowing command
eagle --vcf lic.vcf.gz --geneticMapFile=genetic_map_hg19_withX.txt.gz --outPrefix test.beagleImputed_phased --vcfOutFormat=z --chrom=1 --numThreads=32
Now it throws the error as
ERROR: Invalid chromosome: Y Chromosome number must be between 1 and --chromX (= 23)
I am not sure how to tackle this issue please help. Should we phase individual chromosomes? whats the exact usage of --chrom (i tried interpreting from manual, but its not that much helpful) It would be great help to me if you guys will help me in solving this issue, i am stuck here for almost 4 days without much luck. Its frustrating :(
Thanks a lot for your kind help
Dave !
• 2,988 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Hello David_emir!
Questions similar to yours can already be found at:
We have closed your question to allow us to keep similar content in the same thread.
If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.
Cheers!
PS: question asked twice
Thanks Pierre, by any chance can you help me with a solution for this?