| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,483,254 | 0 | G | A | 54.5% | 0.1 / 10.3 | 11 | S29L (TCA→TTA) | emrK | tripartite efflux pump membrane fusion protein EmrK |
| Reads supporting (aligned to +/- strand): ref base G (4/1); new base A (0/6); total (4/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.52e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-03 | |||||||||||
| Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACC‑TGAGACTTGTGCAGAAATTGGATCTGCATTCCCC‑GTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTG > NC_000913/2483048‑2483468 | tAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACC‑TGAGACTTGTGCAGAAATTGGATCTGCATTCCCC‑GTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTAtttt > 1:99097/1‑251 (MQ=255) tttattgagggggttggTGGCTTCATTTTTATCCAGTGAAACTAAAAGGCCCCTTTTTCGAACGTATTTCTTATCTTTATGATTAACGACAGTGACACTACC‑TGAGATTTGTGAAGAAATTGGATCTGCATTCCCC‑GTGAAATAGGCGTCATCTGTACAAATAATGTCTTCAAATCCCATTGACCAATAGACATAGGCACCTAAAAACGAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTttctt < 2:203993/234‑1 (MQ=255) agggaaacaaaattttCACTTTGTCGAAGGGAGTTGGTATTTTTATGATTAACGCCGGGGACACTCCCGTG‑GATGTGTGAAAAAATGGGATGTGCATGCCCC‑GTGACATTGGCGTCATGTGTACGAATCATTTTTCCTAATCCCATGGACAAATGGGCATAGGCACCTAAAAACGCAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTc < 2:206620/236‑1 (MQ=255) ccTTGTAGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTCCC‑TGAGATTTGTGCAGAAATTGGATCTGTATTCCCC‑GTAACATAGGCGTCATCTGTACTAATCATTTCTTCAAATCCCATTGACCAATAGGCATAGGCACCTAAAAACGCAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCAttt < 1:69157/251‑1 (MQ=255) caaaggtggTTCGTATCTTTAGAATAAACGACAGTGACACTACC‑TGAGACTTGGTCATAAAGTGGATCT‑CATTCCCCCGTGACAAAGCCGTCATCTGTACAAATCATGTCTTCTAATCCCATGAACCAATAGGCATAGGCACTTAAAAACGTAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATa < 2:135929/243‑1 (MQ=255) tCTTTATGATTAACGACAGTGACACTCCC‑TGAGACTTGTGAAGAAATTGGATCTGCATTCCCC‑GGGACATAGGCGTCATCTGTACTAATAATGTCTTCAAATTCCATTGACCAATAGGCATAGGCACCTAAAAACACAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtga < 1:1835/250‑1 (MQ=255) tCTTTATGATTAACGACAGTGACACTACC‑TGAGACTTGTGCAGAAATTGGATCTGCATTCCCC‑GTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtga > 2:200564/1‑250 (MQ=255) aaTTAACGACAGGCACACGACA‑TGAGACTTGCGCAAAAATTGGCTCTGCATGCCCC‑GTAACATAGGCGTCATCTGTACTAATTATTTCTTCTAATTCCATTACCCAAGAGAAAAAGACACCTGAAAACACAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTCAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATGTATGAATtgaatg < 2:121238/250‑1 (MQ=255) aTAGGCGTCATCTGTCATAATCATGTCTTCAAATCCCATTGACCAATAGGCATAGGCACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACg < 1:54189/251‑1 (MQ=255) cATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATg > 1:34379/1‑250 (MQ=255) tcttctAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTg > 1:7910/1‑251 (MQ=255) | TAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACC‑TGAGACTTGTGCAGAAATTGGATCTGCATTCCCC‑GTGACATAGGCGTCATCTGTACTAATCATGTCTTCTAATTCCATTGACCAATAGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTG > NC_000913/2483048‑2483468 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |