Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 2,173,361 | Δ2 bp | 100% | pseudogene (917‑918/1358 nt) | gatC ← | galactitol‑specific PTS enzyme IIC component |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,173,361 | 0 | C | . | 88.5% | 91.4 / 2.9 | 26 | pseudogene (918/1358 nt) | gatC | galactitol‑specific PTS enzyme IIC component |
Reads supporting (aligned to +/- strand): ref base C (1/1); new base . (9/14); total (11/15) | |||||||||||
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: Variant not supported by required number of reads on each strand. | |||||||||||
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch. | |||||||||||
* | NC_000913 | 2,173,362 | 0 | C | . | 88.5% | 86.9 / 4.2 | 26 | pseudogene (917/1358 nt) | gatC | galactitol‑specific PTS enzyme IIC component |
Reads supporting (aligned to +/- strand): ref base C (0/0); major base . (9/14); minor base A (2/0); total (11/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.83e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 | |||||||||||
Rejected as polymorphism: Variant not supported by required number of reads on each strand. |
GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT > NC_000913/2173239‑2173496 || gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTGCACGGCACACAAACAGCAATTAAAATGGTGCGTGGGATACCAGTCAGGCGTGCCGCTACGA‑‑CCGCCGGATGTCCCCgcagca > 1:114188/1‑143 (MQ=255) gACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGAGCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAgcagca < 1:273137/143‑1 (MQ=255) aCCGCCATCGCCACGAAGACGCCCAACATAGCATGTTGAGCAAACTCCAGCACCTGATTACCCGCCACCCAAACAGCAATTAAAATGGTGACTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCaa < 1:587918/143‑1 (MQ=255) tCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGa < 1:439316/143‑1 (MQ=255) cGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGAt < 1:497171/143‑1 (MQ=255) gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATc > 1:78715/1‑143 (MQ=255) gCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAAATAAAATGGCGAGTGGGACAAAAATCAGGCCTGCCGATCCCA‑‑CCGCCGCATGTCCCCGCAGCAACACCGGagc > 2:619139/1‑141 (MQ=255) gcGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAg < 1:403886/142‑1 (MQ=255) gCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGCACAAGGCCAATCa > 2:566549/1‑143 (MQ=255) gtggCAACATCGCCACACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTc < 2:49301/143‑1 (MQ=255) gCCAAACGGCAGCACCTGATTACACGGCACACAAACAGCAATTAAAATGGTGAGTGGGATACAAATCAGGCTTGCAGATACCA‑‑GCGGCGTATGTACCACCAGCAACGCCGGATCCAGGCCAACCAGCAACTCCTGACCGCCGa > 2:274861/1‑143 (MQ=255) cAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCACTCAGGAACTCCTGACCGCCGAAc > 2:468639/1‑143 (MQ=255) cGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATA‑‑ACACCGCCGTATGTCCCAGCCCCAACGCCGGCTCAAGGCCAATCAGGACCACCTCACCGCCGAACGTCg > 1:66469/1‑143 (MQ=255) cagcaCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCt < 2:531986/143‑1 (MQ=255) caCCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAACCAGGCTTGCCGATA‑‑AGACAGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTa > 1:425253/1‑143 (MQ=255) gttccccGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGAc < 2:159962/139‑1 (MQ=255) gcctccaAGCATCAATTAACATGGTGAGTGGGATAAAAATCAGGCTGGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGcc < 2:315001/138‑1 (MQ=255) aaaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTAGGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTa < 2:306213/143‑1 (MQ=255) aaCAGCAATTAACATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGACACGCCGGAGCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg < 1:391229/143‑1 (MQ=255) aaCAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAg < 1:544032/143‑1 (MQ=255) aTTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATggg > 1:336434/1‑143 (MQ=255) gAGTGGGATAAAAATCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccat < 1:15036/143‑1 (MQ=255) aGTGGGATAAAAATCAGGCTTGCCGATAC‑‑CGCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAAccatc < 2:288380/143‑1 (MQ=255) aaaTCAGGCTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGCTAAACCATCCATGATgggc > 2:44551/1‑142 (MQ=255) atcTTGCCGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGGCCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATc < 2:471374/141‑1 (MQ=255) ccGATACCA‑‑CCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGGGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGt > 1:130327/1‑143 (MQ=255) || GACCGCCATCGCCACGAAGAAGCCGATGGTGGCAAGATCGCCAAACGGCAGCACCTGATTACCCGGCACACAAACAGCAATTAAAATGGTGAGTGGGATAAAAATCAGGCTTGCCGATACCACCCCGCCGTATGTCCCAGCAGCAACGCCGGATCAAGGCCAATCAGGAACTCCTGACCGCCGAACTTCGCCTGTAAACGACTACGAGCCTGCTTAGCGATGGGCGTTAAACCATCCATGATGGGTTTAATCACCCGT > NC_000913/2173239‑2173496 |
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: ‑ |