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.4% | 25.2 / 28.5 | 21 | intergenic (+24/‑140) | BW25113_RS25560/acnA | hypothetical protein/aconitate hydratase AcnA |
Reads supporting (aligned to +/- strand): ref base T (3/7); new base G (11/0); total (14/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.03e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCT > NZ_CP009273/1329865‑1330030 | ggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAt > 2:272349/1‑90 (MQ=255) ggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTcct > 1:360748/1‑88 (MQ=255) ctctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTTTTTGcc > 2:483794/1‑90 (MQ=255) ctctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGcc > 2:307942/1‑90 (MQ=255) gggCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTTATTTTCTTTTCCACCAGGTTTTTCCTCttt > 1:664409/1‑90 (MQ=255) gCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTa < 1:272349/90‑1 (MQ=255) gAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGGTGTTATTTGCTTTGCCACAAGGGTTCTCCTCTTTTATCa > 2:523001/1‑90 (MQ=255) cACTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCCCAAGGGTTTTCCTCTTTTTTTAATTTTGGGTGTTTTTAAAttg > 1:683800/1‑88 (MQ=255) cTGATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTTCTTTTCCACAAGGGTTTTCCTCTTTTTTCAATTTTGGGTGTTTTCAAATCGtt > 1:54893/1‑90 (MQ=255) gATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTTCTTTGCCACAAGGTTTCTCCTCTTTTTTCAATTTGGGGTGTTATCAAATCGTTAc > 2:576868/1‑90 (MQ=255) gATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAc < 2:54893/90‑1 (MQ=255) gATCCTTGACCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAc > 1:38869/1‑90 (MQ=255) tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgc < 2:683800/90‑1 (MQ=255) tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgc < 1:149924/90‑1 (MQ=255) tCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAcgc < 2:178775/90‑1 (MQ=255) cTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGa < 1:866040/90‑1 (MQ=255) gACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGt < 1:364609/90‑1 (MQ=255) aCCCCGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGtt > 2:87160/1‑90 (MQ=255) cccGCTGCGGCGGGGGTGTCATTTTCTTTTCCCCAAGGTTTTTCCTCTTTTTTCAATTTGGGGTGTTTTCAAATCGTTACGCGGTGtttt > 1:536939/1‑89 (MQ=255) cGCTGCGGCGGGGGTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgtg > 1:731408/1‑90 (MQ=255) ggcggGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCt > 2:26683/1‑90 (MQ=255) | GGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCT > NZ_CP009273/1329865‑1330030 |
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: ‑ |