Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 4,779,454:1 | +G | 11.3% | intergenic (‑143/+18) | pvdD ← / ← pvdJ | non‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 4,779,454 | 1 | . | G | 11.3% | 190.1 / 21.0 | 53 | intergenic (‑143/+18) | pvdD/pvdJ | non‑ribosomal peptide synthetase/non‑ribosomal peptide synthetase |
Reads supporting (aligned to +/- strand): ref base . (22/25); new base G (3/3); total (25/28) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01 |
ACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGC > NC_002947/4779310‑4779594 | aCTTGACCATCCCCTACCAGGCCACCGGTCTTTGCCTGCGGCACGTTCCAGCCCGCGACAGGCTCATGAGGCTGCCGGACAAATGAAAGCGCCGGACGGCTGCCGGCTCTTGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCgg < 2:662789/150‑1 (MQ=255) aCTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCgg > 2:110078/1‑150 (MQ=255) aCATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGT‑ATGACCGGGCCCgg < 2:1486050/150‑2 (MQ=255) cATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCa > 2:1530753/1‑150 (MQ=255) tctcCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGGGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑AACGGACCCGCaga > 2:234347/1‑150 (MQ=255) tcCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAgc > 1:1282667/1‑150 (MQ=255) gggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTgg < 1:1037990/150‑1 (MQ=255) gggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTgg < 2:384095/150‑1 (MQ=255) gggCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTgg < 1:110078/150‑1 (MQ=255) aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCg > 1:195098/1‑114 (MQ=255) aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCg < 2:195098/114‑1 (MQ=255) aaCGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTg < 1:1530753/150‑1 (MQ=255) cGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCa < 2:970390/150‑1 (MQ=255) cGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCa < 1:739628/150‑1 (MQ=255) ccAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTaa < 2:129264/150‑1 (MQ=255) ccAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTaa < 2:1508783/150‑1 (MQ=255) aGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAAc < 2:921031/116‑1 (MQ=255) aGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGACAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAAc > 1:921031/1‑116 (MQ=255) tgcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAttt > 1:939433/1‑150 (MQ=255) gcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa > 1:1607799/1‑150 (MQ=255) gcCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa < 2:1607799/150‑1 (MQ=255) cGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGGAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGc > 1:772264/1‑150 (MQ=255) aaCAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACAACCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCa > 1:116595/1‑150 (MQ=255) aaaTGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcgg > 2:307202/1‑132 (MQ=255) aaaTGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcgg < 1:307202/132‑1 (MQ=255) gAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGAGGCCAAGGATTTAGCCAACTCGTc < 2:540831/150‑1 (MQ=255) aaaGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACtcgtc < 2:441693/95‑1 (MQ=255) aaaGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACtcgtc > 1:441693/1‑95 (MQ=255) ggACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAtt > 2:237493/1‑128 (MQ=255) ggACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGAtt < 1:237493/128‑1 (MQ=255) ggCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAAt > 2:1076245/1‑143 (MQ=255) ggCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAAt < 1:1076245/143‑1 (MQ=255) tGCCGGCTCTGTGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccc > 1:987142/1‑150 (MQ=255) cGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTcc > 2:611653/1‑132 (MQ=255) cGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTcc < 1:611653/132‑1 (MQ=255) gCTCTGGTTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgccc < 1:123970/150‑1 (MQ=255) gCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCgccctgccc > 1:1009111/1‑150 (MQ=255) tGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcg < 2:1495357/93‑1 (MQ=255) tGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAggcg > 1:1495357/1‑93 (MQ=255) tGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCg > 2:77898/1‑150 (MQ=255) aGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACg < 2:772264/150‑1 (MQ=255) gTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACgg > 2:584801/1‑150 (MQ=255) gACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCt < 2:939433/150‑1 (MQ=255) gACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAAATCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCt > 1:1559992/1‑150 (MQ=255) gctctgcGGGCCCGGT‑ATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc > 2:1044686/7‑122 (MQ=255) gctctgcGGGCCCGGT‑ATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACgcgc < 1:1044686/116‑1 (MQ=255) gcGGGCCCGGT‑ATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa > 1:777346/2‑92 (MQ=39) gcGGGCCCGGT‑ATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTa < 2:777346/91‑1 (MQ=39) cGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGaa < 2:1282667/150‑1 (MQ=255) cGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAAATCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGaa < 2:1559992/150‑1 (MQ=255) cGGGCCCGGT‑ATGACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGaa > 1:713718/1‑150 (MQ=255) ggCCCGGTATG‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGc < 2:996091/139‑1 (MQ=255) gCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGc > 1:1508472/1‑150 (MQ=255) | ACTTGAACATCTCCTACCAGGGCACCGGTCTTTGCCTGCGGAACGTGCCAGCCCGCGACAGGCTCATGATGCTGCCGAACAAATGAAAGAGCCGGACGGCTGCCGGCTCTGGGTGTACTGAAGTTGACGAGCCGGGCCCGGTCAT‑ACCGGGCCCGCAGAGCGCTACTGGGTCAACTCGTCGATCTCCTCTGCAGTAAGACGTTTCAAGGCGGCCAAGGATTTAGCCAACTCGTCCTGAATCGACGCGCCCTGCCCTTGCAACACTGCGACGGCCTCGCTGAACGC > NC_002947/4779310‑4779594 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |