| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,413,735 | 0 | A | T | 47.2% | 32.4 / 22.2 | 36 | Y2* (TAT→TAA) | ydaC | DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein |
| Reads supporting (aligned to +/- strand): ref base A (4/15); new base T (17/0); total (21/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.21e-07 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT > NC_000913/1413586‑1413883 | aTCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGGAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTta < 2:511126/150‑1 (MQ=255) tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTGAATTTtatgatgaaaa > 1:145229/1‑141 (MQ=255) tttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGC‑AATTTTTTTTATGaaaa > 1:12317/1‑147 (MQ=255) ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAacac < 1:92173/149‑1 (MQ=255) ttCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AAATTTTTAAATGaaaaac > 2:603421/1‑146 (MQ=255) tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACAccg < 2:297961/150‑2 (MQ=255) cTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTACATGAAAACCTc > 2:230108/1‑150 (MQ=255) gAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTAttt < 1:464429/151‑1 (MQ=255) aCGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTTAATGAAAACCCCCATtttttttt > 1:377159/1‑148 (MQ=255) ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTg < 2:432288/151‑1 (MQ=255) gCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACCCCTCCATTATTATTTCCAGTGGtt < 1:113580/151‑1 (MQ=255) tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTAAATTAAAACCCCCCTTTTTTTTTTCAGTgg > 1:544837/1‑146 (MQ=255) tcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTc < 2:380114/149‑1 (MQ=255) tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AAATTTATAAATTAAAACCCCCCTTATTATTTTCAGTGGTTCGtttttt > 1:618588/1‑149 (MQ=255) ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTTTATTAACAACCCCCTTATTTTTTTCAGTGGTTCTtttttt > 2:472115/1‑147 (MQ=255) ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCAt < 2:274862/151‑1 (MQ=255) agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCg < 2:301764/151‑1 (MQ=255) tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTTAATGAAAAACTCCCTTTTTATTTTCAGTGGGTGGTTTTTTTCTTCTTTTGAGTg > 1:469362/1‑150 (MQ=255) tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AAATTTTTAAATTAACACCCCCATTATTATTTTCAGTGGGTTGTTTTTTTCATTTTTCgcg > 1:603074/1‑146 (MQ=255) ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTATATAGAAACCTCCCTTTTTTTTTTTCGTGGGTCGGTTTTTTCCTCTTTTGAg > 1:624266/1‑147 (MQ=255) tttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTAAATTAACAACCCCATTATTATTTTCAGTGGGTTGTTTTTTTCCTCTTTTTAGTGgtt > 1:599734/1‑151 (MQ=255) gTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATAAATGAACACCCCCATTATTATTTTCAGTGGTTTGTTTATTTCATCTTTCGAGGGGTTCTTTTTTACTTCCCCCa‑aaa > 1:315388/1‑150 (MQ=255) ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAAcc < 1:583782/150‑1 (MQ=255) tatCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTAAATGAAAAACCCCCTTTTTTTTTTCAGTGGGTCGTTTTTTTCCTCTTTCGAGTGCTTTTTTTTTTCTTTCCCCA‑TAAACCGGTCggggt > 1:82806/1‑146 (MQ=255) tGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTAAATTAAAACCCCCCTTATTTTTTTCCGTGGTTTGTTTTTTTCATCTTTTTAGTGCTTCTTTTTTTCTTTCACCC‑TAACCGGGTCGGGGTTTGCATc > 1:217622/1‑150 (MQ=255) gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACACGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATcc < 2:146978/150‑1 (MQ=255) ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTATATTAAAACCCCCCCTTTTTTTTTCCGTGGGTTTTTTTTTTTTTCTTTTGAGGGCTTCTTTTTTTCTTTCCCCA‑TAACCCGGTCGGGGTTCGCCTTCGt > 2:380352/1‑150 (MQ=255) ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTTTAAATGAACACCCCCATTATTTTTTTCAGTGGGTTGTTTATTTCATCTTTTGCGTGCTTCTTTTTTTCTTCCTCCA‑TAACCGGTTCTGGATTCGCCTCCGt > 2:151887/1‑150 (MQ=255) ccAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGtt < 1:128878/151‑1 (MQ=255) tagtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGA‑AATTTTTTTTAATAACAACCCCCCTTTTTTTTTTAGTGGTTTTTTTTTTTCCTTTTTTGAGTGCTTCTTTTTTTCTTCCTCCATTAACCGGTT‑GGGGTTCGCATCCCtt > 2:514910/1‑146 (MQ=255) gtagGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑TATTTTTTACATAAAAACCTCCCTTATTATTTTCAGGGGGTCGGTTTTTTTTTCTTTTCTGTGCTTCTTTTTTTCTTTCACCA‑TAAACGGGTCTGGGGTCGCCTCCCtt > 2:381031/1‑144 (MQ=255) ccttccttTTCAATAGTGGCGGG‑AATTTTTTATATGAACACCTCCCTTTTTTTTTTCCGTGGTTCGTTTTTTTCCTCTTTTGCGTGGTTCTTTTTTTCTTCCACCC‑TAAACGGGTCGGGGTTTGCATCCGTTACACTTCGCTccg > 1:402823/1‑143 (MQ=255) tttCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaa < 1:465331/150‑1 (MQ=255) ttCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaaa < 2:266617/150‑1 (MQ=255) ttCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGaaa < 2:507967/150‑1 (MQ=255) aTAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCt < 1:503580/151‑1 (MQ=255) gT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGtt > 1:502261/1‑151 (MQ=255) ttATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAt > 1:607236/1‑151 (MQ=255) | ATCAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGT‑AATTTTATACATGAACACCTCCATTATTATTTCCAGTGGTTCGTTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCA‑TAACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGAT > NC_000913/1413586‑1413883 |
| 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: ‑ |