I used 23andMe to download my raw genome. I have it in a .txt file but you can't use the format for real bio programs. I want to make my own library for further analysis. Does anyone know how I can convert .TXT to FASTA, GenBank, or any other usable file type?
•
written
by
someashole •
The commenters are correct re: first figuring out something you want to do with the data, and checking what input formats work with that.
With that said, here are two conversions that might come in handy.
23andMe to VCF: this is now supported by PLINK.
plink --23file [name of your file] --snps-only no-DI --recode vcf
What's --snps-only no-DI, you ask? Well, 23andMe files contain mostly SNP calls, but there are a few indel calls as well. Unfortunately, the actual bases involved in the indels are NOT saved; instead, there's just 'D' for deletion and 'I' for insertion, and you'd need an indel database to determine a valid VCF representation of the call. So we just punt here and filter out all markers with 'D' or 'I' allele codes.
VCF to FASTA:
This can be done with a combination of VCFTools and a Perl script. See here for details.
•
written
by
chrchang523
Once you have the VCF, check the answers in this discussion: New Fasta Sequence From Reference Fasta And Variant Calls File?
This issue I have with the 23file is that I get invalid chromosome code 85280 on line 585542 of .bim file.
(Use --allow-extra-chr to force it to be accepted.)
But when I add that it tells me you cant use --allow-extra-chr as it cannot currently be with --23file.
Anybody now how to fix this?
Jeff
Jeff, try this:
plink --23file 23andmefile.txt Surname Firstname Sex --snps-only no-DI --make-bed --out plink_genome
Log in to answer this question.
Can you provide an example of what your data looks like? Then if you could provide an example of what you want the output to look like that would also help.
there are about 960k lines
I'm not sure how the data should look for a usable format
The output that you get above is already the most compact form that you can get your data in. It represents the differences relative to the reference genome.
You could for example transform this to two diploid genomes in FASTA format but do you realize that your files would then be gigantic ones of many gigabytes and these files would not show you what the the changes were.
The right way to go about this is to formalize what do you want to do next with your data. Then depending on that aim people here can advise what to transform it to.
Could you convert this to FASTA though? Are the genotype alleles listed so that one sister chromosome is always first and the other is always second? Or is it random? If it's random, there's no way to construct FASTA because we don't know if, for example at 752566 and 752721 we have A-A and G-G or A-G and G-A.
If I was going to do anything with it, I think I'd want a VCF file.