Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,502,021 | 0 | T | A | 47.6% | 12.1 / 13.1 | 21 | E384D (GAA→GAT) | yhfW | phosphopentomutase |
Reads supporting (aligned to +/- strand): ref base T (1/10); new base A (10/0); total (11/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.40e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA > NZ_CP009273/3501934‑3502103 | 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) | TCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA > NZ_CP009273/3501934‑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: ‑ |