Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,329,948 | 0 | T | G | 52.2% | 19.9 / 25.4 | 23 | intergenic (+24/‑140) | BW25113_RS25560/acnA | hypothetical protein/aconitate hydratase AcnA |
Reads supporting (aligned to +/- strand): ref base T (4/7); new base G (12/0); total (16/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.35e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTA > NZ_CP009273/1329864‑1330033 | gggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCa > 2:68190/1‑90 (MQ=255) ggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCAt > 2:352233/1‑90 (MQ=255) ttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGTTTctcct > 2:842754/1‑90 (MQ=255) tGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgtt < 1:352233/90‑1 (MQ=255) ggATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCCTTTGTTTTTCCACCAGGGTTCTCCTCTTTTTTCAATTTGGgttg > 1:845176/1‑90 (MQ=255) gATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgttgt < 2:845176/90‑1 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGt > 1:698182/1‑36 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGt < 2:698182/36‑1 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCCTTTTCTTTTCCACAAGGTTTTTCCTCTTTTATCAATTTGGGGTGGTATc > 1:480514/1‑90 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAAGGTTCTCCTCTTTTTTCAATTTTGGGTGTTTTc > 1:951916/1‑90 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAAGTTTCTCCTCTTTTATCAATTTTGGTTGTTATc > 2:832188/1‑90 (MQ=255) tCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctt < 2:425411/43‑1 (MQ=255) tCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCAtttgctt > 1:425411/1‑43 (MQ=255) cTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAAGGTTCTCCTCTTTTTTCAATTTGGGTTGTTATCAAATCGtt > 1:738501/1‑90 (MQ=255) cTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCCCAAGGTTTTTCCTCTTTTTTCAATTTTGGGTGTTTTCAAATCGtt > 1:317797/1‑90 (MQ=255) aTCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCCCAAAGGTTCTCCCCTTTTTTCAATTTTGGGTGTTATGAAATATTTAct > 1:798036/1‑89 (MQ=255) tGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATg > 2:606330/1‑90 (MQ=255) tGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTGCCACAAGGTTTTTCCTCTTTTTTCAATTTGGGTTGTTATCAAATCGTTACGCGATg > 1:686433/1‑90 (MQ=255) cGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTTTCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgtg > 1:256641/1‑90 (MQ=255) cGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGGTGTTATCAAATCGTTACGCGATGTTtgtg > 2:111744/1‑90 (MQ=255) tGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa < 2:738501/90‑1 (MQ=255) tGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa < 1:606330/90‑1 (MQ=255) ggcggGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa < 1:271979/87‑1 (MQ=255) ggcggGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa > 2:271979/1‑87 (MQ=255) ggggTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTa < 2:951916/90‑1 (MQ=255) | GGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTA > NZ_CP009273/1329864‑1330033 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |