| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | NZ_CP009273 | 4,475,282 | 0 | T | G | 55.9% | 27.5 / 40.2 | 34 | T163P (ACC→CCC) | pepA | leucyl aminopeptidase |
| Reads supporting (aligned to +/- strand): ref base T (0/15); new base G (19/0); total (19/15) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 5.39e-10 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-07 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
GCCGAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATCGAAACTGTAGAGCGTC > NZ_CP009273/4475196‑4475364 | gCCGAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACggcggggcgg > 2:1155004/1‑90 (MQ=255) cGAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCa < 2:251334/90‑1 (MQ=255) gagaTCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCGCGGCGGGGCGGCAc > 2:748182/1‑90 (MQ=255) ttttGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGaa > 2:921203/1‑90 (MQ=255) tttGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAac < 1:951003/90‑1 (MQ=255) tgctgcTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAacac < 1:804728/90‑1 (MQ=255) cggcggCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCa < 2:964949/66‑1 (MQ=255) cggcggCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCa > 1:964949/1‑66 (MQ=255) cggcggCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGGCGGCCCGGTGGACACCACCTTTCCCAGCg > 1:37204/1‑90 (MQ=255) gTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCtt < 1:126489/90‑1 (MQ=255) gTGCTGGATCGCGCCCTCACCGCTGGTCAGTTCACGGCGGGGCGGCCCGTTGAAAACCATCTTTCCCCGCGGACGACGCGGGTCGCTCtt > 1:792171/1‑90 (MQ=255) cTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGtt < 2:108503/90‑1 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGGGGGGCGGTGGAAACCCCCTTTCCCCGCGGGCGAGGCGGGTTGGTCTTTTTT‑Gttttt > 1:575113/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGGCCGGTGGACACCCCCTTTCCCCGCGGGCGACGCGGGTTGGTCTTTTTCTGtttt > 1:67315/1‑90 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGGACGTTTTACACCAACTTACGCCGCGGGCGGCGCGGGTCGCTCTTGTTT‑Gttttt > 2:1006948/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCGCGTTGGAAACCATCTTACGCAGCGGGCGACGCGGGTCGCTCTTTTTT‑Ggtttt > 2:910393/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCCCGGTGGACACCCCCTTTCCCCGCGGGGGGCGCGGGTCGCTCTTTTTC‑Gttttt > 1:829047/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGGAAACCATCTTTCGCCGCGGGCGACGCGGGTCGCTCTTGTTT‑Gttttt > 2:478840/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGGTCGCTCTTGTTC‑Gttttt > 2:908920/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGTTGAAAACCATCTTACGCAGCGGGCGACGCGGGTCGCTCTTTTTT‑Gttttt > 2:117538/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGGTTAACACCCCCCTTCCCCCCGGCCGCCGCGGTTCGCCCCTCTTC‑Ggtttt > 1:643494/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGGTTAACACCATCTTTCGCAGCGGGCGACGCGGGTCGCTCTTGTTT‑Gttttt > 2:186922/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGGTGAACACCATCTTACGCAGCGGACGACGCGGGTCGCTCTTGTTT‑Gttttt > 2:654467/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGCTGAACACCATCTTACGCCGCGGGCGACGCGGGTTGCTCCTTTTC‑Gttttt > 2:1083092/1‑89 (MQ=255) gcgcgcTCACCGCTGGTCAGTTCACGGCGGGGCGGCACGATGAACACCACCTTTCACAGCGGGCCACCCGGTTCGCTCTTTTTC‑Ctttct > 2:1066025/1‑88 (MQ=255) cgcTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGc < 1:1155004/90‑1 (MQ=255) aCCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATc < 1:788266/90‑1 (MQ=255) aCCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATc < 2:14454/90‑1 (MQ=255) aCCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATc < 1:1065834/90‑1 (MQ=255) gCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATCGaa < 1:488419/90‑1 (MQ=255) cTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATCGaaa < 1:738067/90‑1 (MQ=255) ggTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATCGAAACt < 1:469012/90‑1 (MQ=255) gTTCACGGCGGGGCGGCACGTTGAACACCATCTTACGCAGCGGGCGACGCGGGTCGCTCTTTTTC‑GTTTTTAGCTGATCGAAACTGTaga > 2:802842/1‑90 (MQ=255) cggcggGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATCGAAACTGTAGAGCGtc < 2:700309/90‑1 (MQ=255) cggcggGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATCGAAACTGTAGAGCGtc < 1:41891/90‑1 (MQ=255) | GCCGAGATCTTTTGCTGCTTTAATCCCGGCGGCAATCGCCAGACCGTGCTGGATCGCGCGCTCACCGCTGGTCAGTTCACGGCGGGTCGGCACGTTGAACACCATCTTACGCAGCGGACGACGCGGTTCGCTCTTGTTC‑GTTTTCAGCTGATCGAAACTGTAGAGCGTC > NZ_CP009273/4475196‑4475364 |
| Alignment Legend |
|---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |