| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 378,843 | T→G | K243N (AAA→AAC) | frmA ← | S‑(hydroxymethyl)glutathione dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 378,843 | 0 | T | G | 95.0% | 65.1 / NA | 20 | K243N (AAA→AAC) | frmA | S‑(hydroxymethyl)glutathione dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base T (0/0); major base G (9/10); minor base C (1/0); total (10/10) | |||||||||||
| 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 | |||||||||||
CCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGATAATCCGACCCGCTTTCGCCTGACGCGCGCCCTGAACCACTG > NC_000913/378719‑378970 | ccGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATGTTCTATCGGCTTGTCGTa > 1:199876/1‑133 (MQ=255) cGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAg < 1:20452/133‑1 (MQ=255) gACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCAtt < 1:9421/133‑1 (MQ=255) gACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCAtt > 1:74809/1‑133 (MQ=255) gACCCCAGCCCCGGTGCGCACTTTCCAGGGCCGCACGCATCTCGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCAtt < 2:266754/133‑1 (MQ=255) aCCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTc > 1:79925/1‑133 (MQ=255) ccccAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCGAACAGGTCATCTTTCATCAGGTTGTCGTCGTCATTCg > 1:268208/1‑133 (MQ=255) cgGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGc > 2:564732/1‑133 (MQ=255) ggTGCGCACTTTCCATCGTCGCCCGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCa < 1:222626/133‑1 (MQ=255) cacgcaTCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGCGTTAATGCAGTGGGAAGCACCGAAGCGAcg > 1:258891/1‑133 (MQ=255) aCGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGAt < 2:74809/133‑1 (MQ=255) gATGCATTTAAAGGTATGGTCGGTACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCggg < 1:286710/133‑1 (MQ=255) ttCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTa < 1:448765/133‑1 (MQ=255) aGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGAt < 1:407243/133‑1 (MQ=255) tGGTCGATCCCCCATTTGCTGATATCCAACCGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGa < 1:183581/133‑1 (MQ=255) tGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGa > 2:76798/1‑133 (MQ=255) cAGGACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGATAATCCGACCCGCTTTCGCCTGAcgcgcg < 2:79925/133‑1 (MQ=255) gACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGATAATCCGACCCGCTTTCGCCTGACGCGCGccc > 1:609791/1‑133 (MQ=255) gACATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGATAATCCGACCCGCTTTCGCCTGACGCGCGccc > 1:56352/1‑133 (MQ=255) aCATCTTTTATCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGATAATCCGACCCGCTTTCGCCTGACGCGCGCCCt > 1:607934/1‑133 (MQ=255) aTCGGGTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGATAATCCGACCCTCTTTCGCCTGACGCGCGCCCTGAACCACTg > 1:206962/1‑133 (MQ=255) | CCGACTGACCCCAGCCGCGGTGCGCACTTTCCAGCGCCGCACGCATCACGTTGACGTTACCGATGCATTCAAAGGTATGGTCGATACCCCATTTGTTGATATCCAACAGGACATCTTTTATCGGTTTGTCGTAGTCATTCGGGTTAATGCAGTCGGTAGCACCGAAGCGACGCGCCAGATCGAATTTCTTCGGGTTGGTATCGATAGCGATAATCCGACCCGCTTTCGCCTGACGCGCGCCCTGAACCACTG > NC_000913/378719‑378970 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |