I am using this package nsdpy to download genome sequences from NCBI nucleotide database.
Specifically I am interested in the whole mitochondrial genome of different species, here I will use a subset of species just as an example.
My Query is:
"{mitochondrion[Title]" -i -t -L subset_species_list.txt
Where :
-i
Adds the taxonomic information to the information lines of the sequences written in the output files (not relevant for this example)-t
option produces a tsv file with the following fields: Taxon name, SequenceID, TaxID, Taxonomic lineage, Sequence length, Sequence.subset_species_list.txt
is a list of species, below you can see its content
Xenopus laevis
Gasterosteus aculeatus
Chironomus riparius
Rasbora heteromorpha
Planktolyngbya limnetica
The output file is quite long and includes the actual sequences. Here I will post the same file where I dropped the column for the sequence and the lineage to keep it compact
My aim is to obtain phylogenetic distances between species.
My questions are:
- What do the sequences for the same species represent?
- Which sequences are best suited for an alignment and how do I know that?
- Should I also look into the orientation of the sequences or are they all on the forward strand?
- Is it better if I filter for Coding sequences (CDS) instead of whole genomes
Also if you noticed any rookie mistake in my process feel free to point that out. It is the first time for me working with DNA sequences and I really want to learn to correct way to do that.