Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,001,832 | 0 | T | G | 60.0% | 7.9 / 18.6 | 20 | G25G (GGT→GGG) | yedE | selenium metabolism membrane protein YedE/FdhT |
Reads supporting (aligned to +/- strand): ref base T (3/5); new base G (12/0); total (15/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.61e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGC > NZ_CP009273/2001746‑2001911 | gagGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGAtt > 1:199064/1‑90 (MQ=255) gCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTa < 2:258531/85‑1 (MQ=255) tACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATtctctc > 2:157022/1‑90 (MQ=255) tcatgtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCtactt < 1:743959/90‑1 (MQ=255) atgtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCtacttac > 2:349053/1‑90 (MQ=255) tgtcatgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGTTTTTTTCttcttact > 2:743959/1‑90 (MQ=255) catgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTctctct < 1:160195/80‑1 (MQ=255) catgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTctctct > 2:160195/1‑80 (MQ=255) catgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTAtt > 1:1240967/1‑90 (MQ=255) tgGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTAtttt > 2:406116/1‑90 (MQ=255) gCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTTTCTTCTTTCTATTTTTGCATTACTGGCACCttt > 1:584125/1‑90 (MQ=255) gCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTTTCTCTACTTTTTATTTTGGCATTACTGGCACCttt > 2:824717/1‑90 (MQ=255) gCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTTTCTACTTACTATTTTGGCATTACTGGCACCttt > 2:1079715/1‑90 (MQ=255) aTTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTc < 2:1240967/90‑1 (MQ=255) aTTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACggg < 2:199064/90‑1 (MQ=255) tGGGCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTTCTTACTATTTTTGCATTACTGGCACCTTTTGGGCTGTCACGGGGGaa > 1:416790/1‑90 (MQ=255) gggCCCCCATCCCTGCGGTCATCGCGGCGGGGATTCTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTTGGCTGTCACGGGTGAAt > 1:384302/1‑90 (MQ=255) ccATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTAccc > 2:257070/1‑90 (MQ=255) aTCGCGGCGGGTATTCTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTGGGCTGTCACGGGGGAATTTACCCGTTGGGGCGGGCAg > 2:867058/1‑90 (MQ=255) tCGCGGCGGGGATTCTCTCTACTTACTATTTTTGCATTACTGGCACCTTTTGGGCTGTCCCGGGGGAATTTACCCCTTGGGGCGGCCAGc > 1:267631/1‑90 (MQ=255) | GAGGCTACATTCATGTCATGGCAGCAATTCAAACACGCCTGGTTGATTAAATTCTGGGCCCCCATCCCTGCGGTCATCGCGGCGGGTATTCTCTCTACTTACTATTTTGGCATTACTGGCACCTTTTGGGCTGTCACGGGTGAATTTACCCGTTGGGGCGGCCAGC > NZ_CP009273/2001746‑2001911 |
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: ‑ |