Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,825,061 | 0 | G | A | 37.2% | 51.4 / 27.1 | 43 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
Reads supporting (aligned to +/- strand): ref base G (25/2); new base A (0/16); total (25/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.77e-10 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAC > NC_000913/1824923‑1825192 | atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑agaaaaga > 1:1117183/1‑140 (MQ=255) atCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑agaaaaga > 2:162189/1‑140 (MQ=255) tATTCCATCCTCGCCCCCGCCCCGTTTTTTTTATTGTTATCGTTCAGGTACCGGTGGATGGTTAGGGGCAACTGGGGTGATTTACGTTAGGGGTGATCGGGTTTTTTTTACAACTTTTGGAGAAG‑AAAAAAAAAAAccgc < 2:1452764/140‑1 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGa > 1:1292370/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGa > 1:546276/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTCATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGa > 2:731617/1‑140 (MQ=255) tCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTc > 1:1403081/1‑140 (MQ=255) gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCa > 1:1256361/1‑140 (MQ=255) gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCa > 1:294962/1‑140 (MQ=255) gCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCa > 1:1028436/1‑140 (MQ=255) ccccccGGCCTGTTTTTTTTATTGTTCTCTTACGGTTCCCGGTGGTTGTTTATCGGAACCCGGGGTTATTGACGTTAGGGGTGTTCAGGCTTTTTTCCCAGCTTTTGGAAAAA‑AGAAAAAAAAACCGCCGATCCTGTCCa < 2:1253948/137‑1 (MQ=255) cccGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCg > 1:1105908/1‑140 (MQ=255) cGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCa > 2:1478039/1‑140 (MQ=255) gcctagTTTTTTTAATTGTTGTCTTACAGTTACCGGTGGTGTTTTATCGGAACCTGGGGTTTTTTGCGTTGGGGGTGACCATTTTTTTTTTAAAATTTTTTGAGAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCAt < 2:1117183/135‑1 (MQ=255) ttCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc > 1:1253948/1‑140 (MQ=255) ttttGTTATCGTTCAGGTAACGGTGGATGGTTATCGGCAACTGGGGTGATTTACGTTAGGGGTGATCGGGTTTTTTTTACAACTGTTGGAGAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg < 2:1105908/138‑1 (MQ=255) tGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTgg > 1:329887/1‑140 (MQ=255) tCAGGTAACGGTGGATGTTTATCGGCAACTGGGGTGTTTTACGTTAGTGGTGATCAGGTTTTTTTTACAATTGTTGGAGAAA‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcg < 2:294962/140‑1 (MQ=255) ggggTAACCGTGGATGTTTATCGGCAACTGGGGTTATTTACTTTAGTGGTGATCAGTTTTTTTTTACAATTGTTTGAAAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg < 2:971767/138‑1 (MQ=255) ggTACCCGTGGATGGTTAGCGGCACCTGGGGTGATTTACGTTAGTGGTGATCAGTTTTTTTTAACAATTGTTTGAAAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg < 2:1292370/140‑1 (MQ=255) ccGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCtt > 2:629134/1‑140 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtt > 2:976432/1‑140 (MQ=255) aTGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGtct > 1:971767/1‑138 (MQ=255) ttATCGGCAACTGGGGTGTTTTACTTTGGGGGTTATCAGGTTTTTTTTACAACTGTTGGAGAAA‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg < 2:329887/140‑1 (MQ=255) ttATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTg > 2:433627/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 1:570518/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 1:977529/1‑140 (MQ=255) cTCGGGTGA‑TTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatt > 2:1531667/1‑140 (MQ=255) ggggTGTTTTAGGTTGGTGGTGATCGGTTTTTTTTTACAATTGTTTGAAAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGatta < 2:1748307/139‑1 (MQ=255) tgttttACTTTAGTGGTGATCAGTTTTTTTTTACAATTTTTTGAGAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTca < 2:977529/137‑1 (MQ=255) gttttCCGTTGGTGGTGATCAGGTTTTTTTTAC‑ACTTTTTGGAAAGAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:1256361/138‑1 (MQ=255) tACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGt > 1:1748307/1‑140 (MQ=255) gtgATCAGTTTTTTTTAACAATTTTTGGAAAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAtt < 2:94950/140‑1 (MQ=255) aTCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAAATTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttttt > 2:981366/1‑139 (MQ=255) gttttttttACAACTGTTTGAGAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt < 1:748483/138‑1 (MQ=255) gCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 1:488567/1‑140 (MQ=255) tttttACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGTAGTTGCAGGCATTTTACCTTTTGcc < 1:629134/140‑1 (MQ=255) ttttACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAATTTTTTTGCTGAAAGAGTATTCAGCAGTTTCAGGCCATTTACCTTTTGCCg > 2:1039437/1‑140 (MQ=255) tttGTTGGAAAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 2:1245666/138‑1 (MQ=255) tcTGTTTGAGAAG‑AGAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 2:546276/139‑1 (MQ=255) tttGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt > 1:1542478/1‑138 (MQ=255) tttGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 2:678387/1‑140 (MQ=255) tttGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:1245666/1‑140 (MQ=255) tttGAGAAG‑AAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGtt < 2:1542478/138‑1 (MQ=255) ag‑agaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc < 1:1531667/140‑1 (MQ=255) | ATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTCATTGTTATCGTTCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAG‑AGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAC > NC_000913/1824923‑1825192 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |