| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,413,735 | 0 | A | T | 45.7% | 33.5 / 21.6 | 35 | Y2* (TAT→TAA) | ydaC | DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein |
| Reads supporting (aligned to +/- strand): ref base A (2/17); new base T (16/0); total (18/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 4.19e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
TTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATG > NC_000913/1413593‑1413884 | ttACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATg > 1:223230/1‑149 (MQ=255) tGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCCTTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGATTTTTTTACAGGAAAAcccc > 1:79228/1‑147 (MQ=255) gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACAGCTc < 2:619254/148‑1 (MQ=255) ctctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAAAACCTCCctt > 1:332365/1‑151 (MQ=255) ctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATTAAAACcccc > 2:2504/1‑146 (MQ=255) ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattat < 2:189022/149‑1 (MQ=255) ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATTAACACCTCCcttattttt > 2:344352/1‑151 (MQ=255) ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattattat < 2:352007/150‑1 (MQ=255) ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCattatt < 2:568539/148‑1 (MQ=255) aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCTCCCTTAtttttt > 2:164843/1‑148 (MQ=255) tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑Tc < 2:171974/149‑1 (MQ=255) tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCg < 2:372212/150‑1 (MQ=255) cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTTATTTTTTTTATGAAAACCTCCCTTATTATTTTTCGTGGT‑TTGtttttttc > 2:209728/1‑146 (MQ=255) gagaGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTTAATGAAAACCCCCATTTTTTTTTTCAGTGGG‑TCGtttttt > 2:32943/1‑144 (MQ=255) gCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGa < 1:499041/151‑1 (MQ=255) cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATGAAAACCCCCATTATTTTTTCTCGTGGG‑TTGTTTTTTTCATCTTTTTTGTGCTtttttt > 2:567807/1‑151 (MQ=255) cATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCtttt < 1:24966/151‑1 (MQ=255) aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTAATTAAAAACCCCCTTTTTTTTTACAATGGTGTCGTTT‑TTTCATCttttaaacgccttttttt > 2:232894/1‑136 (MQ=255) tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTTCATTAAAACCCCCCCTATTTTTTTCCTTGGG‑TCGTTTTTTTCATTTTTTGAGTGCTTCtttttttctttccc > 2:635736/1‑142 (MQ=255) tAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTATATGAACACCCCCCTTTTTTTTTTCCGTGGG‑TTGTTTTTTTCCTCTTTTGAGTGCTTTTTTTTTTC‑TTccacc > 1:17784/1‑148 (MQ=255) ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGCGCCACCTTCCTTTTCAATAGTGGCGGTGAGTTTATTAATGGAGACCTCCCGTAAAATTTCCACCTGG‑TCGTTTTTTTCACCTTTCCTGTGCTTTTTTTTCAC‑TTTCACCCTaaatc > 2:474383/1‑148 (MQ=255) tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCg < 2:332271/150‑1 (MQ=255) tAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGGAATTTTTTTAATGAAAACCTCCCTTTTTTTTTTTAATGGG‑TCGTTTATTTCTCCTTTTTAGTGCTTTTTTTTTTC‑TTCCTCCATAACCgg > 2:309095/1‑151 (MQ=255) aaTATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGAtt < 1:473563/151‑1 (MQ=255) ttGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTACATTAACACCCCCATTATTTTTTTCAGTGGG‑TCGTTTTTTTCCTCTTTTCAGTGGTTCTTTTTTTC‑TTTCCCCATTACCGGGTCCGGGTTCGCCTc > 1:433664/1‑151 (MQ=255) ttGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAAATTTTTAAATAAACACCTCCATTATTTTTTTCAGTGGG‑TTGTTTTTTTCATCTTTTTAGTGGTTCTTTTTTACTTTCCCCCAAAACCGGTTGGGGtttgcat > 1:429758/1‑146 (MQ=255) tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTTTAAATTAAAACCCCCCTTTTTTTTTTCAATGGG‑TTGTTTTTTTCCTCTTTTGAGTGCTTTTTTTTTTC‑TTTCCCCCTAAACCGTTTGGGGTTCGCATcc > 1:7370/1‑149 (MQ=255) tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTacac < 2:498574/150‑1 (MQ=255) agtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACtt < 1:690690/151‑1 (MQ=255) gtatGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑CTCCACCATAGCCGGTTCGGGATTCGCATCCGTTACACTTc < 2:238756/151‑1 (MQ=255) gCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTa < 2:493477/151‑1 (MQ=255) tCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAAt < 1:543558/150‑1 (MQ=255) aTTTTATACATGAACACCGCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTATCGAGTGCTTCTTCTTCGC‑TTCCTTCATAACCGGTTCGGGATTCGCGTCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCTCCGTTa < 2:496746/149‑1 (MQ=255) tttATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTg < 1:124662/150‑1 (MQ=255) ttATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt > 2:67019/1‑151 (MQ=255) tataCATGACCACCTCCATTATCATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATCACCGGTTCGTGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGa > 2:444601/1‑149 (MQ=255) ataCATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg < 1:606995/150‑1 (MQ=255) | TTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTAATTTTATACATGAACACCTCCATTATTATTTCCAGTGGT‑TCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCAC‑TTCCACCATAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATG > NC_000913/1413593‑1413884 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |