Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,118,197 | A→G | 12.5% | intergenic (+29/‑30) | PP_4504 → / → PP_4505 | hypothetical protein/sensor histidine kinase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,118,197 | 0 | A | G | 12.5% | 90.7 / 5.7 | 40 | intergenic (+29/‑30) | PP_4504/PP_4505 | hypothetical protein/sensor histidine kinase |
Reads supporting (aligned to +/- strand): ref base A (15/20); new base G (3/2); total (18/22) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.42e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.00e-01 |
CTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGCCTGGA > NC_002947/5118052‑5118338 | cTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTg > 8:262515/1‑149 (MQ=255) aCGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgg < 2:253151/149‑1 (MQ=255) aCGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgg > 4:168735/1‑149 (MQ=255) gCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgggcc < 3:339826/149‑1 (MQ=255) aCCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGcgggccggg < 7:262515/148‑1 (MQ=255) cagcagAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAGTGATGcgggccgggcat < 5:234278/149‑3 (MQ=255) agaagaAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTcc < 2:318367/149‑1 (MQ=255) gaagTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGtt > 2:110650/1‑149 (MQ=255) agTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGttt < 6:141513/148‑1 (MQ=255) gTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTag < 4:212703/149‑1 (MQ=255) tGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTagga < 1:110650/148‑1 (MQ=255) gACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTaggag < 4:183981/148‑1 (MQ=255) gAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTg > 7:44508/1‑148 (MQ=255) cAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGc < 8:44508/149‑1 (MQ=255) aaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTagga < 5:331388/136‑1 (MQ=255) aaaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTagga > 6:331388/1‑136 (MQ=255) aaaGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGtt < 1:21538/149‑1 (MQ=255) ccgccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTg < 8:168728/149‑1 (MQ=255) ccgccAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCAGTGATGcgggccgg > 8:53511/1‑108 (MQ=255) gCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTc > 5:190909/1‑149 (MQ=255) aGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctg > 8:146142/1‑149 (MQ=255) aGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTctggct < 6:190909/149‑1 (MQ=255) ccAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCt < 2:319513/120‑1 (MQ=255) ccAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCt > 1:319513/1‑120 (MQ=255) cAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTTGGAGGTGCACTTGCGTTCACTGTTGTGGCGCATCCTGGGCAGGTTCTGGCTGGCCa > 8:315745/1‑149 (MQ=255) aGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCAt > 4:47258/1‑149 (MQ=255) tGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCg > 8:259528/1‑147 (MQ=255) aaaaaGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAg < 3:47258/148‑1 (MQ=255) ggCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGt < 7:259528/149‑1 (MQ=255) ggCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCctg < 3:168735/149‑1 (MQ=255) ccAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTc > 2:356418/1‑128 (MQ=255) ccAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTc < 1:356418/128‑1 (MQ=255) cTGTCCGCCACATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGAACTGGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGAAGGCGTGTCGATCATGCTGGGGcagct > 2:236961/1‑146 (MQ=255) tCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCAcc < 5:293985/111‑1 (MQ=255) tCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCAcc > 6:293985/1‑111 (MQ=255) tGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGa < 7:146142/149‑1 (MQ=255) ccggcccgCAGTGATGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGc > 7:308278/1‑149 (MQ=255) gcccgCATCGCTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCaa < 7:76807/137‑1 (MQ=255) gcccgCATCGCTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCaa > 8:76807/1‑137 (MQ=255) ccgCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACACTGGTCGCAGGCCTGTCTCTCATGCTGGGGCACATGCTCAACCAGGACGCCTGGa > 4:234463/1‑149 (MQ=255) | CTCACGGCTGACCAGCAGAAGAAGTTTGACGAACTGCAAAAAGAACGCGCCGCCAAGAAAGCCGAGTGGCAGGAGTTCCAGGCCTGGAAAGCTGAAAAAGGCGCCAAGGCCCAGTAAGCTGTCCGCCCCATGCCCGGCCCGCATCACTGCGGGCCGGGCTTTTCCCGTTTAGGAGGTGCACTTGCGTTCACTGTTCTGGCGCATCCTGGCCAGTTTCTGGCTGGCCATCACCCTGGTCGCAGGCCTGTCGATCCTGCTGGGGCACATGCTCAACCAGGACGCCTGGA > NC_002947/5118052‑5118338 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |