Kuidas leida DNA järjestuse pöördkomplementi

DNA järjestuse pöördkomplement tähistab DNA molekuli vastasahela sisu. DNA molekulid konstrueeritakse sellistena, kuna igal nukleotiidil on komplementaarne nukleotiid teises ahelas, millega eksisteerib mittekovalentne side.

1
Jälgige järjestust tagurpidi, alustades järjestuse viimasest nukleotiidist.

2
Igast nukleotiidist üle minnes lisage selle komplementaarne nukleotiid järgmisele reale, alustades täiendatud stringi lehe vasakust servast. Pidage meeles, et guaniin (G) seostub tsütosiiniga (C) ja adeniin (A) tümiiniga (T).

3
Looge või aktsepteerige sisendfail. See artikkel eeldab, et sisend on FASTA-vormingus ja faili kohta on üks jada. Järgmised sammud eeldavad samuti, et kõik nukleotiidid on ATGC alused.

4
Lugege failist. FASTA vormingu puhul: loobuge faili esimesest reast.Eemaldage kõik ülejäänud reavahetused ja muud tühikud.def init(sequence): sisendiks open(argv[1]): järjestus = “”.join([line.strip( ) rea jaoks input.readlines()[1:]]) tagastusjada

5
Looge räsitabel, mis vastendab iga nukleotiidi selle komplementiga.komplement = {‘A’: ‘T’, ‘C’: ‘G’, ‘G’: ‘C’, ‘T’: ‘A’}

6
Korrake jada läbi ja kasutage komplementaarse jada koostamiseks räsitabeli otsingut. Pöörake saadud vektor.def reverse_complement(seq): alused = [täiendus[alus] järgus oleva baasi jaoks] alused = vastupidised(alused) tagastavad alused

7
Trüki vektori sisu.=result = reverse_complement(seq)print ”.join(result)