| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_000913 | 2,978,164 | T→A | 28.6% | Y245F (TAT→TTT) | lysA ← | diaminopimelate decarboxylase, PLP‑binding |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_000913 | 2,978,164 | 0 | T | A | 28.6% | 41.1 / 2.5 | 33 | Y245F (TAT→TTT) | lysA | diaminopimelate decarboxylase, PLP‑binding |
| Reads supporting (aligned to +/- strand): ref base T (9/14); new base A (3/7); total (12/21) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 7.10e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
CCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACC‑TGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACC > NC_000913/2977927‑2978398 | ccATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACC‑TGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCata < 1:129069/239‑1 (MQ=255) gtgtggggggtAACCCCCGTAAAGTGCGGGGCGCACCAGATTGTGTACCCCGGCACCAACCAGACAAAAGGGGCGCCTCCCGATTTGTTTGACGCTCCGCACC‑TGAGTAATAAAGACGCAAGACTGCGCTACCAGGAAGCGACCGGGTTCAATTTCCAGTTCCACAGGGTGGCCAAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTACCAACCg < 2:210138/231‑1 (MQ=255) gtAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACCCCCCCCACC‑TGCGGCATTAATACGCCCGACTGCGCTACCAAGAAGCGACCCCGGTCAATTTCCAGTTTCACAGGGTGGCCCCAATGGCGGGCGAGCTGCTCACACCCGGCATCCCACA‑CACCATAAAAAAATTTGGGGTCAACCGCCCCCTCa > 2:540102/1‑238 (MQ=255) gtcaatggCCGGGCGCATCAGATCGTTGAACCCGGAATCAACCAGCAAAAAGTGGCGGCCCCCCATTTGTTTGACGCTCCGCACC‑TGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTCCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGAt < 1:610182/232‑1 (MQ=255) tGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCGCCCCATTTGTTTGCCGCTCCCCACC‑TGAGGAATTAATAAGCCCGACTGCGCTACCCGGCAGCGGCCCGGGTCACTTTTCAGCTTCACAGAGGGGCCCCACTGGCGGGCGACTTGCTCACGCCCGGCCTTTCCAA‑GAACCACATAATAATCGTGATCCACAGCCTCTTTACCCCCTCGCTAAGGag > 2:701423/1‑238 (MQ=255) ttgtaccccGGCAACACCCATCAAAAAGGGGCGCCCCCCTTTTTGTTTGGCCCTCCGCACT‑GGAGAAATTAATACGCGAGACTGCGCTACAAGGAAGCGCCCCGGTTCATTTCCCAGTTTCACGGGGTGGCCAAAAGGGCGGGCGATTTGCCCACGCGCGGCATCCCACA‑GACCAAAAAAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCAccgcccg < 2:210380/233‑1 (MQ=255) ggCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTTACGCCCCGCACC‑TGAGTAATTAAAACGCCAGACTGCGCTACCAGGAAGCGACCCGGGTCAATTTCCAGTTTCACAGGGGGGCCCAAAAGGCGGGCGCGTTTCTCACCCCCGGCATTCCCCAGGCCCAAAAAAATGTT‑GGGATCAACCGCCCCCTCAACCCCTTTAAAAAGAACAaaaaaaccaacaccccccaaaaaaa > 2:331006/1‑218 (MQ=255) ggCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGGCGCTCCGCACC‑TGAGTAATTAATACGCCCGACTGCGCTACCAGGAAGCGACCCCGTTCAATTTTCAGTTTCACAGGGTGGCCCCAAAGGCGGGGGATTTGCTCACCCCCCGCATTCCCCA‑ACCCAAAATAAAATTCGGTGTCAACCGCCCCTTGAACCTCCTGTTAAGAAAAAGAAAACCCACCCCCCGCAaaaaaa > 2:506604/1‑236 (MQ=255) gcccccacccccccaAAGTGGGGCCCCCCCCTTTTTTTTGCGCCCCCCCCCGTGAGTAATTAACCGGCCAGACTGGCCCCCCGGGAAGCGCCCCGGTTCAATTCCCATTTTCAAGGGGTGGCCAAAATGGCGGGCGATTGGCCCACGCCCGGCATCCCACA‑GACAAAAAAAATGTCCGGTATCACCCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAg < 2:674181/226‑1 (MQ=255) cTCCCCATTTGTTTGACGCTCCCCACT‑GGAGTAATGAATCCGCCAGCTGGCGTTACCAGGAAGCGACCCGGTTCAATTTCCAGTTCCAAAGGGTGGCCAAAAGGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GCCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGa < 1:540102/238‑1 (MQ=255) tttttgcgtccccGCACG‑GGTGTAATTAACGCGCCAGACTGCGCTCCCGGGAAGCGCCCCGGTTAATTTCCCATTTCCACGGGGTGGCCAAAATGGGGGGAGATTTGCTCCCGCGCGCCATCCCACA‑GACAATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACg < 2:914258/228‑1 (MQ=255) tgtgtgACGCCCCGCACC‑TGAGTAATTAATCCGCCAGACTGCGCTACCGGGAAGCCCCCCGTTTCAATTCCCAGTTTCACAGGGTGCCCAAAAGGGCGGGCGATTTGCCCACGCGCGGCTTTCCACA‑GACCATAATAATGTCCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACg < 1:234772/239‑1 (MQ=255) tgACGCTCCGCACC‑TGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACCCGCGGCATTCCACA‑GACCATAATACTGTGCGGTATCAACCACCTCTTAACCCTGTTGATAATGAACGAAAAGCCCAACCCCCCCAGAAATAGCCTGTAAATCCTGACCGCATTCGGAGACCCGACGCAcc > 2:779469/1‑239 (MQ=255) cccccGCACC‑TTAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTCCCAGTTTCACGGGGTGGCCAAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAAGGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg < 1:665663/236‑1 (MQ=255) ccTCCGCCCC‑TGAGTAATTAATACGCCAGACTGCGCTACAAGGAAGCGACCCGGTTAAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGCCGATTTGCTCACGCGCGGCTTTCCACA‑GACAATAATAAGGTTCGGTACCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAg < 1:335456/238‑1 (MQ=255) cTCCGCACC‑TGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACC‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAACCCCCCCGCCCGCCGACATATCCTGTAAATCCTGACCGAATTCGATGACCGTACGCCCCAAAg > 1:802178/1‑238 (MQ=255) ccagcgccccccccgggaccgcaCCCGGTTCAATTCCGATTTGCACAGGGTGCCCAAAATGGGGGGCGTTTTGCCCCCCCGGGGCATCCCAAA‑GACCAAAAAAATGTTCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCa < 2:426189/217‑1 (MQ=255) cgcTACCAGGAAGCGACCCGGTTTAATTTCCAGTTTCCCAGGGAGGCCCAAATGGCGGGCGATTTGCTCACGCCCGGCATTCCTCC‑GAACATAATAATGTTCGGTATCATCCGCCCCTTCACCCTGTTGATCAGGAACAGAAAGGCCACCGCCCGCAGGACTAGCCTGTCACTCATGACCGCATTCGCTGACCCGAGCAACCCAAGCCCCACAAACCCCTTCCAAGTGGGCATataca > 2:273922/1‑234 (MQ=255) cTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCCGACATAGCCTGTAAATCCTGACCGAATTCGATAACCTCACGCACCATAACACCCAACACCCGTTCCACATGGGAGAAATCaaa > 2:873170/1‑237 (MQ=255) cTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCCTAGCACCACACACCTGTTCCAGATGGGCATAATCAAc > 1:426189/1‑238 (MQ=255) tACAAGGAAGCGCCCCGTTTCAATTCCCAGTTTCACAGGGTGGCCCAATTGGCGGGAGATTTCCCCCCGCGCGGAATTCAACA‑TACAATAATAAGGTTCGTTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGc < 1:723380/239‑1 (MQ=255) caaaTTCCCAGTTCCCATGGGTGACACAAAGGGCGGGCGATTGGCTCACGCGCGTCTTGCCACA‑GACCAAAATAATTTCCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCAtgt < 1:771139/236‑1 (MQ=255) acggggtaccccAAAGGGCGGGCGATTTCCTCACGCGCGGCATCCCACA‑GACCATAATAATGTCCGTTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTg < 1:273922/230‑1 (MQ=255) tGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATCTGCATGTAAATGCCGACCAGCTTCAgatg > 1:855590/1‑239 (MQ=255) gCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGACTTCGATGACCTGACGCACCATAGCACCCCACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAAAGTGCATGTGAATGCCGCCCAACTGCAgacgaa > 2:310886/1‑238 (MQ=255) gCCCAAATGGCGGGCAATTTGCCCACGGGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAgatga < 1:506604/238‑1 (MQ=255) aaggggggggcgTTTTTCCCCCGCGCGGCATTCCCCA‑GCACATAAAAATGTTCGGAACCCACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACg < 1:701423/232‑1 (MQ=255) atcccACA‑GACCATAATAATGTGCGGTACCAACCGCCTCTTCACCCTTTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt < 1:875537/236‑1 (MQ=255) tccccaggACCAAAAAAATGGTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGAAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGAt < 1:880770/231‑1 (MQ=255) caca‑GACCATAAATATGCCTGTTATCAACCGCCTCCTCACCCTGGTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt < 1:493843/239‑1 (MQ=255) caca‑GACCAAAAAAATTTTCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt < 1:331006/239‑1 (MQ=255) tacggACCATAAAAATGGCCGGTACAAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCg < 1:288135/233‑1 (MQ=255) ccaagcccAAAATAATGTTCGGTAACAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGgt < 1:684450/233‑1 (MQ=255) aCCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGAACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACAAGCGGCAGATGATGACGTTGTATCCCGCCCAGTGCGGCGGGCAGAACGGGGTAcc > 2:732303/1‑239 (MQ=255) | CCATCAGCTGCCAGGGCACTGATATGGTGGTAACTACCGTACATTGCCGGGCGCATCAGATCGTTGAACCCGGCATCAACCAGCACAAAGTGGCGGCTCCCCATTTGTTTGACGCTCCGCACC‑TGAGTAATTAATACGCCAGACTGCGCTACCAGGAAGCGACCCGGTTCAATTTCCAGTTTCACAGGGTGGCCCAAATGGCGGGCGATTTGCTCACGCGCGGCATTCCACA‑GACCATAATAATGTTCGGTATCAACCGCCTCTTCACCCTGTTGATAAGGAACAGAAAGCCCACCGCCCGCAGAAATAGCCTGTAAATCCTGACCGAATTCGATGACCTGACGCACCATAGCACCACACACCTGTTCCAGATGGGCATAATCAACGCCAGAACCAATGTGCATGTGAATGCCGACCAGCTGCAGATGATGACGTTGTATCACGTCCAGTGCGGCGGGCAGATCGGTGTACC > NC_000913/2977927‑2978398 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |