| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,485,162 | T→C | 20.5% | intergenic (+68/‑79) | purH → / → purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,485,162 | 0 | T | C | 20.5% | 54.0 / 12.5 | 36 | intergenic (+68/‑79) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
| Reads supporting (aligned to +/- strand): ref base T (12/15); new base C (3/4); total (16/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-02 | |||||||||||
GTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGG > NC_002947/5485015‑5485304 | gTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCATACGATGGGGGGtt > 3:168250/1‑148 (MQ=255) tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAAACGATGGTGGCtta > 6:255019/1‑147 (MQ=255) aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgta < 8:78533/149‑7 (MQ=37) aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCAccatcg < 4:165926/147‑5 (MQ=38) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAAAAAGATGGTGGATTAACCG‑CCGCATTagcg > 4:83977/1‑149 (MQ=255) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCA‑TCGTATTcgcg < 2:133067/149‑1 (MQ=37) gaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATACACCG‑CCGCATTCGCGGGCaa < 8:59113/140‑1 (MQ=255) gaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATACACCG‑CCGCATTCGCGGGCaa > 7:59113/1‑140 (MQ=255) gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCATTCCCCCCAAAAAGATGGTGGATTAAACG‑CCCCATTAGAGGGCAAccccccaaccc > 8:48811/1‑141 (MQ=255) cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTTGCCGCGACTACGATGGTGGGTTCAACG‑CCGCAGTAGCGGGGCAGGCCGCGCCCACAg > 8:123107/1‑149 (MQ=255) taTGGTCTTCACGGGCCTGCCCCACTCCCCCCACAAATTCCCCGGACCCCGCGCCACGCCGAGACCCGGGGGGCGCGGGCTGGCCCGCTAATACGATGGTGGATTCCCCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgcg < 4:168250/148‑1 (MQ=255) ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAATACGATGGTGGATTAACCG‑CCGCCTTAGCGGGCAAGCCCGCTGCCACAGTGTTAGGCgccagcca > 7:249422/1‑149 (MQ=255) aCTGGCATGCGCCACTTCCCCCACTAATTCCGCGGATCCGGAGCCAAGTCTAGATCGTGGGGGCGCGGGCTGGCCCGCTAATACTTTGGTGGATTCCCCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg < 8:137824/149‑1 (MQ=255) gcgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCgg > 2:50879/1‑104 (MQ=255) gcgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCgg < 1:50879/104‑1 (MQ=255) cctcTAATTACGCGGATCCCGAGGCAAGTCGACATCCTGTGGGAGCGGGCTTGCCCGCTAATAGGATGGTGGATTCCCCG‑CCGCATTCGCGGGCAATCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGttt < 5:255019/146‑1 (MQ=39) cACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGtttt < 3:178322/149‑1 (MQ=255) aCGCGGATCCCGAGGCAAGTCGTGATCCTGTGGTAGCGGGCTTGCCCGCTAATACGATGGTGTATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGa < 7:48811/148‑1 (MQ=255) gtgGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt < 5:131275/99‑1 (MQ=255) gtgGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTAGCGGGCAAGCCCGCTCCCCCAGTGTTAGGCGCAAGCCCCCGGATCAcgt > 6:131275/1‑97 (MQ=25) ggggAGGGGGGTTGCCGGCTAATACTTTGG‑GGCTTCCCCGCCCGCTTGCGGGGCAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 5:212933/148‑1 (MQ=255) gggAGCGGGATTGCCCGCTAATACTATGGTGGATTCACGG‑CCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 1:162330/148‑1 (MQ=255) gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat > 2:103424/1‑126 (MQ=255) gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat < 1:103424/126‑1 (MQ=255) cGGGCTTGCCCGCTAATACGATGGTGTATTCACCG‑CCGCTTTCGCGGGCATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGc < 1:88832/149‑1 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg > 8:152136/1‑149 (MQ=255) ccGCGAATGCGGCGGTGAATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 4:26135/1‑149 (MQ=37) ccGCGAATGCGGCGGTGAATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:211659/1‑149 (MQ=37) ccGCGAATGCGGCGGTGAATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:241923/1‑149 (MQ=37) ccGCGAATACGATGGTGGATTCACTG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATTAAAGTTTTGATCATAGGCAGCGGTGGCCGTGAGCCCGCcctggccgg > 7:37673/1‑147 (MQ=255) aTACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa < 3:83977/147‑1 (MQ=255) gATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCGGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc > 2:48229/1‑149 (MQ=255) gATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa < 4:246346/141‑1 (MQ=255) gATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa > 3:246346/1‑141 (MQ=255) ggtggATTCCCCG‑CCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa < 3:99754/148‑1 (MQ=255) gtgaATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAgg < 5:241923/149‑1 (MQ=39) aTTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg > 2:153054/1‑144 (MQ=255) aTTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg < 1:153054/144‑1 (MQ=255) | GTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGG > NC_002947/5485015‑5485304 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |