Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,150 | A→G | 13.5% | intergenic (+56/‑91) | 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,150 | 0 | A | G | 13.5% | 82.8 / 6.4 | 37 | intergenic (+56/‑91) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base A (14/18); new base G (3/2); total (17/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.44e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01 |
GGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAA > NC_002947/5485002‑5485294 | ggtggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATa < 6:244524/149‑1 (MQ=255) gtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACggc > 3:165926/1‑147 (MQ=255) cGATCCGGGATCCGGAAGCCATCGCCGCTGCCGACGAGCCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt < 4:21493/149‑1 (MQ=255) gATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg < 7:95689/149‑1 (MQ=255) gATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt < 2:165038/148‑1 (MQ=255) gTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCATACGATGGGGGGtt > 3:168250/1‑148 (MQ=255) tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAAAACGATGGTGGCtta > 6:255019/1‑147 (MQ=255) gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt > 7:159487/1‑138 (MQ=255) gCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt < 8:159487/138‑1 (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 > 7:59113/1‑140 (MQ=255) gaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATACACCG‑CCGCATTCGCGGGCaa < 8:59113/140‑1 (MQ=255) gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCATTCCCCCCAAAAAGATGGTGGATTAAACG‑CCCCATTAGAGGGCAAccccccaaccc > 8:48811/1‑141 (MQ=255) ccGGCATCGCGATGGTCTTCACTTGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATa < 7:13191/103‑1 (MQ=255) ccGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATa > 8:13191/1‑103 (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 < 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:241923/1‑149 (MQ=37) ccGCGAATGCGGCGGTGAATCCACCA‑TCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:211659/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) | GGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAA > NC_002947/5485002‑5485294 |
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: ‑ |