Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,159 | G→A | 18.5% | intergenic (+65/‑82) | 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,159 | 0 | G | A | 18.5% | 59.6 / 8.9 | 34 | intergenic (+65/‑82) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base G (14/12); new base A (3/3); total (17/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 = 2.49e-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 < 1:50879/104‑1 (MQ=255) gcgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCgg > 2:50879/1‑104 (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 < 1:103424/126‑1 (MQ=255) gCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcat > 2:103424/1‑126 (MQ=255) cGGGCTTGCCCGCTAATACGATGGTGTATTCACCG‑CCGCTTTCGCGGGCATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGc < 1:88832/149‑1 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg > 8:152136/1‑149 (MQ=255) ccGCGAATGCGGCGGTGAATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:211659/1‑149 (MQ=37) ccGCGAATGCGGCGGTGAATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:241923/1‑149 (MQ=37) ccGCGAATGCGGCGGTGAATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 4:26135/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 > 3:246346/1‑141 (MQ=255) gATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa < 4:246346/141‑1 (MQ=255) ggtggATTCCCCG‑CCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCa < 3:99754/148‑1 (MQ=255) gtgaATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAgg < 5:241923/149‑1 (MQ=39) | 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: ‑ |