Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,023,324 | 0 | A | C | 53.8% | 9.0 / 16.8 | 13 | L323V (TTA→GTA) | ygfT | formate‑dependent uric acid utilization protein YgfT |
Reads supporting (aligned to +/- strand): ref base A (5/1); new base C (0/7); total (5/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.66e-03 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03 | |||||||||||
Rejected as consensus: E-value score below prediction cutoff. | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGC > NZ_CP009273/3023244‑3023409 | cAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAgg > 1:237029/1‑90 (MQ=255) aGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGc > 2:533396/1‑90 (MQ=255) gCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAGCGcc < 2:369951/90‑1 (MQ=255) aaTTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAGCGCCAAt > 2:234312/1‑90 (MQ=255) aCATCGCCCTGACCCCCTGCGCGCGCCAGAATATCAGCACCCCCTACCCCTGCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACgg < 2:85718/90‑1 (MQ=255) tCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGa > 1:791009/1‑90 (MQ=255) ccgcccTGACCTCCTGCGCGCGCAAGAATATCAGCCCCCCCTACCCCTGCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGa < 2:237029/86‑1 (MQ=255) gcccTGACCTCCTGCGCGCGCAAGAATATCAGCCCCCCCTACCCCTGCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGaac < 2:268491/88‑1 (MQ=255) cTCCTGCGCGCGCAAGAATACCACCACCCCCTACCCCTGCAGGCCCAGCGCCAACCACCGCCACTTTTTCGCTACGGGGAACAACTTTGc < 2:864444/90‑1 (MQ=255) ccggcgcgcgcAAGAATATCACCCCCCCCTACCCCTGCAGCCCCAGCCCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTg < 1:341158/87‑1 (MQ=255) ccTGCGCGCGCAAAAATACCACCACCCCCAACCCCGGCAGCCCCGGCGCAAACCCCCGCCACTTTTTCGCTACGGGGAACAACTTTGCTg < 2:222670/90‑1 (MQ=255) gcgcgcCAGAATATCAGCACACCCTAACCCTGCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATc > 1:66077/1‑90 (MQ=255) ccTACCCCTGCAGCCCCAGCCCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAgc < 2:733766/90‑1 (MQ=255) | CAGCATACCGCCAATTTCTGGATGGCGATCAAAGACATCGACCTGAACTCCTGCGCGCGCCAGAATATCAGCACACCCTAACCCTGCAGGCCCAGCGCCAATCACCGCCACTTTTTCGCTACGGGGAACAACTTTGCTGACATCGGGACGCCAGCCCATCGCCAGC > NZ_CP009273/3023244‑3023409 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |