Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,783,688 | A→G | 7.3% | L1209P (CTG→CCG) | 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,688 | 0 | A | G | 7.3% | 158.1 / 2.6 | 55 | L1209P (CTG→CCG) | pvdJ | non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base A (27/24); new base G (2/2); total (29/26) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01 |
TGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG > NC_002947/4783541‑4783836 | tGTCGCGCCCATTGGCGCAGGGGTGCCGACGTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAg > 1:21930/1‑149 (MQ=255) gcgcCCATTGGCTCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGt < 7:97120/149‑1 (MQ=255) gcgcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGt > 2:30336/1‑149 (MQ=255) gcgcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGt > 1:34715/1‑149 (MQ=255) gcCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTg > 8:3420/1‑149 (MQ=255) ccATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCt > 6:132436/1‑149 (MQ=255) ccATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCt < 2:190433/149‑1 (MQ=255) ttGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGAc < 1:104878/149‑1 (MQ=255) ttGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGAc < 2:34715/149‑1 (MQ=255) gcgcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCACTAGGCTAGCCTGATGCTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCAt > 2:38024/1‑149 (MQ=255) gcAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCAACg > 7:201659/1‑149 (MQ=255) ccGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCg < 2:7940/149‑1 (MQ=255) cGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCg > 3:106189/1‑148 (MQ=255) cTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACaa > 1:234696/1‑149 (MQ=255) aGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCt < 2:234696/149‑1 (MQ=255) cTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa < 2:40512/149‑1 (MQ=255) aCCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAAc < 2:21930/147‑1 (MQ=255) cccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCAc > 5:265000/1‑149 (MQ=255) ccTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATCTGGCCCCCAGTTCATGCACCGCCAGCAA‑CACCGGTTGATCACAAACGGCCATTTCGAACAATGACCCCCTGAACACCTCGcacc > 1:80990/1‑146 (MQ=255) gCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGc > 5:157339/1‑147 (MQ=255) ctgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGa < 7:3420/149‑1 (MQ=255) tgctgGTAGGCCGTCTGCATGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCATCAACCAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa < 1:191824/126‑1 (MQ=255) tgctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGaa > 2:191824/1‑125 (MQ=255) ctgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAgg < 6:157339/148‑1 (MQ=255) tgGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt > 1:72870/1‑149 (MQ=255) gTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCt < 6:265000/147‑1 (MQ=255) gTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTg < 2:223524/148‑1 (MQ=255) gCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTc < 4:106189/148‑1 (MQ=255) gTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTgg < 5:132436/149‑1 (MQ=255) gTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGAAccaagtttc > 8:94141/1‑142 (MQ=255) tCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTggc < 2:72870/149‑1 (MQ=255) aCAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg > 3:210553/1‑115 (MQ=255) aCAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAg < 4:210553/115‑1 (MQ=255) aTGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAg < 1:30336/149‑1 (MQ=255) gcCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCt > 1:182962/1‑149 (MQ=255) gcCAAGACACCCCATC‑CACCACCAGGAGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCt > 7:155879/1‑149 (MQ=255) ccAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTg > 2:64434/1‑149 (MQ=255) acacCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCAACGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAgg < 8:201659/149‑1 (MQ=255) acCCCATC‑CACGACCAGGTGATGCACCGCCAGCAA‑CAGCCGGTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCGACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGGGGTTTCTAGCGCTGCAAGGCt < 8:155879/149‑1 (MQ=255) acCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCt > 6:52809/1‑149 (MQ=255) cccATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGt > 7:153937/1‑149 (MQ=255) aTC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCa < 1:64434/149‑1 (MQ=255) caccaccaGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCagag > 3:163096/1‑149 (MQ=255) ccaGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGAc > 1:36465/1‑148 (MQ=255) ccaGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACt < 2:36465/149‑1 (MQ=255) aGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTAc > 2:159890/1‑149 (MQ=255) ggTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAgc < 3:128257/122‑1 (MQ=255) ggTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAgc > 4:128257/1‑122 (MQ=255) gATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTg > 4:244273/1‑148 (MQ=255) tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGgcggcg > 6:243603/2‑101 (MQ=255) tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGgcggcg < 5:243603/100‑1 (MQ=255) tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg < 3:143722/145‑1 (MQ=255) tcAGCAA‑CGGCTGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCaggagg > 4:143722/2‑146 (MQ=255) aGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTg < 7:94141/149‑1 (MQ=255) gCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTg > 8:18872/1‑148 (MQ=255) aGCCGTTGCTGCCCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGGTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCgg > 6:141024/1‑149 (MQ=255) | TGTCGCGCCCATTGGCGCAGGGGTGCCGACCTGGACGGCAGGCTCACCGCCTGCCCTAGGCTAGCCTGCTGGTAGGCCGTCTGCAGGTCTTCGAACAGGATGCGCCAAGACACCCCATC‑CACCACCAGGTGATGCACCGCCAGCAA‑CAGCCGTTGCTGACCATCGGCGAGTTCGAACAATGTCGCCCTGAACAGCTCGCCACGGGCCAGGTCGAGGCTGGCCTGCACCTCGTTGGCGGCGGTTTCTAGCGCTGCAAGGCTGTCCAGAGTGACTACCTGCAGGAGGGCCTGCACCGG > NC_002947/4783541‑4783836 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |