Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 3,769,856 | T→A | 33.3% | intergenic (+74/‑91) | yibV → / → yibU | pseudogene, putative cryptic Rhs immunity protein/pseudogene, H repeat‑associated protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 3,769,856 | 0 | T | A | 33.3% | 9.6 / 8.5 | 12 | intergenic (+74/‑91) | yibV/yibU | pseudogene, putative cryptic Rhs immunity protein/pseudogene, H repeat‑associated protein |
Reads supporting (aligned to +/- strand): ref base T (6/2); new base A (2/2); total (8/4) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.47e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01 |
TATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAAATATATTAATAACTCACTCTCCAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACA > NC_000913/3769590‑3770142 | tatCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAAATATATTAATAACTCACTCTCCAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTgaagaag > 1:819907/1‑289 (MQ=255) gatgatAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAAATATATTAATAACTCACTCTCCAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGta < 1:182882/186‑1 (MQ=255) gatgatAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAAATATATTAATAACTCACTCTCCAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGta > 2:182882/1‑186 (MQ=255) tgatgaTTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAAATATATTAATAACTCACTCTCCAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGaacaac > 2:101741/1‑288 (MQ=255) tataAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAAATATATTAATAACTCACTCTCCAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATatcatc > 1:652660/1‑289 (MQ=255) tcCAGGAGTTTTTTTAACTTAGTTTTATAATGAAAACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCt < 1:616950/260‑1 (MQ=255) tcCAGGAGTTTTTTTAACTTAGTTTTATAATGAAAACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCt < 1:616954/260‑1 (MQ=255) tcCAGGAGTTTTTTTAACTTAGTTTTATAATGAAAACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCt > 2:616950/1‑260 (MQ=255) tcCAGGAGTTTTTTTAACTTAGTTTTATAATGAAAACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCt > 2:616954/1‑260 (MQ=255) cAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAAt > 2:290117/1‑289 (MQ=255) aGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAg < 2:819907/289‑1 (MQ=255) aaTACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa > 1:616042/1‑289 (MQ=255) | TATCTGAAAAATCTAAATCCACTCCAACTTGCCGCTGATTGCATTGAAACAGTAAATAAAATAAAAGATAAGAATAAAAAAAATCATTGATATAAATTCAGTTAATATTTGTAGTGATGATAAAAATATAAAATTGAGAGTTAATAGTACTATTATGGCCATTGATGATTCTATAAAGTGCATCGATGAATGATAACTTAATTTTATTTAAAAATATATTAATAACTCACTCTCCAGGAGTTTTTTTAACTTAGTTTTATAATGAATACGATAATTAAAAGTGAAGAAGAAAATAGTTGTATACATATATTCCATGACCCTACATTGTGATCTATAGATTAAATATTTTTTAAATGTGACTGTCATCCTTCAGATAATAAAGCCATCATCGCAATAGATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACA > NC_000913/3769590‑3770142 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |