| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,651,319 | G→C | 12.1% | intergenic (+41/+29) | PP_2323 → / ← aroF‑I | hypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,651,319 | 0 | G | C | 12.1% | 95.5 / 6.0 | 33 | intergenic (+41/+29) | PP_2323/aroF‑I | hypothetical protein/tyr‑sensitive phospho‑2‑dehydro‑3‑deoxyheptonate aldolase |
| Reads supporting (aligned to +/- strand): ref base G (14/15); new base C (2/2); total (16/17) | |||||||||||
| 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.75e-01 | |||||||||||
TCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATA > NC_002947/2651179‑2651458 | tctGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCg > 1:574517/1‑150 (MQ=255) gAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGtttt > 2:1500/1‑150 (MQ=255) ccGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt > 2:557172/1‑142 (MQ=255) ccGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt < 1:557172/142‑1 (MQ=255) gAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCg > 1:785849/1‑150 (MQ=255) cAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCaa < 1:263023/150‑1 (MQ=255) gTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGc > 2:127785/1‑150 (MQ=255) tACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGcc < 2:200328/150‑1 (MQ=255) ccGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTg < 2:613384/150‑1 (MQ=255) gCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttt < 1:335275/109‑1 (MQ=255) gCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGttt > 2:335275/1‑109 (MQ=255) cTCCTATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGTGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGgcagca < 1:195530/150‑1 (MQ=255) tCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt > 1:432918/1‑107 (MQ=255) tCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCACTGCCCGGCGttt < 2:432918/107‑1 (MQ=255) atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAg < 2:334719/125‑1 (MQ=255) atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAg > 1:334719/1‑125 (MQ=255) atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa < 1:745023/134‑1 (MQ=255) atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTa > 2:745023/1‑134 (MQ=255) atatGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAACGCCTGAGCTGACGACAGGAACCAGACACCACAAAAAACGCCGGGCAGTGAACGGCTTTTTTCATTCCCGAAATGTTCAGCCGCGCTTACGCTGCGGCAGCACAt > 1:244528/1‑150 (MQ=255) gtgGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCtt < 2:785849/150‑1 (MQ=255) gAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAg < 1:1500/150‑1 (MQ=255) cGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGAAAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCtt < 1:175419/150‑1 (MQ=255) cATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcat < 2:103205/150‑1 (MQ=255) aTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatg > 1:364404/1‑150 (MQ=255) aTGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGgcatg < 1:127785/150‑1 (MQ=255) gcagcGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCa < 1:98799/150‑1 (MQ=255) gCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGa > 1:626926/1‑150 (MQ=255) gacAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCa > 2:710519/1‑150 (MQ=255) aacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACCTACTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGac > 1:254476/1‑150 (MQ=255) aacaaAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGac > 2:158853/1‑150 (MQ=255) aaaaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGc < 2:626926/150‑1 (MQ=255) aaaCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGTATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCa > 2:725544/1‑150 (MQ=255) aCGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATa < 1:710519/150‑1 (MQ=255) | TCTGGAGTACGCCGAGCAGATGGGGTACACGGTGATCCCGTCTTGCTCATATGTGGAGCGTTACATGGAGCGCCAGCAGCGGCATTCGAGCAAGGCCTGAGCTGACGACAGGCACCAGACACAACAAAAAACGCCGGGCAGTGCCCGGCGTTTTTCATTCCCGCAATGTTCAGCCGCGCTTACGCTGCGGCAGCACATCCTTGAGCTTGGCATGCATGCTGCGCAGGGTTTTCTCGGTAGCCGACCAGTCGATGCAGGCGTCCGTGATCGACACGCCATA > NC_002947/2651179‑2651458 |
| 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: ‑ |