| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,485,207 | G→A | 9.8% | intergenic (+65/‑82) | 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,207 | 0 | G | A | 9.8% | 140.1 / 4.5 | 51 | intergenic (+65/‑82) | purH/purD | bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase |
| Reads supporting (aligned to +/- strand): ref base G (20/26); new base A (3/2); total (23/28) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.47e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01 | |||||||||||
GATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGA > NC_002947/5485058‑5485353 | gATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTATGGGAGCGGGCTTACCCGCGAATACGAtcgtgg > 2:312150/1‑150 (MQ=255) gAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAAGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCa < 2:1070079/150‑1 (MQ=255) tCATCGCTGCTGCCGACGAGGCCGGCATTGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcg < 2:114411/150‑1 (MQ=255) tCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcg > 1:542778/1‑150 (MQ=255) cATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATTCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAAAACTATGGTGGATTCACCGCCGCATTcgc < 1:32230/150‑1 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg > 1:1086044/1‑150 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg < 2:286764/150‑1 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg > 1:286764/1‑150 (MQ=255) aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTcgcg < 1:365745/150‑1 (MQ=255) tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg > 1:757543/1‑150 (MQ=255) ctgctgCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTG‑CGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGc > 2:638434/1‑150 (MQ=2) tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGcc > 2:590224/1‑150 (MQ=255) gcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CCGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGccc > 1:1057357/1‑150 (MQ=25) gCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCcaca < 1:590224/150‑1 (MQ=255) gcagcagcgATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAgtg < 1:893601/144‑1 (MQ=255) cATCGCGATGGTCTTCACTGGCATGCGACACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CCGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAggatc < 2:1057357/150‑5 (MQ=18) tCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCg < 2:1198633/150‑1 (MQ=255) ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑Ac > 1:529629/1‑150 (MQ=255) tCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACg > 1:757937/1‑150 (MQ=255) cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg < 1:584223/114‑1 (MQ=255) cACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg > 2:584223/1‑114 (MQ=255) cATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCAc > 2:736744/1‑150 (MQ=255) ccGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCCCCGCCGCATTTGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGt < 2:851942/150‑1 (MQ=255) gCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGcc < 2:714984/100‑1 (MQ=11) gCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGcc > 1:714984/1‑100 (MQ=11) aCTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGa < 2:757543/150‑1 (MQ=255) aCGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGaaa < 2:443598/150‑1 (MQ=255) ggATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTTGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt > 1:451624/1‑150 (MQ=255) ggATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGttt < 2:1086044/150‑1 (MQ=255) gATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt < 1:1161429/150‑1 (MQ=255) aTCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg < 1:736744/150‑1 (MQ=255) aTCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg < 2:542778/150‑1 (MQ=255) gCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCgg > 1:430085/1‑150 (MQ=255) g‑agaTC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGt < 1:641029/150‑1 (MQ=255) g‑agaTC‑CTGTGGGAGCGGGCTTGCCCGCGAATAAGATGGTGGGTTAACCACCGCATTCGCGGGCAAAGCCGCTGACACAGAGAGCGGCGCGAAGG‑GTGGGAGCACTTGTAGCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGt > 2:136569/1‑150 (MQ=255) cTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGt < 2:757937/150‑1 (MQ=255) cGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc < 2:529629/150‑1 (MQ=255) gCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctg > 2:418353/1‑150 (MQ=255) gCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct < 2:1013606/150‑1 (MQ=255) gCCCGCGAATAAGATGGTCGATTCTCCTCCGCATTCGCGGTCAAGCCAGCCCCCACAGTTTTC‑GTTCTGGCCGGCGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct < 1:625316/150‑1 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt < 2:1093937/110‑1 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt > 1:1093937/1‑110 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa < 1:490780/134‑1 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa > 2:490780/1‑134 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg > 1:1170199/1‑150 (MQ=255) ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctgg > 2:602273/1‑150 (MQ=255) cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGAGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATAATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGa > 2:194160/1‑150 (MQ=255) gAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg > 1:36194/1‑150 (MQ=255) aCGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGc < 2:36194/150‑1 (MQ=255) cGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGcc < 1:136569/150‑1 (MQ=255) gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGAc > 2:366422/1‑89 (MQ=255) gATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGAc < 1:366422/89‑1 (MQ=255) gtggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTGGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTAGCCTGGAAAGTCGCCCAGGa > 2:75551/1‑150 (MQ=255) tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc > 2:534588/1‑126 (MQ=255) tggATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc < 1:534588/126‑1 (MQ=255) | GATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCG‑AGATC‑CTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCC‑ACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCGCCCAGGA > NC_002947/5485058‑5485353 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |