| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,485,383 | C→G | 14.1% | P32R (CCG→CGG) | purD → | phosphoribosylamine‑‑glycine ligase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,485,383 | 0 | C | G | 14.1% | 81.8 / 6.4 | 34 | P32R (CCG→CGG) | purD | phosphoribosylamine‑‑glycine ligase |
| Reads supporting (aligned to +/- strand): ref base C (13/16); new base G (2/3); total (15/19) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01 | |||||||||||
GCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGAAGCACCGCTGGTTATCGGTGTG > NC_002947/5485236‑5485525 | gCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCgg < 1:58685/150‑1 (MQ=255) cTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTTTAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCggg < 1:81596/150‑1 (MQ=255) aGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCGGGGCAACgaa < 1:45929/150‑3 (MQ=255) gtgtTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCGGGGCAACgaag < 1:199444/150‑4 (MQ=255) ttCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTa > 1:156578/1‑150 (MQ=255) ttCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTGTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTa > 1:428798/1‑150 (MQ=255) gcgcAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAGCGCCGGTACCGc > 2:213500/1‑150 (MQ=255) aaGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCAtt < 2:200556/150‑1 (MQ=255) aCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGC > 2:183882‑M1/1‑150 (MQ=255) gggATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCA > 2:393003‑M1/1‑150 (MQ=255) ggATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGACGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAA < 2:273068‑M1/150‑1 (MQ=255) aCTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCAGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCg < 1:38530/150‑1 (MQ=255) tCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCaa > 2:111835/1‑114 (MQ=255) tCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCaa < 1:111835/114‑1 (MQ=255) aTGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCg > 2:239726/1‑139 (MQ=255) aTGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCg < 1:239726/139‑1 (MQ=255) aaaGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTgc > 2:128680/1‑150 (MQ=255) cAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAAttggc > 2:152785‑M2/1‑145 (MQ=255) aGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAttgctccagggcgcacacgtcgatggcgacgttctcgcacttggct < 1:183882‑M1/150‑47 (MQ=255) aGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGcc < 2:156578/150‑1 (MQ=255) gTGAGCACCCCCTGGCATGGAAAGTCGCGCAGGACCGACGCGTGGAGAAAGTTGTCGTTTTCCGGGGTAACGCCGGTAAGGCCATTGAAGCCAAGTACGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCg > 2:117929/1‑150 (MQ=255) tGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCga < 1:47852/150‑1 (MQ=255) tggAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACTTCGACCTg > 1:332250/1‑150 (MQ=255) ccAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCgg > 1:293678/1‑150 (MQ=255) cAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGc > 1:367632/1‑150 (MQ=255) cccACGCGTCGAGAAAGGCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGa > 2:143405/1‑150 (MQ=255) ccACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGaa < 1:128680/150‑1 (MQ=255) cACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGAAg < 2:293678/150‑1 (MQ=255) gcgcgGCGAGAAAGTCTTCGTTGCGCCGGGCAACGCCGGTAGCGCCATTGAAGCCAAGTGCGAGAACGTTGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTGCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGAAGc < 2:137653/149‑1 (MQ=255) tCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGAAGCACCGCTGGTTATCg < 1:393003/150‑1 (MQ=255) ttCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGAAGCACCGCTGGTTATCGgt < 2:367632/150‑1 (MQ=255) ggcGTTGCCCGGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGcc > 1:411926/3‑96 (MQ=255) ggcGTTGCCCGGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGcc < 2:411926/94‑1 (MQ=255) gTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGCAAAACGTCGACCTGACCATCGTCGGCCCGGAAGCACCGCTGGTTATCGGtgtg < 1:213500/150‑1 (MQ=255) ccccGGGCAACGCCGGTACCGCCATTGAAGCCAAttgctccagggcgcacacgtcgatggcgacgttctcgcacttggcttcaa < 1:282584‑M1/84‑51 (MQ=255) ccccGGGCAACGCCGGTACCGCCATTGAAGCCAAttgctccagggcgcacacgtcgatggcgacgttctcgcacttggcttcaa > 2:282584‑M1/1‑34 (MQ=255) | GCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGTCGAGAAAGTCTTCGTTGCCCCGGGCAACGCCGGTACCGCCATTGAAGCCAAGTGCGAGAACGTCGCCATCGACGTGTGCGCCCTGGAGCAATTGGCCGACTTCGCCGAGAAAAACGTCGACCTGACCATCGTCGGCCCGGAAGCACCGCTGGTTATCGGTGTG > NC_002947/5485236‑5485525 |
| 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: ‑ |