Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,026,452 | 0 | A | C | 52.4% | 24.0 / 26.7 | 21 | G141G (GGT→GGG) | yccA | FtsH protease modulator YccA |
Reads supporting (aligned to +/- strand): ref base A (7/3); new base C (0/11); total (7/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAATGAAGACCGCGCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAAC > NZ_CP009273/1026367‑1026533 | gAATGAAGACCGCGCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCg < 2:214189/90‑1 (MQ=255) aTGAAGCCCGCGCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATCCCAATCAGCCCCACCACAATCCCCGcc < 1:394162/90‑1 (MQ=255) cgcgCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCAt > 2:249529/1‑90 (MQ=255) gcaggcagCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCt < 2:875968/90‑1 (MQ=255) ggcagCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTg > 1:749365/1‑90 (MQ=255) tGCAGGAAGATATTCGCACCCATCCCAACCACCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGggtg < 1:515153/90‑1 (MQ=255) gCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTgtgggtgg > 1:561898/1‑90 (MQ=255) ggAAGATATTCGCAACCATACCAATCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCa < 2:737130/90‑1 (MQ=255) ggAAGATATTCGAACCCATACCAATCACCCCCACCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCa < 1:845678/90‑1 (MQ=255) ggAAGATATGCGCACCCATACAAACCAGCACCCCCAAAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCa < 1:493269/90‑1 (MQ=255) aaGATATTCGCACCCATACCAACCACCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGc < 1:598227/90‑1 (MQ=255) aaGATATTCGCAACCATCCCAATCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGc < 1:319917/90‑1 (MQ=255) aaGATATTCGCAACCATCCCAACCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGc < 2:684237/90‑1 (MQ=255) aaGATATTCGCAACCATACCAATCAGCCCCACCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGc < 2:8338/90‑1 (MQ=255) acccATACCAACCAGCCCCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCaga < 1:787552/88‑1 (MQ=255) aaTCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCagaa > 2:645435/1‑90 (MQ=255) tCAGCACCCCCACAATCCCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCagaaga < 2:491919/90‑1 (MQ=255) caccaccacAATACCCGCCATCAGCATACCGCCGAGGAAc < 2:58782/40‑1 (MQ=255) caccaccacAATACCCGCCATCAGCATACCGCCGAGGAAc > 1:58782/1‑40 (MQ=255) accacAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAAc > 2:129883/1‑90 (MQ=255) accacAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAAc > 2:100487/1‑90 (MQ=255) | GAATGAAGACCGCGCTGATCGCCAGATGCAGAGCAGGCAGCTGCAGGAAGATATTCGCAACCATACCAATCAGCACCACCACAATACCCGCCATCAGCATACCGCCGAGGAACGACATATCTTTGCGGGTGGTCAGCACATATGCAGAGCAGCAGAAGAACACTAAC > NZ_CP009273/1026367‑1026533 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |