| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,808,173 | G→A | 20.0% | intergenic (‑282/‑273) | PP_4234 ← / → dsbD‑II | bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2 |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,808,173 | 0 | G | A | 20.0% | 76.1 / 14.2 | 40 | intergenic (‑282/‑273) | PP_4234/dsbD‑II | bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2 |
| Reads supporting (aligned to +/- strand): ref base G (15/17); new base A (5/3); total (20/20) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 6.95e-01 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01 | |||||||||||
TATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAA > NC_002947/4808029‑4808320 | tataGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGACTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCG > 3:112359‑M1/1‑149 (MQ=255) cgcgCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTGAGCGCCTCGGAGATCGAACGCCGACCGCTcggagcttcgcggctaaac > 1:37809‑M1/1‑130 (MQ=255) gggTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacg > 5:62520‑M1/1‑119 (MQ=255) gTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGACTGGGAGATCGAGCGCCGCACGCGaggagcttcgcgcgtaaaccagctcctacgtttgtttcgggccaatta > 1:210107‑M1/1‑99 (MQ=255) gTTTTCACCGCATGGGACGTGTTCGCCTGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCTGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaatta < 2:210107‑M1/147‑49 (MQ=255) cggctgtGAGGTGTTCGCCGGGTGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcg < 4:112359‑M1/144‑57 (MQ=255) tcaggTTTTCGCCGGTCGATTTTCAGGGGCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGTCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggattt < 4:133013‑M1/145‑63 (MQ=255) tCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgg < 7:161566‑M1/90‑14 (MQ=255) tCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgg > 8:161566‑M1/1‑77 (MQ=255) acaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAACCTCCGGCAAACACCTCAAAATTCCAACGCCg < 5:215920‑M1/82‑1 (MQ=255) cGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcgg > 8:240217‑M1/1‑73 (MQ=255) cGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcgg < 7:240217‑M1/76‑4 (MQ=255) ggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAtcctcc > 2:262891‑M1/57‑111 (MQ=255) ggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAtcctcc < 1:262891‑M1/55‑1 (MQ=255) cAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtc > 2:376134‑M1/1‑62 (MQ=255) tcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAAACCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAAtt > 6:215920‑M1/55‑149 (MQ=255) ccTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgtgcgaccgccttggggcatgtctcgat < 2:109221‑M1/149‑93 (MQ=255) gggAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccaccttggggcatgtctcg < 7:356567‑M1/144‑91 (MQ=255) gggAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccaccttggggcatgtctcg > 8:356567‑M1/1‑54 (MQ=255) ggAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcg < 6:62520‑M1/149‑97 (MQ=255) gtgtTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtc < 8:272975‑M1/149‑101 (MQ=255) tgtTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgg > 7:241284‑M1/1‑48 (MQ=255) tgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGccccc > 6:323725‑M1/41‑149 (MQ=255) cGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctggatatcgcgtcgtac > 5:41106‑M1/1‑44 (MQ=255) cGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtac > 7:178361‑M1/1‑44 (MQ=255) ggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaa > 8:356525‑M1/37‑149 (MQ=255) gCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtaca < 1:376134‑M1/148‑106 (MQ=255) agTTTCAGCGCCTGGGAGTTCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcg < 7:195373‑M1/149‑115 (MQ=255) tttCAGCGCCTGGGAGATCGGGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggc < 1:227727‑M1/146‑114 (MQ=255) tttCAGCGCCTGGGAGATCGGGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggc > 2:227727‑M1/1‑33 (MQ=255) tCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcg > 1:80611‑M1/1‑31 (MQ=255) tCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtc < 8:349028‑M1/148‑118 (MQ=255) cgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcgcc < 5:347529‑M1/125‑1 (MQ=255) cgcCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcg < 2:80611‑M1/149‑123 (MQ=255) gcCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgc < 4:84309‑M1/149‑124 (MQ=255) ggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGa > 5:153984‑M1/12‑148 (MQ=255) gcgaagcgcCGCGCGGGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcac > 5:203433‑M1/5‑17 (MQ=255) gcgcCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggag < 8:282535‑M1/149‑138 (MQ=255) gcgcCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgc > 6:291016‑M1/1‑12 (MQ=255) gcgcCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttggttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgc < 5:291016‑M1/137‑126 (MQ=255) cgccgcCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagg < 7:92494‑M1/149‑139 (MQ=255) gccgccCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggaggg < 2:86571‑M1/149‑140 (MQ=255) cgcccgcGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgggcgcctggcacaggagggat < 7:223311‑M1/149‑142 (MQ=255) | TATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAA > NC_002947/4808029‑4808320 |
| 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: ‑ |