Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,963 | 0 | T | G | 55.6% | 9.4 / 24.3 | 18 | D492A (GAC→GCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (1/7); new base G (10/0); total (11/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.51e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-04 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAA > NZ_CP009273/822881‑823051 | aaaGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGcc < 1:189235/90‑1 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAg > 2:223747/1‑90 (MQ=255) gggtggTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGaaa > 2:711097/1‑89 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCGGGAGTCGGGtc > 1:745367/1‑90 (MQ=255) gACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGGCGGTtc > 1:238464/1‑90 (MQ=255) gTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCGGAAGACGGGTCGTCGGGaaaa > 2:491453/1‑89 (MQ=255) cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa < 1:813369/90‑1 (MQ=255) tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTGGGGTCGTCGAGAAACAGAATGtg > 2:314264/1‑89 (MQ=255) tACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAAAGTCCGGGTc > 2:29977/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAAAAGAAAGTCCGGGTCATGc > 2:339921/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 1:339921/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGGGAAAAAAAAAGGCCGGGTCCTGCATCaa > 1:827887/1‑89 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 2:383657/90‑1 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAAc < 2:577651/90‑1 (MQ=255) aGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCCGGAGGCGGGTTGTCGGGAAACAGAAAGTCCGGGTTATGCATCAGCGGACAGgcc > 1:383657/1‑90 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 2:827887/90‑1 (MQ=255) gAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTaaa < 2:289906/90‑1 (MQ=255) gTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTaa > 1:358287/1‑90 (MQ=255) | AAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGGCCAGCGCCAGCCGCTGTTTAAAACCTAA > NZ_CP009273/822881‑823051 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |