| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,364,396 | 0 | A | T | 36.4% | 29.2 / 11.9 | 23 | P530P (CCA→CCT) | yhcD | putative outer membrane fimbrial subunit usher protein |
| Reads supporting (aligned to +/- strand): ref base A (13/1); new base T (0/8); total (13/10) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.81e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGA > NC_000913/3364183‑3364611 | ttaaaaCTGTTCTACCCCCTGTCGGGGTGGTTCAAACATAAGTGGCTATAAAACATATAGTCACGAGGGTGACAATGACACAAACAACCTTAATATTAACTGGAAAAAA‑GCTTTAAAACACC‑CC‑AGGTTTTCGGTAAACTGGCAGCATAAACTAAGGGCTAAAAAAAAAAAGGAAGACGAGGGTGATCTGTTCTCCGCAAACATCAGTATTCCTTTTGGCAGATCaaa < 2:4738/225‑1 (MQ=255) aTTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACt > 1:450738/1‑228 (MQ=255) tttCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGa > 1:529699/1‑228 (MQ=255) tGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATaaaa > 1:467496/1‑224 (MQ=255) gagcaaaaaaCATAGAGGTGAGACGGTGAACATGGCTCAAACAACCATAAAAATTACGGGAAAAAA‑GCGTTTAAACCCGCCC‑AGGGTTCCCGCAAACGGGCAGCATCAATTTGGTGCTTAAAAAAAAAATGAAGAGGAGGGTGATTTTTTCTACTAAAACATCAGTTTTCCTTTTGGCAAATCAAA‑‑CCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCa < 2:623059/222‑1 (MQ=255) ggTTACAAAGACCTAAGCCACCCTAAATATAACCGGAACAAA‑GGCTATAAACCCC‑CCAACGGTTCCCGTCAACTGGAAAAAACATTTGGGTCTTTAAAAAAATATAGAAGAGGATGGTAACCTTTTCTACGTAAACATCAGTTTTCCTTTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgt < 2:390957/228‑1 (MQ=255) cTCAAGCAACCTTAATATTAACTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACt > 1:402871/1‑194 (MQ=255) tCAAGCAACCTAAATTTTAACTGCAAAAAA‑GCTTTCAACCCCG‑CC‑AGGTTTCCCGTAAACGGGCAGCATCAATTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAA‑‑CCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTa < 2:402213/194‑1 (MQ=255) aaagcaCCTTAAATATAACCGGAACAAAGGCGTTAAAACACG‑CC‑AGGGTTCCCGTAAACTGGCAGCATCACTTTAGTGCTTAAAAAAAAAAGGAAGACGATGGTGTTCTGTTCTACGTAAACACCAGTTTTCCTTTTGGCAAAAAAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgag < 2:305197/222‑1 (MQ=255) aaGCAACCTTAATATTAACTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGcc > 1:308492/1‑154 (MQ=255) aaCTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatga > 1:124698/1‑209 (MQ=255) aCTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCaa > 1:116893/1‑129 (MQ=255) cTGAAACAAC‑GCTTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAATTTAGTGCTTAAAAAAAAAATGAAGAGGATGGTGATCTTTTCTCCGTAAACATCATTTTTCCTTTTGGCAGATCaaa < 2:116718/129‑1 (MQ=255) cGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAAc > 2:112860/1‑157 (MQ=255) tttaaaaCACG‑CA‑ACGGTTTCCGTCAATTGGAATCATCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATTTGTTCTACGTCAACATCAGTATTCCTTTTGGAAAATCAAC‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAAcc < 1:113024/153‑1 (MQ=255) ccGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGt > 1:373823/1‑180 (MQ=255) cAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGa > 1:105290/1‑94 (MQ=255) cAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATAACGACGAACGTGAAACGAGCTTGAACGGCAGTATAAGTGCCAATGTCCAtt > 1:472685/1‑228 (MQ=255) aaaaaaTAATGAAGACGATGGTTATCTGTTCTACGTAAACATCAGTTTTCCTTTTGGAAGAACAAA‑‑CACCGCCACACTGTATACTCGCCATGAc < 2:105136/92‑1 (MQ=255) aCATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACt > 1:599508/1‑72 (MQ=255) cATCAGTTTTCCATTTGCAAGACAAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTa < 2:598554/72‑1 (MQ=255) acAGTTTTCC‑TTTTGCACATCAAAACCACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGa < 2:444914/226‑1 (MQ=255) cATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCcta > 2:32586/1‑103 (MQ=255) aTTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCctac < 1:32637/103‑1 (MQ=255) | TTTAAACTGGTCTAACTCCATTCTGGGTGGTTTCAACATCAGTGGCTATAAAACATATAGTTACGACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAA‑GCGTTCAAACACG‑CC‑ACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAA‑‑CACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGAGCTTGAACGGCAGTATCAGTTCCAATCTCCATTACACCCAAGTCAGCCTTGCCGCAGGAGCAAGCGGCAGTGA > NC_000913/3364183‑3364611 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |