| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 5,777,377 | C→G | 17.4% | Y546* (TAC→TAG) | betA‑II → | choline dehydrogenase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 5,777,377 | 0 | C | G | 17.4% | 55.7 / 3.9 | 23 | Y546* (TAC→TAG) | betA‑II | choline dehydrogenase |
| Reads supporting (aligned to +/- strand): ref base C (9/10); new base G (2/2); total (11/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
GGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAA > NC_002947/5777231‑5777516 | ggCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGtactac > 2:200710/1‑150 (MQ=255) aGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCACGTAGTACTTGGCgg < 1:320781/150‑1 (MQ=255) gggCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCggg < 2:218315/150‑1 (MQ=255) gCCGGTGATGATCAGCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGAt > 2:575725/1‑124 (MQ=255) gCCGGTGATGATCAGCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGAt < 1:575725/124‑1 (MQ=255) cTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCAAACAAGATCCGTGGTCGCAACCCACGGCCGCGCAGGACCGCCAATGACTACTTGGCGGGGGCTGCACCAGTGAAAGGCGAGCCGATGCgtgaa > 2:376260/1‑150 (MQ=255) tcatcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCgtgaagtga > 2:414837/1‑150 (MQ=255) atcaCCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGaaa < 2:351883/150‑1 (MQ=255) cTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGc > 2:598843/1‑150 (MQ=255) cGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGCTTGCAAGCCACTACCGCGCAGCTCCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGa > 1:12258/1‑150 (MQ=255) caccacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCa < 1:200710/150‑1 (MQ=255) accacGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAg > 2:458045/1‑150 (MQ=255) gCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGc < 1:414837/150‑1 (MQ=255) aagaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCTGC‑GCGGCCCCtga < 1:598843/150‑3 (MQ=255) aagaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCttt > 2:381286/1‑150 (MQ=255) atgcGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGc < 2:13798/147‑1 (MQ=255) gTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAgc < 1:458045/150‑1 (MQ=255) gCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAggg < 1:589784/106‑1 (MQ=255) gCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAggg > 2:589784/1‑106 (MQ=255) cACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAt < 1:381286/150‑1 (MQ=255) ccGCCACGTAGTACTTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgtg > 2:330849/1‑138 (MQ=255) ccGCCACGTAGTACTTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtgtg < 1:330849/138‑1 (MQ=255) ccGCCACGTAGTACTTGGCGGGGGATGCAACAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGaa > 2:229045/1‑150 (MQ=255) ccGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGaa > 2:93465/1‑150 (MQ=255) | GGCATGAAGGGCCTGCGCGTGGTCGATGCCTCGATCATGCCGCTGATCATCACCGGCAACCTCAACGCCACCACGATCATGATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGC‑GCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAA > NC_002947/5777231‑5777516 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |