| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,862,661 | 0 | T | G | 47.2% | 8.5 / 15.0 | 36 | 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 (6/13); new base G (16/1); total (22/14) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.37e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01 | |||||||||||
CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTTGCGGG‑AATCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGG > NC_000913/3862513‑3862810 | ccGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTTg < 2:298247/151‑1 (MQ=255) gCCCACCAAACCAAAGAGATACATTACGGTCGACATAGAGGTCGCCAGCACGGCGTGTAACGCAACCA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCAAGCGGGTCGGTAATTCCCACCAGCATGGC‑GGATAA‑AGtagtcggga > 1:813186/1‑138 (MQ=255) ccaccaCACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑AGGTGCGGG‑GATCa > 1:586147/1‑147 (MQ=255) accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACGGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑AGGTGCGGG‑GAAcaac‑aa > 2:290416/1‑144 (MQ=255) accaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGTGGGGG‑GATCAGC‑AAA‑cccccc > 2:57000/1‑151 (MQ=255) ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTTA‑GGGTTCGGG‑GATAAGC‑AAG‑ACcccc > 1:478252/1‑150 (MQ=255) agaTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑AGGTGGGGG‑GAACAGC‑AAG‑CCCGCCCC‑cttt > 2:481179/1‑150 (MQ=255) aTTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGGGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑AGGTGCGGG‑GAACAGC‑AAA‑ACCCCCAC‑CTTTTccc > 1:529362/1‑148 (MQ=255) actttGGGGCCGCCAGCACGGCGTGTACCGCAAAAA‑GAAAGGGAGAAATGAACAGAAAGTA‑AAATTCCGGCGGTTGGGTAATCCCCCCCAGCAGGGG‑GGTTAA‑GGTTGGGGA‑AATCAGA‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTg < 2:400249/147‑1 (MQ=255) cATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTTA‑GGGTGCGGG‑GAACCGCAAAC‑CCCCCCCC‑TTTTACCCCTTTTCCCgggc > 1:387899/1‑148 (MQ=255) ttGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGTGCGGG‑AATCAGC‑AAG‑CCCGCCCC‑TTTTTCCCCATTTTTCGGTGCgg > 1:216216/1‑150 (MQ=255) tgcggCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAGCAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAA‑GGTT‑CGGG‑GAGCAGC‑AAGCCCCCCCCC‑TTTTCCCCGATTTTcggtgggg > 1:278474/4‑142 (MQ=255) gccgccAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTTA‑AGTGGCGGG‑AATCAGC‑AAG‑CCCCCCCCTTTTTTCCCCTTTTTCCGGGGGGGcgggg > 1:156729/1‑147 (MQ=255) ccgccA‑CACGGGGGGTACCCGAAACAGGAAAGGGTAAAATAAACA‑GAAGGTAAAATCCCGGGGGTTGGGAAATCCCCCCCACCAGGGG‑TTTTAA‑GGTGGGGGA‑AATCAGC‑AAC‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGaa < 1:150700/151‑1 (MQ=255) ccAGCCGGGCGTGCCCCGAAAACA‑GCAGGGGTGAAATAAAAAGAAGGTT‑AATTTCCGGGGGTTCGGTATTTCCCCCCGGGATGGC‑GGTAAG‑GGTGGCGGA‑AATCGGC‑AGG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGa < 2:426511/147‑1 (MQ=255) cAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTGA‑AGGTGCGGG‑AAACAGC‑AAG‑ACCGCCAC‑TTTTATCCCATTTTTCCGGGCgggggtg > 1:792543/1‑140 (MQ=255) tgtACCGAAAACA‑GCAAGGGAAAAATAAACAGGAAGGT‑AATTCCCAGCGTTCGGGAAATTCCCCCCAGCATGG‑‑GGTTAAGGGTTGGGGA‑AATCAGA‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaaa < 2:156751/150‑1 (MQ=255) tgtACCCCAAACA‑GCAAGGGTAAAAGAAACGGGAAGGA‑AAATCCCAGCGTTTCGGTAATCCCCCCCAGCGTGGG‑GGTTAG‑GTTTGGGGA‑AATAAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTaa < 1:171668/149‑1 (MQ=255) ccGAAAACA‑GCAACGGTAAAATGAACAGGAAGGT‑AATTCCCAGCGGTTCGGTAATCCCCACCAGCAGGGC‑GGTTAA‑GGTTGGGGG‑AATCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGaa < 2:386922/149‑1 (MQ=255) cGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGTGCGGG‑AATCAGC‑AAA‑CCCGCCCC‑TTTTACCCCATTTTCCGGTGGGGCCGTGAAGTACATCGCTAAAGaaa > 2:509544/1‑149 (MQ=255) aCGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGTTCGGG‑GATCAGC‑AAG‑CCCGCCAC‑TTTTTCCCGATTTTCCGGTGCGGCAGTGAAGTACAACGCTAAAGAAATGCCCACGGCAc > 2:530461/1‑150 (MQ=255) aaTGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGTGCGGG‑AATCAGC‑AAA‑CCCGCCAC‑TTTTACCCGATTTTTCGGTGCGGCAGGGAAGGACATCGCTAAAGAAAAGCCCACGGCACCAAAGaat > 1:495924/1‑149 (MQ=255) gaaaGGT‑AAATCCCGGCGGTT‑GGGAATCCCCCACGGGATGGG‑AGTTAA‑GGGTGGGGGAAATCA‑C‑AAGACCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAttttg < 1:321710/144‑2 (MQ=255) aGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGTGCGGG‑AATCAGC‑CAG‑CCCCCCAC‑TTTTACCCCATTTTTCGGTGCGGGCGTGAAGTACATCGCTAAAGAAAAGCCCACCGCACCAAAAATTTTTTAGTTTcc > 1:579947/1‑151 (MQ=255) aaTTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATAA‑GGGTGCGGG‑AATTCGC‑AAG‑CCCGCCCC‑TTTTACCCGACTTTCCCGTGGGGCAGAGAAATACACCGCTAAAGAAAACCCCACGGCCCCCAAAAttttttaaatacccgg > 2:434025/1‑138 (MQ=255) ttCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGGTGCGGG‑GATCAGC‑AAG‑CCCGCCCC‑TTTTACCCGATTTTCCGGGGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCCCGGCACCAAAGATTTTTTAGTTACCGTGGAgg > 2:263220/1‑151 (MQ=255) tCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTTTCGGG‑AAACAGC‑CAG‑CCCCCCCC‑CTTTTCCCCATTTTCCGGTGCGGCAGGGAAGTACATCGCCAAAAAAAAGCCCCCCGCACCAAAGATTTTTTAGTTAccgggccggg > 1:537261/1‑151 (MQ=255) ccAGGGGTTCGGTATTCCCCCCCAGCAGGGC‑GTTTAA‑GGTTGCGGG‑AACCAGC‑AAC‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg < 2:579826/149‑1 (MQ=255) cGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAA‑GGGTGCGGG‑AATCAGC‑AAG‑CCCGCCCC‑TTTTACCCGATTTTCCGGGGCGGGAGGGAAGAAAATCGCTAAAGAAATGCCCACCGCACCAAAAAATTTTTAGTTACCATGCAGGGCaa > 2:104066/1‑149 (MQ=255) gttAATCCCCCCCAGCAGGGG‑GTTTAG‑GGTGGCGGA‑AACCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGcc < 2:9554/148‑1 (MQ=255) ccccccAGCGTGGG‑GGTAAA‑GGTTGGGGG‑AATAAGA‑AGG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCg < 2:30095/148‑1 (MQ=255) cccACCAGCAGGGC‑GGTTAA‑GGTTGGGGG‑AATCAGA‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGc < 2:218444/142‑1 (MQ=255) gCATGGC‑GGTTAA‑GGTTGCGGG‑AATCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACg > 2:282938/1‑150 (MQ=255) gtttAG‑GGTTGGGGG‑AATCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGc < 1:314269/149‑1 (MQ=255) ggTTAA‑GGTTGCGGG‑AATCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGc > 2:663410/1‑151 (MQ=255) ttAA‑GGTTGCGGA‑AACCACC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCa < 2:52484/147‑1 (MQ=255) gtttGCGGG‑AAGCAGC‑ATG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc < 1:37000/148‑1 (MQ=255) gTTGCGGG‑AATCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCgg > 1:484195/1‑151 (MQ=255) | CCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACA‑GCAACGGTGAAATGAACAGGAAGGT‑AAATTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAA‑GGTTGCGGG‑AATCAGC‑AAG‑CCCGCCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGG > NC_000913/3862513‑3862810 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |