| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,808,417 | (A)8→9 | intergenic (‑319/‑123) | waaQ ← / → waaA | lipopolysaccharide core heptosyltransferase 3/KDO transferase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,808,409 | 1 | . | A | 95.8% | 81.7 / ‑3.3 | 24 | intergenic (‑311/‑131) | waaQ/waaA | lipopolysaccharide core heptosyltransferase 3/KDO transferase |
| Reads supporting (aligned to +/- strand): ref base . (0/1); new base A (17/6); total (17/7) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.92e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
AAATAGCATTTCGATAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑AAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACGGGAAAAGT > NC_000913/3808330‑3808477 | aaaTAGCATTTCGATAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑a < 2:681546/81‑1 (MQ=255) gCATTTCGATAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑aaaaaa < 1:681548/81‑1 (MQ=255) gCATTTCGATAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑aaaaa > 2:681576/1‑80 (MQ=255) tttCGATAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑aaaaaaaaa > 1:681579/1‑80 (MQ=255) gATAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑aaaaaaaaatg > 1:681582/1‑76 (MQ=255) aTAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑aaaaaaaaatg > 1:681567/1‑75 (MQ=255) tAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑aaaaaaaaatg > 1:681588/1‑74 (MQ=255) aTGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATg > 1:681574/1‑81 (MQ=255) tGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGt < 2:681548/81‑1 (MQ=255) gTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATg > 2:681572/1‑79 (MQ=255) gTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATg > 1:681585/1‑79 (MQ=255) cAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGc > 1:681580/1‑81 (MQ=255) aGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGc < 2:681549/80‑1 (MQ=255) aCGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCgg > 1:681581/1‑79 (MQ=255) aTTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAAc < 1:681551/81‑1 (MQ=255) tttGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACt > 2:681596/1‑81 (MQ=255) gggCCGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTAttt < 1:681552/81‑1 (MQ=255) tccGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTAttt > 2:681571/2‑79 (MQ=255) ccGAGTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAg > 2:681586/1‑81 (MQ=255) gTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCaa < 1:681553/81‑1 (MQ=255) gTTTCAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCa > 1:681598/1‑80 (MQ=255) cAATGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAAc > 2:681590/1‑81 (MQ=255) cAATGAATCGCATAAAAATAC‑AAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAAc < 2:681553/80‑1 (MQ=255) aaTGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACg > 2:681594/1‑81 (MQ=255) tGAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACgg > 2:681605/1‑80 (MQ=255) gAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACggg > 2:681591/1‑80 (MQ=255) gAATCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACggg > 2:681610/1‑80 (MQ=255) aTCGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACGGGaa > 1:681613/1‑80 (MQ=255) cGCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACGGGAAAAg > 2:681600/1‑81 (MQ=255) gCATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACGGGAAAAGt < 2:681554/81‑1 (MQ=255) cATAAAAATACAAAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACGGGAAAAGt > 1:681584/1‑80 (MQ=255) | AAATAGCATTTCGATAATGTCAGTACGGCCTATGTAGATTTGAAAATGGGGCCGAGTTTCAATGAATCGCATAAAAATAC‑AAAAAAAATGAAGATGTAGCGGCAGATCAGACCAAACTGAAGCTATTTAAGTCAAAAACGGGAAAAGT > NC_000913/3808330‑3808477 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |