Parsing Coding Seq w grep

Question

[maximilian.mayerhofer]

Hello,

I have downloaded the E.coli ref strain K-12 MG1655 coding sequence (https://www.ncbi.nlm.nih.gov/nuccore/NC_000913.3). I want to extract certain sequences of genes like gyrA, 23SrRNA,... etc

when using grep like

grep "[gene=gyrA]" EColiK12MG1655_cds.fasta it simply prints everything starting from gyrA (and the letters turn red and White for unknown reasons).

Best regards

Answer 1

[christopher medway]

I assume you want to extract sequence from the GenBank file. If this is true, there are dedicated libraries that will make this a whole lot easier, See this post