Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 977,118 | T→A | D405E (GAT→GAA) | aroA → | 3‑phosphoshikimate 1‑carboxyvinyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 977,118 | 0 | T | A | 97.7% | 167.5 / ‑3.4 | 44 | D405E (GAT→GAA) | aroA | 3‑phosphoshikimate 1‑carboxyvinyltransferase |
Reads supporting (aligned to +/- strand): ref base T (0/1); new base A (24/19); total (24/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.55e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01 |
ATACAATGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGATCCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGCCTGAA > REL606/977046‑977188 | aTACAATGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAAccc > 1:1446178/1‑76 (MQ=255) tACAATGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCa > 1:1067545/1‑76 (MQ=255) tACAATGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCa > 1:1067557/1‑76 (MQ=255) cAATGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCaaa > 2:868017/1‑76 (MQ=255) aTGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATg > 1:1008990/1‑76 (MQ=255) tGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGc > 2:237766/1‑76 (MQ=255) gATCACCGGATGTCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCa < 1:790970/76‑1 (MQ=255) aTCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCAc > 2:1139563/1‑76 (MQ=255) cACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACgg < 2:173722/76‑1 (MQ=255) cGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCa < 1:1705099/76‑1 (MQ=255) aTGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCaaaa < 2:86344/76‑1 (MQ=255) tGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAAc < 1:850009/76‑1 (MQ=255) ggCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACa < 1:639252/76‑1 (MQ=255) ggCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACa < 1:908386/76‑1 (MQ=255) ggCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACa < 1:1019528/76‑1 (MQ=255) gCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACAt < 1:1160337/76‑1 (MQ=255) cGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACAtt > 1:511113/1‑76 (MQ=255) cGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTAttt < 1:868017/76‑1 (MQ=255) cTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCAGGGCCAAAACATTTCCGGATTATTTCg < 1:1371816/76‑1 (MQ=255) cTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCg > 2:893373/1‑76 (MQ=255) tggtggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCCCGGCCAAAACATTTCCGGATTATTTCGa > 2:1512460/1‑76 (MQ=255) tggtggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGatta > 1:481025/1‑70 (MQ=255) tggtggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGatta < 2:481025/70‑1 (MQ=255) tggtggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGa > 1:1557337/1‑76 (MQ=255) ggtggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGag < 2:1008990/76‑1 (MQ=255) tggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGagca > 1:604914/1‑76 (MQ=255) tggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGagca > 1:604908/1‑76 (MQ=255) ggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGagcag < 1:1796719/76‑1 (MQ=255) ggCGTTGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGagcag > 1:1697372/1‑76 (MQ=255) gTTGTCAGATGCGCCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTg < 1:1696288/76‑1 (MQ=255) tGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGgc > 1:707629/1‑76 (MQ=255) tGTCAGATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGgc > 1:707648/1‑76 (MQ=255) gATACACCAGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGa > 2:743669/1‑76 (MQ=255) acCAGTGACGATTCTTGATCCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGc < 1:858191/76‑1 (MQ=255) aGTGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAg > 1:1168736/1‑76 (MQ=255) tGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGc < 2:1335709/76‑1 (MQ=255) tGACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGc < 2:1238377/76‑1 (MQ=255) gACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCa > 2:559283/1‑76 (MQ=255) gACGATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCa > 2:843535/1‑76 (MQ=255) gATTCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGcc > 2:1880702/1‑76 (MQ=255) ttCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGCCtg < 2:231570/76‑1 (MQ=255) ttCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGCCtg > 1:58759/1‑76 (MQ=255) tCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGCCtga < 1:1891215/76‑1 (MQ=255) tCTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGCCtga < 1:61363/76‑1 (MQ=255) cTTGAACCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGCCtgaa > 1:1283979/1‑76 (MQ=255) | ATACAATGATCACCGGATGGCGATGTGTTTCTCGCTGGTGGCGTTGTCAGATACACCAGTGACGATTCTTGATCCCAAATGCACGGCCAAAACATTTCCGGATTATTTCGAGCAGCTGGCGCGGATTAGCCAGGCAGCCTGAA > REL606/977046‑977188 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |