| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,483,254 | 0 | G | A | 31.2% | 38.2 / 10.3 | 32 | S29L (TCA→TTA) | emrK | tripartite efflux pump membrane fusion protein EmrK |
| Reads supporting (aligned to +/- strand): ref base G (10/12); new base A (0/10); total (10/22) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.31e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
AGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAC > NC_000913/2483008‑2483493 | aGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGaaa < 1:506512/250‑1 (MQ=255) gTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTACTTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCaataa > 2:532116/1‑250 (MQ=255) gAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACt > 2:311788/1‑251 (MQ=255) ggTGCGATTGGGGCATCAGTTTTATCCAGTGAATCTAAAATGTCCCCTGGTCGAACGTGGTTCGTTTCTTTATTATTAACGACAGTGACACTCCCTGAGATTTGTGAAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATTTCTTCTAATTCCATTGACAAAT‑AGAAATAGCCACCTGAAACCGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGa < 2:516614/250‑1 (MQ=255) aGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGtt > 2:70290/1‑248 (MQ=255) gTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCaat < 1:327620/205‑1 (MQ=255) gTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCaat > 2:327620/1‑205 (MQ=255) gggATAGTGGCATAAGTTTAATCCAGTGAAACTAAAATGTCACGTGGTCGAACGTAGTTCGTATTTTTATAATAAACAAAAGTGACACTACCGGAGACTTGTACAGAAATTGAACCTGCATCCCCCGTGAAATAGGCGTA‑ATCTGTACT‑AATCATGTCTTATAATTCAAGGGACAAAT‑AGGAAAAGGCACCTAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATg < 1:311788/249‑1 (MQ=255) gaggggcatatgTTTTATCCAGTGAAACTAAAATGTCCCTTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTTCCTGAGACTTGTGCAGAAATTGGATCTGCTTTCCCCGGAACATAGGCGTC‑ATCTGTACT‑AAACATGTCTTCTATTTCCATTGCCAAAT‑AGGCATAGGCACTTGAAAAGGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGtttt < 1:230678/240‑1 (MQ=255) aTCAGTTTTACCCAGTAAAACTAAAATGTCACCTTGTCGAACGTGGTTCGTATTTTTATGATTAACGACAGTGACACTCCCTGAGATTTGTGCAGAAATGGGATCTGCATCCCCGGTGACATAGGCGTC‑ATCTTTACT‑AATCATGTCTTCTATTCCCATTGACAAAT‑AGGCATAGGCACCGAAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTg < 1:518098/251‑1 (MQ=255) cAGTTTTATCCAGTGAAACTAAAATGTCACTTTGTCAAACGTAGTTCTTATCTTTATGATTAACGACAGTGACCCTCCCTGATACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATGGGCTTC‑TTCTGTACT‑AATCATGTCTTCTAATTCCATTGCCCAAT‑AGGCATAGGCCCCTGAAAACGCAATAAATAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGa < 1:533854/250‑1 (MQ=255) gTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑TTCTGTACA‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTAAAAACGCAATAAAAAAAATTCCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAAtt < 1:329094/251‑1 (MQ=255) tATCCAGAGAAACTAAAAGGTCACTTGGTAGAACGTATTTCGTATCTTTATGATTAAGAACAGTGTCACTCCCGGAGACTTTAGCAAAAATTGTATCGGCATGCCCCGTGAAATAGGCTAC‑ATCGTTACT‑AACCAGGTCTTAAAATTCAATGCACCAAT‑AGGCAAAGGCTCCTAAAACCAAAAAAAAGAATACTACCGCCAATAAAGAAAAGTATTTTCTTCGGTTAGAATGTTTTTTATTTGAATTaag < 2:397821/250‑2 (MQ=255) gtccctttGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACCCTACCTGAGATTTGTGCAAAAATTGGATCTCAATTCCCGGTGACAAGGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCCCCTGAAAACGCAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTAtctctc < 1:263983/245‑1 (MQ=255) tttgcgaaggtggTTCTTATCTTTATGATTACCGACGTTGACCCTCCCGGGGATTTGTGCAGAATTTGGATCTCCATTCCCCGTAACAAAGGCTTC‑ATCTG‑ACTTAATAATGTCTTCTAATCCAATTGACCAAA‑AGGAATAGGCACTTGAAACCGCAATAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTtc < 2:405488/238‑1 (MQ=255) cGTATCTTTATGATTACCGACAGTGACACTCCCTGAGACTTGTGCAAAAATTGGTTCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACA‑AATCATTTCTTCTAATTCCAGTGACCAAT‑AGACAAAGACACCTAAAACCGCAATAAAAAAAATTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAAttt < 2:218526/250‑1 (MQ=255) ttATGATTAACGACAGTGACACTACCTGAGACTTGTGCCGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTAtgaattg > 2:424955/1‑251 (MQ=255) tagaataaaCGACAGTGACACTACATGAGGCTTGTCAAAAAATTGGATCAGCATCCCCGGGGCCATAGG‑GTCCATCTTTACA‑AATCATGTTTTCTAATTCCATTGACAAAT‑AGGCATAGGCCCCTGAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATCTAtgaattg < 1:541911/243‑1 (MQ=255) aaCGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGaaa > 1:377867/1‑250 (MQ=255) gacatgaCCTGAGACTTGTGTAGAATTTGTATCTGCATCCCCCGTGACAGAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATACAATTAACCAAT‑AGGAATAGAAACATAAAAACGCAATAAAAAAAACTTCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGcattc < 2:28517/244‑1 (MQ=255) gagaCTTGTGCAGAAATTGGATCTGCATTCCCCGTAACATAGGCGTC‑TTCTGTACA‑AATCATGTCTTCAAATTCCATTGACCAAA‑AGGCATAGGCACCTGAAAACGCAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCat < 1:71553/249‑1 (MQ=255) tgggCAGAATTTGGATCTGTATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATGCCATTGCCCAAT‑AGGAATAGACACCTAAAAACGCAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTGCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCAGTCATTCTTCATATCaaa < 1:52432/246‑1 (MQ=255) aGAAATTGGATCTGCTTCCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTTTTCTAATTCCATTGACAAAT‑AGGCATAGGCACCTAAAAACGAAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTg < 1:34572/251‑1 (MQ=255) gATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGtt > 2:500322/1‑251 (MQ=255) cGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATccc > 2:394712/1‑250 (MQ=255) gTC‑ATCTTTACA‑AATCATGTTTTCGATTCCCATTG‑CCAATTGGGCAGAGGCCCCTAAAAAGGCAAAAAAAAAAACTGCCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGTATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTagagcgggggggtgggag < 1:115070/251‑18 (MQ=255) gTACT‑AATCATTTCTTCAAATACAATTGACAAAT‑AGGAAAAGCCACATAAAAACGAAATAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAt < 2:138672/251‑1 (MQ=255) gTACT‑AATAATGTCTTCTAATTCCATTGACAAAT‑AGGCATAGGAACCTGAAAACGAAAAAAAAAAACCTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAAt < 1:64416/251‑1 (MQ=255) gTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGt > 2:243282/1‑251 (MQ=255) tctAATTCCATTGACCAAT‑AGGCAGAGGCACCTGAAAACCAAAAAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTg < 2:347524/248‑1 (MQ=255) aggcaaagccaCCAGAAAACGAAAAAAAAAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTa < 1:243282/251‑1 (MQ=255) cataggcaCCTGTAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAc > 2:228800/1‑251 (MQ=255) | AGTTTATTCGTTTGCCGAACAATATTTGCCAGATTATTTTTAGCTTTATTGAGTGCGATAGTGGCATCAGTTTTATCCAGTGAAACTAAAATGTCACCTTGTCGAACGTAGTTCGTATCTTTATGATTAACGACAGTGACACTACCTGAGACTTGTGCAGAAATTGGATCTGCATTCCCCGTGACATAGGCGTC‑ATCTGTACT‑AATCATGTCTTCTAATTCCATTGACCAAT‑AGGCATAGGCACCTGAAAACGCAATAAATAAAACTACCGCCAATAAAGAAAAGTATTTTCTTCTGTTAGAATGTTTTTTATTTGAATTAATCTGTTCCACTATTATCTCTCATTTCTCATAGATGAAATTTATGAATTGAATGAAAAGAGCATTCATTCTTCATATCAAACACAGTGCATGTGTTGCACATTATCCCAGGAACGTTGGAGGCCTAATTTAAATGAGGTGCAAAAAAAACAACAGAGGACTAAAC > NC_000913/2483008‑2483493 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |