Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,108,495 | T→G | E3A (GAA→GCA) | metJ ← | transcriptional repressor protein MetJ |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,108,495 | 0 | T | G | 96.8% | 256.9 / ‑1.2 | 63 | E3A (GAA→GCA) | metJ | transcriptional repressor protein MetJ |
Reads supporting (aligned to +/- strand): ref base T (2/0); new base G (35/26); total (37/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.07e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01 |
AACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGA > REL606/4108422‑4108565 | aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa > 2:997057/1‑76 (MQ=255) aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa > 2:469033/1‑76 (MQ=255) aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa > 2:469022/1‑76 (MQ=255) aaCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCa < 1:489574/76‑1 (MQ=255) aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAg > 1:1047480/1‑76 (MQ=255) aCCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAg > 1:1093357/1‑76 (MQ=255) cGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGcc < 1:107650/76‑1 (MQ=255) aaTCTTTTTGTCTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg < 1:603677/76‑1 (MQ=255) aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg > 1:26419/1‑76 (MQ=255) aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg > 1:1717577/1‑76 (MQ=255) aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg < 1:1648424/76‑1 (MQ=255) aaTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATg < 1:1624418/76‑1 (MQ=255) aTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATga < 1:1245658/76‑1 (MQ=255) aTCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATga < 1:63630/76‑1 (MQ=255) tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGAt < 1:514876/76‑1 (MQ=255) tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGAt < 1:869877/76‑1 (MQ=255) tttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGAt < 1:1606566/76‑1 (MQ=255) ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATa < 1:90512/76‑1 (MQ=255) ttttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATa < 1:550073/76‑1 (MQ=255) tttGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATAc > 1:293712/1‑76 (MQ=255) tGACTTTTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACtt < 2:3623/76‑1 (MQ=255) tGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACtt < 1:1621345/76‑1 (MQ=255) cTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAAt < 2:1434825/76‑1 (MQ=255) tGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATcc > 1:1572549/1‑76 (MQ=255) gTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCct > 2:1254314/1‑76 (MQ=255) ttCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCctc > 1:251769/1‑76 (MQ=255) ttCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCctc > 2:1177463/1‑76 (MQ=255) cACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCtt > 1:1363289/1‑76 (MQ=255) cACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCtt > 2:228305/1‑76 (MQ=255) cACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCtt > 1:1226759/1‑76 (MQ=255) cACTCTTCTTGCAGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCctctc > 2:460068/1‑75 (MQ=255) cttcttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTc > 2:857376/1‑76 (MQ=255) cttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAAt > 1:1163274/1‑76 (MQ=255) ttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATa > 2:826/1‑76 (MQ=255) ttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATa < 1:1411301/76‑1 (MQ=255) ttGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAAGCCTCTTCGTCAATa > 2:1657600/1‑76 (MQ=255) tGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATaa > 2:1639444/1‑76 (MQ=255) gCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATaaa > 2:1163529/1‑76 (MQ=255) gCCGTGCTCAGCGAATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTAAATCCTCCTCACCaaaaaa > 1:534312/1‑76 (MQ=255) ccGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTTAATAAAt < 1:1144723/76‑1 (MQ=255) gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg < 2:1411538/76‑1 (MQ=255) gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg > 1:1738211/1‑76 (MQ=255) gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTg > 2:1652530/1‑76 (MQ=255) gTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTNCGTCAATAAATTg < 1:1190900/76‑1 (MQ=255) gCTCAGCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTgag > 2:704178/1‑76 (MQ=255) cAGCGTATGGGCTGATATATTCGCCGCTACATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagacc > 2:839066/1‑76 (MQ=255) gCGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccag < 1:835373/76‑1 (MQ=255) cGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga > 2:1496725/1‑76 (MQ=255) cGTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccaga > 2:456702/1‑76 (MQ=255) gTATGGGCTGATATATTCGCCGCTCCGTGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac > 1:478177/1‑76 (MQ=255) gTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac > 2:784258/1‑76 (MQ=255) gTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagcc > 2:314863/1‑74 (MQ=255) gTATGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGagaccagac < 2:594227/76‑1 (MQ=255) tGGGCTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATCAATTGAGACCAGACcac > 1:787401/1‑76 (MQ=255) cTGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAg < 1:1279691/76‑1 (MQ=255) tGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGt > 1:530254/1‑76 (MQ=255) tGATATATTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGt > 1:1591030/1‑76 (MQ=255) gATATATACGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGAACACACAgat > 2:775468/1‑74 (MQ=255) atatatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGCCCACACAGTTg < 2:28739/76‑1 (MQ=255) atatatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTg < 1:1242379/76‑1 (MQ=255) atatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAt > 1:1674321/1‑76 (MQ=255) tatTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtg > 1:37674/1‑76 (MQ=255) atTCGCCGCTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGAtgt < 1:1163529/76‑1 (MQ=255) cgccgcTCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTggg < 2:55890/76‑1 (MQ=255) gcTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACt > 2:24903/1‑76 (MQ=255) tCCATGCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGa > 1:31578/1‑76 (MQ=255) | AACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATATATTCGCCGCTCCATTCAGCCATGAGATACTTAATCCTCTTCGTCAATAAATTGAGACCAGACCACACAGTTGATGTGGGTACTGA > REL606/4108422‑4108565 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |