| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 501,564 | 0 | A | C | 32.4% | 62.1 / 11.4 | 37 | *559E (TAA→GAA) | ybaL | putative transporter YbaL |
| Reads supporting (aligned to +/- strand): ref base A (24/1); new base C (0/12); total (24/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.02e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCA > NC_000913/501322‑501796 | tatGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcc > 2:148786/1‑250 (MQ=255) tGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacc > 1:85028/1‑251 (MQ=255) gCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTACCGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGGAGAGTCTTACTGGGCCCGTCAACTTATCGTCGGCTCGCAGGGC‑AGATAAGGCGTTCACGCCGCATCTGGCATATGCCTCTCGATGCCTGATGCGACGCTGGCGCGCCTTATCATGCCTACCTATTTTTGATATTTTACATCCGGCAACCACCGTTTACCCAGTcacc > 1:439493/1‑250 (MQ=255) gCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTcacc > 1:436549/1‑250 (MQ=255) gtgAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccg < 1:307610/251‑1 (MQ=255) tgAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACccgc > 2:345014/1‑251 (MQ=255) aGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTAGCGCGTCTAATCATGCCTACATATTTTTCATATTTTACATCCTGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCgg > 1:39610/1‑250 (MQ=255) tccccctcgtgaaaCGCGCGGCTTGCGGGGGCGTA‑AAGACAGCGTGAAAGATTTTTACTGGGATCGTTAATTTTTCTTCG‑TTCCTTGGCCCAAATAGGGCTTTCACGCCGCATCTGGCATTTGGCTCGCGATGCCTGATGCGGCGCTGGCGCGTCTTAGCATCCCTTCATTTTTTTCATTTTTTACACCCGGCAACCACCGTTTCCCCCGTCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAg < 1:66670/238‑1 (MQ=255) tttaccgggcgGCTTGCCGG‑GGGTGGAAGACAGCGTGGAGGAGTTTTACGGGGTTCTTTAATTTACCGGCGGTTCGTAGCCC‑AGATAGGGCGTTCACGCCGCATCTGGCTTTTGGCTCTCGAGGCCTGATGCGACGGTGGCGCGTCTTATCTGGCCTACATATTTTTCATATTTTACACCCGGCAACCCCCGTTTCCCCCGTCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACg < 1:62615/242‑1 (MQ=255) gcgGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTc > 2:437657/1‑251 (MQ=255) g‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCCTGGTACGGGCGATTTCACGCTCGCCCATCACtc > 2:92580/1‑250 (MQ=255) gcttttccgttcttCTTACGGGGATCTTTAAGTTATCGTCGGTTCGTAGGCC‑AAATAAGGCGTCCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGCCGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGTCACCCCCCCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg < 1:400840/237‑1 (MQ=255) agTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACg > 1:93208/1‑251 (MQ=255) tCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTg > 2:400840/1‑237 (MQ=255) cTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGACGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGACTAAATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCACGTCACCACCTCACCCGCCGGGGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCgg > 1:73522/1‑249 (MQ=255) gTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGc > 1:255013/1‑251 (MQ=255) cgGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACCTATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTc > 1:207062/1‑250 (MQ=255) gTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATATGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCAt > 1:179507/1‑250 (MQ=255) gATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGc > 1:424261/1‑251 (MQ=255) tAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCa > 2:359006/1‑250 (MQ=255) gttccccgcggcttttggttttgggCTCTCGATGCTGGATGCGACGCGGGGGGGTTTTATCATGCCTACATATTTTTCATATTTTCCACCCGGCACCCCCCGTTCCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtc < 1:146527/228‑1 (MQ=255) tcccgccgcATCTGGCGTTTGGCTTTCGATGCGTTATGCGCCGCTGGCGCTTTTTCTCATGCCTTCTTTTTTTTCTTTTTTGACATCCGGCACCCACCGTTTCCCCCCCCCCCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAAtctc < 2:268098/247‑1 (MQ=255) aTTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTcg > 1:51277/1‑251 (MQ=255) ggCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTcgcgcg > 1:40054/1‑251 (MQ=255) gCTCTCGTTGCATGATGCGCCTCTGGCGTGTTTTACCATTCCTAAATATTTTTAATATTTAACACCCGGCACCCACCGTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCgg < 1:359006/251‑1 (MQ=255) cctcgcGATGCCTGAGGCGCCGCTGGCGCTTTTTCTCTTGCCTACATTTTTTTCATATTTTCCACCCGGCACCCCCCTTTTCCCCCGCCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCGCAATATCCGGATTTTTCGCGCgg < 2:51277/246‑1 (MQ=255) ctctCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCggg > 1:232097/1‑251 (MQ=255) tGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATCTTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGc > 1:208334/1‑251 (MQ=255) gACGCTGGCGCGTCTTATCATGCCTACATATTCTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGACGGAGGAGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAAtctc > 1:51977/1‑250 (MQ=255) ttATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATCATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAAcc > 1:298040/1‑251 (MQ=255) tttCATCCCTACATATTTTTCATATTTTCCATCCGGCACCCCCCTTTTCCCCCGTCCCCACCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATaa < 1:251347/246‑1 (MQ=255) tATCATGCCTACATTTTTTTCATTTTTTCCCCCCGGCAACCACCTTTTCCCCCCTCCCCCCCCCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAAcc < 1:254083/250‑1 (MQ=255) ccATGCCTACATTTTTTTCATATTTTACACCCGGCACCCCCCTTTTCCCCCGCCACCCCCTCCCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGtt < 1:430735/250‑1 (MQ=255) tatTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTc > 2:329000/1‑250 (MQ=255) ttttCATTTTTTCCACCCGCCACCCCCCTTTTCCCCCGTCCCCACCCCTCCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTGTCGCGCGGGCAGATGCCACAATCTCACCCGCTTTATAACCGTTGGGAATCGTCAgg < 1:148786/250‑1 (MQ=255) cccggacAACCCCCTTTTCCCCCGGCCCCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATTGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATTGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACATTCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGca < 1:305779/245‑1 (MQ=255) ccaccGTTTACCCCGTCACCACATCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGCTTCGCACCACGTTCGGTGATATACGCAACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATcc > 1:386341/1‑164 (MQ=255) | TATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTG‑AAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCC‑AGATAAGGCGTTCACGCCGCATCTGGCATTTGGCTCTCGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACATATTTTTCATATTTTACATCCGGCAACCACCGTTTACCCCGTCACCACCTCACCCGCCGGTGGCGTTTCCAGCAGTTCCAGCATGGTACGGGCGATTTCACGCTCGCCCATCACTACCTGATTCGCACCACGTTCGGTGATATACGCCACTTCATCGTCATAATGGGCGCGGGCAATAATCTCAATATCCGGATTTTTCGCGCGGGCAGATGCCACAATCTCACCCGCTTCATAACCGTTGGGAATCGTCAGGATCAGCCATTTTGCA > NC_000913/501322‑501796 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |