\begin{Verbatim}[commandchars=\\\{\}]
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{k+kn}{from} \PYG{n+nn}{Bio} \PYG{k+kn}{import} \PYG{n}{SeqIO}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{my\PYGZus{}snp} \PYG{o}{=} \PYG{l+m+mi}{4350}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{n}{record} \PYG{o}{=} \PYG{n}{SeqIO}\PYG{o}{.}\PYG{n}{read}\PYG{p}{(}\PYG{l+s+s2}{\PYGZdq{}NC\PYGZus{}005816.gb\PYGZdq{}}\PYG{p}{,} \PYG{l+s+s2}{\PYGZdq{}genbank\PYGZdq{}}\PYG{p}{)}
\PYG{g+gp}{\PYGZgt{}\PYGZgt{}\PYGZgt{} }\PYG{k}{for} \PYG{n}{feature} \PYG{o+ow}{in} \PYG{n}{record}\PYG{o}{.}\PYG{n}{features}\PYG{p}{:}
\PYG{g+gp}{... }    \PYG{k}{if} \PYG{n}{my\PYGZus{}snp} \PYG{o+ow}{in} \PYG{n}{feature}\PYG{p}{:}
\PYG{g+gp}{... }        \PYG{k}{print}\PYG{p}{(}\PYG{l+s+s2}{\PYGZdq{}}\PYG{l+s+si}{\PYGZpc{}s}\PYG{l+s+s2}{ }\PYG{l+s+si}{\PYGZpc{}s}\PYG{l+s+s2}{\PYGZdq{}} \PYG{o}{\PYGZpc{}} \PYG{p}{(}\PYG{n}{feature}\PYG{o}{.}\PYG{n}{type}\PYG{p}{,} \PYG{n}{feature}\PYG{o}{.}\PYG{n}{qualifiers}\PYG{o}{.}\PYG{n}{get}\PYG{p}{(}\PYG{l+s+s2}{\PYGZdq{}db\PYGZus{}xref\PYGZdq{}}\PYG{p}{)))}
\PYG{g+gp}{...}
\PYG{g+go}{source [\PYGZsq{}taxon:229193\PYGZsq{}]}
\PYG{g+go}{gene [\PYGZsq{}GeneID:2767712\PYGZsq{}]}
\PYG{g+go}{CDS [\PYGZsq{}GI:45478716\PYGZsq{}, \PYGZsq{}GeneID:2767712\PYGZsq{}]}
\end{Verbatim}
