Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,906,011 | 0 | T | G | 54.5% | 39.7 / 50.5 | 33 | T275P (ACC→CCC) | htpX | protease HtpX |
Reads supporting (aligned to +/- strand): ref base T (1/14); new base G (18/0); total (19/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.83e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-08 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG > NZ_CP009273/1905924‑1906093 | ggACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa < 1:443401/90‑1 (MQ=255) gcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGaaca > 1:521042/1‑90 (MQ=255) gcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGaaca > 1:45055/1‑90 (MQ=255) tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAAATCACTGAGCGa > 1:183889/1‑90 (MQ=255) ttCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGAttt < 1:637214/90‑1 (MQ=255) cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGAttttg > 2:35720/1‑88 (MQ=255) cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATTACTTAGCGATTTTGAC‑TTTcc > 1:716997/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGt < 1:261524/90‑1 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCGTGAAAAAACCACAGAGCGGTTTTGAC‑TTTCCGt > 2:741135/1‑90 (MQ=255) tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGGTTTCGAC‑TTACCGt > 1:706268/1‑90 (MQ=255) aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTa > 1:129631/1‑90 (MQ=255) ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTaa < 2:729381/90‑1 (MQ=255) aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGc < 2:128004/90‑1 (MQ=255) cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCagaga < 2:496168/90‑1 (MQ=255) gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGcc > 1:251875/1‑90 (MQ=255) aCGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACCGAGCGGTTTTGAC‑TTTACGTTTATGCAGAGAGGcat > 2:298708/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGc < 2:98634/90‑1 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGc > 1:209366/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGc > 1:433429/1‑90 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAAAAAAACACTGAGCGATTTTGACTTTACCG‑TAATGCAGAGAAGCAACATgg > 1:446617/1‑89 (MQ=255) gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACCCACTGAGCGATTTCGAC‑TTACCCTTAATGGAGAGAGGCAACATGc > 2:36103/1‑90 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGGTg > 1:101367/1‑90 (MQ=255) gCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGCC‑TTACCGTTAATGCAGAGAGCCATCATGCTg > 2:673008/1‑90 (MQ=255) cTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTgg < 1:505824/90‑1 (MQ=255) cTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTgg < 1:86504/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:251875/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 1:619760/90‑1 (MQ=255) ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc < 2:736536/90‑1 (MQ=255) ccAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGGTGCTTCTTTCgggt > 1:663780/1‑88 (MQ=255) gcggcgGGGGGGGCATGAACAACTCACTGAGCGATTTTGAC‑TTACCGTTAATGCAGAGAGCCATCATTCTGGTTGCTTTTTGCGGTTCAt > 2:217453/1‑90 (MQ=255) gcggcgGGGGGGGCATGAAAAACTCACTGAGCGATTTTGAC‑TTTCCGTTTATGCAGAGAGCCATCATGCTGGTTGGTTTTTTCGGTTTAt > 2:217452/1‑90 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 2:407643/90‑1 (MQ=255) cggGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt < 1:723763/90‑1 (MQ=255) ggggggggCATGAACAACTCACTGAGCGGTTTTGAC‑TTTCCCTTAATGCAGAGAGCCATCATTCTGGTTGCTTCTTGCGGTTCATAGCTg > 2:637317/1‑90 (MQ=255) | GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGAC‑TTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG > NZ_CP009273/1905924‑1906093 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |