Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,315,936 | 0 | T | A | 23.6% | 131.3 / 13.0 | 55 | K34* (AAA→TAA) | yciG | protein YciG |
Reads supporting (aligned to +/- strand): ref base T (37/5); new base A (0/13); total (37/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.90e-09 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATA > NC_000913/1315802‑1316072 | aCGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑t > 2:435662/1‑140 (MQ=255) ccccGAAGTTGAATTTATTTTTATCAGCTAATTTTAATAATTAACGGTAAATAAGG‑TTTTGCTTTTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGC‑TTCAGATGCGCGTTGCGGATCTTTTTTAAAA‑TTAc < 2:469932/139‑1 (MQ=255) cGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTAccgc > 1:577776/1‑140 (MQ=255) gCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTAccgcc > 1:844605/1‑140 (MQ=255) cAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTAccgccg > 2:774269/1‑140 (MQ=255) aTCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCC‑GATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCg > 1:943110/1‑140 (MQ=255) cTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGc > 1:975635/1‑140 (MQ=255) tCAGCTGATTTTAATATTAAACGGTAAATCAGG‑TTTTGCTTGTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGC‑TTCAGATGCGCGTTGCGGATCTTTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTa < 2:844605/140‑1 (MQ=255) cTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCACCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGc > 2:791710/1‑140 (MQ=255) ttAATATTTAGGGGTGAA‑CAGGTTTTTCCCTTTTTTTTTCCCCCCGTTTTTTTGCCCCCTTTTTTTCCCCGC‑TCCAGAGGCGGGTTGCGGTTTTTTTTAAAAA‑TTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGt < 2:975635/140‑1 (MQ=255) ttAATAATTACCGGTGAATCAGG‑TTTTGCCTGTTTTATTCCCACCGTTTTGTTGCCCGCTTTTTTTCCCCGC‑TTCAGATGCGCGTTGGGGATCTTTTTAAAAA‑TTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGt < 2:622896/140‑1 (MQ=255) taatTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa > 1:316008/1‑140 (MQ=255) taatTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa > 2:522348/1‑140 (MQ=255) ttAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGc > 1:622896/1‑140 (MQ=255) tAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGcc > 2:726079/1‑140 (MQ=255) aaTCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACg > 2:410855/1‑140 (MQ=255) ccTGTTTTATTCCCCCCGTTTTGTTGCCCGCTTTTTTTCCCCGC‑TTCAGATGCGCGTTGCGGATCATTTTAAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGaa < 2:577776/140‑1 (MQ=255) ccTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGaa > 1:584532/1‑140 (MQ=255) aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAttt > 2:340765/1‑140 (MQ=255) aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAAttt > 1:449999/1‑140 (MQ=255) taCCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGaa > 2:564355/1‑140 (MQ=255) accGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccac > 1:680419/1‑140 (MQ=255) accGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccac > 1:437104/1‑140 (MQ=255) ccGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccacc > 1:518799/1‑140 (MQ=255) ccGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCATTCTCACGGTCTTCGGCGAAATTTCCTGAaccacc > 2:702311/1‑140 (MQ=255) cGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAAccacca > 2:772648/1‑140 (MQ=255) ttttGTTGCCCCCTTTTTTCCCCCCC‑TTCAGAGGCGCGTTGCGGATCTTTTTAAAAA‑TTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACg < 2:210517/139‑1 (MQ=255) cTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACg > 1:698715/1‑140 (MQ=255) cTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACg > 1:684134/1‑140 (MQ=255) tttGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGa > 1:606006/1‑140 (MQ=255) gcccccTTTTTTCCCCCTTTTCAGAGGCGGTTTGGGGATTTTTTTAAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCg < 2:698715/136‑1 (MQ=255) aCCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATgctgtc > 1:961541/1‑135 (MQ=255) cccctttttttCCCCGC‑TTCAGATGCGCGTTGCGGATCTTTTTAAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGc < 1:774269/135‑1 (MQ=255) cccccTTTTTTCCCCGC‑TTCAGAGGCGCTTGGCGGATAATTTTAAAAT‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGc < 1:591439/137‑1 (MQ=255) cGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGcc > 2:8226/1‑140 (MQ=255) cGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGcc > 2:977553/1‑140 (MQ=255) cGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGCCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGcc > 2:583710/1‑140 (MQ=255) gCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCa > 1:155491/1‑140 (MQ=255) gCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCa > 1:56237/1‑140 (MQ=255) c‑TTCAGAGGCGGGTGGGGGTTCTTTTTAAAAT‑TTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg < 1:466402/140‑1 (MQ=255) ctccAGATGCGCTTTGCGGATCTTTTTAAAAA‑TTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcg < 2:606006/137‑1 (MQ=255) ttCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgt > 1:923869/1‑140 (MQ=255) cgctTTGCGAATCTTTTTAAAAT‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCgttgcgttgc < 2:584532/140‑1 (MQ=255) gcgTTGGGGATTTTTTTAAAAT‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCa < 2:437104/140‑1 (MQ=255) cgTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgttgcgttg > 2:930164/1‑137 (MQ=255) cgTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt > 2:362376/1‑140 (MQ=255) cgTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt > 1:301568/1‑140 (MQ=255) cgTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt > 2:878691/1‑140 (MQ=255) cgTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCAt > 1:339924/1‑140 (MQ=255) cgTTGCGGATCATTTTAAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCcgttgcgttg < 1:930164/137‑1 (MQ=255) ttGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTg < 1:878691/140‑1 (MQ=255) tGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGt > 1:891581/1‑140 (MQ=255) gCGGATCTTTTTTAAAT‑TTCCCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGtt < 2:684134/140‑1 (MQ=255) gCGAATCTTTTTAAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGtt < 2:923869/140‑1 (MQ=255) gattttttttAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCa < 1:791710/135‑1 (MQ=255) tCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCAtt < 1:772648/140‑1 (MQ=255) ttttttAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGCCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTaa < 1:583710/139‑1 (MQ=255) tttttAAAATTTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTaa < 1:435662/140‑1 (MQ=255) ttttAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAata < 1:362376/140‑1 (MQ=255) | ACGGCCCGCAGATGCAATTATCTTTATCAGCTGATATTAATAATTAACGGTGAATCAGG‑ATTTGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGC‑TTCAGATGCGCGTTGCGGATCATTTTTAAAA‑TTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCACCACGATGTTCGGCCATGTTATTTCTCCCGTTGCGTTGCATTGTTTCATTAATA > NC_000913/1315802‑1316072 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |