Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,873,040 | C→A | R332L (CGC→CTC) | glmS ← | D‑fructose‑6‑phosphate amidotransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,873,040 | 0 | C | A | 92.9% | 140.1 / 3.6 | 42 | R332L (CGC→CTC) | glmS | D‑fructose‑6‑phosphate amidotransferase |
Reads supporting (aligned to +/- strand): ref base C (0/3); new base A (18/21); total (18/24) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01 |
GGTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAACCAT > REL606/3872967‑3873110 | ggTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAgg > 2:501030/1‑76 (MQ=255) ggTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAgg < 2:783566/76‑1 (MQ=255) cgcgGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCa < 1:1106994/76‑1 (MQ=255) gcgGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCag < 1:583672/76‑1 (MQ=255) ttCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCg < 1:226841/76‑1 (MQ=255) tgtgACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGt > 1:1073547/1‑76 (MQ=255) cAAGGTGATTATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCAc > 2:253431/1‑76 (MQ=255) cAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCAc < 1:652926/76‑1 (MQ=255) cAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCAc > 1:265182/1‑76 (MQ=255) atcaGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTg < 2:13661/76‑1 (MQ=255) tcaGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGc < 2:595559/76‑1 (MQ=255) caGGCTGTTACGACGCACGGCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCt < 1:240464/76‑1 (MQ=255) caGGCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCt > 1:1245650/1‑76 (MQ=255) ggCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAg < 2:472560/76‑1 (MQ=255) gCTGTTACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGc > 1:665337/1‑76 (MQ=255) tGTTACGACGCACNGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGa > 1:853225/1‑76 (MQ=255) tACGACGCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTc > 2:788376/1‑76 (MQ=255) acgacgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCa < 1:323297/76‑1 (MQ=255) acgacgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCa < 1:1149873/76‑1 (MQ=255) cgacgCACGGCAGGTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGCAGTACTGGCTAGCGATTCaa > 2:1060945/1‑76 (MQ=255) gacgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCaaa > 1:1155967/1‑76 (MQ=255) cgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAAcc > 2:57363/1‑76 (MQ=255) cgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAAcc > 1:1258503/1‑76 (MQ=255) cgCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCAACGTCGCACGGAATACCTGCTAGCGATTCAAAcc < 1:99790/76‑1 (MQ=255) gCACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCa < 2:842698/76‑1 (MQ=255) cACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAg < 1:629517/76‑1 (MQ=255) cACGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAg > 2:481481/1‑76 (MQ=255) aCGGCAGCTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGt > 1:810413/1‑76 (MQ=255) cGGCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTa < 2:446387/76‑1 (MQ=255) ggCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAg < 2:364836/76‑1 (MQ=255) ggCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGNTTCAAACCAGTAg < 1:857314/76‑1 (MQ=255) ggCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAg < 2:788882/76‑1 (MQ=255) ggCAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAg > 1:1053908/1‑76 (MQ=255) cAGATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCg < 2:209751/76‑1 (MQ=255) gATTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCggg < 2:865714/76‑1 (MQ=255) aTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGa > 1:737579/1‑76 (MQ=255) aTTTGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGa > 1:565847/1‑76 (MQ=255) tttGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGan < 1:337507/76‑2 (MQ=255) ttGCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGaaa < 1:107204/76‑1 (MQ=255) tGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAAc < 2:565025/76‑1 (MQ=255) gCGATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAAcc > 2:145696/1‑76 (MQ=255) gATAGAGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAACCAt < 2:1258503/76‑1 (MQ=255) | GGTATCCGCGGTTTCGCCAGACTGTGACAAGGTGATCATCAGGCTGTTACGACGCACGGCAGATTTGCGATAGCGGAATTCAGAGGCGATTTCGACGTCGCACGGAATACCTGCTAGCGATTCAAACCAGTAGCGGGAAACCAT > REL606/3872967‑3873110 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |