I am looking for some script which provide the difference from r-CRS to the mtDNA human genome, directly from the GenBank!
e.g. I have a genbank id GQ337588.1
now I search in NCBI and get the fasta file there, after that I need some software which tells me the summary of the difference of this sequence with r-crs.
in above case the differences are-
A73G A263G 315.1C T489C T711C A750G A1041G A1438G A2706G T3394C
G4491A A4769G 5899.1C C7028T T7142C G7697A A8701G A8860G A9242G T9540C
A10289G A10398G C10400T T10873C G11719A C12705T T14308C A14417G C14766T T14783C G15043A G15301A A15326G G15927A C16223T C16234T A16316G T16362C
Thanks for any help.