Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,485,150 | A→G | 19.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 | 19.5% | 76.8 / 15.2 | 41 | intergenic (+56/‑91) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
Reads supporting (aligned to +/- strand): ref base A (10/23); new base G (2/6); total (12/29) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-03 |
GGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC > NC_002947/5485002‑5485297 | ggtggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATa < 8:248429/149‑1 (MQ=255) ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtg < 3:5090/149‑1 (MQ=255) ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGAtg < 7:298673/149‑1 (MQ=255) ttCGATGCGTGATGCTGAAGTCATCGCTGCTACCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCggcg > 5:96456/1‑146 (MQ=255) gTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAtt > 5:331213/1‑148 (MQ=255) cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc < 4:323935/149‑1 (MQ=255) cTGAAGTCATCGCTGCTACCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc < 6:96456/149‑4 (MQ=38) aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt < 4:219963/138‑1 (MQ=255) aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt > 3:219963/1‑138 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc < 1:75130/149‑1 (MQ=37) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCTTAAAAAGATGCTCGAAGTTC‑CTGTGGGAGCGGGCTTTCCCGCGAATAAGATGGTGGATTAACcgccgcct > 1:311442/1‑146 (MQ=255) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCT‑T‑G‑CGC‑CGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc < 3:11084/149‑1 (MQ=14) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg < 6:83360/149‑1 (MQ=37) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCTCGAATACGATGGTGGATTCACCGCCGCATTCGCgg < 2:180772/148‑1 (MQ=255) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg < 7:301339/148‑1 (MQ=37) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg < 1:349425/148‑1 (MQ=37) gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg < 7:53711/148‑1 (MQ=255) tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCACCGCATTCGCGGGCAAGc < 6:319115/149‑1 (MQ=255) cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCg < 8:307241/149‑1 (MQ=255) cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc < 6:331213/149‑1 (MQ=255) aTGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGGTTGCCCGCCCCATTCGCGGGCCAGCCCGCTCCCCCAGTGTTCGgcgc > 7:183718/1‑149 (MQ=255) ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTTCCCCCGGATACGATGGGGGGTTGACCGCCGCATTTGCGGGGAAGCCCGCGGCCACCGGGTTCGGCgccagcca > 3:130080/1‑149 (MQ=255) agtGGCTTGCGCCATTTCCGCCACTAATTCGGCGTATCCCGTGGAA‑T‑G‑TCT‑AGATC‑TTGGGGGCGGGGGTTTGCCCGCTAATACTTTGTTGTATTCCCCGCCGCATTGGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg < 8:33656/147‑1 (MQ=255) aggggCATGCGCAACTTCCGCCATAAATTCCGCGTACCCCCTGCAA‑C‑G‑TCG‑ATATC‑GTGGGGGCGGGGGCTTGCCGGCTAATACTTGGGGGTATTCACCGCCGCTTTCGGGGGAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg < 2:21589/145‑1 (MQ=255) gcCACTTCCGCCACTAATTACGCGTATCCCTTGGCA‑C‑G‑TCG‑AGATC‑CTGTGGGAGCGGGGTGGCCGGCTAATACGATGGTGGATTCACCGCCGCATTGGCGGGCAAGCCCGCTCCcaca < 6:311975/119‑1 (MQ=37) gcCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTTCCCGCGCATAAGCTGGTGGGTTGACCGCCGCCTTAGCGGGCCAGCCCGCGCCccca > 5:311975/1‑119 (MQ=37) ccGCCACTAATTCCGCGGATCCCGCGGCA‑A‑G‑TCG‑AGATC‑CTGGGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg < 2:311442/148‑1 (MQ=255) ccACTAATTACGCGGATCCCGAGCCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACt < 7:86744/136‑1 (MQ=255) gtcA‑A‑G‑TCG‑ATATC‑CTGTGGGAGGGGGCTTCCGCGCAAATACTATTGTGGATTCCCCGCCTCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc < 3:241221/147‑1 (MQ=255) tCG‑AGATC‑CTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCACATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGc < 8:183718/149‑1 (MQ=255) tgGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGTTCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 4:284040/149‑1 (MQ=255) gAGCGGGATGGCCCGCTAATACGATGGTGGTTTCACCGCCGCATTGGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCa < 7:212335/149‑1 (MQ=255) gggCTTGCCCGCTAAAACGATGGGGGATTCCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc < 7:156053/149‑1 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa > 8:348054/1‑116 (MQ=255) ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa < 7:348054/116‑1 (MQ=255) cgacagAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct < 1:324562/144‑1 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 4:174702/1‑149 (MQ=255) cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg > 7:323718/1‑148 (MQ=37) aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg < 3:174702/149‑1 (MQ=255) tACGATGGTTGATTCACCGCCTCTTTCGCGGGCAAGCCCGCTCCCACATTGTTCGGCGCAAGCCACGGGATCACTTGAATCTAGGTTTTGACATTAAAGTTTTGATCATCg < 4:16913/111‑1 (MQ=255) tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg > 3:16913/1‑111 (MQ=255) tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc > 7:29377/1‑149 (MQ=255) | GGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCA‑A‑G‑TCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC > NC_002947/5485002‑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: ‑ |