[oudacontrol]

I am trying to view a portion of a sorted, indexed bam file and I am getting the error:
[sam_header_read2] 26 sequences loaded.
[main_samview] random alignment retrieval only works for indexed BAM files.

It is already sorted and indexed.... I have no idea how to resolve this.

[Strand SI]

I converted some ELAND export files to SAM and now tried to make BAM files out of those, but it is complaining about @SQ headers not being there?!
Does the export2sam.pl tool not create a SAM file format that is correct??

[ElMichael]

Thank you, Heng!

[lh3]

This is fast. If present, the ordering of tags are fixed. You can

if /X0:i\d+).*XA:Z\S+)/

to grep out everything with a single pass. You need to read and split after all. Regex matching should be of similar speed, if not faster. The bad implementation is to loop through the optional tags and then try to find matching.

[ElMichael]

For perl, very easy and efficient:

perl -ne 'print "$1\n" if /XA:i\d+)/'

It's exactly what I want to avoid - checking each line with IF condition.

[lh3]

For perl, very easy and efficient:

perl -ne 'print "$1\n" if /XA:i\d+)/'

[ElMichael]

optional fields in SAM format

Could anybody advise me the best way to parse optional fields in SAM file? I'm interested in tags X0 and XA, and the problem is that they do not have a permanent place in the output string (e.g., sometimes 13th or 14th, 18th, 19th or 20th). Certainly, I can do it myself in perl, but such parsing would take too much time and resources.
What I'm looking for is:
1) either a way to control output tags (their number and order)
2) or (better) some bioperl module, like blast parser, that could parse it effectively.
thanks!

[maubp]

can any 1 tell me please how can i install SAMtools on my windows operating system...Help Please...thanks

Have you tried the provided Windows binaries on sourceforge? i.e. samtools-0.1.12a_i386-win32.zip

[I'd suggest getting hold of a Linux (or Mac OS X) machine if you can, since not that many of the NGS tools will work on Windows]

[seqsyd]

can any 1 tell me please how can i install SAMtools on my windows operating system...Help Please...thanks

[andrewm]

samtools index and alternative alignments

I would like to index/search the alternative alignments found by bwa with samtools.

Currently it seems like samtools index does not index the alternative alignments from the XA flag. It would be great if this worked (maybe as an option).

Another possibility would be for bwa to output alternative alignments as another sam line.

Are either of these possible or planned for future releases?

Andrew

[gaffa]

Does anyone know what has happened to indel lines in the new samtools mpileup format? In the old pileup, each indel is followed by an additional line carrying some further info - when I run mpileup (samtools 0.1.12a) I don't get any such lines. The old pileup also has a flag ("-i") for outputting only indel variants - there does not seem to be such an option for mpileup.

[ashokrags]

perhaps try using reheader option in samtools. i think you can filter out those reads yo want using the view option and then use the reheader option

[Zigster]

What is the safest way to reheader a BAM file generated by an alignment to the human_g1k_v37.fasta genome, i.e.

cat GRCh37/human_g1k_v37.fasta.fai
1 249250621 52 60 61
2 243199373 253404903 60 61
3 198022430 500657651 60 61
4 191154276 701980507 60 61
5 180915260 896320740 60 61
6 171115067 1080251307 60 61
7 159138663 1254218344 60 61
8 146364022 1416009371 60 61
9 141213431 1564812846 60 61
10 135534747 1708379889 60 61
11 135006516 1846173603 60 61
12 133851895 1983430282 60 61
13 115169878 2119513096 60 61
14 107349540 2236602526 60 61
15 102531392 2345741279 60 61
16 90354753 2449981581 60 61
17 81195210 2541842300 60 61
18 78077248 2624390817 60 61
19 59128983 2703769406 60 61
20 63025520 2763883926 60 61
21 48129895 2827959925 60 61
22 51304566 2876892038 60 61
X 155270560 2929051733 60 61
Y 59373566 3086910193 60 61
MT 16569 3147273397 70 71
GL000207.1 4262 3147290265 60 61
GL000226.1 15008 3147294661 60 61
GL000229.1 19913 3147309982 60 61
GL000231.1 27386 3147330289 60 61
GL000210.1 27682 3147358194 60 61
GL000239.1 33824 3147386400 60 61
GL000235.1 34474 3147420850 60 61
GL000201.1 36148 3147455961 60 61
GL000247.1 36422 3147492774 60 61
GL000245.1 36651 3147529866 60 61
GL000197.1 37175 3147567190 60 61
GL000203.1 37498 3147605047 60 61
GL000246.1 38154 3147643232 60 61
GL000249.1 38502 3147682084 60 61
GL000196.1 38914 3147721290 60 61
GL000248.1 39786 3147760915 60 61
GL000244.1 39929 3147801427 60 61
GL000238.1 39939 3147842084 60 61
GL000202.1 40103 3147882751 60 61
GL000234.1 40531 3147923585 60 61
GL000232.1 40652 3147964854 60 61
GL000206.1 41001 3148006246 60 61
GL000240.1 41933 3148047993 60 61
GL000236.1 41934 3148090687 60 61
GL000241.1 42152 3148133382 60 61
GL000243.1 43341 3148176299 60 61
GL000242.1 43523 3148220425 60 61
GL000230.1 43691 3148264736 60 61
GL000237.1 45867 3148309218 60 61
GL000233.1 45941 3148355912 60 61
GL000204.1 81310 3148402681 60 61
GL000198.1 90085 3148485409 60 61
GL000208.1 92689 3148577058 60 61
GL000191.1 106433 3148671355 60 61
GL000227.1 128374 3148779625 60 61
GL000228.1 129120 3148910202 60 61
GL000214.1 137718 3149041537 60 61
GL000221.1 155397 3149181614 60 61
GL000209.1 159169 3149339664 60 61
GL000218.1 161147 3149501549 60 61
GL000220.1 161802 3149665445 60 61
GL000213.1 164239 3149830007 60 61
GL000211.1 166566 3149997047 60 61
GL000199.1 169874 3150166453 60 61
GL000217.1 172149 3150339222 60 61
GL000216.1 172294 3150514304 60 61
GL000215.1 172545 3150689533 60 61
GL000205.1 174588 3150865017 60 61
GL000219.1 179198 3151042578 60 61
GL000224.1 179693 3151224826 60 61
GL000223.1 180455 3151407577 60 61
GL000195.1 182896 3151591103 60 61
GL000212.1 186858 3151777111 60 61
GL000222.1 186861 3151967147 60 61
GL000200.1 187035 3152157186 60 61
GL000193.1 189789 3152347402 60 61
GL000194.1 191469 3152540418 60 61
GL000225.1 211173 3152735142 60 61
GL000192.1 547496 3152949898 60 61

to something that would be more compatible with the UCSC-centric Bioconductor stack, i.e.

cat hg19.fasta.fai
chr1 249250621 6 50 51
chr2 243199373 254235646 50 51
chr3 198022430 502299013 50 51
chr4 191154276 704281898 50 51
chr5 180915260 899259266 50 51
chr6 171115067 1083792838 50 51
chr7 159138663 1258330213 50 51
chr8 146364022 1420651656 50 51
chr9 141213431 1569942965 50 51
chr10 135534747 1713980672 50 51
chr11 135006516 1852226121 50 51
chr12 133851895 1989932775 50 51
chr13 115169878 2126461715 50 51
chr14 107349540 2243934998 50 51
chr15 102531392 2353431536 50 51
chr16 90354753 2458013563 50 51
chr17 81195210 2550175419 50 51
chr18 78077248 2632994541 50 51
chr19 59128983 2712633341 50 51
chr20 63025520 2772944911 50 51
chr21 48129895 2837230949 50 51
chr22 51304566 2886323449 50 51
chrX 155270560 2938654113 50 51
chrY 59373566 3097030091 50 51
chrM 16571 3157591135 50 51

So I want to toss the unscaffolded and haplotyped sequences and rename the rest.

[maximilianh]

The MACS distribution includes two scripts: elandexport2bed.py elandmulti2bed.py elandresult2bed.py

These might be more up-to-date than your scripts and easier to use in some cases.

[xmluo]

I used maq2sam-long to convert maq output to sam format, but all pairing information is missed in the results: MRNM is "*", and MPOS and ISIZE are 0. Can you recommend how to get these information? Thanks, Mei