Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,862,661 | 0 | T | G | 39.5% | 24.0 / 16.0 | 44 | 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 (10/16); new base G (15/2); total (26/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.55e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCG > NC_000913/3862513‑3862804 | ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGtt < 1:175226/150‑1 (MQ=255) ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGt < 1:704252/149‑1 (MQ=255) ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGt < 2:309005/149‑1 (MQ=255) cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGGAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGATTAGGtt > 1:254239/1‑149 (MQ=255) gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATc < 1:813908/151‑1 (MQ=255) ccaccaCACCAAAGAGATACATTACGGTCGACATTAAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGTTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCggggaaaa > 1:142197/1‑142 (MQ=255) cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGccc < 2:721764/151‑1 (MQ=255) cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAAccc > 1:77294/1‑151 (MQ=255) accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAGGTTGCGGGGATCAGCaaaaccccc > 2:682711/1‑144 (MQ=255) ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAACAGCAAACCCGcc > 1:37525/1‑150 (MQ=255) caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGGGGGGAACAGCAAAcccccccc > 1:429229/1‑149 (MQ=255) caAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGATCAGCAAGCCCGCCAc > 1:94278/1‑151 (MQ=255) agaTAAATTACGATCGACATGGGGGCCCCCAGCAGGGCGTGTACCGAAAACAGAAGCGGGGAAAGGAACAGGAAGGTAAATTCCGGCGGTTCGGTAATCCCCACCAGCATGGCGGTTAAGGTTGCGGAAACCAGCAAGCCCGCCACtttt < 2:9302/150‑1 (MQ=255) tACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTTCGGGGATTAGCAAGGCCGCCACTTTTACCCGCttt > 2:616208/1‑151 (MQ=255) aCGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGGAACAGCAAGCCCGCCACCTTTACCCGAattt > 1:800555/1‑151 (MQ=255) gACATTCAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGCGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGTTTCCGGGCCCAGCCAACCCCCCcccttttcccccttttcc > 2:598701/1‑147 (MQ=255) ttGCGGCCGCCAGCACGGCGCGTACCGCAAACAGCAACGGTGTAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGCTGCGGGAATTTGCAAGTGCGCCCATTTTTCCCGACTTTc > 2:57063/1‑142 (MQ=255) tggggCCGCCAGCACGGCGTTTCCCGCAAACAGCAAGGTTGAAAAGAAAAGGAAGGTAAATTCCGGCGGTTCGGTAATTCCCCCCAGCATGGCGGTTAGGGTGGCGGAAATCCGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGc < 2:507419/147‑1 (MQ=255) ccgcaACCACGGCGTGTACCGCAAACAGCAAGGGTAAAAGAAACGGGAAGTAAAATCCCAGGGGTTAGGTAACCCCCACCAGCTGGGCGGTTAAGGGGGCGGAAGTAAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGt < 2:642013/148‑1 (MQ=255) cgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGTGGGGATCAGCAAACCCCCCCCTTTTAACCCATTTTTCGGTGCGGCAggg > 1:120420/1‑146 (MQ=255) cGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAAACAGCAAGCCCCCCACTTTTACCCGATTTTCCGGCGCGGCAGGGGAGGACAACGc > 2:389297/1‑151 (MQ=255) cGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAAGGTGCGGGAATCAGCAAACCCGCCACTTTTTCCCCCTTTTCCGGTGCGGGAGGGGAGTa > 2:426750/1‑145 (MQ=255) ggggtgtCCCGAAAACAGAAACGGTGAAATAAACAGGAAGTTAAATCCCGGCGTTTCGGTAATCCCCCCCAGCATGGCGGTTAGGGTTGCGGGAACCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGt < 1:952968/140‑1 (MQ=255) gtACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAACCCGCCCCTTTTAACCGATTTTTCGGAGCGGCAGTGAAGTACa > 1:210803/1‑141 (MQ=255) cGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCCCTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGaaa > 1:734888/1‑149 (MQ=255) gCAAACAGCAACGGTGAAATGAAAAGGAAGATAAATTCCAGGGGTTCGGTAATCCCCACCAGCAGGGCGTTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAAt < 1:487197/149‑1 (MQ=255) aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGGACATCGCTAAAGAAATAccc > 1:173193/1‑151 (MQ=255) aaaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCCGGAAACAGGAAGCCCGCCACTTTTACCCGCTTTTCCGGTGCGGCAGGGAAGGACATCGCTAAAGAAAAGccc > 2:865768/1‑151 (MQ=255) aaaCAGAACCGGTAAAATGAACGGAAAGAAAAATCCCAGCGGTTCGGTAATCCCCACCAGCAGGGCGATTATGGTTGGGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATg < 1:198589/148‑1 (MQ=255) aaCAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTTCGGGAATCAGCAAGCCCGCCACTTTTTCCCGATTTTCCGGTGCGGCAGTGAAGGACATCGCTAAAGAAATGccc > 1:937104/1‑150 (MQ=255) gaaccGGCGAATTGAACAGAGAGGTAAATCCCGCCTGCCCGTGATTCCCCCCCTGCTTGGCAGTTAAGGGGGCGGAAAGCAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACg < 2:496530/145‑1 (MQ=255) cAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAAACAGCAAGCCCCCCACTTTTACCCGATTTTCCGGGGCGGCAGTGAAGTACATCGCTAAAAAAAAACCCACGGc > 2:144075/1‑150 (MQ=255) ggTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAGGTTGCGGGAAGCAGCAAACCCGCCCCTTTTTCCCCACTTTCCGGTGAGGCAGTGAAGTAAAACGCTAAAAAAATGACCCCGGCAcc > 2:951788/1‑149 (MQ=255) cAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGAACAGCAAACCCGCCACTTTTTCCCGATTTTTCGGTTGGGCAGCGAAGAACATCGCTAAAAGAAACCCCCCGGCACCAAAGAttttt > 1:599032/1‑148 (MQ=255) aGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTTAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCCGGGGAGTAAATCCCCTAAGAAATGCCCACCGCACCCAAGATTTTTTAGt > 1:626467/1‑151 (MQ=255) aaGGTAAATCCCAGGGGTTCGGTAATCCCCCCCAGCATGGGGGTTAAGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTg < 2:3470/145‑1 (MQ=255) ggTAAATCCCGGCGGTTCGGTAATCCCCCCCAGCTTGGCGTTTAGGGTGGCGGGAATCAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTa < 1:266592/148‑1 (MQ=255) aaaTTCCGGCGGTTCGTTAATCCCCACCAGCATGGCGTTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGc < 1:539026/151‑1 (MQ=255) ttCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGGATCAGCAAACCCGCCACTTTTACCCGATTTTTCGGGGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg > 1:786742/1‑151 (MQ=255) cgggggtttCGTTAATCCCCCCCAGCAGGGGGGTAAGGGTGGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg < 2:473063/142‑1 (MQ=255) ccAGCGGTCGGGTAATCCCCACAAGCATGGGGGTTATGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGc < 2:43448/151‑1 (MQ=255) cAGCGGTTCGGTAATCCCCACCAGCAGGGGGTTTAGGGTTGGGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg < 1:138400/149‑1 (MQ=255) ggTTCGGTAATTCCCACCAGCATGGCGGATAAGGGTGTGGGAATCAGCAAGACCCCCACTTTTAACCGATTTTCCGGGGCGGGAGTGTAGTACATCGCTAAAAAAATTGCCACCGCACCAAAGATTTTTGAATTACCGTGCAGGGCaaa > 1:363402/1‑149 (MQ=255) ttCGGTAATCCCCCCCAGCATGGCGGTTAGGGTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGCGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAAcc < 1:682987/150‑1 (MQ=255) tCGGTAATTCCCACCAGCATGGCGGTTAAGGGTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTTCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCCCGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc > 2:991542/1‑151 (MQ=255) caccaGGATGGGGGTTAAGGTTGGGGAAACCACAACGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCg < 2:975423/146‑1 (MQ=255) ccaGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGaaca > 1:708114/1‑150 (MQ=255) ttAAGTTTGCGGAAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCg < 2:37517/150‑1 (MQ=255) | CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGTTAAGGTTGCGGGAATCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCG > NC_000913/3862513‑3862804 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |