Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 1,101,027 | G→T | intergenic (+95/‑250) | ycdB → / → phoH | hypothetical protein/conserved protein with nucleoside triphosphate hydrolase domain |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 1,101,027 | 0 | G | T | 98.5% | 254.1 / ‑5.6 | 68 | intergenic (+95/‑250) | ycdB/phoH | hypothetical protein/conserved protein with nucleoside triphosphate hydrolase domain |
Reads supporting (aligned to +/- strand): ref base G (1/0); major base T (33/33); minor base C (1/0); total (35/33) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01 |
TTTGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACAC > REL606/1100954‑1101095 | tttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTAt > 2:111549/1‑76 (MQ=255) tttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTAt < 1:1167056/76‑1 (MQ=255) ttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGtata < 1:925718/76‑1 (MQ=255) ttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGtata > 1:1582629/1‑76 (MQ=255) ttGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGtata < 2:1487842/76‑1 (MQ=255) gtgCGAAGATGCCTAATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAg < 2:324243/76‑1 (MQ=25) tgCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAgg > 2:841247/1‑76 (MQ=255) tgCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAgg < 1:267874/76‑1 (MQ=255) tgCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGAGTGGGTTTGGTAGGCTGTATAAgg > 2:841243/1‑76 (MQ=21) gcgcAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGc < 2:1230991/76‑1 (MQ=255) cAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGtt < 1:1242773/76‑1 (MQ=255) cAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGtt < 2:510013/76‑1 (MQ=255) cAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGtt < 2:633944/76‑1 (MQ=255) aGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTc > 2:1544932/1‑76 (MQ=255) aGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTc > 2:663421/1‑76 (MQ=255) gATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCa < 1:411851/76‑1 (MQ=255) gATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCa < 2:1367165/76‑1 (MQ=255) aTGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAc < 1:29269/76‑1 (MQ=255) tGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcg > 2:552964/1‑76 (MQ=255) gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc < 1:441312/76‑1 (MQ=255) gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc < 1:485268/76‑1 (MQ=255) gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc < 1:1503357/76‑1 (MQ=255) gCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgc < 1:1409588/76‑1 (MQ=255) ccTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgcc > 2:1446213/1‑76 (MQ=255) cTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccg > 2:943010/1‑76 (MQ=255) cTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccg > 1:391392/1‑76 (MQ=255) cTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCAcgccg > 1:49727/1‑76 (MQ=255) tGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccgc > 2:57718/1‑76 (MQ=255) tGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCAcgccgc > 2:1628073/1‑76 (MQ=255) gATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCa > 2:1551395/1‑76 (MQ=255) gCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATcc < 2:1394964/76‑1 (MQ=255) cGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCg > 2:13186/1‑76 (MQ=255) gCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGa < 1:1543430/76‑1 (MQ=255) cTTGCGCGTCTTATCATGCCTACAATCGGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGAt > 1:705460/1‑76 (MQ=25) gcgTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGc < 2:1333391/76‑1 (MQ=255) cgTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCa > 1:804929/1‑76 (MQ=255) cgTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGCATAAGGCGTTCACGCCGCATCCGGCGATCGTGcc > 1:817302/1‑75 (MQ=255) tCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCAct > 1:866996/1‑76 (MQ=255) tCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCAct > 2:891367/1‑76 (MQ=255) cTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActg > 2:1614263/1‑76 (MQ=255) cTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActg > 2:1287886/1‑76 (MQ=255) cTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActg > 2:1619019/1‑76 (MQ=255) ttATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActga > 1:931031/1‑76 (MQ=255) tATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgat > 2:1131657/1‑76 (MQ=255) aTCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgatg > 2:446449/1‑76 (MQ=255) ccTACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgatgcctgat < 2:605604/76‑1 (MQ=255) tACAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCActgatgcctgatgc < 2:144274/76‑1 (MQ=255) cAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaa > 1:411584/1‑76 (MQ=255) cAATCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaa > 2:431051/1‑76 (MQ=255) aaTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaaa < 2:1106671/76‑1 (MQ=255) aaTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCaaa < 2:581703/76‑1 (MQ=255) aTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAAt < 2:6267/76‑1 (MQ=255) aTCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAAt < 2:1484499/76‑1 (MQ=255) tCAGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATc < 1:634225/76‑1 (MQ=255) aGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCct > 2:128064/1‑76 (MQ=255) aGTGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCct > 1:1436209/1‑76 (MQ=255) tGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgc < 1:1612117/76‑1 (MQ=255) tGCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgc < 2:1153682/76‑1 (MQ=255) gCGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgct < 2:1567498/76‑1 (MQ=255) cGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgctg < 1:311236/76‑1 (MQ=255) cGGGTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCctgctg < 1:907362/76‑1 (MQ=255) ggTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGaa < 2:828216/76‑1 (MQ=255) gTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGaaa > 2:544718/1‑76 (MQ=255) gTTTGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGaaa > 2:1582348/1‑76 (MQ=255) tttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAg > 2:1169529/1‑76 (MQ=255) ttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGc < 1:327158/76‑1 (MQ=255) ttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGc < 1:307483/76‑1 (MQ=255) ttGGTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGc < 1:278423/76‑1 (MQ=255) ggTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcac > 2:852885/1‑76 (MQ=255) gTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcaca > 2:311856/1‑76 (MQ=255) gTAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcaca > 2:1279296/1‑76 (MQ=255) tAGGCTGTATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGcacac > 1:489692/1‑76 (MQ=255) | TTTGTGCGCAGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAATCAGTGCGGGTTTGGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCGATCGTGCACTGATGCCTGATGCAAATCCTGCTGAAAGCACAC > REL606/1100954‑1101095 |
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: ‑ |