| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | REL606 | 3,914,929 | C→T | S215F (TCC→TTC) | ilvE → | branched‑chain amino acid aminotransferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | REL606 | 3,914,929 | 0 | C | T | 97.8% | 155.4 / ‑3.2 | 45 | S215F (TCC→TTC) | ilvE | branched‑chain amino acid aminotransferase |
| Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (25/19); total (26/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01 | |||||||||||
TATATCTCTGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTCCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGATGTGCGTG > REL606/3914856‑3915000 | tataTCTCTGAAGGCGCAGGCGAAAACGTGCTTGNAGTGCAANTTNGTNTGCTCTTCACCCCACCGTTCAccttct > 1:1590207/1‑76 (MQ=255) atatCTCTGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACcttctc > 1:407821/1‑76 (MQ=255) tatCTCTGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCttctcc < 1:1102583/76‑1 (MQ=255) atCTCTGAAGGCGCAGGCGAGAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCcg > 1:94327/1‑76 (MQ=255) atCTCTGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCcg > 1:1307598/1‑76 (MQ=255) atCTCTGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCcg > 1:595154/1‑76 (MQ=255) ctctGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCcgcg < 2:1671318/76‑1 (MQ=255) ctctGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGATCACCTTCTCcgcg > 1:935125/1‑76 (MQ=255) ctctGAAGGCGCAGGCGAAAACCTGTTNGAAGTGANNGATGGAGTGCTGTTCACCCCGCCGTTCATCTTCTCcgcg > 1:1806796/1‑76 (MQ=255) ctGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCt > 1:1311523/1‑76 (MQ=255) ctGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGAGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCt > 1:804307/1‑76 (MQ=255) tggaGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATAGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTg < 2:250158/73‑1 (MQ=255) gAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGc > 2:1269565/1‑76 (MQ=255) gcAGGCGNCAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTa > 2:1470203/1‑76 (MQ=255) cAGGCGAAAACCTGTTTGAGGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTAt < 2:25042/76‑1 (MQ=255) cAGGCGAAAACCTGTTTGAAGTGAAAGNGGNTCNGCCGTTCNCCCCACAGTTCACCTTCTCCGCGCTGCCGGGTAt < 1:1587227/76‑1 (MQ=255) aGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTAtt > 2:180511/1‑76 (MQ=255) aGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTAtt > 2:1520437/1‑76 (MQ=255) ggCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTa < 2:1268465/76‑1 (MQ=255) gAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTAccc < 2:588090/76‑1 (MQ=255) aaaaCCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCg > 2:1313647/1‑76 (MQ=255) aaCCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTg < 1:865592/76‑1 (MQ=255) tttGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCca > 1:387376/1‑76 (MQ=255) tgaaaGTGAAAGATGGTGTGCTGTCCCCCCCCCCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCcat < 2:1175841/73‑1 (MQ=255) tGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCcatc > 2:169064/1‑76 (MQ=255) tGAAGGGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCcatc < 1:1399717/76‑1 (MQ=255) gTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCAtcatca < 1:248058/76‑1 (MQ=255) gAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCaaa < 1:1103481/76‑1 (MQ=255) gAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCaaa > 2:1588330/1‑76 (MQ=255) gAAAGATGGTGTGCTGTTAACCACACCNTTNACNTTTTCCGNGCTGCCGGGTATTACCCGTGATGCCATCATCaaa < 1:1593371/76‑1 (MQ=255) aaaGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAAc > 2:701303/1‑76 (MQ=255) aaGATGGTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACt > 2:1145010/1‑76 (MQ=255) ggTGTGCTGTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGa > 1:182938/1‑76 (MQ=255) gtgCTGTTCCCCCCACCGGTCCGCTTCNCCNCGNTGGCGGGNATGACCAGTGATGCCATCATCAAACTGGCGAAag < 1:1578324/76‑1 (MQ=255) gtgCTGTTCACCCCACCGTGCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAag > 1:1556607/1‑76 (MQ=255) tgCTGTTCACCCCACCGTTCACCTCCTCCGCGCTACCGGGTNTTACCCGTGATGCCATCANCANACTGGCGAAaga > 1:1490589/1‑76 (MQ=255) gTTCACCCCCCCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTg < 2:1350199/76‑1 (MQ=255) gTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTg > 2:1569113/1‑76 (MQ=255) gTTCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTg > 2:1126306/1‑76 (MQ=255) tCACCCCACCGTTTACCTTCNCCGCGCTGCAGGGTANTACCCGTGATGCCATTATCAATCTGGCGAAAGAGCTggg < 2:1825618/76‑1 (MQ=255) tCACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAATNAGCTggg < 2:1476052/76‑1 (MQ=255) cACCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGa > 1:597359/1‑76 (MQ=255) aCCCCACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGaa < 1:444753/76‑1 (MQ=255) ccccACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAAt > 2:749791/1‑76 (MQ=255) ccccACCGTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAAt < 2:1695513/76‑1 (MQ=255) gTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGAtgtg < 2:1421895/76‑1 (MQ=255) gTTCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGAtgtg > 1:1126916/1‑76 (MQ=255) tCACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTCGGAATTGATGTGCg > 2:337041/1‑76 (MQ=255) cACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGATGTGCGt > 2:1360615/1‑76 (MQ=255) cACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGATGTGCGt > 2:1360611/1‑76 (MQ=255) cACCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGATGTGCGt > 2:1018690/1‑76 (MQ=255) aCCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGATGTGCGTg < 1:536111/76‑1 (MQ=255) aCCTTCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGATGTGCGTg < 1:1365348/76‑1 (MQ=255) | TATATCTCTGAAGGCGCAGGCGAAAACCTGTTTGAAGTGAAAGATGGTGTGCTGTTCACCCCACCGTTCACCTCCTCCGCGCTGCCGGGTATTACCCGTGATGCCATCATCAAACTGGCGAAAGAGCTGGGAATTGATGTGCGTG > REL606/3914856‑3915000 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |