|
|
Line 6: |
Line 6: |
| | colspan="2"| | | | colspan="2"| |
| <!-- ##### DO NOT edit above this line unless you know what you are doing. ##### --> | | <!-- ##### DO NOT edit above this line unless you know what you are doing. ##### --> |
| ==Mapping alignments to barcodes using command line BLAST==
| |
|
| |
| '''BLAST command line manual:'''
| |
| http://www.ncbi.nlm.nih.gov/books/NBK1763/
| |
|
| |
| '''BLAST command line dmg for Mac:'''
| |
| ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/LATEST/ncbi-blast-2.2.25+.dmg
| |
|
| |
| '''samtools call for viewing .bam:'''
| |
| samtools view filename.bam > filename.sam
| |
| note: this can also be piped to awk, etc. for parsing
| |
|
| |
| '''UNIX command to convert sam to fasta using awk'''
| |
| awk '{OFS="\t"; print ">"$1"\n"$10}' filename.sam > filename.fa
| |
|
| |
| '''command line call to generate indexed reference files for mapping to:'''
| |
| makeblastdb -in reference.fa -dbtype nucl -parse_seqids -out reference.fa -title "Reference_Title"
| |
|
| |
| '''command line blastn call:'''
| |
| blastn -task 'blastn-short' -query query.fa -db reference_index.fa -outfmt '6 qseqid sseqid pident length mismatch gapopen qstart qend sstart send evalue bitscore qseq sseq' -max_target_seqs 1 -reward 2 -gapopen 4 > output.txt
| |
|
| |
| '''R script for summarizing data:'''
| |
| [~/Work/BLAST/BLAST_parsing.R BLAST_parsing.R]
| |
|
| |
| <!-- ##### DO NOT edit below this line unless you know what you are doing. ##### -->
| |
| |}
| |
|
| |
| __NOTOC__
| |