Biostar Beta. Not for public use.
extracting segments from .qual file
Entering edit mode
3.7 years ago
Stephen • 10

Is there anything similar to _getfasta _(bedtools) but will also work on an associated .qual file in the format:

    50 50 50 50 50 40 40 40 30 30 30 34 34 34 34

I basically need to extract a list of specific segments from the .qual file, and wondering if anyone knows of a suitable tool already available.

Entering edit mode
14 months ago
Cambridge, MA

You can use pyfaidx for this:

pip install (--user) pyfaidx
faidx scaffolds.qual Scaff3:1-8 Scaff3:4-11
50 50 50
50 50 50

Since a "qual" file has 3 characters per position, just make sure you count in threes to preserve your formatting.


Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3