MERS virus sequences dataset.

[xingjianleng]

• [x] Get the MERS virus sequence dataset.

[GavinHuttley]

from cogent3.parse.genbank import RichGenbankParser

path = "/Users/gavin/Desktop/Inbox/mers.gb"
with open(path) as infile:
    n, seq = list(RichGenbankParser(infile, moltype="dna"))[0]

len(seq) # 30111

mers.gb.zip

[GavinHuttley]

I just downloaded this via a search for mers at NCBI Viruses, clicked the "search by virus name", entered MERS, selected the "MERS-CoV", clicked on the "Ref genome" tab and selected the most recent. Downloaded that directly onto my machine in GenBank format.