Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,906,007 | 0 | T | G | 56.7% | 29.9 / 37.5 | 30 | H276P (CAC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (1/12); new base G (17/0); total (18/12) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.50e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG > NZ_CP009273/1905919‑1906093 | gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 1:1145949/90‑1 (MQ=255) gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 2:1589894/90‑1 (MQ=255) gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 1:1350545/90‑1 (MQ=255) cAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtgggt < 1:140889/90‑1 (MQ=255) cgcgCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa < 1:517314/90‑1 (MQ=255) cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaac > 1:871603/1‑90 (MQ=255) tAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCAc < 1:871885/90‑1 (MQ=255) tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACCGa > 1:542034/1‑90 (MQ=255) tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACACACTGAGc > 2:1200866/1‑90 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGa > 1:378906/1‑90 (MQ=255) tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAAAAACTCACTGAGCGGTTTTGAc > 2:652534/1‑90 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTAcc < 2:1041858/90‑1 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGGTTTTGACTTACCGt > 2:133508/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGGGGGATTTCGACTTTCCCTTa > 1:1199844/1‑90 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGTATTAAAAACCCACTTAGCGATTTCGACTTTCCGtgta > 2:1416977/1‑87 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAACAAATCACTTAGCGGTTTTGACTTACCGGTaa > 2:616923/1‑90 (MQ=255) ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACCTACCCTTTATGCa > 1:638529/1‑90 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag < 1:1200866/90‑1 (MQ=255) tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGGTTTTGACTTACCGTTAATGCag > 2:188146/1‑90 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 2:680662/90‑1 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga < 2:1547857/90‑1 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg < 1:1505001/39‑1 (MQ=255) cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg > 2:1505001/1‑39 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa < 2:1054465/72‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTa > 1:1054465/1‑72 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGAGGAACAAATCACCTTGCGGTTTTGACTTACCGTTTATGCagagaggc > 2:871639/1‑88 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc > 1:528813/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc < 2:526162/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAAAAACCCACCGAACGATTTCCACTTACCGTTAATGCAGAGAGGcaacaatc > 1:1547857/1‑87 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGATTTTGACTTTCCGTTAATGGAGAGAGCcaacattc > 2:795260/1‑88 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACCCCCTCAGCGCTTTCCACTTTCCCCTAATTCAGAGAGCCCTCCTGc > 2:1281352/1‑90 (MQ=255) cTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGGTTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTgg > 2:1570607/1‑90 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:528813/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:538037/90‑1 (MQ=255) aTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCgg < 2:638529/90‑1 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 1:92914/90‑1 (MQ=255) gggtggggCTTGGACAAATCACAGAGCGATTTTGCCCTACCGTTAATTCAGGGAGCCATCACCCTGGTTGCTTCTTTTGGGTAATAGCTg > 2:1180413/1‑90 (MQ=255) | GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG > NZ_CP009273/1905919‑1906093 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |