Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,153 | A→G | 17.7% | intergenic (+59/‑88) | 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,153 | 0 | A | G | 17.7% | 79.7 / 10.9 | 42 | intergenic (+59/‑88) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base A (17/16); new base G (5/2); total (22/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.27e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.26e-01 |
ATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC > NC_002947/5485011‑5485297 | aTGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtgg < 1:341626/149‑1 (MQ=255) tGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgccg > 8:209413/1‑149 (MQ=255) gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCCAATACGATGGTGGATTCACcgccgcct > 8:164206/1‑146 (MQ=255) cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc < 2:179155/147‑1 (MQ=37) cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg < 3:83286/128‑1 (MQ=255) cGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtg > 4:83286/1‑128 (MQ=255) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGTCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg < 6:14046/148‑1 (MQ=255) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcgt < 6:219773/148‑2 (MQ=255) tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGGATTCCAAGGTGGATTCACCGCCGCCTTAGCCGGCaacc > 7:178683/1‑147 (MQ=255) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCATGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCGCAAACGAAGGTGGCTTGACCCCCGCCTTAGCGGGGAAGccc > 5:107900/1‑149 (MQ=255) ccGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGc < 7:209413/147‑1 (MQ=255) acgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGAATAACATTGTGGATTAACCGCCGCAATAGCGGGCAAGCCCGca > 6:100055/1‑148 (MQ=255) ggCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTGCCGGGAATACGGTGGTGGTTTGACCGGCGCCTTAGCGGGCCAGCCCGcggcca > 7:281984/1‑144 (MQ=255) ggCATCGCGTTGGTCTTCACTGGAATGCGCCACTTCCGCCAATAATTACGGGGATCCCGTGGCTTGTCGATATCGTGTGGGAGGGGGCTTGCCCGCTAATACTTTGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAgt < 6:107900/149‑1 (MQ=37) ggCATCGCGATGGGCTTCACTGGCAGGCGCCACTTCCGCAACTAATTACGCGGATCCCGAGGCAAGTCGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAgt < 8:260343/149‑1 (MQ=255) aTCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGtt < 3:239251/149‑1 (MQ=255) cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc < 2:88524/125‑1 (MQ=255) cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCAAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc > 1:88524/1‑125 (MQ=255) tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgt > 5:64094/1‑97 (MQ=21) tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatcgt < 6:64094/102‑6 (MQ=21) ccGCCACTAATTCCGCGGAGCCCGCGGCCAGTGGAGCTCCGGGGGGCGCGGGCTTGCCCGCGAATACGATGGTGTCTTCACCGCCGCATTCGCGGGCAAGCCCGCGCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg < 5:83945/148‑1 (MQ=255) ccGCCACAAATTACGCGGACCCCGAGCAATGGCGAGCCCCGGGGGGCGCGGGCTTGCCCGCTAATACGATGGTGGATTCCCCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg < 5:316131/148‑1 (MQ=37) gCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt > 4:112906/1‑75 (MQ=255) gcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCGAAAAACAAGGTTGATTCACCGCCCCAATAGCGGGCGAACCCGCGACCCCAATGTTACCCGCAAGGCACCGGATAACGTTAATCGAGGTTTTTACATAAAAg > 6:71563/1‑149 (MQ=255) gAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCCAATAAGATGGTGGATTAACCCCCCCCATAGCGGGGCACCCCCCGCCCACACGGTTAGGCGCAACCCCCCGGGTCCCATGAATAGAGGGTTTTACAAGAAAATTTTTatg > 6:292055/1‑147 (MQ=255) acTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg < 7:159265/57‑1 (MQ=255) acTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg > 8:159265/2‑58 (MQ=255) tgGGAGCGGGCTTGCCCGGTAATACGTTGGTGGATTCACGGCGGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGCCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 5:71563/149‑1 (MQ=39) tgGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 5:100055/149‑1 (MQ=255) gggCTTGCCCGCGAATACGATGGTGGATTCACCGCTGAATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCACGCCACGGGATAACTTGAATCGAGGTTATGCCATGAAAGTTTTGATCATCGGCCGCGGTGGCCGTGAGCACgac > 2:99448/1‑147 (MQ=38) ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg < 5:343111/95‑1 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg > 6:343111/1‑95 (MQ=255) ccGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGCCAGCGCTGGCCGTGAGCACGCcctggcctg > 5:4005/1‑149 (MQ=25) ccGCGAATGCGGCGGTGAATCCACCATCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 2:244143/1‑149 (MQ=37) cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgcgc > 1:80706/1‑65 (MQ=255) cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGgcgc < 2:80706/65‑1 (MQ=255) cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 6:281792/1‑148 (MQ=255) gcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg > 6:223233/1‑148 (MQ=37) gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGaaa < 5:292055/149‑1 (MQ=255) tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCa < 1:208832/130‑1 (MQ=255) tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCa > 2:208832/1‑130 (MQ=255) tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc > 5:236529/1‑149 (MQ=255) | ATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC > NC_002947/5485011‑5485297 |
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: ‑ |