repeatmasker file analysis
Hi I am new to bioinformatics. I have a command that I don't know what the output file will be? Can anyone help me understand the command? I don't know what the system does in the command line.
awk -F "\t" '{printf $0 "\t"} {system("samtools view -F 0x100 chr1.bam " $6 ":" $7"-"$8 "|" "wc -l")}' chr1RPTS > final.txt
the chr1RPTS is a file in repeatmasker format.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
RNA seq strandedness (reverse-strand)
written by ZuelTech •
-
Seeking advice on methods to calculate ENCODE ChIP-seq data standard M_DISTINCT
written by kalavattam •
-
Adding custom annotation to VEP output VCF
written by avelarbio46 •
-
Multiple commands in a for loop
written by lamia_203
-
Chimeric read insert length in SAM $9
written by Yijun Tian •
-
Number of discovered genes in genome does not correspond with published annotation
written by Jan Röslein •
-
Forum: Awk in Bioinformatics
written by Shicheng Guo
-
what does this formula represent?
written by AAK85 •
-
Extracting discordant Pairs
written by ARB •
-
discordant results between bedTools and awk
written by bioguy24
Output file will be called
final.txt, if the command works as written (and as advertised). Are you getting anything in thefinal.txtwhen you run this command?Hi Yes. In my ch1RPTS file I have 17 columns and after this command my final.txt file has 18th column. I don't understand what the 18th column represents? Generally, in the repeatmasker file there are 17 columns, executing the {system("samtools view -F 0x100 chr.bam .......")} command adds one more column to my output file and I can't figure it out what it is representing? my complete command is this: awk -F "\t" '{printf $0 "\t"} {system("samtools view -F 0x100 chr1.bam " $6 ":" $7"-"$8 "|" "wc -l")}' chr1RPTS > final.txt