Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,329,948 | 0 | T | G | 51.7% | 30.8 / 38.3 | 29 | intergenic (+24/‑140) | BW25113_RS25560/acnA | hypothetical protein/aconitate hydratase AcnA |
Reads supporting (aligned to +/- strand): ref base T (4/10); new base G (15/0); total (19/10) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.00e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTA > NZ_CP009273/1329865‑1330033 | ggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCAt > 1:77858/1‑90 (MQ=255) ttttCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctt < 2:77858/90‑1 (MQ=255) tCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAg < 2:797040/90‑1 (MQ=255) aTGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCCTTTTCTTTTCCCCAAGGtt > 2:70673/1‑90 (MQ=255) tGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCGTTGTTACAAggggt > 2:151022/1‑87 (MQ=255) cttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTctcc < 2:484499/90‑1 (MQ=255) gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTTCTTTTCCCCAAGGTTTTTCCTCTTTTTTCa > 2:263563/1‑90 (MQ=255) gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGGTTCTCCTCTTTTTTCa > 1:12371/1‑90 (MQ=255) gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGGCATTTTTTTTTCCACAAGGTTTCTCCTCttttttta > 2:835905/1‑86 (MQ=255) aCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAt < 2:768114/90‑1 (MQ=255) cctcTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCCTTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTg < 1:899613/90‑1 (MQ=255) ggATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCTCTTTTTTCAATTTTGggtg > 2:658935/1‑90 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCTCTTTTTTCAATTTTGggtgttttt > 2:899613/1‑87 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTGCCCCAAGGTTTCTCCTCTTTTTTTAATTTTGGTTGTTATc > 1:941334/1‑90 (MQ=255) aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTGCCACAAGGTTTCTCCTCTTTTTTCAATTTTGGGTGTTTTc > 1:951962/1‑90 (MQ=255) aCTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTTTCAATTTGGGTTGTTATCaaa > 2:147837/1‑90 (MQ=255) tCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGGTGTTATCAAATc > 1:110856/1‑90 (MQ=255) cACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTTCCACAAGGTTTTTCCTCTTTTTTCAATTTGGGGTGTTTTCAAATCg > 2:47415/1‑90 (MQ=255) cTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCCCTTTTTTCCATTTTGGGTGTTATCAAATCggt > 2:473832/1‑88 (MQ=255) tGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTTTCAATTTGGGGTGTTATCAAATCGTTa > 1:966761/1‑90 (MQ=255) tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgc < 2:888717/90‑1 (MQ=255) ccTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgcg > 1:784971/1‑90 (MQ=255) ttGCCCCCGCTGCGGCGGGGGGGTCAGTTGCTTTGCCCCAAGGGTTCTCCCCTGTTATCAACATGTGTAATTGTCAAAGAGTTAAGCGag > 2:13063/1‑89 (MQ=255) aCCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGtt < 1:889650/90‑1 (MQ=255) cGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgtg > 1:407766/1‑90 (MQ=255) gCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTgtgt < 2:142800/90‑1 (MQ=255) tGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTa < 2:941334/90‑1 (MQ=255) ggggTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTa < 1:147837/90‑1 (MQ=255) ggggGTGTCATTTGCTTTTCCACAAGGGTTTTCCTCTTTTTTCAATTTGGGTTGTTATCAAATCGTTACGCGGTGTTTGTGTTTTCtttt > 1:142800/1‑89 (MQ=255) | GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTA > NZ_CP009273/1329865‑1330033 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |