Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 3,285,696 | C→T | F73F (TTC→TTT) | agaD → | galactosamine‑specific PTS enzyme IID component |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,285,696 | 0 | C | T | 85.7% | 26.5 / ‑1.6 | 14 | F73F (TTC→TTT) | agaD | galactosamine‑specific PTS enzyme IID component |
Reads supporting (aligned to +/- strand): ref base C (1/1); new base T (5/7); total (6/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01 |
CGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGG‑AATTTGCTCAT‑CATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCT‑‑ATTTTTCGCCGTTTACCTGCTTA > NC_000913/3285527‑3285907 | cgtcgCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGTCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCTTATTCCTGATGTTTTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTgg < 1:506438/250‑1 (MQ=255) cAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCGGGTCGGATTTCTGATGGGGTTATTAATGTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCgg < 1:423458/251‑1 (MQ=255) cGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACATGGGCGATGTTGACGATCCTGAAAAAGCTTTATAAGGAAGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGCATTTATTAATACCCACCCGACTCTGGTCGGATTTCTGATGGGGGTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACCCCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGAtt > 1:337285/1‑251 (MQ=255) ccTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAAGTTAGTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTc < 1:228904/136‑1 (MQ=255) ttGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATTGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGc < 1:434680/176‑1 (MQ=255) cTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGcc > 1:216273/1‑178 (MQ=255) ataAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTAGTAATACCCACCCGAAGCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCTCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTGCTGGTTTACTTTGTTGCCGAGTATGGCTGG‑AATTTGCTCAT‑CATTCGCCAGCCAggg < 1:167447/250‑1 (MQ=255) gacAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGG‑GATTTGCTCAT‑CATTTTCCAGCCCggggaaaacgct > 1:220978/1‑242 (MQ=255) cgATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTAc > 1:443522/1‑172 (MQ=255) aaaGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGG‑AATTTGCTCAT‑CATTTGccagccag < 1:294021/213‑1 (MQ=255) ccTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGATATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGGGAATTT‑TTCATTCCTTTTCC‑GCCAGGGGACCCTGCCGGGGGCGATTCTTTATTTTTC‑CCGTTT‑CCTGCttt > 1:198124/1‑250 (MQ=255) cGAATTTCTTAATACCCCCCCGAATCTGGTCGGATTTCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGAACTGTTTGGccc > 1:153844/1‑120 (MQ=255) cGAATTTATTAATACCCACCCGAATCTGGTCGGATTTCTGATGGGGTTATTAAGTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTTTTTGGCCCAATCGCCGGTATTGGCGATGCGATTTTCTGGt < 1:426793/152‑1 (MQ=255) gTCGGATTTCTGATGGGTTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGgcgatgc < 1:188540/113‑1 (MQ=255) | CGTCGCTGCTGCAAGCGAGCTTTAACTACGAAAGGATGCAGGCGGGCGGTTTTACCTGGGCGATGTTGCCGATCCTGAAAAAGATTTATAAGGACGACAAACCGGGCTTAAGCGCGGCGATGAAAGATAACCTCGAATTTATTAATACCCACCCGAATCTGGTCGGATTCCTGATGGGGTTATTAATTTCGATGGAAGAAAAAGGAGAAAACCGCGACACCATTAAAGGCCTCAAAGTGGCACTGTTTGGCCCAATCGCCGGGATTGGCGATGCGATTTTCTGGTTTACTTTGTTGCCGATTATGGCGGG‑AATTTGCTCAT‑CATTTGCCAGCCAGGGAAACCTGCTGGGGCCGATTCT‑‑ATTTTTCGCCGTTTACCTGCTTA > NC_000913/3285527‑3285907 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |