Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,016,669 | 0 | T | A | 41.4% | 17.9 / 38.3 | 29 | S26T (TCA→ACA) | sxy | CRP‑S regulon transcriptional coactivator Sxy |
Reads supporting (aligned to +/- strand): ref base T (8/9); new base A (7/5); total (15/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.10e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.35e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGCACAGTACTGTGTAAAACATCCGC > NZ_CP009273/1016555‑1016786 | tgAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGa < 2:97676/139‑1 (MQ=255) aaTCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGAcc > 2:588687/1‑139 (MQ=255) aaTCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCAAAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGAcc > 1:432451/1‑139 (MQ=255) tatgATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGgtgt > 2:100277/1‑139 (MQ=255) ttAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGAt > 2:483949/1‑139 (MQ=255) aTCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAg > 2:68830/1‑138 (MQ=255) aTCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAg > 2:199602/1‑138 (MQ=255) aTCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATgg < 1:223577/113‑1 (MQ=255) aTCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATgg > 2:223577/1‑113 (MQ=255) aTCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAg > 2:233303/1‑138 (MQ=255) aTCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAg > 2:223584/1‑138 (MQ=255) taAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCaaa < 2:177627/137‑1 (MQ=255) taAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCaaa < 1:199602/137‑1 (MQ=255) taAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGc < 1:588687/134‑1 (MQ=255) aTCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca < 1:483949/139‑1 (MQ=255) aaGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGCACAGTa < 1:556709/139‑1 (MQ=255) aaGAATACCTGGCAACGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGCACAGTa < 1:223584/139‑1 (MQ=255) ggCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCg < 2:328064/105‑1 (MQ=255) ggCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCg > 1:328064/1‑105 (MQ=255) aCGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca > 1:484736/1‑120 (MQ=255) aCGTTGGGCACAATTGAATACCGAACATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca < 2:484736/120‑1 (MQ=255) ggCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca > 2:210456/1‑114 (MQ=255) ggCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca < 1:210456/114‑1 (MQ=255) aTCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca > 1:583195/1‑97 (MQ=255) aTCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca > 1:534443/1‑97 (MQ=255) aTCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca < 2:534443/97‑1 (MQ=255) aTCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGca < 2:583195/97‑1 (MQ=255) aTCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGCACAGTACTGTGTAAAACATCCGc < 2:311364/119‑1 (MQ=255) aTCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGCACAGTACTGTGTAAAACATCCGc > 1:311364/1‑119 (MQ=255) | TGAATCATTCTTTTATGTTATGATTTTAAAAGGAATTTTATGAAAAGCCTCTCCTATAAGCGGATCTATAAATCACAAGAATACCTGGCAACGTTGGGCACAATTGAATACCGATCATTGTTTGGCAGTTACAGCCTGACCGTTGACGACACGGTGTTTGCGATGGTTTCTGATGGTGAGTTGTATCTTCGGGCTTGTGAGCAAAGTGCACAGTACTGTGTAAAACATCCGC > NZ_CP009273/1016555‑1016786 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |