| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 2,666,921 | G→A | 7.5% | L29L (CTG→CTA) | prpF → | aconitate isomerase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 2,666,921 | 0 | G | A | 7.5% | 140.0 / 2.5 | 53 | L29L (CTG→CTA) | prpF | aconitate isomerase |
| Reads supporting (aligned to +/- strand): ref base G (20/29); new base A (2/2); total (22/31) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01 | |||||||||||
ATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGATCCTGTCGCAAAGCATCAA > NC_002947/2666772‑2667070 | aTGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTg > 1:990062/1‑150 (MQ=255) aTCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCg > 2:1001069/1‑150 (MQ=255) tCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGa < 1:515925/150‑1 (MQ=255) tCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGa < 1:1067917/150‑1 (MQ=255) ggggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTACCCGAGCa > 1:576559/1‑150 (MQ=255) ggCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAgg < 2:1118360/150‑1 (MQ=255) aaaGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCaaa < 2:990062/150‑1 (MQ=255) gCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATACGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGCGCAGGCGCAAATTc > 2:24656/1‑150 (MQ=255) ggCCCCAGCGTATTTCACCGGACTATCTGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCa < 1:1100429/150‑1 (MQ=255) gCCCCAGGGTATTTCAACGGGCTTACAAGGTGCCGATGAAAGTGTCACATGTACTCCAAATCATAATTCCCGCCAGCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGAAGGCGCAAATTCCAg < 2:403471/150‑1 (MQ=255) ccccAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCTCCTGCAAGACCTGCCCTAGCAGGCGCAAATTCCAgg < 2:13585/150‑1 (MQ=255) cTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCgctgc < 1:290383/150‑1 (MQ=255) aTCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAg < 1:748553/114‑1 (MQ=255) aTCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAg > 2:748553/1‑114 (MQ=255) tCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCgctgctgc > 2:1092120/1‑150 (MQ=255) gCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGt > 2:1123390/1‑150 (MQ=255) gCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGt > 2:102096/1‑150 (MQ=255) ccAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAg > 2:386034/1‑104 (MQ=255) ccAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAg < 1:386034/104‑1 (MQ=255) ccAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTaa > 1:591899/1‑150 (MQ=255) ccAAGAAAATGTCACATGTACCCCACATCAAAATCCCCGCCACCTACATCCGCGGCCGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGACCAGGCGCCAATTCCAGGTCCCGCCCGCGATGCGCCGCTGCTGCGGGTaa > 2:99080/1‑150 (MQ=255) cAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAgg > 1:1255471/1‑117 (MQ=255) cAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAgg < 2:1255471/117‑1 (MQ=255) aaGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATc > 2:980910/1‑150 (MQ=255) aaGAAAATGTCACATGTACCCCAGAGCAAAATCCCCACAACCTACATCCGCGGCGTCACCGGCAAAGGCGTGTTCTTCCGCCTCCAAGACCTGCCCGAGCAGGCGCACATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATc > 2:538745/1‑150 (MQ=255) aGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCg < 1:490221/150‑1 (MQ=255) gTCACATGTACCCCAGATCAAAATCCCCCCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCCGGCCTGTGCCCCTTTGTTCGGCCGCTATGCGCTGCTGCTGCGGGTAATCGGCAGccc < 1:24656/150‑1 (MQ=255) acaTGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGa < 1:1123390/150‑1 (MQ=255) tGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTACCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGAccc > 1:718502/1‑150 (MQ=255) tGTACCACAGATCAAAAACCCCGCCACCTACATCCGCCGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGAccc < 2:160794/150‑1 (MQ=255) gTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATc > 1:83957/1‑135 (MQ=255) gTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATc < 2:83957/135‑1 (MQ=255) aGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGgc > 1:1013790/1‑82 (MQ=255) aGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGgc < 2:1013790/82‑1 (MQ=255) aTCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATc < 1:1001069/150‑1 (MQ=255) ccGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACg < 1:1092120/150‑1 (MQ=255) aCCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATg < 2:9457/150‑1 (MQ=255) ccTACATCCGCGGCGGCACCAGCAAAGGCGTGTTGTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATgg > 1:103191/1‑150 (MQ=255) ctCCGCGGCGGCACAAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCACGTTCTGGCGCAAATTCCAGGTCCCCCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGggcggc < 2:359886/149‑1 (MQ=255) ccGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCa < 2:382709/107‑1 (MQ=255) ccGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCa > 1:382709/1‑107 (MQ=255) gcggcgGCACCAGCAAAGGCGTGTTCTTCCGCCTGCTACACGTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCa < 2:103191/150‑1 (MQ=255) gcggcgGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGTTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCa > 2:1063865/1‑150 (MQ=255) cggcACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTACCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTc < 2:718502/150‑1 (MQ=255) aCCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGc > 1:1116965/1‑150 (MQ=255) aGCAAAGGCGTGTTCTACAGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCAcc < 1:1028997/150‑1 (MQ=255) gCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCa < 1:980910/150‑1 (MQ=255) gCAAAGGCGTGTTCTTCCGCCTGCAAGACCTACCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCa < 2:576559/150‑1 (MQ=255) gCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGa < 1:102096/150‑1 (MQ=255) tcttcCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGTTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGa < 1:1063865/150‑1 (MQ=255) tcttcCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGa < 2:591899/150‑1 (MQ=255) gCCTGCAAGACCTGCCCGAGCAGGCGCAAATGCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGATCCTGt > 2:789330/1‑150 (MQ=255) aGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGATCCTGTCGCAAAg > 1:463619/1‑150 (MQ=255) gACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGATCCTGTCGCAAAGc < 2:463619/150‑1 (MQ=255) gCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGATCCTGTCGCAAAGCATCaa < 1:789330/150‑1 (MQ=255) | ATGAATCGGGGGCCGCAAAGCGGCCCCAGGGTATTTCACCGGACTATCAGGATGCCCAAGAAAATGTCACATGTACCCCAGATCAAAATCCCCGCCACCTACATCCGCGGCGGCACCAGCAAAGGCGTGTTCTTCCGCCTGCAAGACCTGCCCGAGCAGGCGCAAATTCCAGGTCCCGCCCGCGATGCGCTGCTGCTGCGGGTAATCGGCAGCCCCGACCCTTACGGCAAGCAGATCGACGGCATGGGCGGCGCCACGTCGAGCACCAGCAAGACCGTGATCCTGTCGCAAAGCATCAA > NC_002947/2666772‑2667070 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |