Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | AE015451 | 5,087,105 | A→T | 12.5% | L73Q (CTG→CAG) | astB ← | N‑succinylarginine dihydrolase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | AE015451 | 5,087,105 | 0 | A | T | 12.5% | 115.8 / 9.3 | 40 | L73Q (CTG→CAG) | astB | N‑succinylarginine dihydrolase |
Reads supporting (aligned to +/- strand): ref base A (10/25); new base T (2/3); total (12/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01 |
GCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTCC > AE015451/5086873‑5087334 | gcCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAgcg > 1:13501/1‑236 (MQ=255) gCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGGCCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCTGcgcc < 2:147892/226‑1 (MQ=255) cccGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAaca < 2:629490/234‑1 (MQ=255) atatTTGCAGTTGAGGTTGGCGGCGGGGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGcgcc < 1:2345/196‑1 (MQ=255) aGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGt > 1:300789/1‑235 (MQ=255) gTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCg < 1:639255/235‑1 (MQ=255) gTTGGCGGCGGGGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCgg < 1:371774/236‑1 (MQ=255) cACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTc < 1:289177/235‑1 (MQ=255) cACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTc < 1:720891/235‑1 (MQ=255) gCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCct < 2:327831/215‑1 (MQ=255) aaCTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCCCGCGGGGTGGACGCctg > 1:232073/1‑202 (MQ=255) aCTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGACCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCa > 2:194060/1‑116 (MQ=255) aCTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCg > 2:358415/1‑218 (MQ=255) gCCACCAGCACCGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACa < 1:228930/235‑1 (MQ=255) caGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGTACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGgcc < 2:388106/234‑1 (MQ=255) caGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGgccg < 1:358416/235‑1 (MQ=255) aaCGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACa < 2:501547/226‑1 (MQ=255) gggcTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGACCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCa < 1:91564/200‑1 (MQ=255) cgcgcTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCTGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCa > 2:698791/1‑189 (MQ=255) cgcgcTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCTGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCa < 1:698792/189‑1 (MQ=255) cgcgcTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACa < 2:263503/204‑1 (MQ=255) cgcgcTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAgg > 1:545688/1‑236 (MQ=255) gcgcTGGATAACTTCGGCATCGCTGCCCCTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGTCCCGCGCTTTCATCTTCGCCACGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTg > 2:40361/1‑180 (MQ=255) gcgcTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCTGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTgg < 1:425020/158‑1 (MQ=255) gcgcTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCTGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTgg > 2:425019/1‑158 (MQ=255) cgcgGGATAACTTCGGCAGCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCct < 2:208329/163‑1 (MQ=255) aTAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCaaa < 1:575572/236‑1 (MQ=255) cctCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCg < 2:660927/169‑1 (MQ=255) ttCGGCATCGCTGCCACTGACCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAg < 1:253262/188‑1 (MQ=255) cATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCAt > 2:160357/1‑195 (MQ=255) aTCGCTGCCACTGAGCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCccacca < 1:353182/215‑1 (MQ=255) tGAACCCGAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACt < 2:647042/221‑1 (MQ=255) aaCCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGCGGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGt < 2:430684/191‑1 (MQ=255) cccAGCCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTGGAAGCCCATCTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCa < 2:167275/154‑1 (MQ=255) gcggcgCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCt < 1:160357/235‑1 (MQ=255) ggcgCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTcc > 1:681499/1‑235 (MQ=255) cgcAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCgg < 2:113582/100‑1 (MQ=255) cgcAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTAAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACa < 1:172619/160‑1 (MQ=255) gcAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGgccg < 2:477729/164‑1 (MQ=255) gcAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGgccg > 1:477730/1‑164 (MQ=255) aGCGCCTCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCGGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTtgt < 1:17393/172‑1 (MQ=255) | GCCCAGCACGCGGCTGGTGGTCGGGTGCTCGATGCTGCGGTGATATTTGCAGTTGAGGTTGGCGGCGGTGAAGTGCACGCGACCGTCAGCAGTGTCGGCACTCGGGCTGACTGTTGCGGCGTTGGCCACCCACATGCTCGAGGCCGAGCAACTGGCCACCAGCAACGGCATGGCCTCGCGGGCGGCGCGCTGGATAACTTCGGCATCGCTGCCACTGAACCCCAGGCGGCGCAGCGCCGCCACGTCCGGGCGCTCTTGTGGCGCCAGAACACCCTGCTTGAAGCCCATGTCGGCCAGCGCTTTCATCTTCGCCAGGCCCTGGCGCGCGGCCTCACGCGGGTTGGACGCCTGCTGGCTATTGCTCTGCGAAGCCACGTTACCGTAAGACAGGCCGCCATAGTTGTGGGTAGGCCCCACCAGGCCATCAAAATTCACTTCATAGGATTTCATCGGCTAGGCTCC > AE015451/5086873‑5087334 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |