Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,455,629 | T→C | D15G (GAC→GGC) | pyrB ← | aspartate carbamoyltransferase catalytic subunit |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,455,629 | 0 | T | C | 90.9% | 106.7 / 0.9 | 33 | D15G (GAC→GGC) | pyrB | aspartate carbamoyltransferase catalytic subunit |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (13/17); total (14/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.56e-01 |
CTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAA > REL606/4455559‑4455700 | ctctGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGttt < 2:438944/76‑1 (MQ=255) ctGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTAt < 2:890127/76‑1 (MQ=255) ctGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTAt > 2:5655/1‑76 (MQ=255) ctGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTAt > 1:364645/1‑76 (MQ=255) tGGTTGCTGTTTTTCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGTCTAAGGTCGTTTATg < 2:634521/76‑1 (MQ=255) ggTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATgg > 1:1019819/1‑76 (MQ=255) gTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGa > 2:330171/1‑76 (MQ=255) ttGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGaa > 1:324537/1‑76 (MQ=255) gggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATgg < 1:709777/76‑2 (MQ=255) ggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGat > 1:1044317/1‑76 (MQ=255) ggTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGat > 1:706727/1‑76 (MQ=255) ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATg < 2:436676/76‑1 (MQ=255) tGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATgn < 1:768041/76‑2 (MQ=255) cTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGttt > 2:85804/1‑76 (MQ=255) tttCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGatat < 2:843677/76‑1 (MQ=255) tCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAg < 1:719291/76‑1 (MQ=255) tCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAg < 1:1087649/76‑1 (MQ=255) cTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGAtt < 2:486870/76‑1 (MQ=255) cTGTCACCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGAtt < 2:1270606/76‑1 (MQ=255) tGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTa < 1:1184157/76‑1 (MQ=255) cGCCAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGcc > 1:907382/1‑76 (MQ=255) cAGCACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATc > 2:231552/1‑76 (MQ=255) cACCAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCttt < 2:370777/76‑1 (MQ=255) ccAGATTAAGGTCATCGCGACTAAGGCCGTTTATGGAAATCATATGTTTCTGATATAGCGGATTAGCCATCTTTTa < 2:1139980/76‑1 (MQ=255) aaGGTCATCGCGACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTg > 2:290039/1‑76 (MQ=255) gTCATCGCGACTAATGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACg < 2:395277/76‑1 (MQ=255) cgcgACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTgg < 2:387051/76‑1 (MQ=255) gACTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCa < 1:398706/76‑1 (MQ=255) cTAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaa < 1:707928/76‑1 (MQ=255) tAAGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaa > 2:382711/1‑76 (MQ=255) aaGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaa > 2:279563/1‑76 (MQ=255) aaGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaa < 2:213439/76‑1 (MQ=255) aaGGCCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCaaaaa > 1:1244828/1‑76 (MQ=255) | CTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAA > REL606/4455559‑4455700 |
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: ‑ |