Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | NC_002947 | 5,903,165 | A→C | 32.3% | intergenic (‑129/‑139) | PP_5175 ← / → PP_5176 | HlyD family secretion protein/hypothetical protein |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,903,165 | 0 | A | C | 32.3% | 18.8 / 23.2 | 31 | intergenic (‑129/‑139) | PP_5175/PP_5176 | HlyD family secretion protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (7/14); new base C (7/3); total (14/17) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.21e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.62e-01 |
GCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAA > NC_002947/5903024‑5903304 | gCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCcgc > 8:67682/1‑148 (MQ=255) aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAaga < 6:223628/148‑1 (MQ=255) gcttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTAt < 7:67682/148‑1 (MQ=255) ttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAg > 7:150174/1‑149 (MQ=255) aaaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGcc < 6:218900/149‑1 (MQ=255) aTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTAGC‑‑GGGCGAGCCCGCGAAGAGGGCGCCACAGCCATGTATAAGTATATGCCat > 3:97091/1‑149 (MQ=255) tGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAagag > 8:51856/1‑117 (MQ=255) tGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAagag < 7:51856/117‑1 (MQ=255) tgtgGGTTGCGCCTGGTTGCGTTGTTCAGGGTTCCAAAAACATCACTGAGGGGGGCAATGCTGGATTTTTGGGTGCCTGG‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCGGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCa < 4:97091/148‑1 (MQ=255) gggTTGCTACTGGGAGGGTTTGTCCGGGTTCGATAAACATCCCTGAGGGGGGCCATGCGGGAATTTTTGGTGCGGGG‑‑‑‑GGGGGCCTCTGCGC‑‑GGGCGACCCGGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGaaa < 6:175023/149‑1 (MQ=255) gCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTAGC‑‑GGGCGCGCCCGCGAAGAGGCCTCCCCAGCCATGTATAAGTATCTGCCATATACCGAAACCAt > 6:224832/1‑148 (MQ=255) cGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGAGGGGGCCAATTAGCAGGGGCGAGCCCGCGAAGAGGCCCCCGCAGGCac > 1:218806/1‑122 (MQ=255) taCTGTGAGCTTAGTTCTGGTTTCTATAAACATCACTGACGGGGGCCTTGCTGGATTTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCTAAGAGGCCTc < 2:6708/106‑1 (MQ=255) aCTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGAGGCCTCTTCGC‑‑GGGCTCGCCCGCGaa < 5:257394/98‑1 (MQ=255) aCTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGAGGCCTCTTCGC‑‑GGGCTCGCCCGCGaa > 6:257394/1‑98 (MQ=255) gAGCGTAGTTCCGGTTTCTTTACAAACCACTTCGGGGGGCCAGGCTGGCATTTTTGGGGCCGGG‑‑‑‑GGGGGCCTCTGCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCg < 6:149270/149‑1 (MQ=255) aGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGcc < 4:202392/124‑1 (MQ=255) aGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGcc > 3:202392/1‑124 (MQ=255) ttCAGGGCTCGATAAACATCACTGACGGGGTCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATccc < 8:128311/149‑1 (MQ=255) aCGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTcac < 5:224832/147‑1 (MQ=255) ggggCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGa < 1:257740/149‑1 (MQ=255) tGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTcac < 5:265439/149‑1 (MQ=255) tGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTATTAGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATATGCCATATCCACAAACCAAAAGGCAACGGCCAGACCCGGTCCAACCCCTACCATTCACAGGGAGGCGcacaa > 6:113525/1‑147 (MQ=255) cATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCaa < 8:150174/149‑1 (MQ=255) tctTCGC‑‑GGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCttt > 5:264773/1‑105 (MQ=255) tctTCGC‑‑GGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCttt < 6:264773/105‑1 (MQ=255) tctTCGC‑‑GGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGcc > 5:221644/1‑149 (MQ=255) ctTCGC‑‑GGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTcac > 1:133846/1‑119 (MQ=255) ttCGC‑‑GGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCg > 4:182002/1‑148 (MQ=255) tCGC‑‑GGGCTCGCCCGCGAAGAGGCCCCCGCAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGaa > 8:196042/1‑149 (MQ=255) cgc‑‑gGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCg < 3:182002/146‑1 (MQ=255) | GCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGC‑‑‑‑GGGGGCCTCTTCGC‑‑GGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGGCCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAA > NC_002947/5903024‑5903304 |
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: ‑ |