| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,764,774 | 0 | G | A | 42.9% | 41.8 / 10.2 | 35 | intergenic (‑388/+160) | sufA/ydiH | Fe‑S cluster assembly protein/uncharacterized protein |
| Reads supporting (aligned to +/- strand): ref base G (13/7); new base A (0/15); total (13/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.87e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGC > NC_000913/1764627‑1764921 | aaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGaat > 1:434190/1‑151 (MQ=255) aTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGaata < 1:706951/151‑1 (MQ=255) gACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGaataa < 2:308860/150‑1 (MQ=255) tttAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCAt > 2:383178/1‑151 (MQ=255) agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt > 2:248561/1‑151 (MQ=255) aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc < 2:25405/150‑1 (MQ=255) cGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGAGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAAcaca > 1:719334/1‑150 (MQ=255) aaaaaataatccccTAAA‑GTCTGTGCGACAC‑GGATTAAGGTGTTT‑CTCCCCCGCCGGGCGTTTATCTCGGCCCTTCTCTG‑GGATTAGCCTTTTGAGGGGGAGAAAAAAAAACCCAAATAAATCCACTCCGCCGTTTTAACACAAACTGCg < 1:592022/139‑1 (MQ=255) aaaaaagaaaCCCTTAAGGTCGGGGCGAAACAGG‑TTAAGG‑GTTTCCTCCCCCACCCGGCTCTTATCTCCGGCACTCTAGGTGGCTTAG‑CTTTTAAGGGGGAAATAAAAATAATCTCAAAAACCTAACCCGCCGTTTTAACACAAACt < 1:664570/140‑1 (MQ=255) cccTTAA‑GTCGGTGCGACACAGGTTAAAGG‑TTTT‑CTCCCCCATCCGGCCTTTATCCCGGGCACTCTCATGGGCTTTC‑CTCTTAAAGGGCTAATAAGAATAATCGTAAAATGCCAACCCGCCGTTTTAACACAAACTGCGATTAGtatta < 1:99306/149‑1 (MQ=255) aa‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTg > 2:313261/1‑150 (MQ=255) a‑GTCTGGGCAACACAGCCTTAAGG‑TTTT‑CTCCCCCACCCGGCCTTTATCCCCGCCACTCTCAGTGGCTTAG‑CTTTTGAAGGGGCAATAAAAATAATTTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa < 1:416154/151‑1 (MQ=255) acacggttttAGGG‑TTTT‑CTCCCCCACCCGGCGTTTATCCCCGCCCCTCTAGGTGTTTTAG‑CTTTTGAAGGGGCAAAAAAAATAATCTAAAAAAACAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAatat < 1:557352/139‑1 (MQ=255) acacAGGTTTAAGG‑TTTC‑CCCCCCCCCCCGGCCCTTACCCCCGCCACTCTCAGTGGCTTAG‑CTTTTGAAGGGGCAATAAAAATAATCTAATAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCa < 1:808549/149‑1 (MQ=255) gCTAAAGG‑TTTT‑CTCCCCCACCCGGCGTTTACCCCCGCCACCCTAAGTGGTTTAT‑CTTTTGAGGGGGCAAAAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt < 1:86796/149‑1 (MQ=255) ttAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt > 2:854023/1‑147 (MQ=255) ttAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt > 1:803856/1‑147 (MQ=255) ttAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt > 1:473029/1‑147 (MQ=255) gg‑gTTT‑CTCCCCCATCCGGCGCTTATCCCCGGCACTCTCAGGGGTTTAG‑CTTTTGAAGGGGCGAAAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATaa < 2:481825/149‑1 (MQ=255) ccccACCGGGCTCTTACCCCCGCCCCCCTCAGTGGTTTAG‑TTTTTAAAGGGGCAAAAAAAATAATCTAAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt < 1:261396/148‑1 (MQ=255) ccATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga > 1:179291/1‑150 (MQ=255) aTCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAgg > 2:600182/1‑151 (MQ=255) ggCGTTTTTCCCCGGCACTCTCATTGGTTTAG‑CTTTTAAAGGGGAGAAAAAAATATTTTAATAAATCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCAt < 1:695345/151‑1 (MQ=255) cTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGa > 1:4118/1‑150 (MQ=255) ttATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGCATTAAGAAGGTCATCGAAc > 2:715493/1‑151 (MQ=255) aTCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt > 1:406333/1‑151 (MQ=255) ggCACTCCGAGTGCCTTGG‑TTCTGGAGGGGGCGAAAAAAATAACCTCAAAAAGCAACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACgg < 2:825023/151‑1 (MQ=255) ctcAGTGGCTTAG‑CTCTTGAGGGGGCGATAAAAATAATCTAAAAAAGCAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGt < 2:734384/149‑1 (MQ=255) ggCTTAG‑CTTTTGAAGGGGCAATAAAAATAATCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAAGTTa < 2:106590/145‑1 (MQ=255) ttAG‑CTTTTGAGGGGGAGAAAAAAAAAATCTCAAAAAGATAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATc < 1:14420/145‑1 (MQ=255) agggggaaaTAAGAAAAACCTCAAAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGAtt < 2:206055/137‑1 (MQ=255) aaGGGGCGATAAGAATAATCTAAAAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTg < 2:467134/147‑1 (MQ=255) ggggCGATAAAAATAATCTAATAAAGTAAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTg < 2:189705/145‑1 (MQ=255) gggCAAAAAAAAAAATCTAATAAAGCTACCCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACt < 1:91041/151‑1 (MQ=255) aaGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGc < 1:360970/150‑1 (MQ=255) | AATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAA‑GTCTGTGCGACACAGGCTTAAGG‑GTTT‑CTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAG‑CTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGC > NC_000913/1764627‑1764921 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |