Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,502,021 | 0 | T | A | 50.0% | 15.7 / 12.2 | 22 | E384D (GAA→GAT) | yhfW | phosphopentomutase |
Reads supporting (aligned to +/- strand): ref base T (0/11); new base A (11/0); total (11/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.84e-06 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA > NZ_CP009273/3501950‑3502103 | aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccccc < 2:332235/90‑1 (MQ=255) aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGccccc < 1:502221/90‑1 (MQ=255) gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGGCCCcccaac > 1:187770/1‑88 (MQ=255) gtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACAGTAGCCCCCACATc > 1:140672/1‑90 (MQ=255) ctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCcccacccc > 2:3194/1‑86 (MQ=255) tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAAATCACACACGGTAGCCCCCCCATCAGaaaa > 1:136816/1‑89 (MQ=255) aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtg < 1:368840/90‑1 (MQ=255) aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCgtgg < 2:228488/90‑1 (MQ=255) ccGGAGGGAGGAAAGAAAAGAGCGACCAATTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAACgtggg > 1:64697/1‑89 (MQ=255) gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGCGGGAAAAATCACACACGGTAGCCCCCCCATCAAAAAACgggggggg > 1:443658/1‑86 (MQ=255) gggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCacac < 2:423251/89‑1 (MQ=255) agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc < 1:525215/90‑1 (MQ=255) agaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc < 1:165089/90‑1 (MQ=255) aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTg < 2:574838/90‑1 (MQ=255) gagCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACGGTAGCCCCCCCATCACAAAGCGGGGGGGGCACACCCCGCTGCGTa > 1:251919/1‑90 (MQ=255) gagCGACCATTTTGTGGCGGGGGCGCGCGGGAAAAATTACACACCGGAGCCCCCCCATCAAAAAGAGTGGGGGGCACACCGAGCTGCGtc > 1:423251/1‑89 (MQ=255) agCGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATTACACACGGGAGCCCCCCCCTCAGAAAAAGGGGGGGGCCCCCCCAGCTGCGGAg > 1:372091/1‑90 (MQ=255) agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg < 2:359052/90‑1 (MQ=255) gCGACCATTTTGTGGCGGTGGCGCGCGGGAAAAATCACACACCGGAGCCCCCCCCTCAAAAAACGGGGGGGGCGCACCCAGCCGCGTAGc > 2:332256/1‑90 (MQ=255) cGACCATTTTGTGGCGGTGGGGCGCGGGAAAAATCACACACCGGAGCCCCCCCATCAAAAAACGTGGGGGGCCCACCCCGCCGCGTAGCg > 2:96680/1‑90 (MQ=255) ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa < 2:16514/90‑1 (MQ=255) ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa < 2:573470/90‑1 (MQ=255) | AGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCA > NZ_CP009273/3501950‑3502103 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |