Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,140,032 | 0 | A | C | 53.3% | 18.0 / 21.3 | 15 | V22G (GTT→GGT) | yghW | DUF2623 domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base A (5/2); new base C (0/8); total (5/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 6.99e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGC > NZ_CP009273/3139961‑3140112 | tttCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCg > 2:253028/1‑90 (MQ=255) ttCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGt < 1:334363/90‑1 (MQ=255) aCATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCtt > 1:269950/1‑90 (MQ=255) gtgAGTACCCTAAAACAAAACCGCTTTTATAATCGGCACAGAATTTTGTCACATTACCCGCACTGTCGGCTTGCGTTGTTTTTAATCCCg < 1:43895/90‑1 (MQ=255) gAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGcc > 1:88614/1‑90 (MQ=255) aaaCCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCa > 2:35995/1‑90 (MQ=255) cccGCTTTTATAACCGCCACAGATTTTTGTCACATTACCCGCACTGTCGGCATGCGTTTTTTTTAACCCCGCCATTAAGCCTTTACCaaa < 1:240794/89‑1 (MQ=255) cgcgTTTATAATCGGCAAAAAATTTTGTCACATTACCCGCACTGTCGGCATGCTTTGTTTTTAATCCCGCCATTAAGCCTTTACCAAAAt < 1:97379/90‑1 (MQ=255) gcgTTTATAATCGGCACAGAATTTTGTCACATTAACCGCAc < 1:18284/41‑1 (MQ=255) gcgTTTATAATCGGCACAGAATTTTGTCACATTAACCGCAc > 2:18284/1‑41 (MQ=255) taATCGCCAAAAATTTTTGTCAAATTACCCGCACTGTCGGCATGCGTTGCTTTTAACCCCGCCATTAAGCCTTTACCAAAATGGTTATTc < 1:145736/90‑1 (MQ=255) aTCGGCACAGAATTTTGCCACATTACCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCAt < 1:140356/90‑1 (MQ=255) acaAAATTTTGTCACATTACCCCCCCTGCCGGCATGCGTTGTTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCa < 1:22039/90‑1 (MQ=255) tttGTCACATTACCCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCAtt < 2:235935/90‑1 (MQ=255) gTCACATTACCCGCCCTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGc < 2:283321/90‑1 (MQ=255) | TTTCGTACATCCGGTGTGAGTAGCCTAATACAAAACCGCGTTTATAATCGGCACAGAATTTTGTCACATTAACCGCACTGTCGGCATGCGTTGCTTTTAATCCCGCCATTAAGCCTTTACCAAAATGGTTATTCATCGTTCACCTCATTCGC > NZ_CP009273/3139961‑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: ‑ |