| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 1,592,373 | C→A | 21.3% | intergenic (‑129/+501) | yneO ← / ← lsrK | pseudogene, AidA homolog/autoinducer‑2 (AI‑2) kinase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 1,592,373 | 0 | C | A | 21.3% | 71.8 / 4.3 | 28 | intergenic (‑129/+501) | yneO/lsrK | pseudogene, AidA homolog/autoinducer‑2 (AI‑2) kinase |
| Reads supporting (aligned to +/- strand): ref base C (12/10); new base A (3/3); total (15/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01 | |||||||||||
AGCAGGCCTGAAATACCTGTAGAGTGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTGTATTGAGTGTTAGTTGGAGGTGGGAAGGTGTCAATTTGGCTAATTTGCTAATGCTTTCGAATTGATAGATTAAAAATTGCATAGTTCTTATTTTATTTAAA > CP009273/1592286‑1592533 | aGCAGGCCTGAAATACCTGTAGAGTGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTGTATTGAGTGTTAGTTGGAGGTGGGAAGGTGACAATTTGGCTAATTTGCTAATGCTTTCGAATTGATAGATTAAAAATTGCATAGTTCTTATTTTATTTaaa > 1:372891/1‑248 (MQ=255) cTGTAGAGTGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAgg > 1:227660/1‑127 (MQ=255) gTAGAGTGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGt > 2:58169/1‑119 (MQ=255) agagTGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTg < 2:227656/127‑1 (MQ=255) agTGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCaga < 1:58169/119‑1 (MQ=255) tGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTGTATTGAGTGTTAGTTGGAGGTGGGAAGGTGTCAATTTGGCTAAt > 2:620676/1‑167 (MQ=255) aTTCCATATCACGCTATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTGTATTGAGTGTTAGTTGGAGGTGGGAAGGTGTCAATTTGGCTAATTTGc < 1:620687/167‑1 (MQ=255) cgcgATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAAtt > 2:418341/1‑96 (MQ=255) cgcgATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGt > 2:272415/1‑98 (MQ=255) cgcgATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCAACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACtta > 1:122113/1‑89 (MQ=255) cgcgATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCAACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACtta > 1:97941/1‑89 (MQ=255) gcgATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTAAGACATTa > 1:37077/1‑75 (MQ=255) aTAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCaga < 1:272420/98‑1 (MQ=255) aTAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCa < 1:418346/96‑1 (MQ=255) aTAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCAACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTaa < 2:97941/89‑1 (MQ=255) aTAGATCCTATTCATATAAACGCTCCATATACAAACAATACAGCAACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTaa < 2:122112/89‑1 (MQ=255) tAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTAAGACATTaaaaa < 2:37077/75‑1 (MQ=255) tCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCATAGGTCTGt > 1:426060/1‑82 (MQ=255) ccATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTaaaaa > 1:278986/1‑52 (MQ=255) ccATATACAAACAATACAGCAACGAAAGAGATGAgtc > 1:92071/1‑37 (MQ=255) tataCAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCATAGGTCTGTATTg < 2:426054/82‑1 (MQ=255) ataCAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGaa < 2:278981/52‑1 (MQ=255) ataCAAACAATACAGCAACGAAAGAGATGAGTCGTCa < 2:92071/37‑1 (MQ=255) aacaaTACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTGt > 1:150280/1‑72 (MQ=255) aTACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAAtt > 2:408693/1‑55 (MQ=255) aTACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTGTATTg < 2:150277/72‑1 (MQ=255) aGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAAtt > 1:566231/1‑51 (MQ=255) aGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCa < 1:408698/55‑1 (MQ=255) | AGCAGGCCTGAAATACCTGTAGAGTGCAATTCCATATCACGCGATAGATTCTATTCATATAAACGCTCCATATACAAACAATACAGCCACGAAAGAGATGAGTCGTCAGACATTAAAAATGAAACTTATTAAATTGTCAGAGGTCTGTATTGAGTGTTAGTTGGAGGTGGGAAGGTGTCAATTTGGCTAATTTGCTAATGCTTTCGAATTGATAGATTAAAAATTGCATAGTTCTTATTTTATTTAAA > CP009273/1592286‑1592533 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |