Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,350,458 | 0 | T | G | 59.3% | 20.8 / 31.8 | 27 | T304P (ACC→CCC) | sapA | peptide ABC transporter substrate‑binding protein SapA |
Reads supporting (aligned to +/- strand): ref base T (0/11); new base G (16/0); total (16/11) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 7.67e-08 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-06 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTG > NZ_CP009273/1350369‑1350543 | ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGgt < 2:232126/90‑1 (MQ=255) aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTTTGCGgggt > 2:983931/1‑90 (MQ=255) taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTa < 1:1378818/90‑1 (MQ=255) taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTa < 1:786817/90‑1 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGCGGGGGTaaaag > 2:492128/1‑88 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTTGGGGGGGTTtaaag > 2:955141/1‑88 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTTGGGGGGGTTaaaag > 1:886211/1‑88 (MQ=255) ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTTGGGGGGGTTaaagc > 1:390825/1‑88 (MQ=255) ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTAAATGc < 1:1133930/90‑1 (MQ=255) gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGG‑TTTGGGGGGGTTAAATGc > 1:75220/1‑90 (MQ=255) gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGGGTTTAAAGcc > 2:1133930/1‑90 (MQ=255) aTTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGGGGGGTTTAATGCCa > 2:870949/1‑90 (MQ=255) aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTTGCGGGGTTTAAAGCCAGATAg > 2:944121/1‑90 (MQ=255) aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGG‑TTTGGGGGGGTTAAA‑GCCAGGTAgggg > 1:1467795/1‑88 (MQ=255) ggCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTAAATGCCAGATAGGCGa < 2:1232397/90‑1 (MQ=255) tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTc < 2:639527/90‑1 (MQ=255) gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGGGGGGTTTAATGCCCGATAGGGGGCGGTCATTCCAGGACGc > 2:992537/1‑90 (MQ=255) gcgcATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑GTTGGGGGGGGTTAATGGCCGAGAGGCGGCGCTCATCCCCGGGCGCCG‑GGGTaaaa > 1:200659/1‑89 (MQ=255) aTGGCGGACAGCGGGATTATTTAGCGGC‑GG‑ATTTGCGGGGGTTAAAGCCCGAAAGGCGGAGCTCATCCCAGGGCGCAGCCGTTaaagcca > 2:32275/1‑86 (MQ=255) tGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGGGTTAAAAGCCAGATAAGCGACGGTCATCCCCGGACGCAG‑CGTTAAAAGCAAg > 2:617932/1‑90 (MQ=255) cAGCGGGATTATTTAGCGGC‑GG‑TTTTGCGGGGTTTAAAGCCAGAAAAGCGGCGTTCATCCCCGGGCGCAG‑CGGTAAAAGCAAACGCgggg > 2:832912/1‑89 (MQ=255) aGCGGGATTATTTAGCGGC‑GGTTTTGGGGGGGTTAATTCCAAAAAAGGGGCGGTTCTTCCCCGGA‑GCAG‑CGGTAAAAAAAAACGCggggg > 1:396826/1‑88 (MQ=255) gCGGGATTATTTAGCGGC‑GG‑TTTGGCGGGGTTAAAAGCCAGAAAAGGGAGGGTTATCCCCGGACGCCG‑CGGTAAAAGCAAGCGCgggggg > 1:145356/1‑87 (MQ=255) cGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGgtcgt < 2:1221341/90‑1 (MQ=255) attTAGCGGC‑GG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGGTCGTCACGCaa < 1:955141/90‑1 (MQ=255) attTAGCGGC‑GG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGGTCGTCACGCaa < 1:450706/90‑1 (MQ=255) g‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgct < 2:562623/90‑1 (MQ=255) g‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgct < 1:324799/90‑1 (MQ=255) gCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggctg < 1:32275/90‑1 (MQ=255) | ATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GG‑TTTGGCGGTGTTAAATGCCAGATAGGCGACGTTCATCCCAGGACGCAG‑CGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTG > NZ_CP009273/1350369‑1350543 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |