Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 443,161 | 0 | T | A | 52.9% | 9.5 / 21.0 | 17 | M1M (ATG→TTG) † | cyoE | heme o synthase |
Reads supporting (aligned to +/- strand): ref base T (7/1); new base A (0/9); total (7/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.11e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATAGCGATGATTAG > NZ_CP009273/443078‑443247 | cagcagGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAAccgc > 2:236580/1‑90 (MQ=255) aaccccccAATCCCCGAGATCAGGTTCCAAAAGATGATGCCTGTTTTCTTTCCTTGCTGGTATTGTTC‑AAACATCATAACCGCCGCTCtt < 1:218581/87‑1 (MQ=255) ccGAGTTCGGTTTCCAAAAGATGTTGCTTGTTTTTGTTCCTTGCGGTTTTTGTTA‑AAAAATCAAAACCGCCGCTCTTAGTGcatcatcat < 2:235042/90‑1 (MQ=255) ccGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAACCGCCGCTCTTAGTGcatcatcat > 2:396890/1‑90 (MQ=255) ccGAGATAGGTTTGAAAAAGTTGTGGCCTGTTTTTTTTACTTGCGGTTTTTGTTA‑AAAAATAAAACCCGCCGCTCTTAGTGcatcatcat < 2:573519/90‑1 (MQ=255) agaTCAGGTTCCAAAAGATGATCCTTGGTTTGGTTCCTTGCAGTTTTTCTTA‑AAACAACAAAACCGCCGCTCTTAGTGCATCATCAtgtt < 1:290068/90‑1 (MQ=255) gaTCAGTTTCCAAAAGTTGATGCCTGTTTTCTTTTCTTGCAGTTTTTGTTA‑AAACATAAAAACCGCCGCTCTTAGTGCATCATCAtgttg < 2:178084/90‑1 (MQ=255) ggTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGtt > 1:107283/1‑90 (MQ=255) ttGCAAAAGATGATCCTTGTTTTCTTTTCTTGCAGGTTTTGCTA‑AAACACAAAAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGa < 1:244445/90‑1 (MQ=255) tGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAg > 2:552457/1‑90 (MQ=255) aaaGAGGTTCCCTGTTTTCTTTTTTTG‑AGGTTTTTCTTAAAACATAAAACCCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTc < 2:153242/90‑1 (MQ=255) atgCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACAtct > 2:304169/1‑88 (MQ=255) atgCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATaa > 1:247125/1‑90 (MQ=255) tctttACTTGCAGTTTTTGCTA‑AAACATCAAAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGa < 1:236580/87‑1 (MQ=255) tcTTGCAGGTTTTGTTA‑AAACATCAAAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTa < 1:89031/89‑1 (MQ=255) cggttATTGTTA‑AAACATCAAACCCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAAcc < 2:15261/86‑1 (MQ=255) gCTT‑AAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATAGc > 1:243312/1‑90 (MQ=255) tcatAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATAGCGATGATTAg > 1:112758/1‑90 (MQ=255) catAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTa < 2:33059/67‑1 (MQ=255) catAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTa > 1:33059/1‑67 (MQ=255) | CAGCAGGAATCCCCCAATCACCGAGATCAGGTTGCCAAAGATGATGCCTGGTTTCGTTACTTGCAGGTATTGCTT‑AAACATCATAACCGCCGCTCTTAGTGCATCATCATGTTGTAGTTGAGGTTCCACATAATCCAGATGGAGCCTACAACCAGGATAGCGATGATTAG > NZ_CP009273/443078‑443247 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |