Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_000913 | 456,490 | T→C | 10.1% | intergenic (+59/‑187) | tig → / → clpP | peptidyl‑prolyl cis/trans isomerase (trigger factor)/proteolytic subunit of ClpA‑ClpP and ClpX‑ClpP ATP‑dependent serine proteases |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_000913 | 456,490 | 0 | T | C | 10.1% | 182.7 / 9.3 | 59 | intergenic (+59/‑187) | tig/clpP | peptidyl‑prolyl cis/trans isomerase (trigger factor)/proteolytic subunit of ClpA‑ClpP and ClpX‑ClpP ATP‑dependent serine proteases |
Reads supporting (aligned to +/- strand): ref base T (26/27); new base C (3/3); total (29/30) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01 |
AAAGTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACTAGCTGATAATCCGTCCATA > NC_000913/456378‑456626 | aaagTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGa > 3:560017/1‑139 (MQ=255) aaagTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGa < 4:560017/139‑1 (MQ=255) gTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGt < 3:1133877/119‑1 (MQ=255) gTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGt > 4:1133877/1‑119 (MQ=255) aaaaGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑ttttt < 2:235673/106‑1 (MQ=255) aaaaGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑ttttt > 1:235673/1‑106 (MQ=255) ttcagtcACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg < 1:476534/142‑1 (MQ=255) ataaaCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTg > 2:359621/3‑107 (MQ=255) aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTg < 1:359621/107‑1 (MQ=255) aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGc < 3:52739/120‑1 (MQ=255) aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGc > 4:52739/1‑120 (MQ=255) aGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg > 2:759506/1‑148 (MQ=255) aaaCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc < 4:462748/144‑1 (MQ=255) aaaCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc > 3:462748/1‑144 (MQ=255) aCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGa < 4:507079/109‑1 (MQ=255) aCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGa > 3:507079/1‑109 (MQ=255) aCCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGtt < 1:759506/148‑1 (MQ=255) gatgaACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTaaca < 1:554150/140‑1 (MQ=255) gatgaACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTaaca > 2:554150/1‑140 (MQ=255) aCCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg < 1:633182/122‑1 (MQ=255) aCCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg > 2:633182/1‑122 (MQ=255) ccAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑ttttt > 1:1106693/1‑76 (MQ=255) ccAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑ttttt < 2:1106693/76‑1 (MQ=255) cagcagGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg < 4:1128964/120‑1 (MQ=255) cagcagGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCgg > 3:1128964/1‑120 (MQ=255) tAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGgtg < 4:636979/113‑1 (MQ=255) tAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGgtg > 3:636979/1‑113 (MQ=255) tAATTTACGCAGCATAACGCGCTAAATCCGCTCAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTc < 3:740738/72‑1 (MQ=255) tAATTTACGCAGCATAACGCGCTAAATCCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTc > 4:740738/1‑72 (MQ=255) aTTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGc < 3:362777/81‑1 (MQ=255) aTTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGc > 4:362777/1‑81 (MQ=255) cagcaTAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGAt > 1:480397/1‑148 (MQ=255) cagcaTAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTc > 2:974921/1‑63 (MQ=255) cagcaTAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTc < 1:974921/63‑1 (MQ=255) cgcgCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAg > 4:426176/1‑84 (MQ=255) cgcgCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAg < 3:426176/84‑1 (MQ=255) aaaTTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGaa > 4:350470/1‑123 (MQ=255) aaaTTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGaa < 3:350470/123‑1 (MQ=255) tCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCAGTGTTAGCGTAACAACAAAAGGTTGTTATGCTTGAAata < 1:328625/122‑1 (MQ=255) tCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATCTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAata > 2:328625/1‑122 (MQ=255) cacaAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAg > 2:581943/1‑144 (MQ=255) cacaAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAg < 1:581943/144‑1 (MQ=255) caAAGGCCCGTCACCGCCAGGTGGTGGGC‑Ttttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttc < 3:924722‑M2/108‑79 (MQ=255) caAAGGCCCGTCACCGCCAGGTGGTGGGC‑Ttttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttc > 4:924722‑M2/1‑30 (MQ=255) caAAGGCCCGTCACCGCCAGGTGGTGGGC‑Ttttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggtt > 3:769524‑M2/1‑30 (MQ=255) caAAGGCCCGTCACCGCCAGGTGGTGGGC‑Ttttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggtt < 4:769524‑M2/78‑49 (MQ=255) gCCCGTCACCGCCAGGTGGTGGGCTTTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACTAGCTg < 2:480397/148‑1 (MQ=255) cccGTCACCGCCAGGTGGTGGGC‑Ctttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgcc > 2:82436‑M2/1‑24 (MQ=255) cccGTCACCGCCAGGTGGTGGGC‑Ctttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgcc < 1:82436‑M2/121‑98 (MQ=255) ccGTCACCGCCAGGTGGTGGGC‑Ctttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgccagtacagcttcaacagcctgttcttcca < 2:218670‑M2/148‑126 (MQ=255) cGTCACCGCCAGGTGGTGGGC‑Ctttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttct < 2:5113‑M2/95‑74 (MQ=255) cGTCACCGCCAGGTGGTGGGC‑Ctttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttct > 1:5113‑M2/1‑22 (MQ=255) ccGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc < 2:879819/64‑1 (MQ=255) ccGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTc > 1:879819/1‑64 (MQ=255) cGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt < 2:667297/133‑1 (MQ=255) cGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt > 3:56768/1‑133 (MQ=255) cGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt > 1:667297/1‑133 (MQ=255) cGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACt < 4:56768/133‑1 (MQ=255) aGGTGGTGGGC‑Ctttgtgcgaatttagcgcgttatgctgcgtaaattacgcctgctggttcatcagctcgttgaaagtggtttctttttcagtcactttcgctttcgccagtacagcttcaacagcctgttcttccagagcaacattg > 1:218670‑M2/1‑12 (MQ=255) ggtggtGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAggg > 1:186930/1‑125 (MQ=255) ggtggtGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAggg < 2:186930/125‑1 (MQ=255) c‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCg > 1:186983/1‑59 (MQ=255) c‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCg < 2:186983/59‑1 (MQ=255) | AAAGTGACTGAAAAAGAAACCACTTTCAACGAGCTGATGAACCAGCAGGCGTAATTTACGCAGCATAACGCGCTAAATTCGCACAAAGGCCCGTCACCGCCAGGTGGTGGGC‑TTTTTTTTGTCATGAATTTTGCATGGAACCGTGCGAAAAGCCTCTTTCGGTGTTAGCGTAACAACAAAAGATTGTTATGCTTGAAATATGGTGATGCCGTACCCATAACACAGGGACTAGCTGATAATCCGTCCATA > NC_000913/456378‑456626 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |