| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 4,043,030 | 0 | A | G | 62.5% | 5.5 / 22.8 | 24 | N206D (AAT→GAT) | srkA | Cpx stress response Thr/Ser protein kinase; MazF antagonist protein |
| Reads supporting (aligned to +/- strand): ref base A (5/4); new base G (15/0); total (20/4) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.19e-02 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-05 | |||||||||||
TGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGA > NC_000913/4042890‑4043172 | tGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGGTACATGGAGACCGCCACCGCGGGGATATtctc > 1:158740/1‑149 (MQ=255) acTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATTGAGAATGCCACGCCGGGGAAATTTTCTGGcgcggagggccca > 1:514288/1‑140 (MQ=255) cTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTTCAAGGAGACCGCCACCCCGGGGATATTTTCTGGGGCGGAGGTTCCAtgt > 2:465510/1‑151 (MQ=255) gggTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCAACGTTCTGCGGGTAAATGGAGACTGCCACGCCGGGGAAATTTTCTGGGGCGAGGGGCCAATGTTTTTTGATATgg > 1:159807/1‑151 (MQ=255) gTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCGTCTGCGGGTACATGGGGACTGCCACGCCGGGGATATTTTCTGGCGCGGTGGGCCAATGGTTTTTGATTTgg > 2:382710/1‑149 (MQ=255) aaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGc < 2:801871/151‑1 (MQ=255) ggCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGGCCAATGTTTGTTGATCTGGATGATGCACg > 2:330689/1‑150 (MQ=255) gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGGTCTTCGGCTACATGGGGAATGCCACGCCGGGGATATTCTCTGGCGCGGTGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTcc > 1:730085/1‑150 (MQ=255) gAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAAATTTTGACGTTCTTTGGGCACATCGTGACTGGCCCGCCGGGGAAATTTTCCGGGGCCACGGTCCACTTGTTTTTGATCTTGATGATGGACGTAATGGTccc > 1:601971/1‑150 (MQ=255) aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCCACATTGGGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTTTTTATCTGGGTGATGCACGGAATGGTCCAGcc > 2:181071/1‑151 (MQ=255) aGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTCAACGTTCTGGGGGGACAAGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGGTGGTCCCATGTTTGTTGGTCTGGGTGATGCCAGGAATGGGCCCGCCGtt > 1:35664/1‑148 (MQ=255) tGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTCCCGTTTTTCGGCTACATGGAGACTGCCACCCCGGGGCTATTCTTTGGCGCGATGGTGCCATGTTTGGTGATCTGGATGATGCACGTAATGGGTCCGCCGTTTAGGGTTtgtg > 1:230264/1‑151 (MQ=255) gccgccGTTACAGCACACTGGCGGGAAGATTTCACCGGTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTTTTTATCTGGATGATGCACGTAATGGTCCAGCCGTTCCGGGTTTTTGGgtgt > 2:469463/1‑151 (MQ=255) ccgTTACAGCACACTGGCGGGAAGATTTTACCGGTCTGCGGCCACATTGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTTTTGATCTGGATGATGCACGGAATGGTCCAGCCGTTCAGGGTTTTTGGATGtt > 1:625016/1‑148 (MQ=255) ccgTTACAGCACACTGGCGGGAAGATTTCACCGTTCTTCGGCTACATGGAGACGGCCACGCCGGGGATATTTTCTGGCGCGGTGGTCCCATGTTTGTTGGTCTGGGTGATGCCCGGAATGGGCCCGCCGGTCAGGATTTTTGGGTGTTg > 1:233744/1‑149 (MQ=255) ggCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGc > 2:838058/1‑151 (MQ=255) aaGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGagca < 1:583302/151‑1 (MQ=255) tttCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTATGGCGCGATGGTCCAATGGTTGTTGGTCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCgg > 2:165726/1‑151 (MQ=255) cACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATg < 1:181043/151‑1 (MQ=255) ccGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGc > 1:400292/1‑150 (MQ=255) ggAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTg < 1:86327/151‑1 (MQ=255) gaCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGa > 2:439814/1‑149 (MQ=255) gCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATga > 2:611090/1‑149 (MQ=255) ccGGGAATATTCTCTGGCGCGATGGTACAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATga > 1:499811/1‑148 (MQ=255) | TGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGA > NC_000913/4042890‑4043172 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |