Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 126,116 | C→T | 9.4% | intergenic (+26/+24) | PP_0120 → / ← thrB | zinc ABC transporter/homoserine kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 126,116 | 0 | C | T | 9.4% | 141.4 / 5.6 | 53 | intergenic (+26/+24) | PP_0120/thrB | zinc ABC transporter/homoserine kinase |
Reads supporting (aligned to +/- strand): ref base C (27/21); new base T (3/2); total (30/23) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01 |
CTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGCGC > NC_002947/125980‑126260 | cTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGaagaag < 2:29251/150‑1 (MQ=255) cGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGAc > 2:109634/1‑150 (MQ=255) gCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGAcgcg < 1:512434/150‑1 (MQ=255) tGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCTGt > 1:360343/1‑150 (MQ=255) ttGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAgcg > 2:116060/1‑150 (MQ=255) gTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGaa > 1:584447/1‑150 (MQ=255) ggCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaa > 2:84415/1‑122 (MQ=255) ggCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCgaa < 1:84415/122‑1 (MQ=255) ccGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCAc > 2:44856/1‑109 (MQ=255) ccGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCAc < 1:44856/109‑1 (MQ=255) ggTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAAt < 2:402407/150‑1 (MQ=255) aaaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCAc < 2:155290/145‑1 (MQ=255) aaaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCAc > 1:155290/1‑145 (MQ=255) aaGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACctgctg < 2:584447/150‑1 (MQ=255) ggTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGgcg > 1:176125/1‑150 (MQ=255) aGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTgcgc > 1:712675/1‑150 (MQ=255) gcTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCACTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGt < 2:310058/102‑1 (MQ=255) gcTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCACTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGt > 1:310058/1‑102 (MQ=255) gggAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCAc < 1:152132/149‑1 (MQ=255) gagaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACtt > 2:757259/1‑150 (MQ=255) agaAGCTGGGCGGTGACCTGTCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTAAAGCGCGACCGGCAGGTCCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTc < 1:381166/150‑1 (MQ=255) gaAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCg < 2:176125/150‑1 (MQ=255) cGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGTCTCCTTCGCGGGCACGCCCGCGAAGAAGCAAACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAActcgctc < 1:53242/150‑1 (MQ=255) cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTTTACAGCGCGAACGGCAGGTGCAATGCCACCTGGTGGCGATGCGACAGACGCACTTCGAACTCGCTCGGATTGt > 1:310611/1‑150 (MQ=255) cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGt > 1:715308/1‑150 (MQ=255) cTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTCCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGt > 2:577751/1‑150 (MQ=255) tGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTg > 1:247854/1‑150 (MQ=255) ggCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTgg < 1:116060/150‑1 (MQ=255) cTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatcatc > 2:153753/1‑150 (MQ=255) cTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcat < 1:577751/150‑1 (MQ=255) tAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCCATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATACATCCCGGCGa > 1:166868/1‑150 (MQ=255) tAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcatcca < 2:380235/141‑1 (MQ=255) tAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCAcatcca > 1:380235/1‑141 (MQ=255) ccTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGa > 2:636130/1‑150 (MQ=255) ccTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGa < 2:247854/150‑1 (MQ=255) cTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc < 1:777056/108‑1 (MQ=255) cTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc > 2:777056/1‑108 (MQ=255) cTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGAc > 1:397445/1‑150 (MQ=255) cTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGAc > 2:704454/1‑150 (MQ=255) ccGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGc > 2:682716/1‑150 (MQ=255) ggCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc > 2:549747/1‑101 (MQ=255) ggCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc > 2:549649/1‑101 (MQ=255) ggCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc < 1:549649/101‑1 (MQ=255) ggCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAAc < 1:549747/101‑1 (MQ=255) gcttcTTCGCGGGCGTGCCCGCGAAGGGGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTg > 2:570067/5‑114 (MQ=38) gcttcTTCGCGGGCGTGCCCGCGAAGGGGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTg < 1:570067/110‑1 (MQ=38) ttcTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCaa > 2:176120/3‑150 (MQ=255) ttcTTCGCGGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCaa > 1:696944/3‑150 (MQ=255) cccTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTAGAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCaa > 1:701262/1‑150 (MQ=255) ttCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGa < 2:544630/104‑1 (MQ=255) ttCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGa > 1:544630/1‑104 (MQ=255) ttCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatc > 1:267250/1‑113 (MQ=255) ttCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGatc < 2:267250/113‑1 (MQ=255) cgcgGGCGTGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAgg < 1:176120/150‑1 (MQ=255) gggCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGgcgc < 1:704454/150‑1 (MQ=255) | CTTGCCGAGCTGGATGCCTTGGGTGGCACCGATCCGGTGGATGGCAAAGGGTATGAGCGCTTGCTGGAGAAGCTGGGCGGTGACCTGGCAGGCTGCCTGGAGCAGTTGTAAGCCTGTGCCGGCCCCTTCGCGGGCACGCCCGCGAAGAAGCAGACGCGGTCTACAGCGCGAACGGCAGGTGCAATGCCACCTGCTGGCGCTGCGCCAGACGCACTTCGAACTCGCTCGGATCGTGGATCATCACATCCATCCCGGCGAACGACTCGGCCGCAATCAGGCGC > NC_002947/125980‑126260 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |