Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,141,214 | 0 | T | G | 73.1% | 4.6 / 17.1 | 26 | G20G (GGT→GGG) | yegD | molecular chaperone |
Reads supporting (aligned to +/- strand): ref base T (2/5); new base G (19/0); total (21/5) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.19e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAA > NZ_CP009273/2141136‑2141293 | ccGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCCTGCGTGGCGGGAAAACCCAttt > 2:262016/1‑90 (MQ=255) gCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGGGACGGGGAAACGCATTTTCt > 2:858439/1‑90 (MQ=255) gCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGGCATGCGGGGCGGGAAAACCCATTTTCt > 2:428070/1‑90 (MQ=255) ttattTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATTCGTGGCGGGAAAACGCCTTTTCTAAAAAAGGa > 2:605979/1‑90 (MQ=255) attTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGaaa < 1:1208774/90‑1 (MQ=255) gtgtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGGGGGAAAACGCCTTTTGTAAAAAAAGAAAAAGAc > 2:1143250/1‑90 (MQ=255) gtTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTGCTAAAAATGGAAAAAGACAg > 1:405637/1‑90 (MQ=255) tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGGCGGGAAAACGCATTTGTTAAAAAAGGAAAAAAACAGCACGctgctg > 1:742548/1‑90 (MQ=255) tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGctgctg < 1:141177/90‑1 (MQ=255) tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGctgctg < 1:262016/90‑1 (MQ=255) tttGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGctgctg < 2:607626/90‑1 (MQ=255) cAGCAAACTGTTCAGTGGCGGTCCAGCGTGTCGGGAAAACCCCTTTTCTTAAAAAAGAAAAAAAAAACACGCCGCTGGCCTCAATGCttt > 2:1362754/1‑90 (MQ=255) gCAAACTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCttttc > 1:274415/1‑88 (MQ=255) aaCTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCCTTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTtcgcg > 2:218073/1‑90 (MQ=255) aaCTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTtcgcg > 1:354057/1‑90 (MQ=255) aaCTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCCTCAATGCTTTtcgcg > 1:1188138/1‑90 (MQ=255) cTGTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAACGAAAGCACGCTTCTGCCTTTAATGGTTTgcgcggc > 1:51847/1‑88 (MQ=255) tGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCa < 2:742548/90‑1 (MQ=255) gTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTTCTAAAAAAGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGccca > 1:1195123/1‑88 (MQ=255) gTTCAGTGGCGGTCATGCGTGACGGGAAACCGCATTTGCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGccca > 1:1405439/1‑88 (MQ=255) gTTCAGTGGCGGTCATGCGGGACGGGAAAACGCATTTGCTAAAAATGGAAAAAGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCaa > 1:381575/1‑90 (MQ=255) cAGTGGCGGTCATGCGTGACGGGAAACCGCATTTGCTAAAAAAGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTTCGCGCCAAcgc > 1:47816/1‑90 (MQ=255) gCGGTCATGCGTGGCGGGGAAACCCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCCATGCTTTTCGCGCCAACCCGCGaa > 2:1029225/1‑90 (MQ=255) ggTCCTGCGTGGCGGGAAACCCCATTTTCTAAAAAAAGAAAAAGACAGCACCCCGCTGCCTTCAATGCTTTTCGCGCCCACCCGTGAAGc > 2:815106/1‑90 (MQ=255) ggTCATGCGTGACGGTAAACCGCATTTTCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGc > 1:757107/1‑90 (MQ=255) tGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTaa > 2:1084634/1‑90 (MQ=255) | CCGGCATGGGGATTATTTCGTGTTTATTGGTTTTGATTACGGTACAGCAAACTGTTCAGTGGCGGTCATGCGTGACGGTAAACCGCATTTGCTAAAAATGGAAAACGACAGCACGCTGCTGCCTTCAATGCTTTGCGCGCCAACGCGTGAAGCGGTAA > NZ_CP009273/2141136‑2141293 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |