Predicted mutation | ||||||
---|---|---|---|---|---|---|
evidence | seq id | position | mutation | annotation | gene | description |
RA | NZ_CP009273 | 2,885,672 | T→C | S33P (TCC→CCC) | queD → | 6‑carboxytetrahydropterin synthase QueD |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,885,672 | 0 | T | C | 92.9% | 120.4 / 1.9 | 42 | S33P (TCC→CCC) | queD | 6‑carboxytetrahydropterin synthase QueD |
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (23/16); total (24/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.67e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01 |
TACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTC > NZ_CP009273/2885542‑2885804 | tACGGCGGTTTTATCGTCCTGTAGAGAAATTTTGGTGGCCACCACGTTATTTAAAGATTTCACCTTCGAAAGCGCTCACCGCTTACCACACGGCCCGGAAGGGCATAAATGGGGGCGCCTGCACGGGCATCCCTTTATg > 1:324193/1‑139 (MQ=255) cggcggTTTTATCGTCCTGTAGAGAAATTATGATGTCCCCCACGGTATTTTAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGt > 2:177539/1‑139 (MQ=255) ttttttCGTCCCGTAGAGAAATTATGATGTCCCCCACGTTTTTTTAAGATTTCACCTTTGAAGCCGCTCAACCCTTACCACCCGTCCCGGAAGGGGATAAATGTGGTCGCCCGCCAGGGGATCCCTTTATGGGGCGAcc > 2:250384/1‑138 (MQ=255) aTTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTc > 1:843037/1‑138 (MQ=255) atgatgTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGAt < 2:192805/138‑1 (MQ=255) atgatgTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATc > 1:560993/1‑139 (MQ=255) caccacGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACg > 2:381249/1‑139 (MQ=255) tATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGACGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGAt > 2:809904/1‑139 (MQ=255) tATTTAAAGATTTCACCTCAGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCCTAAATGTGGTCGCCTGCACGATCTTCCCTTTATGGTGCGACTGGTAATTACAGGGGACGACGATCCGCATACGGGCTGGAt > 2:175940/1‑139 (MQ=255) tttAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTa > 2:409358/1‑139 (MQ=255) tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGa > 2:76013/1‑91 (MQ=255) tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGa < 1:76013/91‑1 (MQ=255) tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGa < 2:239832/133‑1 (MQ=255) tCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGa > 1:239832/1‑133 (MQ=255) ccGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAgc < 1:409358/139‑1 (MQ=255) cggtCACCGCTTACCCCACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAgcg < 2:643445/136‑1 (MQ=255) aCCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaacc > 1:510672/1‑138 (MQ=255) aCCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACggg > 1:87060/1‑96 (MQ=255) aCCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACggg < 2:87060/96‑1 (MQ=255) aCCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAaac < 1:358885/137‑1 (MQ=255) gTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTAc > 1:493547/1‑138 (MQ=255) gTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTAc > 2:38995/1‑138 (MQ=255) ggCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCg < 1:630610/66‑1 (MQ=255) ggCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCg > 2:630610/1‑66 (MQ=255) ggCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGa < 1:687795/139‑1 (MQ=255) ggCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGa < 2:193069/139‑1 (MQ=255) ggCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGa < 2:493547/139‑1 (MQ=255) gCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt < 2:59857/139‑1 (MQ=255) gCATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt < 1:177539/139‑1 (MQ=255) cATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCAt < 1:40040/72‑1 (MQ=255) cATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCAt > 2:40040/1‑72 (MQ=255) cATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt > 1:153490/1‑138 (MQ=255) cATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt > 1:251184/1‑138 (MQ=255) cATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt > 2:354299/1‑138 (MQ=255) cATAAATGTGGTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt > 2:338551/1‑138 (MQ=255) ggTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGc < 2:72547/86‑1 (MQ=255) ggTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGc > 1:72547/1‑86 (MQ=255) ggTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt > 2:170580/1‑129 (MQ=255) ggTCGCCTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGAt < 1:170580/129‑1 (MQ=255) ccTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGCACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAAt > 2:78046/1‑139 (MQ=255) cTGCACGGGCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATg < 1:32245/139‑1 (MQ=255) gggCATCCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTc < 1:338551/139‑1 (MQ=255) | TACGGCGGTTTTATCGTCCTGTAGAGAAATTATGATGTCCACCACGTTATTTAAAGATTTCACCTTCGAAGCCGCTCACCGCTTACCACACGTCCCGGAAGGGCATAAATGTGGTCGCCTGCACGGGCATTCCTTTATGGTGCGACTGGAAATTACCGGGGAAGTCGATCCGCATACGGGCTGGATTATCGATTTCGCTGAACTAAAAGCGGCGTTTAAACCAACCTACGAGCGCCTCGATCACCATTATCTCAATGATATTC > NZ_CP009273/2885542‑2885804 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |