Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NC_000913 | 2,131,667 | T→C | D135G (GAC→GGC) | wcaC ← | putative glycosyl transferase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 2,131,667 | 0 | T | C | 92.5% | 122.8 / ‑0.1 | 40 | D135G (GAC→GGC) | wcaC | putative glycosyl transferase |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (20/17); total (21/19) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.96e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01 |
CCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCT > NC_000913/2131535‑2131795 | cccGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCa < 2:64426/139‑1 (MQ=255) gcgcGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGa > 1:180579/1‑139 (MQ=255) cTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTa < 1:1274292/139‑1 (MQ=255) cggATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTc > 2:377569/1‑139 (MQ=255) cggATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTc > 1:376108/1‑139 (MQ=255) gATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCAc > 1:209851/1‑139 (MQ=255) ggTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAg > 1:271842/1‑139 (MQ=255) tCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAgcg < 2:180579/139‑1 (MQ=255) tCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAgcg > 1:308966/1‑139 (MQ=255) ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac > 2:142113/1‑139 (MQ=255) ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCaa > 2:86671/1‑131 (MQ=255) ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCaa < 1:86671/131‑1 (MQ=255) ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac > 2:749788/1‑139 (MQ=255) ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac > 1:1204076/1‑139 (MQ=255) cATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca < 2:1204076/139‑1 (MQ=255) ccTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac > 1:53503/1‑128 (MQ=255) ccTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac < 2:53503/128‑1 (MQ=255) cTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGAtt < 1:377569/139‑1 (MQ=255) cTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGAtt > 1:855737/1‑139 (MQ=255) aGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTa > 2:463112/1‑139 (MQ=255) ttCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac < 2:485271/114‑1 (MQ=255) ttCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac > 1:485271/1‑114 (MQ=255) cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTc > 1:1200960/1‑73 (MQ=255) cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTc < 2:1200960/73‑1 (MQ=255) cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg > 1:475888/1‑113 (MQ=255) cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg < 2:475888/113‑1 (MQ=255) cgcAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAg < 1:142113/139‑1 (MQ=255) cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAg < 2:300036/139‑1 (MQ=255) cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAg < 1:462979/139‑1 (MQ=255) gcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTgcag > 2:1264682/1‑123 (MQ=255) gcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTgcag < 1:1264682/123‑1 (MQ=255) cAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGa > 2:1311458/1‑139 (MQ=255) aGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGac > 1:1023318/1‑139 (MQ=255) acacTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCg < 1:1311458/139‑1 (MQ=255) acTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCt < 2:303550/139‑1 (MQ=255) acTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCt < 2:917613/139‑1 (MQ=255) tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGc < 1:786231/78‑1 (MQ=255) tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGc > 2:786231/1‑78 (MQ=255) gTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac < 2:942248/72‑1 (MQ=255) gTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcaccac > 1:942248/1‑72 (MQ=255) | CCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCT > NC_000913/2131535‑2131795 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |