| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,152,494 | G→T | A87S (GCG→TCG) | metC → | cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,152,494 | 0 | G | T | 93.1% | 87.0 / ‑5.5 | 29 | A87S (GCG→TCG) | metC | cystathionine beta‑lyase/L‑cysteine desulfhydrase/alanine racemase |
| Reads supporting (aligned to +/- strand): ref base G (0/2); new base T (21/6); total (21/8) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.90e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTC > NC_000913/3152274‑3152734 | aGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCtt > 1:142099/1‑251 (MQ=255) gACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCGTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCtttt < 2:264625/251‑1 (MQ=255) cagaaaaaaaTACACTCTCGGCGCGGTAAATAGCGTGAGTCGGCTCGCTTCTTCGCTGGTTGTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTGGGTCTAGGGACGGCGCGGAACGTTAACCCATTGCTCCTTACAACAAGCGATGTTTGAACTGG‑AAGGTGGCGCAGGCGGCTTGCTAGTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCg < 1:299666/247‑1 (MQ=255) ccGGCGCGGTAAATAGCGTGAGTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGGAGAAGGCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCAtgtg < 1:98259/250‑1 (MQ=255) ggTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTtgatg > 1:12031/1‑250 (MQ=255) tGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGcc > 1:186585/1‑251 (MQ=255) ttCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGcc < 2:182625/248‑1 (MQ=255) ttCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGcc > 1:182625/1‑248 (MQ=255) tGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTAt > 1:157490/1‑233 (MQ=255) tGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGc > 1:189839/1‑250 (MQ=255) aCAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCaa > 2:166358/1‑250 (MQ=255) ggtAGAATCCAAGAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTATATGGACGGCGCGGAACGTACACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACAAACACAGCCTATGAACCGAGTCAGGATTTCTGTAGCAACATcc > 2:34150/2‑251 (MQ=255) aaCACGCGACACGTAATCGCGCCAATGGAAAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTTCTATTTCCCTGCGGGTCGGCAGCGGTGGCTAATTCCATTATTCCTTTTATTGAACCGGGCGATCATTTGTTGATTACCCTAACCGCCTATGAACCGAATCCGGATTACTGTCGCAAAATCCtacggaacatgggc > 2:1793/1‑238 (MQ=255) gTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc > 1:177290/1‑250 (MQ=255) aaTGGAGCGTTGTTCTCTGGACGGAGCGAAACGTCAACCCATTTCTCCTTCCAACAAGCGATGTGTGAACAGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGATAGTTCCATTATTGCTTTTATCGCACAGGACGATCATGAGTTGATGACCAACACCGCCTATGAACAGAGTCAGGATTTCTGTAACAAACTCCTCAGCACACTGTGCGTAAAGCCAGCATGGATTGATcc > 2:28160/1‑251 (MQ=255) ggAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACGAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt > 1:217461/1‑251 (MQ=255) ggAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCt > 1:175715/1‑251 (MQ=255) ggAGAGTTATTCTATGGACGGCGCGGAACGTTAACCAATTTCTCCATACAACCAGCGATGTGTGAACTGG‑ATGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTCTGAACCGAGGCAGGATTTCTGTAGCAAAATCCTCAGACAACTGGGCGTAACGACATCTTGGATTGATCCGCt > 2:137720/1‑251 (MQ=255) tCTATGGAAGGAGCGGAACGTTAACCAATTTCTCCTTACAACACGCGATGTGTGCACTGG‑CAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc > 2:31667/1‑251 (MQ=255) cGGCGCGGACCGTTACCCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGTTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTTTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAg < 1:229228/185‑1 (MQ=255) cGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAg > 2:229228/1‑185 (MQ=255) ggAACTTTAACCCTTTCCCCCTTCCAACAAGCGTTTGGTGACCTGG‑AGGGGGGCGCAGGCGGCGTGCTATTCCCCTGCGGGTCGGCAGCGGTTGCTAATCCCATTCTTGCTTTTATCGAAGAGGGCGTGCTTGGGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCa < 1:188565/251‑1 (MQ=255) ggAACGTGAACCCATTCCCCCGTACAACAAGCGATGTGTGAACGGG‑AGGGGGGCGCGGGCTGCGTGTTTTTCCCCGGCGGGTGGGCAGCGGTTGCTAATTCCTTTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACTCCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGTCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATAGCGTTAAGCa < 2:175715/251‑1 (MQ=255) gAACGTTACCCCATTCCCCCTTAAAACAAGCGATGTGTGA‑CTGGGAGGGGGGCGGGGGTTGCGTGTTATTCCCCGGCGGGGCGGCAGCGTTTGCTATTCCCATCCTTGCTTTTACCGAACTGGGCGAACGTGTTTTGAGGCCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGTGCCGCTGATTGGTGCCGATATCGTTAAGCAt < 1:278269/251‑1 (MQ=255) tAACCCATTTCTACTTACAACAAGCGCTGTGTGAACTGG‑AAGGTGGCGCAGGCTACGTGCTATTTCCCAGCGGGTCGGCAGCGGCTGCTAATTCCCTTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGCTGACCAACACCGCCTATGAACCGAGTCAGGCTTTCTGTAGCAAAATCCTCAGCCAACTGGGCGTAACGACATCATGGTTTGATCCGCTGTTTGACGCCGCTATCGTTAAGCATCTGc > 2:30546/1‑249 (MQ=255) tgtgAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATc > 1:179955/1‑251 (MQ=255) gtgAACTGG‑CAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACCGGGCGATCATGTGTTGATGCCCAACACCGCCTATGAACCGCGTCAGGATTTCTGTAGCCAAATACTCAGCAAACTGGGCGTAACGACATAATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGCGTTTCTGGAATCg > 2:114071/1‑251 (MQ=255) tCCCTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTc > 1:78361/1‑251 (MQ=255) cTGCGGGTCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg > 1:12208/1‑141 (MQ=255) | AGGACGCAGCAAAAAATACACTCTCGGCGCGGTAAATAGCGTGATTCAGCGCGCTTCTTCGCTGGTCTTTGACAGTGTAGAAGCCAAAAAACACGCGACACGTAATCGCGCCAATGGAGAGTTGTTCTATGGACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGG‑AAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAGGCTCCATCACCATGGAAGTCCACGACGTTC > NC_000913/3152274‑3152734 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |