Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 997,447 | 0 | T | G | 51.5% | 26.2 / 48.3 | 33 | V62G (GTA→GGA) | elfG | fimbrial protein |
Reads supporting (aligned to +/- strand): ref base T (3/13); new base G (17/0); total (20/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 9.77e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-05 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
GGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGT > NZ_CP009273/997359‑997529 | gggCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTa < 2:574926/90‑1 (MQ=255) ggCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAg > 1:906780/1‑90 (MQ=255) aTACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACa > 2:236620/1‑90 (MQ=255) aTACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACa > 2:502323/1‑90 (MQ=255) tactacTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGaaaa > 1:31887/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTTAAGAAAAAatt > 1:511519/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTTAAGAAAAAatt > 1:573189/1‑90 (MQ=255) actCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAttt > 2:1146842/1‑90 (MQ=255) tCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTa > 2:161701/1‑90 (MQ=255) tCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAAatttt > 1:849894/1‑89 (MQ=255) gggAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAg < 2:1021616/90‑1 (MQ=255) gAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTTTCAGTg > 1:384122/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGa > 2:779447/1‑90 (MQ=255) aTCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGACAAATTTTATCAGTGGa > 2:1101401/1‑90 (MQ=255) ttGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgagtcga < 1:360455/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:373218/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAAcgaatcgaat < 2:629425/90‑1 (MQ=255) ggCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAGAAAAATTTTTTCAGTGGAAagaatcgaat > 1:1082311/1‑90 (MQ=255) cAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTTAAAAAAAATTTTATCAAGGGAACCAATTGAAAGattt > 2:1071383/1‑89 (MQ=255) ccTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTg < 2:511519/90‑1 (MQ=255) cTGGAGGAAGAAGCAAATACGGCAGGGGGAGTTAAAAACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGc > 1:401523/1‑90 (MQ=255) gaggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGAt < 1:779447/90‑1 (MQ=255) aggaAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATa < 2:433687/90‑1 (MQ=255) gaagaaGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATaaa < 1:1034181/90‑1 (MQ=255) aagaagCAAATACGGCAGGGGGAGTTAAAGACAAAATTTATCAGTGGAACGAATCGAATGATTATTTTGTTTCCTGTGATTGCGATAAAg > 1:1107648/1‑90 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa < 1:745157/90‑1 (MQ=255) aagCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACa < 2:832250/90‑1 (MQ=255) gCAAATACGGCAGGGGGAGTTAAAGACAAAATTTATCAGGGGAACGAAACGAATGATTATTATGTTTTCTGTGATTTCGATAAAAACaaa > 1:533203/1‑89 (MQ=255) aaTACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTc < 1:606633/90‑1 (MQ=255) ggCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAg > 2:616373/1‑90 (MQ=255) gCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGt < 1:1071383/90‑1 (MQ=255) gCAGGGGGAGTTAAAGACAAATTTTTTCAGTGGAACGAATCGAATGATTATTATTTTTCCCGTGATTGGGATAAAAACAAAGTCAGAAGt > 1:1029313/1‑90 (MQ=255) gCAGGGGGAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAgg > 1:614092/1‑89 (MQ=255) | GGGCAGTGTCATACTACTCAGGGGAATCCGTATATTGGCGTCAATTTTGGCGTTAAAACCCTGGAGGAAGAAGCAAATACGGCAGGGGTAGTTAAAGACAAATTTTATCAGTGGAACGAATCGAATGATTATTATGTTTCCTGTGATTGCGATAAAGACAATGTCAGAAGT > NZ_CP009273/997359‑997529 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |