Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 822,955 | 0 | T | G | 61.1% | 13.2 / 19.4 | 18 | T495P (ACC→CCC) | ybhF | ATP‑binding cassette domain‑containing protein |
Reads supporting (aligned to +/- strand): ref base T (0/7); new base G (11/0); total (11/7) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 3.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGG > NZ_CP009273/822870‑823023 | gCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAgg < 2:524352/90‑1 (MQ=255) tCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCa < 1:772972/90‑1 (MQ=255) ggtggtGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGAAgg > 1:851783/1‑89 (MQ=255) tgACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACGCCGGGAGTCgggt > 2:649666/1‑88 (MQ=255) aCCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCCGGAGGCGGTtcg > 1:590380/1‑90 (MQ=255) cATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGAAGTCGGTtcgtc > 2:364294/1‑90 (MQ=255) cGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGaaa < 2:50578/90‑1 (MQ=255) cGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGaa < 1:389456/90‑1 (MQ=255) tttCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAACGCCGGGAGTCGGTTCGTCGGGAAAAAAAAAGTc > 1:45412/1‑90 (MQ=255) gtctcTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTcc < 2:699753/89‑1 (MQ=255) ctctACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCgg < 2:471667/90‑1 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGTTCGTCGGGAAACAGAAAGTCCGGGTCATGc > 2:242609/1‑90 (MQ=255) aTGCTGTTGATGTGCAGCCAAAATTCACGGCGGGGGAGGGGGGCAAAGCCGGAAGTCGGTTCGTCGAGAAAAAGAATGTCCGGGTCATGc > 2:235176/1‑90 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAg < 2:183729/90‑1 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGGAACCCCGGAAGGTGGGTTGTCGGGAAAAAGAAAGTCCGGGTTATGCCTCaa > 2:824447/1‑89 (MQ=255) gTTGATGTGCAGCCAAAATTCACGGCGGGGGGGGGGGGCAAAGCCGGGAGGCGGGTCGTCGAGAAAAAAAAAGTCCGGGTCCTGCATCAg > 1:183729/1‑90 (MQ=255) tgtgCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGGCGGGTCGTCGAGAAAAAAAAAGTGCGGGTCATGGATCAGCGAAc > 2:853240/1‑90 (MQ=255) gCAGCCAAAATTCACGGCGGGGGAGGGGGGCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCAAGCATCAGCGAACAgg > 2:264829/1‑90 (MQ=255) | GCTTCATCCATAAAGTGGGTGGTGACCATCACCGTGACGCCTTTCTCTACCATGCTGTTGATGTGCAGCCAAAATTCACGGCGGGTGAGGGGGTCAACGCCGGAAGTCGGTTCGTCGAGAAACAGAATGTCCGGTTCATGCATCAGCGAACAGG > NZ_CP009273/822870‑823023 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |