| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 4,586,393 | A→G | T123T (ACT→ACC) | hsdR ← | endonuclease R Type I restriction enzyme |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,586,393 | 0 | A | G | 92.0% | 78.4 / ‑1.3 | 25 | T123T (ACT→ACC) | hsdR | endonuclease R Type I restriction enzyme |
| Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (12/11); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01 | |||||||||||
CGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGA > NC_000913/4586257‑4586523 | cGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGt > 1:215767/1‑139 (MQ=255) cTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAg < 2:123198/139‑1 (MQ=255) cgcTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAg < 1:395596/126‑1 (MQ=255) cgcTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAg > 2:395596/1‑126 (MQ=255) tgttTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCa > 2:701546/1‑94 (MQ=255) tgttTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCa < 1:701546/94‑1 (MQ=255) gttTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACAt < 1:714623/139‑1 (MQ=255) gcgcagccctttCTTGGTGATAGA‑GTTTTGCCCTCGGTCCCGGAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCa < 1:501166/129‑1 (MQ=255) gCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGc < 2:717492/139‑1 (MQ=255) acGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGc > 2:447001/2‑139 (MQ=255) aGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGt > 2:407153/1‑139 (MQ=255) tCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGAt > 1:536169/1‑139 (MQ=255) ccTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATc < 1:535579/139‑1 (MQ=255) tAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGat > 2:309665/1‑139 (MQ=255) ggTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACgg > 1:395983/1‑68 (MQ=255) ggTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACgg < 2:395983/68‑1 (MQ=255) accacGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGt < 1:306199/122‑1 (MQ=255) accacGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGt > 2:306199/1‑122 (MQ=255) cGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTggg > 1:606119/1‑103 (MQ=255) cGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTggg < 2:606119/103‑1 (MQ=255) cacaAACACCGGCACCGGGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATg > 2:646592/1‑139 (MQ=255) ggcaccggGAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAAtt > 1:568283/1‑139 (MQ=255) ggAAGTCATAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGaaa > 1:141026/1‑139 (MQ=255) tAATCTTTGGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACagaga < 2:141026/139‑1 (MQ=255) aTCTTTAGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGca < 2:232881/57‑1 (MQ=255) aTCTTTAGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGca > 1:232881/1‑57 (MQ=255) | CGACTTCTGCCTGAGTCTGCGCTTTTTCTCGCACCTGCTGTTCAAGCTGTTGTTTTAGCGTCAGCACTTCCTGGTGATAGAGGTTTTCACCACGTTCCGGCAACACAAACACCGGCACCGGGAAGTCATAATCTTTAGTGACCAGACGGTAGTACCAGACAGCCAGGCGGAACCCGAGTCGCAGGCACATCTGGGCATCGTTGAGATCGTTATGATATTCGTGCACCGCCTGGTTACCAATGCGGCGTAATTTGTGAAATACAGAGA > NC_000913/4586257‑4586523 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |