Open subwaystation opened 4 years ago
Oh, just discovered there is fpa keep
:
[heumos@wave ~]$ fpa keep --help
fpa-keep
fpa keep only mapping match this constraints
USAGE:
fpa keep [FLAGS] [OPTIONS] [SUBCOMMAND]
FLAGS:
-c, --containment Keep only containment mapping
-d, --dovetail Keep only dovetail mapping
-h, --help Prints help information
-i, --internalmatch Keep only internal mapping
-m, --same-name Keep only mapping where reads have same name
-V, --version Prints version information
OPTIONS:
-l, --length-lower <length_lower> Keep only mapping with length lower than value
-L, --length-upper <length_upper> Keep only mapping with length upper than value
-n, --name-match <name_match> Keep only mapping where one reads match with regex
-s, --sequence-length-lower <sequence_length_lower>
Keep only mapping where one reads have length lower than value
-S, --sequence-length-upper <sequence_length_upper>
Keep only mapping where one reads have length upper than value
The working commands:
zcat Sc+Sp.pan.paf.gz | fpa keep -L 10000 | pigz >Sc+Sp.pan.fpal10k.paf.gz
zcat Sc+Sp.pan.paf.gz | fpa keep -L 2000 | pigz >Sc+Sp.pan.fpal2k.paf.gz
Running the major script I get:
Actually
fpa
does not have this option: