Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,150 | A→G | 20.0% | 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 | 20.0% | 59.7 / 10.2 | 30 | intergenic (+56/‑91) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base A (9/15); new base G (3/3); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.60e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01 |
GTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAA > NC_002947/5485003‑5485294 | gtggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATAc < 6:93272/149‑1 (MQ=255) gtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAt < 6:179956/149‑1 (MQ=255) ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtg > 5:258263/1‑149 (MQ=255) gATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg < 5:249098/149‑1 (MQ=255) gATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg < 2:403156/149‑1 (MQ=255) tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa < 5:35740/149‑1 (MQ=255) tGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGa < 7:140823/149‑1 (MQ=255) gATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCAc < 4:200484/149‑1 (MQ=255) gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcg > 4:372847/1‑145 (MQ=17) aaGTCATGGGTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTAACTGGCATGCGCCACTTCCGCCAATAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGGCGGTGAATCCAccatc < 2:201717/146‑4 (MQ=37) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc < 3:372847/149‑1 (MQ=14) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTTGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg < 6:285029/148‑1 (MQ=25) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg < 3:10703/148‑1 (MQ=37) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg < 2:242928/148‑1 (MQ=255) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg < 6:258263/148‑1 (MQ=255) ctgctgCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGAGC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGc < 1:8377/149‑1 (MQ=255) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtg < 2:19036/117‑1 (MQ=255) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtg > 1:19036/1‑117 (MQ=255) gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAcc < 4:175726/120‑1 (MQ=255) gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAcc > 3:175726/1‑120 (MQ=255) cATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATa > 2:302943/1‑78 (MQ=255) cATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATa < 1:302943/78‑1 (MQ=255) ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGccc > 4:80301/1‑149 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt < 6:111911/113‑1 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt > 5:111911/1‑113 (MQ=255) ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:409081/1‑149 (MQ=37) ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGGTTTGATCATCGGCAGCGGTAGCCGTGAGCACGCcctggcctg > 4:76769/1‑149 (MQ=25) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCTCTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 4:232609/1‑149 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:92006/1‑149 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 8:370374/1‑149 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATAGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 7:255683/1‑149 (MQ=255) gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa < 7:370374/149‑1 (MQ=255) | GTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAA > NC_002947/5485003‑5485294 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |