Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,862,664 | 0 | C | G | 42.0% | 27.3 / 22.0 | 50 | pseudogene (940/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 C (15/14); new base G (4/17); total (19/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.74e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01 |
CCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG > NC_000913/3862516‑3862809 | cccATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCg < 1:784174/150‑1 (MQ=255) cATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑a < 2:627102/151‑1 (MQ=255) gTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑AAT‑c < 1:709203/151‑1 (MQ=255) cccACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGCGGG‑TAAA‑GGTTGgggggacag > 1:449533/1‑144 (MQ=255) caccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑GTAA‑AGGTGGGGG‑GAT‑AAACaa > 1:692395/1‑150 (MQ=255) ccacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑GTAA‑GGGTGCGGG‑GAT‑CAGCAAGcc > 2:343203/1‑151 (MQ=255) cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑CTAA‑GGGTGCGGG‑AAA‑CAGCAAGccc > 2:1098991/1‑151 (MQ=255) caccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑AAT‑CAGCAAGCCCGc < 1:534791/151‑1 (MQ=255) ccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑AAT‑CAGCAAGCCCGCCa < 2:1043178/151‑1 (MQ=255) aTACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTTA‑AGTTGGGGG‑GAT‑CAACAAGCCCCCCACTTTTTcc > 2:67330/1‑151 (MQ=255) tACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGGTGGGGG‑AAT‑CAGCAAGCCCCCCACtttttacc > 1:865111/1‑147 (MQ=255) aCATTACGGTCGACATGGAGGCCGCCAGCAGGGCGTTTCCCGAAAACAGCAACGGTAAAAGGAACAGGAAGGTAAATC‑CCAGCGGTTCGGTAATCCCCACCA‑GCATGGG‑GG‑TTAG‑GGTGGCGGG‑AAT‑CAGCAAGCCCGCCACTTTTACCCg < 1:557138/151‑1 (MQ=255) ggccgacttggggcccgccAGCGCGGCGTGTCCCGAAAACAGAAAGGGAAAAATGAACAGAAAGGTAAATC‑CCGGCGGTTGGTTAATCCCCCCCA‑GCAGGGG‑GT‑TTAG‑GGTGGGGGA‑AAT‑CAGAAAGCCCGCCACTTTTACCCGAtttt < 1:1081950/135‑1 (MQ=255) gTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑AGGTGCGGG‑GAT‑CAGCAAACCCCCCCCTTTTTCCCCAtttttc > 2:923365/1‑148 (MQ=255) tggggcccgccAGCGGGGCGTTTCCCCCAACCAGAACGGGTAAAAGAAACAGAAAGGAAAATC‑CGGGCGTTTCGGTATTCCCCCCCA‑GCAGGGG‑GT‑TTAG‑GGTGGGGGA‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCgg < 1:1105189/139‑1 (MQ=255) gccgccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑AGGTGCGGG‑GAA‑CAGCAAACCCCCCCCTTTTTCCCGATTTTTCGGGGCGGGAggg > 2:1027256/1‑148 (MQ=255) ccgccAGCCCGGGGTGCCCCAAAAACAGAACCGGTAAAATAAACAGGAAGGAAAATC‑CCGGCGGTTCGGTAATCCCCACCA‑GCATGGC‑GG‑TTAA‑GGTGGCGGA‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGa < 2:186612/150‑1 (MQ=255) gccAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑GTAA‑AGTTTCGGG‑GAA‑CAACAAAACCGCCCCTTTTACCCGATTTTCCCGGGCGGCCGTg > 2:209899/1‑147 (MQ=255) caaGCAGGGCGTGACCCAAAACCAGCAAGGGTAAAATGAAAGGGAAGTTAAATT‑CCAGCGTTCCGATAATCCCCCCCA‑GCGTGGG‑GG‑TGAG‑GGTGGGGGA‑AAT‑CAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCgg < 1:1038736/139‑1 (MQ=255) aGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGGTGGGGG‑GAA‑CAACAAACCCCCCCCCTTTATCCCATTTTTCCGTGCGGCAgag > 2:390423/1‑142 (MQ=255) cGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑AGGTGCGGG‑GAT‑CAACAAGCCCCCCCCCTTTTCCCCATTTTTCCGTGCGGGAGGGGAGTACa > 2:1102695/1‑147 (MQ=255) ggggtgtCCCACAAAAAGCAACGGAAAAAGGAACAGGAAGAAAAATC‑CCGGCGGTTCGGAAATCCCCCCCA‑GCAGGGC‑GG‑TTAG‑GGTGGCGGG‑AAT‑CAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGt < 1:407948/135‑1 (MQ=255) ggggtgtACCGAAAACAGCAACGGTGAAATGAACAGGAAGTTAAATT‑CCAGCGGTTCGGTATTCCCCACCA‑GCAGGGG‑GT‑TTAG‑GGTTGGGGG‑AAT‑CAGCAGGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTa < 2:964679/141‑1 (MQ=255) ggggtaacccacAAAAAGAAAGGGTAAAATAAACAGAGAGGTAACTG‑CCCGCGTTTCGGT‑ATTCCCCCCACGTAGGGG‑GG‑TAAA‑GGTTGGGGC‑AAG‑CAGCAACTCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCg < 1:226979/138‑1 (MQ=255) cAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑AGTTGCGGG‑GAT‑CAGCAAACCCCCCCCTTTTTCCCCATTTTTCGGTGCGGGCGGGAAGGACATCGCTAAAGAAATTcc > 2:598556/1‑151 (MQ=255) aaggggaaaaTAAACAGGAAGGAAAATC‑CCGGGGTTGCGGTATTCCCCCCCA‑GCATGGC‑GG‑TTAG‑GGTGGGGGA‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGc < 2:449547/142‑1 (MQ=255) aagggTAAAATAAACAGGAGGGAAAATC‑CCGGCGTTTCGGTAATTCCCCCCA‑GCATGGC‑GT‑TTAG‑GGTTGCGGA‑AAT‑CAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGc < 2:748959/146‑1 (MQ=255) aCGGTAAAAGAAACAGAAAGGTAATTC‑CCAGGGGTTGGTAAATCCCCCCCG‑GCAGGGC‑GG‑TTAG‑GGTTGGGGG‑AAT‑CAGAAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGc < 1:622765/148‑1 (MQ=255) gggtaaaaaaaaagggAAGAAAAATT‑CCGGCGGTTCGTTAATCCCCCCCA‑GCAGGGC‑GG‑TTAG‑GGTGGGGGA‑AAG‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCAc < 1:1118047/135‑1 (MQ=255) aaaaGACCAGGAAGATAACTT‑CCGGTGTTTTGGTACTCCCCACCG‑GTGTGGC‑GG‑TAAA‑GTGTGGGGC‑AAT‑CACCATGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAt < 2:99586/151‑1 (MQ=255) aaaTGAACAGGAAGGTAAATC‑CCAGCGGTTCGTTAATCCCCACCA‑GCAGGGC‑GG‑TTAA‑GGTGGCGGA‑AAT‑CAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCaaa < 1:499886/148‑1 (MQ=255) aaaGGGAGGTTAATTCCCAGCGTTCCGTTAATCCCCACCA‑GC‑GGGG‑GG‑TTAAGGGTGGGGGG‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGAtttt < 1:68153/147‑1 (MQ=255) ggAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑AGGTGCGGG‑GAT‑CAGCAAGCCCGCCACTTTTACCCCATTTTTCGGTGCGGCAGGGAAGGACATCGCTAAAAAAAAGCCCACCGCACCCAAGATTTTTTAGtt > 2:826107/1‑151 (MQ=255) aGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑GAT‑CAGCAAGCCCCCCACTTTTAACCGAATTTTCCGTGCGGCAGTGGAGAACATCCCTAAAGAAAAGCCCACCGCACCAAAGAATTTTTAGGTa > 2:652247/1‑149 (MQ=255) aaaTT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑AAT‑CAGCAAACCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTTAGTTACCGt > 2:933342/1‑149 (MQ=255) aTT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTTA‑GGTTGCGGG‑AAT‑CAGCAAGCCCCCCACTTTTTCCCGATTTTTCGGGGCGGCAGTGAAGGACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTTAGTTACCGTg > 2:581191/1‑148 (MQ=255) t‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑GAT‑CAGCAAGCCCCCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAAATCGCTAAAGAAAAACCCACGGCCCCAAAGATTTTTTAGTTACCCTGCAagg > 2:1163898/1‑151 (MQ=255) cccagcgataggttAATTCCCACCA‑GCATGGG‑GG‑TAAA‑GGTAGCGGG‑AAA‑CAGAAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAgg < 2:934339/137‑1 (MQ=255) aGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑AAA‑CAGCAAACCCGCCACTTTTACCCGATTTTCCGGGGCGGCAGGGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg > 2:486575/1‑148 (MQ=255) ggTTCGTTAATCCCCCCCA‑GC‑TGGG‑GG‑TTTA‑GGGTGGGGGAAAC‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAggg < 1:923309/145‑1 (MQ=255) aacccccACCA‑ACAGGGC‑GT‑TAAA‑GGTTGGGGA‑AAA‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCtt < 1:1188196/145‑1 (MQ=255) ccccccA‑GCAGGGG‑GT‑TTAG‑GGTTGGGGA‑AAC‑CAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGa < 1:511458/151‑1 (MQ=255) ccccccA‑GCA‑GGC‑GGTTTAG‑GGTTGGGGA‑ACC‑CAACACCCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCtt < 1:67325/145‑1 (MQ=255) ccccccA‑GAATGGG‑GG‑TTAA‑GGTTGCGGG‑AAA‑CACGAGATCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCgg < 2:204951/149‑1 (MQ=255) gCAGGGG‑GG‑TTAG‑GGTGGCGGG‑AAC‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACg < 2:528215/150‑1 (MQ=255) ggC‑GG‑TTAA‑GGTTGCGGG‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAg > 1:696385/1‑120 (MQ=255) ggggg‑TTAG‑GGTGGGGGA‑AAT‑AAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGAttt < 2:630462/147‑1 (MQ=255) gggggtAA‑GGTTGGGGG‑AATCCAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTc < 1:812280/146‑1 (MQ=255) gcg‑TTGA‑GGGTGGGGGAAAT‑CAACAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTc < 2:3888/148‑1 (MQ=255) aa‑GGTTGCGGG‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGc > 2:200975/1‑150 (MQ=255) ggTTGCGGA‑AAT‑CAGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCg < 1:821480/151‑1 (MQ=255) cGGG‑AAG‑CCGCAACCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc < 1:978606/145‑1 (MQ=255) | CCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAATT‑CCAGCGGTTCGGTAATTCCCACCA‑GCATGGC‑GG‑TTAA‑GGTTGCGGG‑AAT‑CAGCAAGCCCGCCACTTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCG > NC_000913/3862516‑3862809 |
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: ‑ |