Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,434,438 | 0 | T | C | 72.7% | 24.3 / 15.5 | 22 | E216G (GAA→GGA) | fabB | beta‑ketoacyl‑ACP synthase I |
Reads supporting (aligned to +/- strand): ref base T (2/4); new base C (7/9); total (9/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTA > NZ_CP009273/2434307‑2434572 | aGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGgt > 1:557259/1‑139 (MQ=255) gCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGgtg < 1:326761/139‑1 (MQ=255) aTTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGt < 1:353646/139‑1 (MQ=255) tCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTAt < 2:557259/138‑1 (MQ=255) tgtgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAg < 2:452355/139‑1 (MQ=255) tgAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAgcg < 2:353009/139‑1 (MQ=255) cacGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCAt > 1:492145/1‑139 (MQ=255) ttCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTcc > 1:345294/1‑139 (MQ=255) cgccgccgcCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTcgcc > 2:112951/1‑139 (MQ=255) ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa < 2:67284/115‑1 (MQ=255) ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa < 2:506135/115‑1 (MQ=255) ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa > 1:67284/1‑115 (MQ=255) ccgccgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCa > 1:506135/1‑115 (MQ=255) cgccAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCgccgccgcc < 1:15737/139‑1 (MQ=255) aaaCCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGt < 1:112951/139‑1 (MQ=255) aaaCCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGt < 1:332427/139‑1 (MQ=255) aaaCCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGt > 1:224453/1‑139 (MQ=255) cGGTGAGCGTCGTAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGc < 2:224453/139‑1 (MQ=255) cgtcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTg > 2:113869/1‑139 (MQ=255) gtcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTgg < 2:345294/139‑1 (MQ=255) tcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGtt > 2:27228/1‑35 (MQ=255) tcgtAAGTACGGGAGGCTTTTCCCGGGGTGTCGtt < 1:27228/35‑1 (MQ=255) ttttcCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTa > 1:479531/1‑139 (MQ=255) | AGCCAACGATTTCAGCATAGATGTGAGCACCACGCGCCAGCGCGTGTTCCAGCTCTTCAACCACTACCATACCGCCGCCGCCAGCGATAACGAAACCGTCACGGTGAGCGTCGTAAGTACGGGAGGCTTTTTCCGGGGTGTCGTTGTATTTAGTAGACAGCGCACCCATTGCGTCGAATTCGCAAGCCATTTCCCAGCACAGCTCTTCGCCGCCGCCAGCAAACACGATGTCCTGTTTGCCCAGTTGGATCTGCTCTACTGCGTTA > NZ_CP009273/2434307‑2434572 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |