Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 510,080 | 0 | A | C | 51.7% | 36.6 / 44.5 | 29 | V79G (GTG→GGG) | ybbJ | NfeD family protein |
Reads supporting (aligned to +/- strand): ref base A (13/1); new base C (0/15); total (13/16) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.06e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCA > NZ_CP009273/509994‑510164 | ccAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCAccc > 2:287331/1‑90 (MQ=255) ggggAGATTCCAGCAAAAACCGCCGCCCAACCCCCGCCCCCCCCCGCTGTTTA‑AAATGACTTTCGCTGTTTTTTTGTCCGCCCCCccgcc < 2:1045239/88‑1 (MQ=255) cccaccaAAAAACGCCGCCCAATCGGCGCCTCCCCGCGCTGTTTT‑AAATGCCTATCGCTGTGTTTTTGTCCGCGCCCCCGCCGCGACAAc < 1:1199744/89‑1 (MQ=255) caaaaaaCGCCGCCCAATCAGCTGCCCCCCGCGTTGTTTA‑AAATGCCTATCGTTGTGTTTTTGTTCGCGCCCCCGCCGCGACAACCAttt < 2:562742/87‑1 (MQ=255) aaaCGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTccac > 1:788357/1‑90 (MQ=255) tCGGCCAATCAGCTGCTGCCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCcaccacct > 1:1018688/1‑88 (MQ=255) gCCAATCAGCTGCTCCCCGCGCTGTTTA‑AAATGCCTATCGCTGTGTTTTTGTTCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCa < 1:635392/90‑1 (MQ=255) gctgctGCCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGa > 1:149699/1‑90 (MQ=255) tgctgcCCGCGTTGTTTA‑AAATGCCTTTCGCTTTTTTTTTGTCCGCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagc < 2:475982/90‑1 (MQ=255) tgccgcCCGCGCTGTTTA‑AAATGACTCCCGTTGTTTTTTTTTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagc < 2:209797/90‑1 (MQ=255) cgccgcCCGCG‑TGGTTTAAAATGACTCCCGTTGTGTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGagc < 2:1054453/86‑1 (MQ=255) gcCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCg > 1:349910/1‑90 (MQ=255) gcCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCg > 2:598435/1‑90 (MQ=255) gcgcTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAg > 1:209797/1‑90 (MQ=255) ccgcTGTTTA‑AAATGCCTCCCGTTGTGTTTTTTTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAg < 2:1002042/89‑1 (MQ=255) ggtGGTTA‑AAATGCCTCCCGCTGTGTTTTTGTCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGa < 2:1206672/88‑1 (MQ=255) gcTGTTTT‑AAATGACTATCGCTGTTTTTTTTTTCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGa < 2:349910/90‑1 (MQ=255) ttaaaaTGCCTCCCGTTGTGTTTTTTTCCCCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCa < 2:149699/87‑1 (MQ=255) taaaaTGCCTCCCGCTGTGTTTTGGTCCCCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCaa < 2:702371/88‑1 (MQ=255) aCTATCGCTGTGTTTTTGTTCGCGCACCCGCCGCGACAACCATTTccacc < 1:342107/50‑1 (MQ=255) aCTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTccacc > 2:342107/1‑50 (MQ=255) aCTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAc > 2:523542/1‑90 (MQ=255) cTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAcc > 2:1071759/1‑90 (MQ=255) cTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAcc > 2:1047360/1‑90 (MQ=255) cTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACAATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTAcc > 1:1028556/1‑90 (MQ=255) cTGTGTTTTTGTTCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGc < 1:363402/90‑1 (MQ=255) tgtgTTTTTGTCCCCGCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGcc < 2:1268935/90‑1 (MQ=255) tttttttCCGCCCCCCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTc < 1:1155072/89‑1 (MQ=255) cgcgCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCa > 2:540111/1‑90 (MQ=255) | CCAGCGGAGATTCCAGCACAAAACGTCGGCCAATCAGCTGCTGCCCGCGCTGGTTT‑AAATGACTATCGCTGTGCTTTTGTTCGCGCACCCGCCGCGACAACCATTTCCACCACAGCCAGGCGGCGAGCAGCGTCAGGATGGCAAACATTACCCCTTGCCACTCCCAACCCA > NZ_CP009273/509994‑510164 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |