This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Parsing Coding Seq w grep

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

genome gene

1 answer

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

Log in to answer this question.