Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,710,374 | 0 | T | G | 55.0% | 11.8 / 26.4 | 20 | G156G (GGT→GGG) | yiaD | OmpA family lipoprotein |
Reads supporting (aligned to +/- strand): ref base T (1/8); new base G (11/0); total (12/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.14e-05 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCAT > NZ_CP009273/3710285‑3710463 | ggcgCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGggt < 2:202087/90‑1 (MQ=255) gACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGGTTTATCCACAGCACGGGGGGGCACGACCCGaa > 2:523397/1‑90 (MQ=255) gCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGGGGTCACGACCTGAACATGc > 1:58084/1‑90 (MQ=255) gTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGt < 1:205634/90‑1 (MQ=255) gAAAGAGTATCCGAAAACGGCGGTTTACGGGATTGGGTTTACCCACAGCACGGGGGGGCACGACCCGAACAAGCGTCTCTCCCCGCaaag > 2:289768/1‑88 (MQ=255) gAAAGAGTATCCGAAAACGGCGGTTAAAGGGATTTGGTTTACCCACAACACGGGGGGGCACGACCCGAACATGCGGCTCTCCCACCAACg > 2:399769/1‑90 (MQ=255) agagTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGGGGGCACGACCTGAAAAAGCGTCTCTCCCCGCAACGGGc > 2:75876/1‑90 (MQ=255) gagTATCCGAAAACGGCGGTTAACGGGATTGGTTTTAACGACAGGACGGGGGGGCGCGACCTGAACAAGCGGCTCTCCCCGCAACGTggg > 1:254923/1‑88 (MQ=255) agTATCCGAAAACGGCGGTTAACGTGATTTGTTATACCCACAGCACGGGGGGGCACGACCCTAAAATGCGTCTCTCCCCGCAACGGGCgg > 1:503590/1‑90 (MQ=255) agTATCCGAAAACGGCGGGTTAAGTGATTGGGTTTAACGGCAAAACGGGGGGGGAGGACCTGAAAAAGTGGGTCTCCCAGCAAAGCGCgg > 2:411031/1‑90 (MQ=255) aTCCGAAAACGGCGGTTAACGTGATTGGTTATACCCACAGCACGGGGGGGCACGACCTGAAAATGCGTCTCTCCCCGCAACGGGCGGGtt > 2:80910/1‑90 (MQ=255) aTCCGAAAACGGCGGGTAAAGTGATTGGTTATACCCACAGCCCGGGGGGGCACGACCCGAACAAGGGGGTCTCCCCGCAACGGGCGGagt > 2:350848/1‑88 (MQ=255) gcggGTAAAGTGATTGGGTTTACCGACAGCACGGGGGGTCACGACCTGAACATGCGGCTCTCCCAGCAAAGTGGGGGTTCCCTTGCCAgc > 1:231798/1‑90 (MQ=255) tAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGtt < 1:393502/90‑1 (MQ=255) cGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGAt < 1:79482/90‑1 (MQ=255) ttGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCAccc < 1:80910/90‑1 (MQ=255) taCCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGt < 2:10002/90‑1 (MQ=255) cAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCg > 2:35690/1‑49 (MQ=255) cAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCg < 1:35690/49‑1 (MQ=255) cGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGcc < 2:231798/90‑1 (MQ=255) cGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGcc < 1:248179/90‑1 (MQ=255) ggtggtCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGc < 2:226919/90‑1 (MQ=255) tggtCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCAt < 1:251736/90‑1 (MQ=255) | GGCGCTAACACCCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCAT > NZ_CP009273/3710285‑3710463 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |