Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,611,050 | 0 | A | C | 53.3% | 13.2 / 17.5 | 15 | T38T (ACA→ACC) | arsC | arsenate reductase (glutaredoxin) |
Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/8); total (6/9) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.40e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGTTCAAAGAGCCGGGAAACGCTGAACCTGCTGAAAGAAAAC‑GGCGTGGAGCCAGAAGTGGTGCTCTACCTTGAGACACCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAAGAGGATCTTTA > NZ_CP009273/2610973‑2611139 | tGTTCAAAGAGCCGGGAAACGCTGAACCTGCTGAAAGAAAAC‑GGCGTGGAGCCAGAAGTGGTGCTCTACCTTGAGACACCCGCCGATGCg > 1:67521/1‑90 (MQ=255) ccGGGAAACGCTGAACCTGCTGAAAGAAAAC‑GGCGTGGAGCCAGAAGTGGTGCTCTACCTTGAGACACCCGCCGATGCGGCAACGCTgcg > 2:361086/1‑90 (MQ=255) ccGCTGACCCTGC‑GAAAAAAAACGGGGGGGGAGCAAGAAGTGGTGCTCTCCTTTGAGACCCCCGCCGATGCGGCAACGCTGCGCGATTTg < 2:122992/89‑1 (MQ=255) acccTGACCCTGCAAAAAAAAACC‑GGGGGGGACCAAAAGGTGGTCCTCCACTTTGAGCCCCCCGCCGATGCGGCAACGCTGCGCGATTTg < 2:537393/87‑1 (MQ=255) gcaaaaaaaaaCC‑GGGGGGGAGCAAGAGGGGGTCCTCTCCCTTGAGACCCCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACt < 2:67521/86‑1 (MQ=255) gcaaaaaaaaaCC‑GGCGGGGACCCAGAAGTGGTGCTCTCCCTTGAGACACCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACt < 1:786906/86‑1 (MQ=255) taaaaaaaaCC‑GGCGGGGAGCCAGAAGTGGTGCTCTCCCTTGAGACCCCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTgg < 1:154425/88‑1 (MQ=255) aaagaaaAC‑GGCGTGGAGCCAGAAGTGGTGCTCTACCTTGAGACACCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTgggg > 1:486372/1‑90 (MQ=255) aaC‑GGCGTGGAGCCAGAAGTGGTGCTCTACCTTGAGACACCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAAt > 1:747676/1‑90 (MQ=255) cGTGGACCGAAAAGTGTTCCTCTCCCTGGAGACCCCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAgcgc < 1:65324/90‑1 (MQ=255) ggAGCCAGAAGTGGTGCTCTCCTTTGAGACCCCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCg < 1:223099/90‑1 (MQ=255) ccAGAAGTGGTGCTCTACCTTGAGACACCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGaa > 1:142434/1‑90 (MQ=255) tgCTCTCCTTTGAGCCCCCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGcc < 1:404115/90‑1 (MQ=255) tGAGACCCCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAAGAGGa < 1:361086/90‑1 (MQ=255) gagaCACCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAAGAGGAt > 1:172455/1‑90 (MQ=255) acacCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAaga < 1:23806/83‑1 (MQ=255) acacCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAaga > 2:23806/1‑83 (MQ=255) acCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAAGAGGATCTTta > 1:697229/1‑90 (MQ=255) acCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAAGAGGATCTTta > 1:79713/1‑90 (MQ=255) acCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAAGAGGATCTTta > 2:353634/1‑90 (MQ=255) | TGTTCAAAGAGCCGGGAAACGCTGAACCTGCTGAAAGAAAAC‑GGCGTGGAGCCAGAAGTGGTGCTCTACCTTGAGACACCCGCCGATGCGGCAACGCTGCGCGATTTGCTGAAAATACTGGGGATGAATAGCGCCCGTGAATTGATGCGCCAGAAAGAGGATCTTTA > NZ_CP009273/2610973‑2611139 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |