Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,061,180 | 0 | A | C | 67.5% | 49.0 / 41.7 | 40 | V486G (GTT→GGT) | yihQ | sulfoquinovosidase |
Reads supporting (aligned to +/- strand): ref base A (6/7); new base C (14/13); total (20/20) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAA > NZ_CP009273/4061044‑4061315 | ggTGAAGGCGCTGAAACCGCCCCAGCGCAGCAGCAGCTTTTTTCTGCGTTTCATCCCAAGCGGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAAcc < 2:579762/138‑1 (MQ=255) cgggCGCAGAAACCGCACCAGCGCGGCAGCACCTTTTTGCGGCGCTTCATCCCAAACGGGGGGGTTAAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACg < 2:762319/137‑1 (MQ=255) cTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGc > 2:426299/1‑139 (MQ=255) gCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCAtcgctg > 1:762319/1‑136 (MQ=255) cgcAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTc > 2:790714/1‑138 (MQ=255) cagcagcTTTTTGCTGCGCTTAATCCAAAAGGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCg < 2:12009/139‑1 (MQ=255) ctgcGTTTCATCACAAACGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgc < 2:719800/139‑1 (MQ=255) cgcTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccg > 2:646416/1‑139 (MQ=255) cTTCATCACAAGCGGGGGGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcccgcc < 2:194598/139‑1 (MQ=255) aaCAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCAtgg < 1:790714/139‑1 (MQ=255) ggtggtGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcc > 1:210571/1‑120 (MQ=255) ggtggtGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcc < 2:210571/120‑1 (MQ=255) ggtggtGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAg > 1:277300/1‑139 (MQ=255) ggtggtGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCAt > 2:672958/1‑89 (MQ=255) ggtggtGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCAt < 1:672958/89‑1 (MQ=255) gtttAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAgg < 2:818401/78‑1 (MQ=255) gtgtAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAgg > 1:818401/1‑81 (MQ=255) ccGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGc > 2:657380/1‑139 (MQ=255) gCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTa < 2:291583/139‑1 (MQ=255) cgaTGTCGCTATGGTGCAGGCCATGTCCGGTCATTGCGAGCGACAGCGCCGCCGGGACACCCGCCGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTAcc < 2:56172/137‑1 (MQ=255) cAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTAcc > 2:124983/1‑139 (MQ=255) aaTGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcg > 2:285464/1‑103 (MQ=255) aaTGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcg < 1:285464/103‑1 (MQ=255) cAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcc > 2:271329/1‑90 (MQ=255) cAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTcgcc < 1:271329/90‑1 (MQ=255) gCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATa < 1:825793/139‑1 (MQ=255) gTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGa > 1:357702/1‑138 (MQ=255) gTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGa > 2:27051/1‑138 (MQ=255) gTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCAtggtgg > 2:262406/1‑94 (MQ=255) gTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCAtggtgg < 1:262406/94‑1 (MQ=255) gTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAggc > 2:849260/1‑138 (MQ=255) gTCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGa > 2:422750/1‑139 (MQ=255) tCATTGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGAt < 1:657380/139‑1 (MQ=255) tGCCAGCGACAGCGCCGCCGGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCg > 1:16295/1‑139 (MQ=255) cGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGc > 2:772942/1‑139 (MQ=255) cgccgcCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCcaca > 2:229932/1‑68 (MQ=255) cgccgcCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCcaca < 1:229932/68‑1 (MQ=255) ccgccgGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGt < 1:820043/139‑1 (MQ=255) ccgccgGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGt < 1:658353/139‑1 (MQ=255) cgccgGGACACCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGtt < 2:846564/139‑1 (MQ=255) acacCCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCaa < 2:853333/139‑1 (MQ=255) | CGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGGTCATTGCCAGCGACAGCGCCGCCGGGACAACCGACGCCAGGCCATCGTCGAGACTCCAGTCGACGTTCTGGTCGCCCGCCCACATCATGGTGGAGTATTTCTGGCTACCGGTAGAACCGGCGCGCATAAAGAAAAGGATCTCGCCGAGCTTGCCCGTTTCTTCAA > NZ_CP009273/4061044‑4061315 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |