Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,506,593 | G→A | 14.8% | intergenic (+52/‑147) | PP_4839 → / → cycA | iron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,506,593 | 0 | G | A | 14.8% | 59.8 / 4.0 | 27 | intergenic (+52/‑147) | PP_4839/cycA | iron‑regulated membrane protein/D‑alanine/beta‑alanine/D‑serine/glycine permease |
Reads supporting (aligned to +/- strand): ref base G (16/7); new base A (2/2); total (18/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.82e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 |
GGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTAT > NC_002947/5506469‑5506741 | ggTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCCCGTTGTAGGAGCGGGTTTa > 4:314983/1‑149 (MQ=255) gATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAAACCCCACCAACGACGGACCGCGTTGTAGGATCGGGTTTACCCCCGAAGAGGCCGGTACAGGCGAAAt > 8:41302/1‑149 (MQ=38) gtggtggTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGGCTTGTGTTGCCTGTGCTGGCCGCTTCGCGGGTAACCCCGCTCGTCCGAGGGACGGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGa < 3:314983/147‑1 (MQ=38) gccgccGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACcgcg > 6:272364/1‑85 (MQ=255) gccgccGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACcgcg < 5:272364/85‑1 (MQ=255) ggCCAGCGCTTGAGCCTTGTGTTGCCGGTGCTGGCCTCTTCGCGGGTAACCCCGCTCCTACGAGGGACCGCGTTGTAGCAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGc < 7:281192/119‑1 (MQ=255) gAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTg > 8:335918/1‑149 (MQ=255) gtgtTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcg > 7:275003/1‑149 (MQ=255) cTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTtagcct < 7:335918/148‑1 (MQ=255) ccTCTTCGCGGGTAAACCCGCTCCTCCGAGGGCCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGcc < 4:90765/118‑1 (MQ=255) ccTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGGGTTGTAGGAGCGGGTTTACCCGCGAAGCGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCTGAGTCGATCTGTCGCGGCCCGCTCTGGCACGCGCGTGTTAGCCTAGGCg > 1:259248/1‑149 (MQ=38) ccTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGcc > 3:90765/1‑118 (MQ=255) ctctTCGCGGGTAAACCCGCTCCTACGAGGGCCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTgg > 5:336524/1‑126 (MQ=255) ttCGCGGGTAAACCCGCTCCTACGAGGGGCCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTa > 4:332329/1‑149 (MQ=255) tAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGa > 2:67351/1‑148 (MQ=255) aaaCCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCGTCGCCGGCTCCCTCTGAc > 6:244817/1‑148 (MQ=255) aCCCGCTCCTACGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGAc > 3:111599/1‑146 (MQ=255) aCCCGCTCCTACGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGAc > 3:106245/1‑146 (MQ=255) ccGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACaag > 4:137901/1‑147 (MQ=255) gCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAagaaga > 3:101736/1‑149 (MQ=255) tCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcgc < 5:169434/113‑1 (MQ=255) tCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAAcgcgc > 6:169434/1‑113 (MQ=255) cTACGAGGGACCGCGTTGTAGGTGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGc < 7:41302/149‑1 (MQ=255) cTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTg > 4:182665/1‑148 (MQ=255) aCGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGTCTACCTCTGACAAGAAGACTGCAg > 7:205273/1‑149 (MQ=255) cGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCa < 4:111599/147‑1 (MQ=255) cGAGAGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCa < 4:106245/147‑1 (MQ=255) gggACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGACTACCTCTGACAAGAAGACTGCAGCCTa > 2:64759/1‑149 (MQ=255) ggACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCtatgggtgaggtttgcgtcataggctgcagtcttcttgtcagaggtagccggctaggctaacac > 2:195769‑M1/1‑85 (MQ=255) | GGTGTCGATGGTGGTGATGTGGGGGGTGGATAGTTTGGTGGTGAGGCGCCGCCGGGTCTTGGCCAGCGCTTGAGCCTTGTGTTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACGAGGGACCGCGTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACAGGCGAAATGCTGAGCCGAAAATGCAGAGTCGATCTGTCGCGCCCCGATCTGGAACGCGCGTGTTAGCCTAGCCGGCTACCTCTGACAAGAAGACTGCAGCCTAT > NC_002947/5506469‑5506741 |
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: ‑ |