Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NC_002947 | 5,903,165 | 0 | A | C | 40.0% | 19.2 / 16.8 | 30 | intergenic (‑129/‑139) | PP_5175/PP_5176 | HlyD family secretion protein/hypothetical protein |
Reads supporting (aligned to +/- strand): ref base A (5/13); new base C (7/5); total (12/18) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.36e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAAT > NC_002947/5903018‑5903305 | gCGGACGCAGCAACGCCAGGCGGGCTCCTTGCTTGTGCACAAATGATTTGTGGGAGGCGACTGTGAGCGTAGTTCAGGGGTCGATAAACATCACTGACGGGGGCCATGCGGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGa < 7:142315/148‑1 (MQ=255) ggACGCATCAACGCCATGCTGGCTCCTTGCTGGTGCACAAATGATTTGTGGGAGGCGACTGTGAGCGTAGGTCAGGGTGCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGc < 5:119307/148‑1 (MQ=255) aCGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCGGGGGCGccc > 4:30947/1‑148 (MQ=255) gCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGCGCCcg > 5:77943/1‑147 (MQ=255) tCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGGTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGCGCCCGCGaa < 7:12804/148‑1 (MQ=255) cAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCTCGCCCGCGAag < 4:154541/148‑1 (MQ=255) aaCGCCATGCTGGCTCCTTGCTTGTGCACAAATTATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGCGCCCGCGAaga > 7:7038/1‑148 (MQ=255) aCGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACCGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCGTTTTTGGTGCCTGCGGGGGCCACTTGGCGGGGGCGCCcgcccagac > 7:62022/1‑142 (MQ=255) aCGCCATGCTGGCTCCTTGCTTGGGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGGTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCGCTTCGCGGGCGCGCCCGCGAaga < 7:30935/147‑1 (MQ=255) aCGCCATGCGGGCTCCTTGCGTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTGCAGGGGTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGGGGAGGCCTCTTCGCGGGCTCGCCCGCGAaga < 7:74224/147‑1 (MQ=255) tCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCGGGCGCGCCCGCGCAGAGGCCGCCACACCCa > 1:147365/1‑148 (MQ=255) ttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCa < 8:30763/148‑1 (MQ=255) ttgTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCa < 7:1585/148‑1 (MQ=255) aaaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTAGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGc > 7:43965/1‑148 (MQ=255) aaaTGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCAATTAGCGGGGGAGCCCGCGAAGCGGCCTCCCCCCCCATGTATCAGTATCGGc > 8:80307/1‑148 (MQ=255) gtgGGAGGCGCCTGTGCGCGTAGTTCAGGGTTCTATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGGGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCa < 1:59512/147‑1 (MQ=255) ttCTGGGTTCGATAAACATCAGTGAGGGGGGCCAGGCTGGCTTTTTTGGTGCCGGGGGGGGCCTCTGCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATc < 1:101842/147‑1 (MQ=255) ggggggCCATGCGGGCATTTTTGGGGCCTGGGGGGGCCTCTGCGCGGGCGAGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAg < 2:39915/147‑1 (MQ=255) gggCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGa < 3:39686/148‑1 (MQ=255) ggCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTACCGGGCGCGCCCGCCCAGAGGCCTCCACAGCCATGTTTAAGTATATGCCATATACAGCAACCAAAAGGCATAGGCCCAACCCCG‑TACAACCGCAACCATTAACAGGGAg > 6:65055/1‑148 (MQ=255) tGCTGGCTTTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTca < 8:62022/148‑1 (MQ=255) tGGCATTTTTGGTGCGTGCGGGGGCCTGTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTcaca < 2:147365/147‑1 (MQ=255) aTTTTTGGGGCGGGGGGGGGCCTCTGCGCGGGCGCGCCCGCTAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCaa < 1:162364/148‑1 (MQ=255) tttttGGTGCCGGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGGTCCAACCGCTCCCTATCACAGGGAGGCTCACATTCAAg < 8:136994/148‑1 (MQ=255) cGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTTCCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATg < 7:116515/148‑1 (MQ=255) gggggCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGc < 6:85455/148‑1 (MQ=255) tctTCGCGGGCTCGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAACCCATCAGGCCTCGG‑CCATCCCCGTTCCACCCGCTACCTTTCACAGGGCGGCTCCCATTCAAGGAAGATCAAAATGCCGGTAGc > 6:155000/1‑148 (MQ=255) tctTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCCTCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGc > 8:56224/1‑148 (MQ=255) tctTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGc > 4:159520/1‑148 (MQ=255) cgGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGA‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAAt < 8:105622/147‑1 (MQ=255) gAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGtt > 8:112483/1‑134 (MQ=255) aGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTa < 7:112483/134‑1 (MQ=255) | GCGGACGCATCAACGCCATGCTGGCTCCTTGCTTGTGCACAAATGATTTGTGGGATGCGACTGTGAGCGTAGTTCAGGGTTCGATAAACATCACTGACGGGGGCCATGCTGGCATTTTTGGTGCCTGCGGGGGCCTCTTCGCGGGCGAGCCCGCGAAGAGGCCTCCACAGCCATGTATCAGTATCTGCCATATCCAGAAACCATCAGGCATCGG‑CCATCCCCGTTCCAACCGCTACCTTTCACAGGGAGGCTCACATTCAAGGAAGATCAAAATGCCGTTAGCCGAAT > NC_002947/5903018‑5903305 |
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: ‑ |