Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,413,735 | 0 | A | T | 42.0% | 56.5 / 37.4 | 50 | Y2* (TAT→TAA) | ydaC | DUF1187 family protein, Rac prophage; putative double‑strand break reduction protein |
Reads supporting (aligned to +/- strand): ref base A (4/25); new base T (21/0); total (25/25) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.76e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG > NC_000913/1413588‑1413885 | cAGTTTTACTTTCTTCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGAGTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTtata < 2:51810/150‑1 (MQ=255) aGTTTTACTTTCTGCGCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTtttat > 2:536762/1‑149 (MQ=255) aCTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAac < 1:1048150/150‑1 (MQ=255) tCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑A‑AAA‑TTTATTCATTAACACCT‑c > 1:266547/1‑151 (MQ=255) gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTT‑AAT‑TTTTTTCATGAA‑ACCT‑cc > 2:645493/1‑149 (MQ=255) gCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCa < 2:649203/150‑1 (MQ=255) tctctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑TAA‑TTTTTTCATTAGAACCT‑CCctt > 2:687367/1‑150 (MQ=255) ctGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCatta < 1:794903/148‑1 (MQ=255) ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCatta‑ttatt < 1:513312/151‑1 (MQ=255) ggAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCatta‑ttat < 2:670075/150‑1 (MQ=255) ccTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTg < 1:421651/151‑1 (MQ=255) cTGCTTTTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑g < 2:874763/148‑1 (MQ=255) ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTACATGAACACCT‑CCCTtt‑tttttt > 1:817548/1‑140 (MQ=255) ttcttTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAAAAAAAACCT‑TCACTA‑TTTTTTTCA‑ATgggt > 1:697499/1‑147 (MQ=255) tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑TAAAT‑TTTTAAAATGAAAACCC‑CCCTTT‑TTTTTTCCG‑GGGGtgggtt > 1:904546/1‑145 (MQ=255) tttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑ttt < 1:731829/150‑1 (MQ=255) ttttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAA‑TTTTTTTATGAAAAACC‑CCCCTA‑TTAttttttgtgggttgtgt > 2:1109947/1‑135 (MQ=255) ttACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTAtt < 2:670095/150‑1 (MQ=255) tACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTATATGAACACCT‑CCCTTA‑TTTTTTTCatcgggttgtttttt > 1:17719/1‑135 (MQ=255) ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTACATGTACACCC‑CCCTTA‑TTTTTTTCA‑ATGGGTTG‑tttttt > 1:380391/1‑147 (MQ=255) ccTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCAt < 1:1102343/151‑1 (MQ=255) cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTATATTAAAAACC‑CCCTTT‑TTTTTTTCA‑GTGGGTCG‑ttttttt > 2:542734/1‑146 (MQ=255) cTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAATTAAAACCT‑CCCTTT‑TTATTTTCC‑GTGGGTCG‑tttttttc > 1:1037128/1‑146 (MQ=255) agCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTTAATAAACACCTACCTTTT‑TTTTTTCAT‑GTGGTTGT‑TTTTTTCcttctt > 2:1073302/1‑149 (MQ=255) cATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAg < 1:2306/151‑1 (MQ=255) tttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTg < 1:444150/150‑1 (MQ=255) ttttCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTTAAAAAACACCC‑CCCTTT‑TTTTTTTCT‑GTGGTTTG‑TTTTTTTCCTTTTTTGAGt > 1:1067449/1‑148 (MQ=255) cGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGcttct < 1:772348/150‑1 (MQ=255) aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGGTT‑AAT‑TTTTTA‑ATTAACACCC‑CCCTTT‑TTTTTTTCT‑GTGGTTTG‑TTTTTTTCATTTTTTTTGTGCtttttttt > 2:249596/1‑151 (MQ=255) aTTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑TAT‑TTTATTAATGAACACAT‑CCCTTATTTTTTTCCT‑GGGGTTGG‑TTTTTTctttcttttggaggtttttttt > 2:616191/1‑137 (MQ=255) tCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAA‑TTTTTATATGTAAAACC‑CCCTTA‑TTTTTTTTATGGGGTTGG‑TTATTTCCTTTTTTC‑AGTGCTTCttttt > 1:364178/1‑149 (MQ=255) gATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AATATTTTTATAAGAAAACAC‑CATTTA‑TTTTTTCTA‑‑TGGTTTGTTTTATTCTA‑CTTTCTAGTGCCTttttttttttt > 2:1046190/1‑146 (MQ=255) aTTTGGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATAAATAAAAACAT‑CCATTT‑TTTTTTTCA‑GTGGGTGG‑GTTTTTTTATCTTTTGAGTGGTTGtttttttctt > 2:379832/1‑145 (MQ=255) tCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑G‑AAT‑TTTTTAAAAGAAAACCC‑CCACTA‑TTTTTTTCC‑ATGGGTCG‑TTTTTTTCTTTTTTTTAGTGCGTttttttttctt > 2:753933/1‑142 (MQ=255) ttatttttcaTTTTGAATATCTTGTCCATTTAGTAGAAGTGCCACCTTCCTTTTCAATAGTGGCGC‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GGGGTTCT‑TTTATTCCATCTTTCGATTGCTTCTTTTTCATTTccac < 2:880934/139‑1 (MQ=255) ttAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AAcc < 1:1022052/150‑1 (MQ=255) aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCgg < 2:904546/150‑1 (MQ=255) aGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGt < 1:458344/151‑1 (MQ=255) aTTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑CAT‑TTTTTATAAATAAACCC‑CCACTT‑TTATTTTCC‑CTAGGTCG‑TTTTTTTCCTCTTTTGTCGGCTTTTTTTTTTCTTCCACCCT‑AAcccg > 1:276474/1‑146 (MQ=255) ttttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCCCTTCCACCAT‑AACCg < 2:747797/146‑1 (MQ=255) ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCgg < 1:981593/150‑1 (MQ=255) ttGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCgg < 1:300861/150‑1 (MQ=255) tGTCCAGTTAGTAGGAGGGCCACCTTCCTTTTCAATAGTGGCGG‑T‑GAT‑TTTTTTATTTAATAACT‑CCCTTC‑TTTTTTTTT‑GTTGTTCG‑GTTTTTTCTTTTTTCCCATGCTTCTTTTTTTCTTTCTCCCT‑AATCGGCTTggggtgaccct > 2:962340/1‑141 (MQ=255) gTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATcc < 2:141344/150‑1 (MQ=255) tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTACATGAACAACC‑CCACTATTTTTTTCAA‑GTGG‑TCG‑TTTTTTTTATCTTTTTTGTTCTTTTTTTTTTCTTTCACCCC‑AACCCGGTCGGGGTTCTCATcc > 2:205122/1‑149 (MQ=255) tCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAATTAACACCC‑CCCTTT‑TTATTTTCA‑GTGGGTCG‑TTTTTTTCCTCTTTCTTGGGCTTCTTTTTTTCTTTCCCCAT‑AACCAGTTTGGGGTTCGCATcc > 1:925962/1‑149 (MQ=255) agGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTTTAAAAGAACACCC‑CCCCTT‑TTATTTTCA‑GTGGGTTG‑TTTTTTTCTTTTTTTGAGGGCTTTTTTTTTTTTTCCCCCCTAAACCGGTTGGGGGTT‑GCATCCCTTACACtt > 1:1170449/1‑148 (MQ=255) gAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCt < 2:328860/150‑1 (MQ=255) gTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCgg < 1:502659/151‑1 (MQ=255) gCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTa < 2:855017/151‑1 (MQ=255) tttCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGa < 2:774537/149‑1 (MQ=255) gCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCAcc > 1:382050/1‑151 (MQ=255) ataCATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATgg < 1:760778/151‑1 (MQ=255) ataCATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg < 2:704889/150‑1 (MQ=255) ataCATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATg < 1:1209361/150‑1 (MQ=255) | CAGTTTTACTTTCTGCTCTCTGGAAACGCCTGCTTCTTTTTTACCTGAGAGCATTTTTTCGCATTCTGATTTCGTTAGTTTAGATTTTGAATATCTTGTCCAGTTAGTAGGAGTGCCACCTTCCTTTTCAATAGTGGCGG‑T‑AAT‑TTTATACATGAACACCT‑CCATTA‑TTATTTCCA‑GTGGTTCG‑TTTATTCCATCTTTCGAGTGCTTCTTTTTCACTTCCACCAT‑AACCGGTTCGGGATTCGCATCCGTTACACTTCGCTCGGTAATATCCTGAAATGGCTTTCACCGTTACTGATGG > NC_000913/1413588‑1413885 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |