| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 3,140,032 | 0 | A | C | 52.2% | 23.1 / 26.1 | 23 | V22G (GTT→GGT) | yghW | DUF2623 domain‑containing protein |
| Reads supporting (aligned to +/- strand): ref base A (7/4); new base C (0/12); total (7/16) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.35e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGC > NZ_CP009273/3139946‑3140112 | ggCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCg > 2:322866/1‑90 (MQ=255) ttttCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGc > 2:541103/1‑90 (MQ=255) ttCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGt < 1:53552/90‑1 (MQ=255) aCATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCtt > 2:363432/1‑90 (MQ=255) aGCCTAATAAAAACCCGCTTTTATAACCGGCAAAAAATTTTGTCACATTACCCGCACTGTCGGAATGCGTTGCTTTTAATCCCGCCATTa < 2:56233/90‑1 (MQ=255) aaTACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCg > 2:309592/1‑79 (MQ=255) aaTAAAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCg < 1:309592/79‑1 (MQ=255) aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCa > 2:393243/1‑90 (MQ=255) cgcgTTTATAATCGGCTTAAAATTTTGTCACTTTACCCGCGCAGTGGGGGTGCGTTGGTTTTAATCCCGCCATTAAGCCTTTACCAAAAt < 2:212553/90‑1 (MQ=255) ttATAATGGGAAAAAATTTTTGCCAAATTACCCCCCCGGCGGGAATGTGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTa < 1:363432/90‑1 (MQ=255) taACCGGCACAAAATTTTGTCACATTACCCCCACTGCCGGCTTGCTTTGTTTTTAACCCCGCCATTAAGCCTTTACCAAAATGGTTATTc < 2:359219/90‑1 (MQ=255) aTCGGCACAAAATTTTGTCACATTACCCGCCCTGCCGGCATGCGTTGTTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCAt < 1:322866/90‑1 (MQ=255) ggCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGt > 2:197070/1‑90 (MQ=255) gCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGTTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCAt < 1:264373/86‑1 (MQ=255) gCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCAt > 2:264373/1‑86 (MQ=255) cacaGAATTTTGTCACATTACCCGCCCTGTCGGCATGCGTTGCTTTTAACCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTc < 1:386702/90‑1 (MQ=255) tttGTCACATTACCCGCCTTGTCGGCATGCTTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCAtt < 1:309325/90‑1 (MQ=255) tttGTCACATTACCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCAtt < 1:228927/90‑1 (MQ=255) tttGTCACATTACCCCCACTGTGGGCATGCTTTGTTTTTAACCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCAtt < 2:485407/90‑1 (MQ=255) gTCACATTACCCGCACTGTCGGCATGCGTTGTTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGc < 1:493040/90‑1 (MQ=255) gTCACATTACCCGCACTGTCGGCATGCGTTGTTTTTAACCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGc < 1:299107/90‑1 (MQ=255) gTCACATTACCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGc < 1:237395/90‑1 (MQ=255) gTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGc < 1:197070/90‑1 (MQ=255) | GGCGATCTCCGGTCTTTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGC > NZ_CP009273/3139946‑3140112 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |