| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 1,429,964 | 0 | C | A | 61.1% | 12.4 / 26.1 | 18 | intergenic (‑88/+53) | uspF/ompN | universal stress protein UspF/porin OmpN |
| Reads supporting (aligned to +/- strand): ref base C (6/1); new base A (0/11); total (6/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 3.77e-04 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-04 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
CCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAA > NZ_CP009273/1429881‑1430040 | cctcctGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAagga < 2:717393/90‑1 (MQ=255) tATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCa > 2:431746/1‑90 (MQ=255) aaGTGTACCGCAAAATGACTTTTTTTGTGAAAGGACGGAGTTCACTTTTTTAATTTAGATAAAAGGAGGAAAAAAAAGAAAAAAAAGccc < 1:295343/90‑1 (MQ=255) aaGTGAAACGCCAGATGACCTTTTTTGTGAAATGAGGGAGTTCACATTTTTAATTTAGATAAAAGGAGAAAAAAAAAGCAAAAAAAGccc < 2:41681/90‑1 (MQ=255) ccAGATGACTTTTTTTGTGAAAGGAGGGATTTCACTTTTTTAATTTAAAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAg < 1:233199/90‑1 (MQ=255) ttttttGTGTAAGGAGGGATTCCATTTTTTTAATTAAAAAAAAAGGAGAAAAAAAAAGAAAAAAACCCCCGCCATAACAGCGGGCaggag < 1:100480/89‑1 (MQ=255) ttttttGTATAATGACGAATTCCATTTTTTAAATTAAAAAAAAAGGAGAAAAAAAAAAAAAAAAAACCCCGCCATAACAGCGGGCaggag < 1:776935/89‑1 (MQ=255) ttttttGTAAAAGGAGGAATTTCAATTTTTAAATTAAAAAAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCaggag < 1:358907/89‑1 (MQ=255) ttttttGGGTAATGGCGGTGTCCATATTTTTAATAAAAAAAAAAGGAGAAAAAAAAAAAAAAAAAAGCCCGCCATAACAGCGGGCaggag < 1:594285/89‑1 (MQ=255) ttGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAtt > 2:610437/1‑90 (MQ=255) tAATGCGGGAGTTCATTTTTTTATTTAAGATAAAAGGAGAAAAAAAAAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:850770/90‑1 (MQ=255) tAATGACGGATTTCACTTTTTTAATTTAGATAAAAGGAGGAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:332424/90‑1 (MQ=255) tAATGACGGAGTTCACATTTTTAATTTAAAAAAAAGGAGAAAAAAAAAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:27063/90‑1 (MQ=255) tAAGGACGGAGTTCCCTTTTTAAATTTAAAAAAAAGGAGGAAAAAAAAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGaa < 2:51612/90‑1 (MQ=255) ggAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaa > 1:332424/1‑90 (MQ=255) ggAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATaa > 1:300015/1‑90 (MQ=255) tCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAg > 2:772396/1‑90 (MQ=255) ttttAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa > 2:9948/1‑90 (MQ=255) | CCTCCTGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGATCAAAGGAGGAAGAATAAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAA > NZ_CP009273/1429881‑1430040 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |