Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | REL606 | 3,535,089 | (A)6→5 | coding (9/597 nt) | yhhF → | predicted methyltransferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | REL606 | 3,535,084 | 0 | A | . | 94.9% | 172.1 / 3.8 | 39 | coding (4/597 nt) | yhhF | predicted methyltransferase |
Reads supporting (aligned to +/- strand): ref base A (1/1); new base . (13/24); total (14/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01 |
AGTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGG > REL606/3535010‑3535156 | aGTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaa > 2:600973/1‑76 (MQ=255) aGTCTATCACTTAACCTAATTCACATCACCGCCAGCAAGTATGTGTTATCTGGCGGATTGANCAATTTATCATGaa > 1:1144387/1‑76 (MQ=255) gTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaa < 2:1224731/76‑1 (MQ=255) gTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaa > 1:937743/1‑76 (MQ=255) gTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaa > 1:352215/1‑76 (MQ=255) tATCACTTAACTTAATTCACATCACTGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaac > 2:382944/1‑75 (MQ=255) tATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaac < 1:441032/76‑2 (MQ=255) aTCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaacc < 2:391716/76‑3 (MQ=255) aTCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaacc > 2:169113/1‑74 (MQ=255) tCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaaccg > 2:532247/1‑73 (MQ=255) tCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGaaaaaccg > 2:795113/1‑73 (MQ=255) cttaacttaaTTCACATCACCGCCTGGAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAAt < 2:77038/76‑1 (MQ=255) taacttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCa < 1:611612/76‑1 (MQ=255) aacttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCAt < 2:1205989/76‑1 (MQ=255) aacttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCAt > 1:563858/1‑76 (MQ=255) cttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTc < 1:1061965/76‑1 (MQ=255) cttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTc > 2:631831/1‑76 (MQ=255) ttaaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTcc < 1:539490/76‑1 (MQ=255) aaTTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCgg > 1:252262/1‑76 (MQ=255) ttCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCa < 1:650592/76‑1 (MQ=255) ttCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCa < 1:650604/76‑1 (MQ=255) tCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAg > 1:430344/1‑76 (MQ=255) caTCACCGCCTGCANGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCgg > 1:3186/1‑76 (MQ=255) aTCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGc < 2:657742/76‑1 (MQ=255) cAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCat > 2:735498/1‑76 (MQ=255) cAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCat > 2:156883/1‑76 (MQ=255) aGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCatta < 1:489840/76‑1 (MQ=255) aGTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCatta < 1:489828/76‑1 (MQ=255) gTATGTGTTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCattat > 2:181750/1‑76 (MQ=255) tgtgTTATCTGGCGGATTGAGCAATTGATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTANTgg < 1:1204118/76‑1 (MQ=255) ntgtTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGc > 1:957496/2‑76 (MQ=255) gtgtTATCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGc < 1:944812/76‑1 (MQ=255) ttATCTGGCGGATTGAGCAATTTATCAT‑NTAAAACCGAAACATTCCGGGANCGTCCAAATCCGCATTATTGGCggg > 1:195049/1‑76 (MQ=255) aTCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCa < 1:344517/76‑1 (MQ=255) tCTGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAg < 1:651686/76‑1 (MQ=255) tGGCGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTg < 2:1144387/76‑1 (MQ=255) tGGCGGANGGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAggg > 2:1075336/1‑74 (MQ=255) ggCGGATTGAGCACTTTATCNTG‑AAAAACCGAATCANTCCGGCAGCGGCCAACTCCGCATTATTGGCGGGCAGTgg < 2:1336303/76‑1 (MQ=255) cGGATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCg > 1:764349/1‑76 (MQ=255) cGGATTGAGCAATTTATCATG‑AAAAACCGAAACATACCGGCAGCGGCCAAATCCGCATTATTGGCGGGAAGTGGCg < 2:783914/76‑1 (MQ=255) ggATTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGa < 2:443742/76‑1 (MQ=255) aTTGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAgg < 2:181987/76‑1 (MQ=255) tGAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGcc < 2:550342/76‑1 (MQ=255) gAGCAATTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCg < 1:709653/76‑1 (MQ=255) aaTTTATCATGAAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTaa < 2:1269921/76‑1 (MQ=255) aaTTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGGGGCGAGGCCGCaaa > 2:876147/1‑76 (MQ=255) aTTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAAc < 1:600896/76‑1 (MQ=255) aTTTATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAAc < 1:511523/76‑1 (MQ=255) ttATCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGGGGCCGCAAAcac > 2:240405/1‑74 (MQ=255) tCATG‑AAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCg < 2:352215/76‑1 (MQ=255) cATGAAAAAA‑CGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGGCCTAAACTccccg > 2:451021/1‑74 (MQ=255) | AGTCTATCACTTAACTTAATTCACATCACCGCCTGCAAGTATGTGTTATCTGGCGGATTGAGCAATTTATCATGAAAAAACCGAATCATTCCGGCAGCGGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGG > REL606/3535010‑3535156 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |