Note that there are some explanatory texts on larger screens.

plurals
  1. POfind unique first top and bottom lines of fastq file from fasta file
    primarykey
    data
    text
    <p>I have 2 files one is <code>fasta</code> file and other one is <code>fastq</code> file. I want to take the <code>fasta</code>, read each line and search each line in the <code>fastq</code> file and print top line and bottom lines. This is what I have</p> <p><strong>fasta file</strong></p> <blockquote> <p>read1</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAG</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGA</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAGG</p> </blockquote> <pre><code>for seq in `cat sequences`;do grep -A 2 -B 1 $seq FAP.1.txt;done </code></pre> <blockquote> <p>@DH1DQQN1:269:C1UKCACXX:1:1107:20386:6577 1:N:0:TTAGGC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC</p> <p>+</p> <p>CCCFFFFFHGHHHJIJHFDDDB173@8815BDDB###############</p> <p>@DH1DQQN1:269:C1UKCACXX:1:1114:5718:53821 1:N:0:TTAGGC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</p> <p>+ ;@?DBD&lt;@@FFHHH&lt; <p>@DH1DQQN1:269:C1UKCACXX:1:1209:10703:35361 1:N:0:TTAGGC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</p> <p>+</p> <p>@@@FFFFFHGHHHGIJHFDDDDDBDD69@6B-707537BDDDB75@@85</p> <p>@DH1DQQN1:269:C1UKCACXX:1:1210:18926:75163 1:N:0:TTAGGC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAG</p> <p>@CCFFFFFHHHHHJJJHFDDD@77BDDDDB077007@B###########</p> </blockquote> <p>From this we can see that <code>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</code> appeared twice, but I want to print only once. how can I do that?</p> <p><strong>Final output file</strong></p> <blockquote> <p>@DH1DQQN1:269:C1UKCACXX:1:1107:20386:6577 1:N:0:TTAGGC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC</p> <p>+</p> <p>CCCFFFFFHGHHHJIJHFDDDB173@8815BDDB###############</p> <p>@DH1DQQN1:269:C1UKCACXX:1:1114:5718:53821 1:N:0:TTAGGC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</p> <p>+</p> <p>;@?DBD&lt;@@FFHHH&lt; <p>@DH1DQQN1:269:C1UKCACXX:1:1210:18926:75163 1:N:0:TTAGGC</p> <p>AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAG</p> <p>+</p> <p>@CCFFFFFHHHHHJJJHFDDD@77BDDDDB077007@B</p> </blockquote>
    singulars
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    plurals
    1. This table or related slice is empty.
    1. This table or related slice is empty.
    1. This table or related slice is empty.
 

Querying!

 
Guidance

SQuiL has stopped working due to an internal error.

If you are curious you may find further information in the browser console, which is accessible through the devtools (F12).

Reload