Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,168 | C→T | 14.7% | intergenic (+74/‑73) | 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,168 | 0 | C | T | 14.7% | 74.8 / 5.9 | 34 | intergenic (+74/‑73) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base C (13/16); new base T (3/2); total (16/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.48e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-02 |
TCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGT > NC_002947/5485030‑5485314 | 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 > 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 < 1:153054/144‑1 (MQ=255) aTTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAg > 2:153054/1‑144 (MQ=255) cg‑ccgcATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCACGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTGTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGCGCCAcgcgg > 7:216619/1‑148 (MQ=255) | TCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCG‑CCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGACCCACGCGT > NC_002947/5485030‑5485314 |
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: ‑ |