Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 4,108,409 | T→C | K32E (AAG→GAG) | metJ ← | transcriptional repressor protein MetJ |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 4,108,409 | 0 | T | C | 95.7% | 274.8 / 0.6 | 69 | K32E (AAG→GAG) | metJ | transcriptional repressor protein MetJ |
Reads supporting (aligned to +/- strand): ref base T (0/3); new base C (31/35); total (31/38) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
TCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATA > REL606/4108336‑4108479 | tCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCaa < 1:1695581/76‑1 (MQ=255) tCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCaa < 1:1390642/76‑1 (MQ=255) cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAag > 2:1458643/1‑76 (MQ=255) cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAag < 2:1819120/76‑1 (MQ=255) cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAag > 2:1740359/1‑76 (MQ=255) cGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTT‑ATCGGTGAGGATTTTTAACACCTCAaga < 1:724212/76‑1 (MQ=255) tGTTGGTATCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAgg > 1:495482/1‑76 (MQ=255) ngTTGGTATCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAgg > 1:110525/2‑76 (MQ=255) gTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGa < 1:1067199/76‑1 (MQ=255) ttGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGaa > 2:1572340/1‑76 (MQ=255) tGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAAt < 1:505196/76‑1 (MQ=255) gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATgg < 1:842911/76‑1 (MQ=255) gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATgg < 1:1212055/76‑1 (MQ=255) gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATgg < 2:213376/76‑1 (MQ=255) gTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATgg < 1:1484735/76‑1 (MQ=255) tAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGa > 1:1019985/1‑76 (MQ=255) aGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGaa > 1:1215601/1‑76 (MQ=255) aGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGaa > 2:308063/1‑76 (MQ=255) cGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAAc > 2:1455336/1‑76 (MQ=255) gTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAAcc < 2:77844/76‑1 (MQ=255) aCGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTa < 2:1095718/76‑1 (MQ=255) cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTaa < 1:659084/76‑1 (MQ=255) cGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTaa < 1:18462/76‑1 (MQ=255) ncAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAAt > 1:896390/2‑76 (MQ=255) aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCt < 2:359982/76‑1 (MQ=255) aGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCt > 2:323533/1‑76 (MQ=255) gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt > 1:1049242/1‑76 (MQ=255) gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt > 1:544903/1‑76 (MQ=255) gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt < 2:747678/76‑1 (MQ=255) gttgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttt > 1:799798/1‑76 (MQ=255) tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttttt < 1:1134293/76‑1 (MQ=255) tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttttt < 1:898696/76‑1 (MQ=255) tgttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCttttt > 1:1407709/1‑76 (MQ=255) gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTg < 2:687399/76‑1 (MQ=255) gttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTg > 1:291112/1‑76 (MQ=255) ttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGa < 1:1454262/76‑1 (MQ=255) ttCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGa < 1:516797/76‑1 (MQ=255) ccTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTg < 2:1407709/76‑1 (MQ=255) tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGtt > 2:564790/1‑76 (MQ=255) tGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGtt > 1:1245366/1‑76 (MQ=255) gacgacGCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTc > 2:578843/1‑76 (MQ=255) cgacgCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCAc < 1:1521117/76‑1 (MQ=255) acgCGTGCGTTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCActc > 1:842197/1‑76 (MQ=255) gcgtgcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtctt < 2:486953/76‑1 (MQ=255) gtgcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct < 2:873698/76‑1 (MQ=255) gtgcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACtcttct > 1:55595/1‑76 (MQ=255) gcgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTg < 2:1115910/76‑1 (MQ=255) cgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc > 2:350500/1‑76 (MQ=255) cgtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGc < 1:242853/76‑1 (MQ=255) gtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc < 2:1111795/76‑1 (MQ=255) gtTCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGcc > 1:1471748/1‑76 (MQ=255) ttCATCGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCg < 1:485690/76‑1 (MQ=255) cGGTGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTc < 1:1045355/76‑1 (MQ=255) tGAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGc < 2:48807/76‑1 (MQ=255) gAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCg > 1:1775946/1‑76 (MQ=255) gAGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCg > 1:1775940/1‑76 (MQ=255) aGGATTTTTAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGt < 2:1562840/76‑1 (MQ=255) tttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATgg < 1:6670/76‑1 (MQ=255) tttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATgg < 1:224196/76‑1 (MQ=255) ttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg > 1:1540483/1‑76 (MQ=255) ttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg > 1:883999/1‑76 (MQ=255) ttttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATggg > 1:1434591/1‑76 (MQ=255) tttAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGc > 1:824396/1‑76 (MQ=255) tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg < 2:470171/76‑1 (MQ=255) tAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg < 2:671992/76‑1 (MQ=255) tAACACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTg > 2:488127/1‑76 (MQ=255) aaCACCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGa > 2:1881324/1‑76 (MQ=255) acacCTCAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGat < 1:1261106/76‑1 (MQ=255) cacCTCAAGAGGAATGGAAACCGTAATCGTTGTGACTTGTTCACTCTGCTTGCCGTGCTCAGCGTATGGGCTGata < 2:329375/76‑1 (MQ=255) | TCGCTGTTGGTAGCGTGACGCAGGTTGTTCACCTGACGACGCGTGCGTTCATCGGTGAGGATTTTTAACACCTTAAGAGGAATGGAAACCGTAATCTTTTTGACTTGTTCACTCTTCTTGCCGTGCTCAGCGTATGGGCTGATA > REL606/4108336‑4108479 |
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: ‑ |