| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,329,948 | 0 | T | G | 55.6% | 19.1 / 21.5 | 18 | intergenic (+24/‑140) | BW25113_RS25560/acnA | hypothetical protein/aconitate hydratase AcnA |
| Reads supporting (aligned to +/- strand): ref base T (2/6); new base G (10/0); total (12/6) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.51e-03 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGT > NZ_CP009273/1329865‑1330025 | ggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCAt > 1:225100/1‑90 (MQ=255) gggCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGGTTTTCCTCttt > 2:239781/1‑90 (MQ=255) ccGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTAt < 1:542198/90‑1 (MQ=255) ctctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGGTTCTCCTCTTTTATCAATTTgg > 1:285015/1‑90 (MQ=255) ctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgt > 1:117707/1‑90 (MQ=255) gATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgttgt < 2:117707/90‑1 (MQ=255) tCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCCTTTTCTTTTCCACAAGGTTTTTCCTCTTTTTTCAATTTTGgttgttt > 2:451583/1‑89 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGTTTCTCCTCTTTTTTCAATTTGGGTTGTTATc > 1:509035/1‑90 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGTTTTTCCACAAGGGTTTTCCTCTTTTTTCAATTTTGGGTGTTATc > 2:774939/1‑90 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATc > 1:238475/1‑90 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGGTGTTATc > 1:310517/1‑90 (MQ=255) tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCcaca > 1:317816/1‑43 (MQ=255) tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCcaca < 2:317816/43‑1 (MQ=255) tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgc < 2:183865/90‑1 (MQ=255) ccTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgcg < 1:239781/90‑1 (MQ=255) cGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGTTTCTCCTCTTTTTTCAATTTTGGGTGTTATCAAATCGTTACGCGATGTTtttg > 2:889593/1‑90 (MQ=255) cGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAAGTTTTTCCTCTTTTTTCAATTTGGGGTGTTATCAAATCGTTACGCGATGTTtttg > 2:81904/1‑90 (MQ=255) gCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTgtgt < 2:227121/90‑1 (MQ=255) | GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGT > NZ_CP009273/1329865‑1330025 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |