| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | CP009273 | 2,087,615 | G→T | 19.3% | D223Y (GAC→TAC) | hisB → | fused histidinol‑phosphatase/imidazoleglycerol‑phosphate dehydratase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | CP009273 | 2,087,615 | 0 | G | T | 19.3% | 65.0 / 3.2 | 30 | D223Y (GAC→TAC) | hisB | fused histidinol‑phosphatase/imidazoleglycerol‑phosphate dehydratase |
| Reads supporting (aligned to +/- strand): ref base G (11/13); new base T (3/3); total (14/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01 | |||||||||||
CGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGTTTTGTGCTACCGATGGACGAATGCCTTGCCCGCTGCGCGCTGGATATCTCTGGTCGCCCGCACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCT > CP009273/2087392‑2087862 | cGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCCCGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATcac > 2:492197/1‑244 (MQ=255) cgcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACcaca < 1:492205/244‑1 (MQ=255) aCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATcgac < 1:321251/228‑1 (MQ=255) tGGCAAATGTTTGGCGAGCAACTCACCACACGTTACCTTTACGCTCACGTAGCGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCTACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCg < 1:65988/274‑1 (MQ=255) gCCAATGATTGGCGAGCAACTCACCAGTCGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATAGACGATcac > 1:404192/1‑223 (MQ=255) gCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGCGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCTACCTCTATATCGCCGATCACCACACCGGCGAAGATACCGGCCTGGCGCTGGGCGAAGCGTNAAAAAGCGcc > 2:65988/1‑274 (MQ=255) gTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCTACCTCTATATCGACGATcacca > 2:551609/1‑179 (MQ=255) tGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCTACCTCTCTATCGACGATCACCACAcc < 1:551620/179‑1 (MQ=255) cGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGGTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCAAACCGTAGAAGATACCGGCCTGGCGCTGGGCGAAGCGCCAAAACTCGCCCTCCTAAACACACGCGGTCATTGCCTGCTtttttttttacctcctttc > 2:351812/1‑259 (MQ=255) cGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCAGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGTTTTGTGCTACCGATGGACGa < 1:351817/275‑1 (MQ=255) gagaCTCAGATTGACGTTCGGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTAGATCGTCGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGaa < 1:261240/195‑1 (MQ=255) aGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTatat > 1:530029/1‑130 (MQ=255) gtgGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATcgac < 2:530018/130‑1 (MQ=255) gCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTatat > 2:637401/1‑88 (MQ=255) gCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATcac > 1:9341/1‑98 (MQ=255) tGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGATATCAACGTCAAAGGCGACCTCTATATCGACGATCACcaca < 2:9341/98‑1 (MQ=255) tGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATcgac < 1:637412/88‑1 (MQ=255) aGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCTACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTaa > 1:97873/1‑111 (MQ=255) gATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGttt > 1:585239/1‑153 (MQ=255) cGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCTACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAAt < 2:97873/111‑1 (MQ=255) gCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGTTTtgtg < 2:585228/153‑1 (MQ=255) ggcggTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTatat > 2:307944/1‑44 (MQ=255) ggcggTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAg > 1:166131/1‑67 (MQ=255) gcggTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCgg > 2:219004/1‑109 (MQ=255) gTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATcgac < 1:307949/44‑1 (MQ=255) gTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATAc < 2:166128/67‑1 (MQ=255) tttCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGAc < 1:219008/109‑1 (MQ=255) tttCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGTTTTGTGCTACCGATGGACGAATGCCTTGCCCGCTg > 1:202623/1‑168 (MQ=255) ttCCGCATGGAAATCAATGTCAAAGGCGACCTCTATAGCGACGATAACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGTTTTGTGCTACCGATGGACGAATGCCGTGCCCGCTGCGCGCTGGATATCTCTGGTCGCCCGCACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGct < 1:324718/275‑1 (MQ=255) cGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGTTTTGTGCTACCGATGGACGAATGCCTTGCCCGCTgcgcg < 2:202619/168‑1 (MQ=255) | CGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCATGGAAATCAACGTCAAAGGCGACCTCTATATCGACGATCACCACACCGTCGAAGATACCGGCCTGGCGCTGGGCGAAGCGCTAAAAATCGCCCTCGGAGACAAACGCGGTATTTGCC‑GCTTTGGTTTTGTGCTACCGATGGACGAATGCCTTGCCCGCTGCGCGCTGGATATCTCTGGTCGCCCGCACCTGGAATATAAAGCCGAGTTTACCTACCAGCGCGTGGGCGATCTCAGCACCGAAATGATCGAGCACTTCTTCCGTTCGCT > CP009273/2087392‑2087862 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |