| Predicted mutation | |||||||
|---|---|---|---|---|---|---|---|
| evidence | seq id | position | mutation | freq | annotation | gene | description |
| RA | NC_002947 | 4,783,736 | C→T | 12.1% | R1208Q (CGG→CAG) | pvdJ ← | non‑ribosomal peptide synthetase |
| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NC_002947 | 4,783,736 | 0 | C | T | 12.1% | 82.3 / 3.7 | 33 | R1208Q (CGG→CAG) | pvdJ | non‑ribosomal peptide synthetase |
| Reads supporting (aligned to +/- strand): ref base C (14/15); new base T (2/2); total (16/17) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01 | |||||||||||
GCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCT > NC_002947/4783590‑4783874 | gcgcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGtt > 1:479168/1‑150 (MQ=255) aTTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGAc > 2:92945/1‑150 (MQ=255) tGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGTCAGCAACGGCTGTTGCTggcgg < 1:111779/150‑5 (MQ=255) cgcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCg < 1:491990/150‑1 (MQ=255) aCCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACa > 1:328169/1‑150 (MQ=255) tCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACa > 1:228995/1‑150 (MQ=255) aCCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGaa > 2:407489/1‑131 (MQ=255) aCCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGaa < 1:407489/131‑1 (MQ=255) ccTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGc < 1:290679/150‑1 (MQ=255) ccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACgg > 2:268924/1‑150 (MQ=255) ccTACGCTAGCTTGCTGGTAGGCCGTCTGCCGGTCTTCGAACAGGATGCGCCAAGACCCCCCATCCCCCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACgg < 2:104604/150‑1 (MQ=255) tAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGc > 2:330384/1‑150 (MQ=255) gctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGc < 2:328169/150‑1 (MQ=255) tgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTg > 1:371693/1‑150 (MQ=255) gCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTg < 1:33018/149‑1 (MQ=255) gCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTg > 2:33018/1‑149 (MQ=255) ccGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGTCAGCAACGGCTGTTGCTggcgg < 1:316695/86‑5 (MQ=38) ccGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGTCAGCAACGGCTGTTGCTggcgg > 2:316695/1‑82 (MQ=38) cTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTggcgg < 1:268924/150‑1 (MQ=255) gcCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTg < 1:31985/55‑1 (MQ=255) gcCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTg > 2:31985/1‑55 (MQ=255) cccATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGAGCTGCAAGGCTGTc > 1:126382/1‑150 (MQ=255) ccATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTggcggc < 2:160124/125‑1 (MQ=255) ccATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTggcggc > 1:160124/1‑125 (MQ=255) caccaccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGt < 2:371693/150‑1 (MQ=255) accaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACt < 2:126260/150‑1 (MQ=255) ccaGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTa > 1:94603/1‑150 (MQ=255) gATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCa < 1:92945/150‑1 (MQ=255) cGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGc < 1:330384/150‑1 (MQ=255) gctaGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGcc < 2:127605/147‑1 (MQ=255) tcAGCAACGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCt > 1:370674/2‑150 (MQ=255) cAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggag < 2:264607/138‑1 (MQ=255) cAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggag > 1:264607/1‑138 (MQ=255) | GCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCCACCACCAGGTGATGCACCGCCAGCAACAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCT > NC_002947/4783590‑4783874 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |