| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 1,829,328 | C→G | 15.6% | G136A (GGC→GCC) | astD ← | succinylglutamic semialdehyde dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,829,328 | 0 | C | G | 15.6% | 87.5 / 11.2 | 32 | G136A (GGC→GCC) | astD | succinylglutamic semialdehyde dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base C (19/8); new base G (3/2); total (22/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.37e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01 | |||||||||||
ATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACC > NC_000913/1829185‑1829473 | aTAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGc > 4:394897/1‑149 (MQ=255) tAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCg > 3:260833/1‑149 (MQ=255) aaCGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCggg > 3:211479/1‑149 (MQ=255) tACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCGgtgtt < 4:334252/149‑2 (MQ=255) tCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCt < 3:394897/149‑1 (MQ=255) gCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcg > 3:200707/1‑149 (MQ=255) cACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgccgcgcc > 3:223479/1‑149 (MQ=255) gTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGc > 1:32769/1‑108 (MQ=255) gTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGc < 2:32769/108‑1 (MQ=255) gggTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTAc < 4:195901/149‑1 (MQ=255) tttAAAGGTAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACgct > 2:261167/1‑149 (MQ=255) aGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCg > 2:168510/1‑149 (MQ=255) aTAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCgg > 4:79120/1‑101 (MQ=255) aTAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCCGCACGGCGTGCTGGCgg < 3:79120/102‑2 (MQ=255) tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgc < 3:17237/112‑1 (MQ=255) tGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTcgc > 4:17237/1‑112 (MQ=255) gtgtTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACg > 2:374619/1‑149 (MQ=255) aCCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGt < 1:374619/149‑1 (MQ=255) gCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGcc < 3:43772/149‑1 (MQ=255) ggAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGa > 4:251419/1‑149 (MQ=255) aaCGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGata > 4:170340/1‑149 (MQ=255) atGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGa < 3:170340/149‑1 (MQ=255) tGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAt < 4:223479/149‑1 (MQ=255) tCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGAttt > 2:163109/1‑149 (MQ=255) ggCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGat > 2:80707/1‑149 (MQ=255) gACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGcc > 4:49800/1‑149 (MQ=255) gggAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCAc > 2:173966/1‑149 (MQ=255) ttATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGt > 1:402350/1‑149 (MQ=255) taCGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGc > 1:435124/1‑149 (MQ=255) gCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCt > 1:244471/1‑149 (MQ=255) cacCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgc < 1:163109/149‑1 (MQ=255) atcGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgcg > 4:195692/3‑149 (MQ=255) atcGCCCGCACGGCGTGCTGGCGGTGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAgcg > 1:282022/3‑149 (MQ=255) aCGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTAcc > 2:148725/1‑149 (MQ=255) | ATAAACGCATTACCGCTTCGCCACTCCACGGTGTCAGTTCGCTGGGTTTAAAGATAATGGTGTTACCTGCCAGCAATGCCGGAACGATATGTCCGTTCGGCAAATGACCAGGGAAATTATACGGCCCAAACACCGCCAGCACGCCGTGCGGGCGATGTCGCAGGCTCGCCGCGCCGTCCGGCATTTCACTACGCTGCTCGCCGGTACGAACGTGATACGCCTTAATTGATATCGCGATTTTATTGATCATCGCCGTCACTTCGGTTGCCGCTTCCCAGCGCGGCTTACC > NC_000913/1829185‑1829473 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |