| Predicted mutation | ||||||
|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | annotation | gene | description |
| RA | NC_000913 | 3,886,834 | G→A | D3N (GAT→AAT) | mnmE → | 5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 3,886,834 | 0 | G | A | 80.0% | 19.7 / 2.0 | 15 | D3N (GAT→AAT) | mnmE | 5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit |
| Reads supporting (aligned to +/- strand): ref base G (0/3); new base A (6/6); total (6/9) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 2.29e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01 | |||||||||||
CCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGC > NC_000913/3886646‑3887071 | ccATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCGGCGAGAAGAAAAAATCCTGATTCGGTGAGATTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCAAAGCGCATCAGCCACCCAAAGCAACAGGAACATCATGAGCAATAATGACCATATCGTAGCCCAGGACACGCCTCCGGAACGTggc > 1:28210/1‑233 (MQ=255) agcagcTGNTNNNCNNTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATNCTGATTCGGTGAGTTNTCNCTAANATANGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCAATAATGACACTATCGTan > 2:223617/1‑194 (MQ=255) agcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCAATAATGACACTATCGTAg < 1:223617/195‑1 (MQ=255) gTGGTCTGGAAAAACGTGGCCTGCATAGCCGAGAGAAGAAAAAATCCTGATTCGGCGAGTTTTCGCTAACCTAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCAACACAAAGCAACAGGAACATCATGAGCAATAATGacac > 1:122261/1‑172 (MQ=255) ggAAAGACGTGGCCTGCATAGCCGCGAGAGGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGACTGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATTAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGGGGCGGCGTTGGCATCCTGCGCATCTCCGGCCTCAAAGCCCGTGAAGTTGCCg < 1:168641/248‑1 (MQ=255) ncTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTACGATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCCCACAAAGCACCAGGAACATCATGAGCAATAATGCCACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTg > 1:61676/2‑247 (MQ=255) ggTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCATTTTTTCTTTTCGTAGGGCGGATAAGCACTGCGCATCCGCCACACAAAGCAACAGGCACATCATGAGCAATAATGACACTATCGTAGACCAGGCAACGCCTCAGGGCCGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGcc > 1:50664/1‑201 (MQ=255) tCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCATTTCGTAGGGCGTATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCAATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCTGCGTTGGCATCCTGCGCATCTCCTGCTTCAAAGCCCGTGAAGTTGCCGAAACCGGGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAg < 1:290077/250‑1 (MQ=255) agggcgttaatttGACCGCCTTTTTTCTTTTCGTAGGGCGAATAAGCACCGCCCATCCGCCACACAAAGCAACAGGACCATCATGAGCAATAATGACACTATCGTAGCCCAGGCCACGCCGCCGGGACGTGGCGGCGTTGGCACGCTGCGCATCTCCGGCTTCAAGCCCCGTGAAGTTGCCGAAACCGTGCTGCGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGCCGGCa < 1:305814/240‑1 (MQ=255) gACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCa < 1:160028/251‑1 (MQ=255) gATAAGCACCGCGCATCCGCCGCACAAAGCAACAGGATCATCATGAGCAATAATGACACGGTCGTAGCCCAGTTCACGCCTCCGGGACGTGGGGGCGTTGGCATCCTGCGCATCTCCGGGTTCAAAGcc < 1:266011/129‑1 (MQ=255) acacaAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCATTCCACGCCCCCTGGACGTGGCGGCGtt < 1:270695/78‑1 (MQ=255) caAAGCAACAGGAACATCATGAGCAATAATGACACTATCGTAGCCCAGGcaac > 1:230567/1‑50 (MQ=255) caAAGCAACAGGAACATCATGAGCAATAATGACACTANNGNNNCNCAGGCCAc < 2:230567/53‑1 (MQ=38) gAACATCATGAGCAATAATGACCCTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCGACGCCGATTATCTTCCGTGTAAAGCCGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGc < 1:233682/251‑1 (MQ=255) | CCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGC > NC_000913/3886646‑3887071 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |