Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,994,529 | 0 | A | C | 58.8% | 4.8 / 14.8 | 17 | T289P (ACC→CCC) | xdhA | xanthine dehydrogenase molybdenum‑binding subunit XdhA |
Reads supporting (aligned to +/- strand): ref base A (6/1); new base C (0/10); total (6/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 5.66e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTC > NZ_CP009273/2994440‑2994607 | agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTa > 1:101630/1‑90 (MQ=255) agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTa > 1:256646/1‑90 (MQ=255) agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTa > 2:108733/1‑90 (MQ=255) agCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTa > 1:536210/1‑90 (MQ=255) aTGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTACCCGc > 2:352233/1‑90 (MQ=255) gggcgccATCCCGGACCAGCAAGCTTGGCGGC‑CTCCGGGTAAATTTT‑CCCTTTAGCC‑GGAAAGGGTTTTCCCCCCCACCCCGTCCCCGcc < 2:101630/89‑1 (MQ=255) gcATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCt > 1:2299/1‑90 (MQ=255) gCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACg > 1:192923/1‑90 (MQ=255) tGGCGGC‑TTTCGGTTAAATTTTCCCC‑TTAGCG‑GTGAAGGGTTTTTCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAAt < 2:256646/90‑1 (MQ=255) gtaaaaTTT‑CCC‑CTAACCGGTGAAGATTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACcgc < 1:112328/88‑1 (MQ=255) tt‑CCCTTTAGCT‑GTAAAGTGTTTTCCCCCCCACCCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGaa < 1:39272/90‑1 (MQ=255) ttAGCC‑GTGAAGAGTTTTCCCCCCCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGa < 1:195293/90‑1 (MQ=255) ggCC‑GGGAAGAGTTTTCCCCCCCCCCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa < 2:536210/89‑1 (MQ=255) cgCC‑GGGATGGGTATTCCCCCGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa < 2:524669/89‑1 (MQ=255) aGCC‑GTGAAGAGTTTTCCCCCGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa < 1:232662/90‑1 (MQ=255) aGCC‑GTAAAGGGTTTTCCCCCCCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGaaa < 1:108733/90‑1 (MQ=255) gTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGtt > 2:157610/1‑90 (MQ=255) tGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTta > 1:288941/1‑90 (MQ=255) gtgtTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCt > 1:456283/1‑90 (MQ=255) ccGCACCCCGTCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc < 1:34501/89‑1 (MQ=255) ccGCACCCCGCCCCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTc < 2:456283/89‑1 (MQ=255) | AGCCAATGGCGGCATTCCTGACCAGCAAGCTTGGCGGCATTCCGGTGAAAGTTT‑CCC‑TTAGCC‑GTGAAGAGTGTTTCCTCGCAACCCGTACCCGCCACGCTTTTACCATTGACGGGCAAATGGGCGTGAACCGCGACGGAACATTGAAAGGTTATAGTCTGGATGTTC > NZ_CP009273/2994440‑2994607 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |