Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 3,454,136 | 0 | A | C | 52.4% | 22.6 / 25.2 | 21 | T321P (ACC→CCC) | gspF | type II secretion system inner membrane protein GspF |
Reads supporting (aligned to +/- strand): ref base A (8/2); new base C (0/11); total (8/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGAT > NZ_CP009273/3454048‑3454225 | ccTCAACAACCCCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAAc < 1:378351/90‑1 (MQ=255) caacaaCCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGc > 1:507666/1‑90 (MQ=255) caacaaCCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTAGCTGGAACAAACCGc > 1:212229/1‑90 (MQ=255) aaTCCGCCGGCGTCGGGCAAATGCGGAAGAGAACGTCCGCCGGGGTAACACCATTCATTTTTCGCTGGAACAACCCGC‑AATTTT‑Cccgcc < 1:338343/90‑1 (MQ=255) cccggCGTCTGGAAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATTTTTCGCTGAAAAAACCCGC‑AATTTC‑CCCGCCgatga < 1:589116/86‑1 (MQ=255) cccAGCGTCTGGAAAAGGCGGAAGAGAACGTTCGCCGGGTAAACAGCATTCATTTTTCGCTGGAAAAACCCGA‑AATTTT‑CCCGCCgatga < 1:346353/89‑1 (MQ=255) gggCAAAGGCGGCAGAAAACTTTCGCCGGGGAAACAGCATTCTTCTTTCGCTGAAACAACCCGC‑ATTTTC‑CCCGCCGATGATGCTCTACa < 2:133350/89‑1 (MQ=255) aaTGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAACCCGC‑AATTTT‑CCCGCCGATGATGCTCTACAtggtg < 1:178051/90‑1 (MQ=255) ggCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCctc > 2:79726/1‑90 (MQ=255) ttCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGaaaaaa > 1:64822/1‑90 (MQ=255) ccGCCAGGGAAACAGCATCCCTTTTCCGC‑GGAAAACCCCCCAAATTTC‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAg < 2:687534/89‑1 (MQ=255) aGGGAAACAGCTTTCCTTTTTCGCGGGAACAACCCGA‑ATTTTC‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGc < 2:515787/90‑1 (MQ=255) aaCAGCATTCTTTTTTCGCTGAAACAACCCGA‑AATTTC‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTc < 2:552397/90‑1 (MQ=255) aaCACCATTCTTTTTTCGGT‑GAAAAAACCCA‑AATTTTCCCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTc < 2:507666/90‑1 (MQ=255) aCAGCATTCATCTTTCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCg > 2:386081/1‑90 (MQ=255) aGCATTCTTTTTTGGCTGGAACAACCCGC‑ATTTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGc < 1:609799/90‑1 (MQ=255) aGCATTCTTTTTTCGCTGAAACAACCCGC‑ATTTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGc < 2:52133/90‑1 (MQ=255) cATCCATTTTCCGCTGAAACAACCCGC‑ATTTTC‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGcac < 1:1859/90‑1 (MQ=255) aTTCATCTTTCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGcaca > 2:471732/1‑90 (MQ=255) ttCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGt > 2:515251/1‑90 (MQ=255) ttCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGt > 2:521633/1‑90 (MQ=255) aaCCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGa > 1:309463/1‑90 (MQ=255) aCCGA‑ATTTTC‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGAt < 2:453073/90‑1 (MQ=255) | CCTCAACAACCTCGAAATTCGCCAGCGTCTGGCAAATGCGGCAGAGAACGTTCGCCAGGGTAACAGCATTCATCTTTCGCTGGAACAAACCGC‑AATTTT‑CCCGCCGATGATGCTCTACATGGTGGCCTCTGGCGAAAAAAGCGGGCAGCTCGGCACATTAATGGTCAGAGCCGCAGAT > NZ_CP009273/3454048‑3454225 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |