Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,819,318 | 0 | G | A | 42.3% | 26.2 / 12.7 | 26 | intergenic (‑124/+79) | ves/spy | environmental stress‑induced protein Ves/ATP‑independent periplasmic protein‑refolding chaperone Spy |
Reads supporting (aligned to +/- strand): ref base G (15/0); new base A (0/11); total (15/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.29e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGC > NZ_CP009273/1819230‑1819405 | tCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGagaaaaga > 2:1871245/1‑90 (MQ=255) cGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAAccgccg > 1:2147757/1‑90 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc > 1:28919/1‑90 (MQ=255) gATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATcc > 2:2013687/1‑90 (MQ=255) gtttATCGCCAACTGGGGTGATTTACTTTAGTGGTGATCAGTTTTTTTTTACAATTTTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGt < 1:898915/88‑1 (MQ=255) aTCGGCAACTGGGGTGTTTTCCGTTAGTGGTGATCAGGTTTTTTTTACAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCAc < 2:1060218/90‑1 (MQ=255) aTCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCAc > 1:564531/1‑90 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 1:1591725/1‑90 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 2:753600/1‑90 (MQ=255) aCTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa > 2:862157/1‑90 (MQ=255) ttattttAGGTTAGTGGTGATCGGGTTTTTTTTCCAATTTTTTGAAAAAAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa < 1:931028/86‑1 (MQ=255) gtgtTTTCCTTTGGTGGTGATCGGTTTTTTTTTACAATTGTTTGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa < 2:2428110/90‑1 (MQ=255) gtgtTTTACTTTGGTGGTGATCAGGTTTTTTTTACAATTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 2:2133659/90‑1 (MQ=255) gtgtTTTACGTTAGTGGTGATCGGGTTTTTTTTACAACTTTTGGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCaa < 1:862157/90‑1 (MQ=255) gTGATTTACGTTAGTGGTGATCAGTTTTTTTTTACAACTGTTTGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCaa < 2:2409759/90‑1 (MQ=255) tattttACGTTAGTGGTGATCGGTTTTTTTTTACATCTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAg < 1:2209934/87‑1 (MQ=255) ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg > 2:931028/1‑90 (MQ=255) ttACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAg > 2:1194441/1‑90 (MQ=255) tACTTTAGGGGTGTTCGGGTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAACCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGt < 1:1965099/90‑1 (MQ=255) tAGTGGTGATCAGTTTTTTTTTACAATTTTTGGAAAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACa < 1:1201984/90‑1 (MQ=255) tttttACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAg > 2:337699/1‑90 (MQ=255) aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCAtt > 1:2187460/1‑50 (MQ=255) aCAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGttttt > 2:2422798/1‑90 (MQ=255) aCAACTGTTTGAGAAGAAAAAAAAAAACCGCCGATCCTGTCCACCGCAtt < 2:2187460/50‑1 (MQ=255) gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTTAGTTTTTTGGCtt > 1:2043864/1‑89 (MQ=255) gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTTGCTg > 1:1938585/1‑90 (MQ=255) gTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGCCTTAAGGTTTTTTGCTg > 2:1209155/1‑90 (MQ=255) agaaaagaaaACCGACGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGGCTTAAGTTTTTTGTCTGAAAGattatt > 1:1142298/1‑90 (MQ=255) aagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc < 1:753600/90‑1 (MQ=255) aagaaaaCCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTcagc > 1:108622/1‑90 (MQ=255) | TCAGGTAACCGTGGATGGTTATCGGCAACTCGGGTGATTTACGTTAGTGGTGATCAGGCTTTTTTTACAACTGTTTGAGAAGAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGTCTTAAGTTTTTTGGCTGAAAGATTATTCAGC > NZ_CP009273/1819230‑1819405 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |