| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,862,667 | 0 | G | A | 42.3% | 24.7 / 10.3 | 27 | pseudogene (937/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 G (14/1); new base A (0/11); total (14/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.55e-06 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGG‑CCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGGAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAG > NC_000913/3862481‑3862899 | aGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGG‑CCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑AGGT‑GCGGGGATTAGCAAACCCGCCACtt > 1:972008/1‑207 (MQ=255) atatttcccccACCAACCAAAAGAGATACATTCGGGTCGACAGGGAGC‑CCGCCAGCCGGGGGTGTCCCAAAAACAGCAAGGGTAAAATGAACAGGAAGATAAA‑TTCCCGCGTTTCGGTAATC‑CCCCCCAGCAGGGG‑GGTTAG‑GGTG‑GGGGAAATCAACAACCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTAc < 2:706168/195‑1 (MQ=255) aTGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGG‑CCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑AGGT‑GTGGGGATCAGCAAGCCCGCCACTTTTAACCCGATTTT‑CGGGG‑AGGCAGGGAAGGACATCGCTAACGAAAAGGCCAcag > 1:148234/1‑221 (MQ=255) cccAACACCCAAAAGGGATACAATACGGTCAGAATTGAGGCCCGCCA‑CCCGGGGTTTCCCCAAAACAAGAAAGGATAAAATAACAGGGAAGTAAATTTCCGGGGGTTCGGTAATC‑CCCCCCGGCAGGGG‑GGTAAG‑GGTG‑GGGGAAACCAGCAACCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCa < 2:273614/215‑1 (MQ=255) caccacaccaAAGAGATACATTACGGTCGACATTGAGG‑CCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑AGGT‑GCGGGGATCAGCAAACCCGCCACCTTT‑ACCCCATTTTTCCGGG‑CGGCAGGGAAGGAAATCGCTAAAg > 1:706168/1‑203 (MQ=255) caAAGAGATACATTACGGTCGACATTGAGG‑CCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑GGGT‑GCGGGGAACAGCAAACCCCCCACTTTT‑TCCCGATTTTCCCGTG‑CGGCAGTGAAGTAAAACGCTAAAGAAAAGCCCACGGCACCcaa > 1:273614/1‑214 (MQ=255) caAAGAGATACATTACGGTCGACATTGAGG‑CCGCCAGCACGGCGTGTACCGCAAACACCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCCGCGGTTCGGTAATT‑CCCACCAGCATGGCGGGTTAG‑GGTT‑CGGGGAATAACAAAGCCCCCCCCTTTT‑ACCC‑ATTTTTCTGTGCCGGCAGTGAGGTAAACCGCGAAAAAAAGCCCCCCGCCACCAAAGATCTTTAAGTTAACG‑C‑CAGGGCa > 2:457275/1‑239 (MQ=255) cGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGGTAA‑GGGT‑GCGGGGAACAGCAAACCCGCCACCTTT‑ATCCCATTTTTCCGGG‑CGGCAGGGAAGGACAACCCTAAAGAAAAACCCACCGCACCCAAGAATTTTTAGGTACCGTGCCAGGGAAAAACC‑CCTTCCCGGGACAACGAATTCAACGGcctcgcatct > 2:823447/1‑230 (MQ=255) tgtACCGAAAACAGAAACGGTGAAATGAACAGGAGGGTAAA‑TCCCAGCGGTTCGTTAATC‑CCCACCAGCATGGG‑GGTTAA‑GGTT‑GCGGGAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑Cgg < 1:1031/130‑1 (MQ=255) aCAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑GGTT‑TCGGGGATTAGCAAACCCGCCACCTTT‑ACCCGATTTTCCGGTG‑GGGCTGGGGAGTa > 2:1031/1‑130 (MQ=255) tacaaaGGTAGGAAAA‑TTCCAGGGGTTCGTAAATT‑CCCACCAGTGTGGC‑GTTTAA‑GGTT‑GCGGAAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGcc < 1:198992/234‑1 (MQ=255) aaaaggaaggaAAA‑TCCCGGGGGTTCGGTAATC‑CCCCCCAGCATGGC‑GGTTAG‑GGTG‑GGGGAAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCa < 1:990553/130‑1 (MQ=255) tAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑GGGT‑GCGGGGATCAGCAAACCCCCCACCTTT‑ACCCCGTTTTTCGGTG‑CGGGAGGGGAGGACAACGCTAAAAAAAATCCCACGGCCCCaaa > 2:990553/1‑133 (MQ=255) ttCCAGCGGTTCGGTAATC‑CCCCCCAGCATGGC‑GGTTAG‑GGTT‑GCGGAAACCAGCAACCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTAc < 2:783078/100‑1 (MQ=255) tcccggggTTTGGTTATTCCCCCCCAGC‑TGGG‑GGTTAAGGGTT‑GGGGAAACCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGAtgct < 1:497493/212‑1 (MQ=255) tCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGGGATCAGCAAGCCCCCCACCTTT‑TCCCGATTTTCCGGGG‑CGGCAGGGAAGTACAACGCCAAAAAAATGCCCACCGCGCCCAAGATTTTTGTGTTTCCGTG‑CAGGGCAAAAACGCCCCCCCCGGAAAAAGATTTTATCGGGTCTGGGCGCACAATGAACACTTTGCA‑ATGCGCGGCCCAGGACATCTGAATtcc > 2:293196/1‑239 (MQ=255) ccgggggTTGGGTACTC‑CCCACCAGCATGCC‑GGTAAG‑GGTT‑GCGGAAATCAGCAACCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGCCCAGTACATCTGAATgccg < 1:352445/234‑1 (MQ=255) tCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGGGATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGGGAAGTACATCGCTAAAg > 1:783078/1‑100 (MQ=255) caccaGCTGGGC‑GGTGAA‑GGTT‑GCGGTAAGCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc < 1:293196/239‑1 (MQ=255) ccagttgggC‑GTTTAA‑GTTT‑GCGGAAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCa < 1:460919/232‑1 (MQ=255) cATGGC‑GGTTAA‑GGTT‑GCGGGAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg > 1:888728/1‑120 (MQ=255) cATGGC‑GGTTAA‑GGGT‑GCGGGAATTAGCAAGCCCGCCCCTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGaa > 1:764527/1‑69 (MQ=255) ggC‑GGTTAA‑GGTT‑GCGGGAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGGG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTg > 1:273898/1‑117 (MQ=255) ggTTAA‑GGTTGGCGGAAATCAACAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa < 1:552331/238‑1 (MQ=255) t‑GCGGAAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAg < 1:7700/238‑1 (MQ=255) cGGGAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTaaa > 1:461597/1‑65 (MQ=255) cGGAAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTaaa < 2:461597/65‑1 (MQ=255) | AGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGG‑CCGCCAGCACGGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATT‑CCCACCAGCATGGC‑GGTTAA‑GGTT‑GCGGGAATCAGCAAGCCCGCCACTTTT‑ACCCGATTTTCCGGTG‑CGGCAGTGAAGTACATCGCTAAAGAAATGCCCACGGCACCAAAGATTTTTGAGTTACCGTG‑CAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTC‑TTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAACTGTCCGTAG > NC_000913/3862481‑3862899 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |