| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 1,825,061 | 0 | G | A | 25.5% | 95.8 / 19.9 | 51 | intergenic (‑124/+79) | ves/spy | HutD family protein Ves/ATP‑independent periplasmic chaperone |
| Reads supporting (aligned to +/- strand): ref base G (35/3); new base A (0/13); total (35/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.18e-09 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TATCAAAGTATTCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTC > NC_000913/1824922‑1825197 | tatCAAAGTATTCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGagaaaag > 2:1454037/1‑140 (MQ=255)tatCAAAGTATTCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGAGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGagaaaag > 1:327441/1‑140 (MQ=255) atCAAAGTATTCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGagaaaaga > 1:1654549/1‑140 (MQ=255) atCAAAGTATTCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGagaaaaga > 1:1613426/1‑140 (MQ=255) gttttccttccTCGCCCCCGCCCT‑TTTTTTTTGT‑TTGTTATCGTACGGTGACCCGGGG‑ATGTTTTCCGGCACCGGGGGTGATTTCCGTGGGTGGTGATCTGGTTTTTTTT‑ACAACTGTTTG‑AAAAAAAAAAAGAAAAccg < 2:1613426/137‑1 (MQ=255) tCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGa > 1:746259/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGa > 1:80353/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGa > 1:852329/1‑140 (MQ=255) tCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGa > 1:1265469/1‑140 (MQ=255) cATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATc > 1:382604/1‑140 (MQ=255) tCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTc > 1:567290/1‑140 (MQ=255) gctcTCCGGGCCC‑TTTCTTTCCT‑TTTTTATC‑TTGGAGTAACCGGGGTATGTTTAGCGACACCTGGGGTGATTTACGTTAGTGGTCATGGTGTTTTTTTT‑ACACCTGTTTG‑AAAAGAAAAAAGAAATCCGCCGATCCTGTc < 2:683791/137‑1 (MQ=255) gCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCa > 2:1814113/1‑140 (MQ=255) ccTCCCGGCCC‑GTTTTTTTCA‑TTTTTATCGTTCAGGTACCCGTGG‑TTGTTTACCGGCAACTGGGGTGATTTACGTTAGGGGGGATCAGGTTTTTTTT‑ACAATTGTTGG‑AGAAAAAAAAAAAAAACCGCCGATCCTGTCCa < 2:723054/139‑1 (MQ=255) cTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:723054/1‑140 (MQ=255) ctccGGCCCTGTTTTTTTCT‑TTGTTATCGTGCGGTTAACGGGGG‑TTGGTTAGCGGAATCTCGGGTTATTTTGTTTATGGCTGACCAGG‑TTTTTTT‑CAAATGGATTGAAGAAAAAAAAAAAAAACCGCCGATCCTGTCCAc < 2:695274/138‑1 (MQ=255) gCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt > 1:695274/1‑140 (MQ=255) agacagtttttCG‑TTGTTATCTTTCGGTTAACGGTGG‑TTGGTTAGCGGCAACTGGGGTGTTTTACGTTGGTGGTGATCAGCTTTTTTTT‑ACAACTTTTTG‑AGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTa < 2:1557912/133‑1 (MQ=255) tttttCATTTGTTTTTGTTCAGCAAACGGGGG‑ATGGTTATCGAAATCTGGGTTTTTTTTCGTTAGTGGT‑CTCGGGCTTTTTTTAACATTTGTTTG‑AAAAGAGAAAAAAAATCCGCCGATCCTGTCCACCGCATTACTGCa < 2:1678324/139‑1 (MQ=255) ttCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt > 1:1443183/1‑140 (MQ=255) ttttttttATCGTTCAGGTACCGGTGG‑GTGTTTATGGGCAACTGGGGTGATTTACGTTAGGGGTGATCGGTTTTTTTTT‑ACAATTTTTGG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTa < 2:1756442/134‑1 (MQ=255) ta‑TTTTTATCGTTCAGGTAACGGGGG‑ATGTTTTTCGGCAACGGGGGTGTTTTACGTTGGTGGTGATCGGGTTTTTTTT‑ACAATTGTTTG‑AAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg < 2:852329/139‑1 (MQ=255) tATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAgctgtct > 1:1557912/1‑134 (MQ=255) cAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcgg > 1:82674/1‑140 (MQ=255) ggTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg > 1:87749/1‑140 (MQ=255) ggTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg > 2:998148/1‑140 (MQ=255) gTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg > 2:1180061/1‑140 (MQ=255) aTGGTTATCGGCAACTCGGGGGTTTTACGTTAGGGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGAAGTGGCCAAGACCGGCGGCCTTAAGttt > 2:510498/1‑140 (MQ=255) tATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTgg > 1:135170/1‑140 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGa > 1:561583/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 2:737103/1‑140 (MQ=255) cTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGat > 2:1391333/1‑140 (MQ=255) gttttACGTTGGTGGTGATCGGGTTTTTTTT‑ACAATTGTTGG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcag < 2:135310/138‑1 (MQ=255) cGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTg > 2:559911/1‑140 (MQ=255) ttggtggtgATCGGGTTTTTTTT‑ACAACTTTTTG‑AGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 1:494405/137‑1 (MQ=255) ttAGTGGTGATCGGGTTTTTTTT‑ACAACTGTTTG‑AAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 2:87749/140‑1 (MQ=255) ttAGTGGTGATCAGGTTTTTTTT‑ACAATTGTTTG‑AGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCa < 2:933604/140‑1 (MQ=255) gggtggtgATCAGTTTTTTTTT‑ACAACTGTTTG‑AAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAg < 2:91812/138‑1 (MQ=255) tggtgATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCa > 2:11838/1‑140 (MQ=255) gtgATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTTAAAAATTATTCAGCAGTTGCAAGCAtt > 1:652552/1‑140 (MQ=255) aTCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTa > 1:1376260/1‑140 (MQ=255) ggCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtc > 1:91812/1‑139 (MQ=255) gCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCttt > 2:1663723/1‑140 (MQ=255) gCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 1:500533/1‑140 (MQ=255) gCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCttt > 1:1684235/1‑140 (MQ=255) gCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTCGCTGAAAGATTATTCAGGAATTTCAGGCATTTTAccctg > 1:1255775/1‑137 (MQ=255) gCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGATTTTTGGCCGGAAGATTATTGCGCAGTTGCAAGCATTTTACCttt > 2:251540/1‑140 (MQ=255) ttttt‑ACATGTGTTGA‑AAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGcc < 2:82674/140‑1 (MQ=255) cAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTGCCGCTggg > 2:416337/1‑139 (MQ=255) aCTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACg < 1:1391333/140‑1 (MQ=255) ttGTTGG‑AAAAAAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGt < 2:724564/139‑1 (MQ=255) tttG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:724564/1‑140 (MQ=255) tttG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:1782491/1‑140 (MQ=255) tttG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCt > 1:222284/1‑140 (MQ=255) aaaaagaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTc < 2:500533/137‑1 (MQ=255) agagaaAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc > 1:685671/1‑140 (MQ=255) aaaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCtt > 1:261539/1‑140 (MQ=255) aaagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTtc < 1:251540/140‑1 (MQ=255) | TATCAAAGTATTCCATACTCGCTCCCGGCCC‑GTTCTTTTCA‑TTGTTATCGTTCAGGTAACCGTGG‑ATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTT‑ACAACTGTTTG‑AGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTC > NC_000913/1824922‑1825197 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |