Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | CP009273 | 2,820,787 | C→A | 11.2% | Q23K (CAG→AAG) | srlB → | glucitol/sorbitol‑specific enzyme IIA component of PTS |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | CP009273 | 2,820,787 | 0 | C | A | 11.2% | 92.5 / 2.1 | 32 | Q23K (CAG→AAG) | srlB | glucitol/sorbitol‑specific enzyme IIA component of PTS |
Reads supporting (aligned to +/- strand): ref base C (13/15); new base A (2/2); total (15/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01 |
CACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTC > CP009273/2820532‑2821038 | cACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑CCCACCATCACCCGTATCGGCGCGCGTCCCATTACCGCCCCCAGTACCCAGATGCTCATCCCCTTTc > 1:605399/1‑274 (MQ=255) gCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGgcg < 1:218720/275‑1 (MQ=255) gcggcCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAAGCTGAAACATGCCCGTTATTTATAAG‑ACCACCATCACCCATATCGGCGCGAGTGCCATTTACGCCCTCAGTGACCAGATGCCCATCACCTTTCTTGAAGGCGCGCCTGCGGACCTCg > 1:223492/1‑275 (MQ=255) gtctgcATCCCGGTCGGGGTGTCGCTGGCGGCAGCCCGTTAGGACACGGGTCGCGTCGTGGGCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAATGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTa < 2:223490/268‑1 (MQ=255) aGCCCGTGAGGACACGGTTCGCGTCGTTGTCCCTTCTGTACTGTTGAGCCGCTTGTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAAGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGt < 1:588668/275‑1 (MQ=255) gcgcgcCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGCCCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTgg > 1:417582/1‑275 (MQ=255) tCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTg > 1:137119/1‑189 (MQ=255) cTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAg < 2:137118/189‑1 (MQ=255) gTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGc > 2:657150/1‑170 (MQ=255) gTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACg < 1:657164/170‑1 (MQ=255) gTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACAAGATGCTCATCACCtt > 2:363727/1‑110 (MQ=255) gTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACAAGATGCTCATCACCtt > 2:268505/1‑110 (MQ=255) tatctatcAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACAAGATGCTCATCACCTTTCGt < 1:363734/110‑1 (MQ=255) tatctatcAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACAAGATGCTCATCACCTTTCGt < 1:268509/110‑1 (MQ=255) cAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGc > 1:418335/1‑148 (MQ=255) atggagGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTAACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCATGTCACCCTgcgc < 1:95796/272‑1 (MQ=255) agagGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACg < 2:418328/148‑1 (MQ=255) agGATGAAAAATGACATTTATTTACAGGAACCACC‑TCCCCATTACCGGAGCGAGTGCCCTTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGTCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCATGTCACCCTGCGCTTCg < 2:417575/275‑1 (MQ=255) gACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCAACGGCGAGTTGGAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTTGCAGCGTGGCGGAAGACAACCCTCGCGAAATGGGTCAGGGCTCCCGGCTCTGCGCTggggtctctg > 2:294948/1‑267 (MQ=255) ttatttatCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTTACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCGTGTCACCCTGCGCTCCGCTGGTTTAAACAAAGCg > 1:375398/1‑275 (MQ=255) ttatttatCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCg > 1:190810/1‑210 (MQ=255) ttatCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGc < 2:190808/210‑1 (MQ=255) caccaTCACCCGTATCGGCGCGAGTACCATTGACGCACTCAGTGACCAGACGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCGTCATTCACTGCCACCGCGAGTTGAAAGGTGCGATCCATCCCGGTTTGCACTTTTCACTCGGGCAGCACAGCTCTCCGGCGACCGCTGTGTGCAACGCAGCGGATGCCAACCTCCGCGAGCTGGGTCATGTCACCCGGGGGTTCGGTGGGTTAAACGAAGGGGACTGCCCGGGCa > 2:699804/1‑275 (MQ=255) tATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGgc > 2:269338/1‑152 (MQ=255) gggCGCGAGTTCCATTGACGCCCTCATTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCa < 1:369756/273‑1 (MQ=255) ggCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGgcagca < 1:269342/152‑1 (MQ=255) cgcgAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAag > 2:716684/1‑74 (MQ=255) cgcgAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGcc < 1:294953/275‑1 (MQ=255) aGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTa < 1:716698/74‑1 (MQ=255) tGCCATTGGCGCCCTCATTGACCAGATGCTAATCACCTGGCGTGACGGCGCGCCTGCGGACATCAAAGAGTATTGCTTAATTCATTGCCACGGGGAGTTGAAAGGTGCACTCCATCCCGGTATACAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTc < 1:688710/274‑1 (MQ=255) aGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGt > 1:35468/1‑82 (MQ=255) aCCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGaa < 2:35467/82‑1 (MQ=255) | CACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTATCAATAGAGGCTGAAACATGACCGTTATTTATCAG‑ACCACCATCACCCGTATCGGCGCGAGTGCCATTGACGCCCTCAGTGACCAGATGCTCATCACCTTTCGTGAAGGCGCGCCTGCGGACCTCGAAGAGTATTGCTTCATTCATTGCCACGGCGAGTTGAAAGGTGCACTCCATCCCGGTTTGCAATTTTCACTCGGGCAGCATCGCTATCCGGTGACCGCTGTTGGCAGCGTGGCGGAAGACAACCTTCGCGAACTGGGTCATGTCACCCTGCGCTTCGATGGTTTAAACGAAGCGGAATTTCCGGGCACTGTCCATGTGGCAGGCCCTGTC > CP009273/2820532‑2821038 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |