Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,862,661 | 0 | T | G | 41.5% | 29.2 / 20.7 | 55 | pseudogene (943/1617 nt) | glvC | pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component |
Reads supporting (aligned to +/- strand): ref base T (5/26); new base G (21/1); total (28/27) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.05e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01 |
CATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGG > NC_000913/3862518‑3862805 | cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCggg > 1:902311/1‑150 (MQ=255) cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑AGGTGCggg > 2:802758/1‑150 (MQ=255) gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATc < 1:991787/151‑1 (MQ=255) gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCgggga > 1:18151/1‑147 (MQ=255) tGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCa < 1:604068/150‑1 (MQ=255) cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑a < 1:633/151‑1 (MQ=255) cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTA‑AGTTGCggg > 2:511030/1‑143 (MQ=255) ccaccaCACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGTggggaaca > 1:47513/1‑142 (MQ=255) caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAg < 2:157271/151‑1 (MQ=255) caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGATCAGc > 2:1105383/1‑148 (MQ=255) caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTA‑GGGTGCGGGGATCAGC‑a > 1:351332/1‑149 (MQ=255) accacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAg < 2:974883/150‑1 (MQ=255) cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGc < 2:487864/149‑1 (MQ=255) aaaGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTTAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGAACAGC‑AAACCCGcccc > 1:1000271/1‑148 (MQ=255) agagATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTTGGGGAAAAAGC‑AAACCCCCCCCttt > 1:935567/1‑151 (MQ=255) ttACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGATAATC‑AAACCCGCCACCTTTTCCCGAtt > 2:328574/1‑151 (MQ=255) tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAA‑GGGTGCGGGGATTAGC‑AAGACCGCCCCTTTTTCCCGCttt > 1:1036633/1‑151 (MQ=255) cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATTAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCTGCATGGCGGCTAA‑AGGTGGGGGAAAAAGG‑AAACCCCCCCCCTTTCCCCCAtgt > 1:444784/1‑147 (MQ=255) cGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGAACAGC‑AAACCCGCCACCTTTACCCCAtttt > 1:1110584/1‑150 (MQ=255) aTTGAGCCCGCCAGCACGGGGTGCACCGCAAACAGCAACGGAAAAATGAACAGAAAGGTAAATCCCAGCGGTTCGGTAATCCCCCCCAGCATGGCGTTTAA‑GGTGGCGGAAATCAGC‑ACGCCCGCCACTTTTACCCGATTTTCCGGTg < 2:419164/148‑1 (MQ=255) gcggCCCCCAGCAGGGCGCCTCCCGAAAAGAGCACCGGTAAAAGAAACAGGAAGGTAAATCCCGGCGGTGCGATAATCCCCACCAGAATGG‑GGTTAAGGGTTGCGGCAAACACC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAg < 2:334836/149‑1 (MQ=255) gccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAA‑GGTTGCGGGAATCAGC‑CAGGCCGCCACTTTTTCCCCACTTTTCGGTGTgg > 2:10241/1‑145 (MQ=255) ccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑AGGTGCGGGAATCAGC‑AAGCCCGCCCCTTTTTCCCCCTTTTCCGGTGCGGcgggg > 2:213463/1‑142 (MQ=255) cAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGAACAGC‑AAGCCCGCCACTTTTACCCGCTTTTTCGGTGCGGGAggg > 2:420850/1‑143 (MQ=255) gCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGGATCAGC‑AAACCCCCCACTTTTACCCGATTTTTCCGTGCGGCAGTGAAGTACATCGc > 1:1003469/1‑149 (MQ=255) tgtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATAA‑GGGTGCGGGAAACAGC‑AAGCCCGCCCCCTTTACCCGCTTTTCCGGTGCGGCAGTGAAGTACAACGCtca > 2:456970/1‑147 (MQ=255) ccGCAAACAGCAAGGGTAAAAGAAACGGGAAGGTAAATTCCGGGGTTTCGTTAATCCCCCCCAGCATGGCGGTTAG‑GGTTGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa < 1:929259/150‑1 (MQ=255) cGCAAACAGCAACGGTAAAATAAACAGGAAGGTAAATCCCGGCGGTTCGGTAATTCCCCCCAGCAGGGGGTTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa < 1:401511/149‑1 (MQ=255) gCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTCA‑AGGGGGGGGAATCAGC‑AAGCCCGCCCATTTTTCCCCGTATTCCCGGGCGGCAGGGCGGAACCTCGGTTAAGAAATg > 1:280253/1‑150 (MQ=255) accAGCACGGGAGAAAGGAACAGAAAGATAAATCCCAGCGTTTCGGTATTCCCCCCCAGCGTGGGGTTTAG‑GGTGGGGGAAAGCAAC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATg < 2:581352/145‑1 (MQ=255) cAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCa > 1:551992/1‑149 (MQ=255) cGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGGATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAAGCCCACGGCAcc > 1:347763/1‑150 (MQ=255) gTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGATGCGGGAAACAGC‑AAGCCCCCCACTTTTATCCGATTTTTCGGTGCGGCAGCGGAGTACATCGCCAAAAAAATTCCCACGGcccc > 2:214607/1‑148 (MQ=255) aaaaGAACAGGAAGGAAAATTCCGGCGTTTGGGTAATCCCCACCAGCAGGGCGGTTAA‑GGTGGCGGAAACCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGa < 1:132012/150‑1 (MQ=255) aaTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGGAGGGAAGTACCTCGCTAAAGAAAATCCCACGGCACCAAAGAtt > 2:299723/1‑151 (MQ=255) aTGAACAGGAAGGTAAATTCCGGGGTTTCGGTATTCCCCACCAGCATGGCGTTTAG‑GGTTGCGGAAATCAGA‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt < 1:1013111/151‑1 (MQ=255) aTGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTA‑TGCTTCCGGGATCAGC‑CAACCCCCCCCTTTTACCCGATTTTCCTGCGCGGCAGTGGAGTACAACGCTAAAAAAAAGCCCACCGGGCCAAAGAttt > 1:544434/1‑151 (MQ=255) taaaCAGGAAGGAAAATTCCAGCGGTTCGGTAATCCCCACCAGCTTGGCGGTTAG‑GGTTGGGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttt < 1:921932/148‑1 (MQ=255) tGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGAATCAGC‑AAGCCCCCCACTTTTACCCCATTTTCCGGTGCGGCAGGGAAGTACATCGCTAAAGAAAAGCCCACGGCACCAAAGAtttt > 2:248559/1‑151 (MQ=255) aacggaaaGGAAATTACCGGCGGTGCGTTATTCCCCACCAGCGGGGCGATTAA‑GGGGGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt < 1:1062072/144‑1 (MQ=255) aaaaGGAAGAAAATTTCCAGCGGTTCGGTAATCCCCACCAGCAGGGGGGTAAG‑GGTGGCGGAATTCAGA‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttt < 2:489789/147‑1 (MQ=255) cAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑AGGTTCGGGAATCAGC‑AAGCCCCCCACTTTTACCCGAATTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAAACCCACCGCACCAAAGAttttt > 2:284695/1‑148 (MQ=255) ggAAGGTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCAGGGGGGTTAA‑GGTGGGGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGa < 2:523712/148‑1 (MQ=255) aaGGAAAATCCCAGCGGTTCGGTAATCCCCACCAGCATGGGGTTTAG‑GGTGGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTAc < 2:551866/151‑1 (MQ=255) gTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGGTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGTACATCGCTAAAAAAATGCCCACGGCCCCAAAGATTTTTGAGTTACCg > 1:629831/1‑150 (MQ=255) aaaTTCCGGCGGTTCGGTAATCCCCACCAGCAGGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc < 1:26347/151‑1 (MQ=255) aaaTCCCAGCGGTTCGTTAATCCCCACCAGCAGGGCGGTTAG‑GGTGGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc < 2:546857/151‑1 (MQ=255) ggTTCGGTAACCCCCACCAGCGTGGCGGTTAA‑GGTGGCGGAAATCAGC‑AACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCaaaa < 2:1116550/150‑1 (MQ=255) gcgatAATTCCCCCAAGCGGGGCGGTAAA‑GGTTGC‑GGAACCAACAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc < 2:265275/147‑1 (MQ=255) ggTATTCCCCCCCAGCGGGGCGGTAAG‑GTTTGGGGAAACCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAcc < 1:93146/147‑1 (MQ=255) ccccccAGCATGGGGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTcc < 2:458776/147‑1 (MQ=255) ccccccAGCATGGCGTTTAA‑GGTGGGGGAAACCAGC‑AGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGa < 1:315070/151‑1 (MQ=255) gCAGGGCGGTTAA‑GGTTGCGGAAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaacaac < 1:802964/149‑1 (MQ=255) cttGGCGGTTAA‑GGTCGCGGAAAACAGC‑ACCGCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACg < 1:1026315/147‑1 (MQ=255) ttAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCgg > 2:341562/1‑151 (MQ=255) | CATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAA‑GGTTGCGGGAATCAGC‑AAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGG > NC_000913/3862518‑3862805 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |