Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_012967 | 3,296,372 | 1 | . | G | 25.0% | 42.3 / 17.8 | 20 | pseudogene (338/405 nt) | ECB_RS16265 | DUF1016 N‑terminal domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base . (2/13); new base G (5/0); total (7/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGAGGCTTTTTTACCTCTTTTTTCAACA‑TGTTGAAATT‑CGCC‑AGACAGTGTC‑TGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑TTTGGCAGACGCTGTCTGCCAAATT > NC_012967/3296302‑3296441 | aTGAGGCTTTTTTACCTCTTTTTTCAACA‑TGTTGAAATT‑CGCT‑AGACAGTGTC‑TGGCG‑AATTAACACCAttt > 1:8571406/1‑71 (MQ=255) aTGAGGCTTTTTTACCTCTTTTTTCAACA‑TGTTGAAATT‑CGCC‑AGACAGTGTC‑TGGCG‑AATTAACACCATT‑g < 2:8448589/72‑1 (MQ=255) tttttttttCCCTTTTTTTTCAAAA‑AATTGAAATTCCCCC‑AAAAAGTTTG‑GGGGGAAATTAACACCATT‑ggg < 2:2297029/69‑1 (MQ=255) tgtttttttACCTTTTTTTTACAAA‑TTTTGAAATTCCCCC‑AAACGTTTTC‑GGGGG‑AATTAACACCATT‑GGGa < 1:2027351/69‑1 (MQ=255) ggtttttttACCTTTTTTTTCAAAA‑TTGTGAAATT‑CCCC‑AAAAAGTTTT‑GGGGGAAATTAACACCATT‑GGGa < 1:568728/69‑1 (MQ=255) ttttttCCTTTTTTTTCAAAA‑TGTTGAAATTCCGCA‑AGACGTTGTC‑GGGCG‑AATTAACACCATT‑GGGAATT‑c < 1:4552586/72‑1 (MQ=255) tttttCCCTTTTTTTACAAAT‑TTTTAAATTT‑CCCC‑AGACGGTGTC‑GGGCG‑AATTAACACCATT‑GGGAATT‑cc < 1:1892808/72‑1 (MQ=255) tttttACCTCTTTTTTCAACA‑TGTTGAAATT‑CGCC‑CGACAGTGTC‑TGGCG‑GAATAACACCCTT‑TGGGAtttc > 2:6659350/1‑70 (MQ=255) ttttttttCAAAATTGTTAAATTC‑CGCC‑AGACGGTGTC‑TGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGa‑ca < 2:3899002/70‑1 (MQ=255) tttttttcCAACA‑TGTTAAAATC‑CGCA‑AGACAGTGTG‑GGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGa‑ca‑c < 1:2209889/70‑1 (MQ=255) tttttttaAAAAA‑TGTTAAATTT‑CCCC‑AGACAGTGTC‑TGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGa‑ca‑c < 2:4059739/70‑1 (MQ=255) cAACA‑TGTTGAAATT‑CGCC‑AGACAGTGTC‑TGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑Ca < 1:8595819/72‑1 (MQ=255) aaCA‑TG‑TGAAATT‑CCCCAAAACAGTGTC‑GGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAg < 2:3719639/72‑1 (MQ=255) cGCC‑CGACAGGGTTTTGGGG‑AATTAAAACCATT‑GGGAATT‑CCCCCAA‑AA‑CCA‑CCT‑C‑CAG‑AATTCCC‑C‑TCGGc > 2:8338375/1‑72 (MQ=255) aGACAGTGTC‑TGGCG‑AATTAACACCATTGGGGAATT‑CC‑CCGA‑CA‑CCA‑CCT‑T‑CCGAAATTCCC‑C‑CCGGCaaaa > 1:7090820/1‑71 (MQ=255) aGACAGTGTC‑TGGCG‑AATTAACACCATT‑GGGGAAT‑TCGCCGA‑CA‑CCA‑ACT‑GCCAA‑AATTCCC‑‑‑TCGGCaaaaa > 1:8089642/1‑70 (MQ=255) gACAGTTTC‑TGGGG‑AATTAACACCATT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑t < 1:6862823/72‑1 (MQ=255) gACAGTGTC‑GGGGA‑AATTAACCCATTT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑t < 2:7457900/72‑1 (MQ=255) aCAGTGTC‑TGGCG‑AATTAACACCATTGGGGAATT‑TC‑CCGG‑CA‑CCACCCT‑G‑CGA‑AATTACC‑A‑TCGGCAAAGA‑t > 1:7000154/1‑72 (MQ=255) aCAGTGTC‑TGGCG‑AATTAACACCATTGGGGAATT‑CC‑CACA‑CA‑CCC‑CCTCC‑CAG‑AATTCCC‑C‑CCGGTCAAGA‑t > 1:3733483/1‑72 (MQ=255) aGTGTC‑TGGCG‑AATTAACACCATTGGGGAATT‑CC‑CCGA‑CACCCA‑CCC‑G‑CAA‑AATTCCC‑C‑TCGGCAAAAA‑ttt > 2:6306532/1‑72 (MQ=255) aGTGTC‑TGGCG‑AATTAACACCATT‑GGGAATT‑TCGCCGA‑CA‑CCC‑TCT‑G‑CAA‑AAATCCC‑A‑TCCGCCAAGA‑tttt > 2:8616839/1‑71 (MQ=255) aGTGTC‑GGGAA‑AATTAACNCCCTT‑TGGAATTCCCGAAAA‑CC‑CCA‑TCT‑G‑AAA‑ATTCCCC‑A‑TCGGCAAAGA‑ttt < 1:3417708/72‑1 (MQ=255) gtC‑GGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑TTTGGCa < 2:8129102/72‑1 (MQ=255) caaaaTAAACA‑CATT‑TGGAATTCCCGCAGA‑CC‑CCA‑TCT‑G‑AAG‑ATTCCCA‑A‑TCGGCAAAGA‑TTTGGCAGACGc < 2:6266259/69‑1 (MQ=255) g‑AATTAACAACATG‑TGGAATT‑CCGGAGA‑AC‑CCC‑TCC‑G‑GAA‑AATTCCCAA‑TCGGCAAAGA‑TTTGGCAGACGCt < 2:1917297/72‑1 (MQ=255) aaTTAACACCATT‑GGAAATC‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑TTTGGCAGACGCTGt < 1:8135161/72‑1 (MQ=255) aTTAACACCATTGGGGAATT‑CC‑CAGA‑CA‑CCC‑ACT‑G‑CAGAAATTCCC‑T‑CCGCCAAAAA‑TTTTGC‑GACGCCGtt > 2:4129530/1‑71 (MQ=255) ttAACACCGTG‑GGAAATT‑CCCCACA‑CC‑CCT‑TCG‑G‑CAA‑ATTCCCC‑A‑TCGGCAAAGA‑TTTGGCAGACGCTGTCt < 1:6778812/72‑1 (MQ=255) tAACCCCATT‑GGGAATTCCCGAAGA‑CA‑CCT‑TCT‑G‑‑AA‑AATTCCCAA‑TCGGCAAAGA‑TTTGGCAGACGCTGTCt < 1:4705771/72‑1 (MQ=255) tAACACCATT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑TTTGGCAGACGCTGTCTg < 1:6773987/72‑1 (MQ=255) aaCACCATG‑GGGATNN‑NCGCAGA‑CA‑CCA‑TCT‑G‑CAN‑AATTCCC‑A‑TCGGCAAAGA‑TTTGGCAGACGCTGTCTGc < 1:6227890/72‑1 (MQ=255) acacCATT‑TGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGATTTTGGCGGA‑GGTGTCTGcc > 1:6603587/1‑72 (MQ=255) cATT‑TGGAATT‑CCAGAGACCC‑CCA‑TCT‑G‑‑AA‑AATTCCCAA‑TCGGCAAAGA‑TTTGGCAGACGCTGTCTGCCaaa < 2:1836039/72‑1 (MQ=255) aTT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑TTTGGCAGACGCTGTCTGCCAAAtt < 2:8580058/72‑1 (MQ=255) tt‑GGGAACT‑CCGCAGA‑CC‑CCC‑TCT‑G‑CAA‑AATTCCC‑ATTCGGCAAAGA‑TTTGGCAGACGCTGTCTGCCAAAtt < 2:5998901/72‑1 (MQ=255) | ATGAGGCTTTTTTACCTCTTTTTTCAACA‑TGTTGAAATT‑CGCC‑AGACAGTGTC‑TGGCG‑AATTAACACCATT‑GGGAATT‑CCGCAGA‑CA‑CCA‑TCT‑G‑CAG‑AATTCCC‑A‑TCGGCAAAGA‑TTTGGCAGACGCTGTCTGCCAAATT > NC_012967/3296302‑3296441 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |