Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,889,008 | 0 | A | C | 55.6% | 10.9 / 18.5 | 18 | V44G (GTG→GGG) | ygcQ | electron transfer flavoprotein subunit alpha/FixB family protein |
Reads supporting (aligned to +/- strand): ref base A (3/5); new base C (0/10); total (3/15) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.86e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.47e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ACAGTTCATCGCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTT > NZ_CP009273/2888925‑2889091 | ccAGTTCATCCCCAACCGTCCCCGGGGGAAAGAAAACACCGCCCGCGGGCGTTCGCTGCCACTGTTCCACCAACGCACCCAGCCCCTGtt < 2:321815/89‑1 (MQ=255) aGTTCACCGCCAACCGTCCCCGGCGGAAAGAAAACAACGCCCGCGGGCGTTCGCTGCCCCTGTTCCACCAACGCACCCAGCCCCTGTTcc < 2:725981/90‑1 (MQ=255) tCGCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACt > 2:331474/1‑90 (MQ=255) gCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTAc < 1:261532/90‑1 (MQ=255) aaaCGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAgg > 1:319268/1‑90 (MQ=255) cGTCCCCGGCGGAAAGAAAACACCGCCCCCGGCCGTCCGCTGCCACTGTTCCACCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTg < 2:730373/90‑1 (MQ=255) cccGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTGAGGt < 1:56313/90‑1 (MQ=255) cccGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGt < 1:1597/90‑1 (MQ=255) cccGGCGGAAAGAAAACACCGCCCGCGGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGt < 2:319268/90‑1 (MQ=255) gcggAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCCCCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCaa < 1:175790/90‑1 (MQ=255) ccgcccGCGGGCGTCCGCTGCCACTGTTCCCCCAACGCACCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc < 2:655792/86‑1 (MQ=255) acctCCGCGGGCTTTCGCTGCCACTGTTCCACCAACGCACCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc < 2:7429/87‑1 (MQ=255) acGTCCGCGGGCGTTCGCTGCCACTGTTCCCCCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc < 1:5686/90‑1 (MQ=255) acGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGcc < 1:331474/90‑1 (MQ=255) cGTTCGCTGCCCCTGTTCCACCAACGCATCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAAcc < 1:112865/90‑1 (MQ=255) gCTGCCACTGTCCCACCAACCCCTCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGa < 1:712625/90‑1 (MQ=255) aaCGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCgccagc > 1:540717/1‑90 (MQ=255) tCCAGCCCCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCtt < 1:485624/90‑1 (MQ=255) | ACAGTTCATCGCCAAACGTTCCCGGCGGAAAGAGAACAACGTCCGCAGGCGTTCGCTGCCACTGTTCCACCAACGCATCCAGCACCTGTTCCGCCACTACAGGCTGAGGTTCAATCTGCCAGTGCGCCAGGGTACAACCGCTGAAATCCTGCGCCGCCAGCCAGCTT > NZ_CP009273/2888925‑2889091 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |