| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913_3_bme_pgi | 2,173,361 | Δ2 bp | 100% | pseudogene (1‑2/442 nt) | gatC ← | pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913_3_bme_pgi | 2,173,361 | 0 | C | . | 90.7% | 157.6 / 8.6 | 43 | pseudogene (2/442 nt) | gatC | pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC |
| Reads supporting (aligned to +/- strand): ref base C (2/2); new base . (21/18); total (23/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
| * | NC_000913_3_bme_pgi | 2,173,362 | 0 | C | . | 88.4% | 149.6 / 9.3 | 43 | pseudogene (1/442 nt) | gatC | pseudogene, galactitol‑specific enzyme IIC component of PTS;transport; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system galactitol‑specific enzyme IIC |
| Reads supporting (aligned to +/- strand): ref base C (3/0); new base . (20/18); total (23/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.43e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01 | |||||||||||
| Rejected as polymorphism: Variant not supported by required number of reads on each strand. | |||||||||||
| Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
CGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTA > NC_000913_3_bme_pgi/2173213‑2173503 || cGGATCAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccac < 1:203670/149‑1 (MQ=255) aCAGATTTCCACGATGCACGGCGACCGCCATCGCCACTAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACCCAAACAGCAATAAAAATGGTGAGTGGGATACAAATAAGACTGGCGGATaccaccgcc > 2:55383/1‑146 (MQ=255) aTTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATaccaccggagaat > 2:898547/1‑142 (MQ=255) tttCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATg > 1:898446/1‑149 (MQ=255) cACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CCCCGCCGTATGTccc > 2:176344/1‑149 (MQ=255) cGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAAc < 2:506020/149‑1 (MQ=255) gaacGCCATCGCCACGAAGAAGCCGATGGTAGCAAAATAGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGAAACTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCTTATGTCCCAGCAGCAACGCCg < 2:69863/146‑1 (MQ=255) gCCATCGCCACGAAGAAGCCGATTGTGGCAAGATCGCCAAACGGAAGCAATTTATTACCCGGCACACAAACAGCAATTAAAATGGTGTGTGGTATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc < 1:752370/149‑1 (MQ=255) ccATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCa > 2:377554/1‑149 (MQ=255) cATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCaa < 2:612430/149‑1 (MQ=255) gtCGCCAGGAAGAAGCCGATGGTGGCAAGATCCCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg < 2:348958/148‑1 (MQ=255) agcGCCACAAAGAAGGCGATGGTGTCAAGATCTCTAAACGGCAGAACCTGATTACCCGGCACACAAACAGCAATGAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg < 1:767774/147‑1 (MQ=255) aTCGCCACGAAGAAGCCTATGGTGGCAAGACCGCCAAACGGCATCACCTGTTTACCCGTCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg < 1:552480/149‑1 (MQ=255) aCGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACTTGATTACCCGGCACACAAACAGCAAGTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAAt > 2:188296/1‑149 (MQ=255) cGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAGGCAGCAGCGCCGGCTCAAGGCCAAtg > 2:734239/1‑148 (MQ=255) aagaagCTGATGGTGGCAAGATTGCCAAATGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg < 1:359922/149‑1 (MQ=255) aagaagCTGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAg < 1:809039/149‑1 (MQ=255) aagCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CAGCCGTAGGTCCCAGCAGCAACGACGGAGCAAGGCCAATCAGGaa > 2:349715/1‑149 (MQ=255) ccGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAATGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCAAGCAGCAACGCCAGATCAAGGCCAATCAGGAAATc > 1:479895/1‑149 (MQ=255) ggCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTTCCAGCAGCGACGTGGGATCAAGGCCAATCAGGAATTCTTGAccgc > 1:137808/1‑149 (MQ=255) gCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAAGAGCACTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCACCGTATGTCCCAGCACCAACGCCTGACCAAGGCCAACCAGCCACTCCTTAccccc > 2:38825/1‑149 (MQ=255) gCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGAAGCAACGACGGATCAAGGCCAATCAAGAACTCCTGACCGCCGAACTTCg > 1:82580/1‑149 (MQ=255) cAAACGGCAGCACCTGATTACCTGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGcc < 1:362630/149‑1 (MQ=255) aaCGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAAACAAGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTg > 2:576897/1‑149 (MQ=255) cTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATA‑CA‑CCCACCGTATGTCCCAGCAGCAACGGCGGGTCAAGGCACATCAGGACCTCCTCACAGCCGAACTACCACTGTCAActtccgc > 2:248026/1‑143 (MQ=255) gATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGa < 1:624617/149‑1 (MQ=255) aTTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAg > 1:747175/1‑149 (MQ=255) aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc > 1:332627/1‑127 (MQ=255) aCCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTc < 2:332627/127‑1 (MQ=255) aaTTAAAATGGTGAGTGGGATAAAAATAAGGCTTGCCGATACCA‑‑CCACCGTATGTCCCAGCAGAAACGCCGGATCAAGGCCAATCCCGAACTCCTGACCGCCGAACTTCGCCTATAAACGACTGCGACCCTGCCTCACGATGGCCCTTa > 2:173958/1‑149 (MQ=255) aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTGCCAGCAGCAACGCCGGATCACGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGACCCTCCTTCGCGACGGGCGTTaa > 1:646874/1‑149 (MQ=255) aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCAAGAAAAAGACCAATCAAGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaa > 1:10416/1‑149 (MQ=255) ggaaaaTGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCC‑‑ACGTATGTCCCAGCAACAACGCCAGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTAGACTGTAAGCGACTACAGGCCTGCTTAGAGATAGGAGTTaga > 2:286009/3‑147 (MQ=255) atagaaTGGTGATTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTaaa < 1:199100/145‑1 (MQ=255) gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatga < 1:898547/149‑1 (MQ=255) gTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCatgatg < 1:188296/149‑1 (MQ=255) gggATCAAAATCCGCCGTGCCGATACCA‑‑CCGCTTTATGACCCACCAGCAACGGCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATggg < 1:841051/149‑1 (MQ=255) ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt > 1:629492/1‑149 (MQ=255) ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt < 2:747175/149‑1 (MQ=255) ggATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGt > 2:760705/1‑149 (MQ=255) gATAAAAATCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGtt < 1:751185/149‑1 (MQ=255) tAAAAATCAGGCTTGCCGATACCA‑‑CCACCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTGCTGACCGCCGAACTTCGCCTTTAAACGACTACAAGCCTGCTTAGCCATGGGCGTTAAACCATCCATAATTGGttga > 2:230197/1‑147 (MQ=255) aaaaTCAGGCTTGCCGATAC‑‑CTCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAAt < 1:687624/149‑1 (MQ=255) aTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCAc > 1:334398/1‑149 (MQ=255) gCTTGGCGATTCCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTg < 2:390527/149‑1 (MQ=255) cGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGACAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGCGTTGAATCACCCGTGGCAATa > 2:880880/1‑149 (MQ=255) || CGGAACAGATTTCCACGATGCACGGCGACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGTGGCATTA > NC_000913_3_bme_pgi/2173213‑2173503 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |