Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 1,018,977 | C→A | intergenic (+40/+36) | matP → / ← ompA | macrodomain Ter protein/outer membrane protein A |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 1,018,977 | 0 | C | A | 90.0% | 51.5 / ‑2.2 | 20 | intergenic (+40/+36) | matP/ompA | macrodomain Ter protein/outer membrane protein A |
Reads supporting (aligned to +/- strand): ref base C (2/0); new base A (10/8); total (12/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.95e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01 |
TGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGC > NC_000913/1018733‑1019222 | tGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCgc > 1:230313/1‑250 (MQ=255) cgcAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTAcc < 2:356560/251‑1 (MQ=255) cTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAACAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTAc > 1:361045/1‑250 (MQ=255) cTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACgagatcggaagacacacgtctgaactcc > 1:187783/1‑227 (MQ=255) cTGGCAACGGCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACg < 2:187783/224‑1 (MQ=255) gggcTTGCGCAGCGTCGCGGTAAAACGCGTGCTGAAACGATTGTACATCTGATGTAAGATGCGGAAAACAAAGAGAAATACGCGAAGAAAATGTCTTCACTGAAGCAAGATCTGCAGGCATTGCTGGGTAAGGAATAACGGACGAAAGTCAGTTCAATTTCCTAAAGGCAAAAAAAAACCCGGAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCtt < 2:172488/248‑1 (MQ=255) aaCGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGcgacg > 2:335216/1‑250 (MQ=255) cGCTTTCTGAAACGAATTTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACGAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTGTTTCTACCAGACGAGAACGTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATc < 2:140700/251‑1 (MQ=255) ttCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCgg < 2:228475/250‑1 (MQ=255) cTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGc > 2:93504/1‑251 (MQ=255) aaCGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGCGAAATACGCGAATAAACTGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAACCAgg > 2:298956/1‑251 (MQ=255) gCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGt > 2:169975/1‑251 (MQ=255) aTGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAGGGAATCACTGACGAAAGTCAGTTCAATTTACTCAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGCCGCGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGCGCAGCACGCTg > 1:414088/1‑251 (MQ=255) cGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGttt < 1:173218/251‑1 (MQ=255) cAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTg > 1:126638/1‑251 (MQ=255) cTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCAc > 1:369246/1‑251 (MQ=255) aTAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGAAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTAGCCCATACCACGTGCGGAGATCTTGTCTGCCggg > 1:388687/1‑251 (MQ=255) cGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCttt > 1:49711/1‑251 (MQ=255) aaaGTCAGTTCAATTTACTAAAGGCAAAAAAAAACCCGCAGCATCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGACAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTgg < 1:88582/251‑1 (MQ=255) tACTAAAGGCAAAAAAAAACCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCGGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATca < 1:217430/250‑1 (MQ=255) aaaaaaCCCGCAGCAGCGGGGTTTTTCTACCAGAAGAGAACTTAAGCATGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCATTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGTAGATCAGGTAATCAACCACAGACTGAGc > 1:116997/1‑251 (MQ=255) | TGCGGAACACCAGCATACGCGCAAAAAATCGATCGATCTGGAATTTATTGTCTGGCAACGTCTGGCTGGTCTTGCGCAGCGTCGCGGTAAAACGCTTTCTGAAACGATTGTTCAGCTGATTGAAGATGCGGAAAACAAAGAGAAATACGCGAATAAAATGTCTTCTCTGAAGCAGGATCTGCAGGCATTGCTGGGTAAGGAATAACTGACGAAAGTCAGTTCAATTTACTAAAGGCAAAAAAAACCCCGCAGCAGCGGGGTTTTTCTACCAGACGAGAACTTAAGCCTGCGGCTGAGTTACAACGTCTTTGATACCTTTAACTTCGATCTCTACGCGACGATCCGGAGCCAGGCAGTCGATCAGTGCAGCACGCTGTTTCACGTTGTCACAGGTGTTGCCAGTAACCGGGTTGGATTCGCCCATACCACGTGCGGAGATCTTGTCTGCCGGGATACCTTTGGAGATCAGGTAATCAACAACAGACTGAGC > NC_000913/1018733‑1019222 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |