Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 1,350,458 | 0 | T | G | 52.4% | 15.3 / 26.1 | 21 | T304P (ACC→CCC) | sapA | peptide ABC transporter substrate‑binding protein SapA |
Reads supporting (aligned to +/- strand): ref base T (2/8); new base G (11/0); total (13/8) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 2.21e-04 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.23e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
ATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAG > NZ_CP009273/1350376‑1350547 | aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAAt < 1:725231/90‑1 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGGTTTGGG‑GGGGTTaaag > 2:39421/1‑89 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTTGG‑GGGGTTTAAtt > 1:1082631/1‑89 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGCGGGGGTTAAAt > 1:724396/1‑90 (MQ=255) tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGGTTGGG‑GGGGTTaaaag > 1:1349299/1‑88 (MQ=255) tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTc < 1:742756/90‑1 (MQ=255) tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTc < 2:675227/90‑1 (MQ=255) cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAAAGCCCGATAGGG‑GGCGGTCCTCCCCAGACg > 2:903547/1‑90 (MQ=255) gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGGGTTAAATGCCAGATAGGC‑GGCGTTCATCCCAGGACGc > 2:504168/1‑90 (MQ=255) aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTTGC‑GGGGTTAAAAGCCAGATAGGC‑GGCGTTTATCCCAGGGCGCAGc > 2:666823/1‑90 (MQ=255) aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGG‑GGGGTTTAAAGCCAGAGAGGCGGACCTCCTTCCC‑GGACGCAGc > 2:716107/1‑90 (MQ=255) aTGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGGGGTTAAAGCCAGAAAAGC‑GGCGGTCATCCCCAGACGCAGCGTTAAAAGCaa > 1:578996/1‑90 (MQ=255) tGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GGCGTTCATCCCAGGACGCCGCGTTAAAAGCAAg > 2:262385/1‑90 (MQ=255) ggCGGACAGCGGGATTATTTAGCGGC‑GGTTTTGG‑GGGGTTAAAAGCCAGATAAGG‑GGCGGTCATCCCAGGACGCAGCGGTTAACGCAAga > 2:923711/1‑89 (MQ=255) cAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGgt < 1:262385/90‑1 (MQ=255) gggATTATTTAGCGGC‑GGGTTGGC‑GGGGGTTAATGCCCGAAAGGG‑GGGGTTTATCCCCGGGAGCAGCGGTAAAAGCAAGCGCGGggcggc > 1:1342180/1‑88 (MQ=255) ggATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCa < 2:34925/90‑1 (MQ=255) ggATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCa < 1:923711/90‑1 (MQ=255) ttAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCaaaa < 2:1349299/90‑1 (MQ=255) gcggc‑gGTTTGGG‑GGGGGTAAATGCCAAAAAAGC‑GGCGTTTATCCCCGGACGCAGCGTTAAACGCAAAAGCGGGGCGGCACGCCAAAAgg > 1:540619/1‑90 (MQ=255) ggC‑GGTGTTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgctggct < 2:1342180/90‑1 (MQ=255) tgtTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag > 1:670951/1‑90 (MQ=255) | ATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGC‑GGTTTGGC‑GGTGTTAAATGCCAGATAGGC‑GACGTTCATCCCAGGACGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAG > NZ_CP009273/1350376‑1350547 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |