Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,502,016 | 0 | A | G | 52.4% | 24.7 / 29.9 | 21 | F386S (TTC→TCC) | yhfW | phosphopentomutase |
Reads supporting (aligned to +/- strand): ref base A (0/10); new base G (11/0); total (11/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.84e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA > NZ_CP009273/3501928‑3502103 | aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGaa < 2:266023/90‑1 (MQ=255) gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGgaaaaat > 1:462642/1‑88 (MQ=255) gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGGGGgaaaaat > 2:112612/1‑88 (MQ=255) gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGgaaaaat > 1:221304/1‑88 (MQ=255) gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGgaaaaat > 1:276090/1‑88 (MQ=255) tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATc > 1:305709/1‑90 (MQ=255) ggATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTcac < 2:153763/90‑1 (MQ=255) tAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGcccc < 1:112612/90‑1 (MQ=255) tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAATTCACACACGGTAGCCCCCCCAt > 1:111926/1‑90 (MQ=255) gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCATc > 1:56619/1‑90 (MQ=255) gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGGCGGGAAAAATTACACACGGGAGGCCCCCCATc > 1:226190/1‑90 (MQ=255) gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCcccaac > 1:145665/1‑88 (MQ=255) gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGGGCGGGAAAAATCACAAACCGGAGCCCCcccacc > 2:376972/1‑88 (MQ=255) cAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg < 2:175561/90‑1 (MQ=255) gaaagaaaAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACAAGGGAGCCCCCCCCACACGAAACGTGGGGCGCCCAccc > 2:242407/1‑89 (MQ=255) aagaaaagaGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATAAAACACAGGAGCCCCCCCCCCAGAAAACGGGGGGGGCCCACCCAg > 2:127048/1‑90 (MQ=255) agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc < 2:260622/90‑1 (MQ=255) agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc < 1:127048/90‑1 (MQ=255) aagaGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATTACACACCGTAGCCCCCCCCACACAAAACGGGGGGGGGGCCCCGAGCTGCg > 1:400248/1‑90 (MQ=255) gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa < 2:350073/90‑1 (MQ=255) agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg < 2:145665/90‑1 (MQ=255) gCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACACCCCGCCGCGGAGc > 2:48373/1‑90 (MQ=255) gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc < 2:346977/90‑1 (MQ=255) aCCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGTAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCggt > 1:314144/1‑88 (MQ=255) ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa < 1:104703/90‑1 (MQ=255) ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa < 2:314144/90‑1 (MQ=255) | AACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA > NZ_CP009273/3501928‑3502103 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |