Thank you very much! Yes, I read somewhere that many people criticize gbff for this sort of thing, so I started playing around with the .gff. However, the next question that raises is: in the case of organisms with multiple chromosomes/plasmids, how do you then get the subseqs out of the gbff (once you have the coordinates from the gff) by pointing for a specific chromosome? I will now look into your code and see how it goes, I hope I can convert it to regular script! Thank you very much!
<quote>Have you printed out the value of $pt to verify that it is either "CDS" , "rRNA" or "tRNA" ?</quote>
Yes, it was retrieving "CDS", "rRNA" and "tRNA". Which is strange because then I can find some "Source" within "@two"