mutiple fasta to single fasta
i have a huge reference genome with a lot of contigs, it looks something like this.
>aalba5_s00000010
TTGTCTGCTTCACAGTACAGCTAGAAAATTATGAATTCATTTCCCCACATCAAGCAACCCCTGCTTATTC
>aalba5_s00000011
ACTTGGAATGGGATCTTGTTGGGGGGCCAACAGAACCATAAGGGCAATGGCTGCAATCTTTGATAAGATC
>aalba5_s00000012
TGTAGCAAACAGCTACGGAAAAATTTTAAAAATTTTCGAAATTTAAATCTGGGGTTCCCTTTCCTGTGTA
GATGTATTCCCTTTTTAAAGGTTTTCCTAGGACTTGCAGTCATTAATGAGACGTCTTCTCATGATATCCT
AATTTTTGGAAGATGCCTCCTACATCAGGAATCTTTGCTGCCACTTGTCTCTTTCATCAGCCAGATGTCT
how can i subset this that i have a file each with the filename of the name of the contig (examplea alba5_s00000010.fa) conatining its sequence?
tanks for the help
• 1,736 views
•
link
updated
by
Renesh
•
written
by
setschmann •
•
written
by
setschmann •
2 answers
• 0 views
•
link
written
by
Pierre Lindenbaum
Log in to answer this question.
More posts like this
-
magpurify errors
written by shevch2009 •
-
Correcting coordinates in gff3 file
written by hulu •
-
Mitochondrial contig identification in fungal genome assemblies
written by VineethVenumadhavan •
-
freebayes only picking up variants in a single contig
written by gt •
-
Tool: Filtering BAM files efficiently
written by ddeemer •
-
How to filter a VCF file with a list of CHR or contig IDs?
written by acoles •
-
Assemble bacterial .fastq files and find differences (SPAdes followed by MUMmer)
written by LRStar
-
Rename fasta file after organism name obtained from NCBI
written by Guillaume.Tahon •
-
Issue using IDBA
written by Chvatil
-
python script to sort blast results and fetch best matches from NCBI
written by skbrimer
You can try below Python code for your file
Replace the n with number of sequences in your file. Read more here https://www.reneshbedre.com/blog/filereaders.html#split-fasta-file-into-multiple-fasta-files