Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 385,080 | 0 | A | C | 69.4% | 59.1 / 30.0 | 36 | I35S (ATC→AGC) | hemB | porphobilinogen synthase |
Reads supporting (aligned to +/- strand): ref base A (5/6); new base C (10/15); total (15/21) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.00e+00 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
CGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAG > NZ_CP009273/384996‑385161 | cgTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGc < 1:268656/90‑1 (MQ=255) tGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAgg < 2:748670/90‑1 (MQ=255) gTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGt < 2:169023/90‑1 (MQ=255) tctGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTa < 2:283110/90‑1 (MQ=255) ctGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTaa < 1:204854/90‑1 (MQ=255) aTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTa < 1:56551/90‑1 (MQ=255) aTGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTa < 2:202227/90‑1 (MQ=255) gcgcATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAg > 2:555628/1‑90 (MQ=255) cATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAgtgt > 1:508221/1‑90 (MQ=255) aCGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtc < 2:415051/90‑1 (MQ=255) cGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtct < 1:461730/90‑1 (MQ=255) cGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtct < 1:761769/90‑1 (MQ=255) ccTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTctct > 2:309032/1‑90 (MQ=255) tGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTaa > 1:575220/1‑77 (MQ=255) tGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTaa < 2:575220/77‑1 (MQ=255) tGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTctct < 2:364926/88‑1 (MQ=255) tGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTctct > 1:364926/1‑88 (MQ=255) ggCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCa > 1:430558/1‑90 (MQ=255) gCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCt > 1:596176/1‑68 (MQ=255) gCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCt < 2:596176/68‑1 (MQ=255) tCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcg > 2:328960/1‑90 (MQ=255) cAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAgg > 1:129272/1‑53 (MQ=255) cAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAgg < 2:129272/53‑1 (MQ=255) cAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgc > 1:784474/1‑90 (MQ=255) aCGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgcgc < 1:621631/90‑1 (MQ=255) aCGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgcgc < 1:705060/90‑1 (MQ=255) gCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtct < 2:546350/70‑1 (MQ=255) gCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGtct > 1:546350/1‑70 (MQ=255) gCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCgcgcagc > 1:33351/1‑90 (MQ=255) gCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCgcgcagc > 2:290973/1‑90 (MQ=255) tttGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgc > 2:501609/1‑83 (MQ=255) tttGTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAgcgc < 1:501609/83‑1 (MQ=255) gTAGTCGTCAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAgg < 2:321042/90‑1 (MQ=255) cAATTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgc > 2:165574/1‑90 (MQ=255) aTTTCTTCTTCAACAAAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTgcgc < 1:478897/90‑1 (MQ=255) caacaaAGCTCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAg < 1:555628/90‑1 (MQ=255) | CGTGCCAGATGTTTCTCTGGAATGCGCATCACGCCTGGCATGGCTTCAACGGCTTTGTAGTCGTCAATTTCTTCTTCAACAAAGATCGGCAACACCAGGTCGTTAAGGCTAAGTGTTGTCTCTTCAAACATAGCGCGCAGCGCAGGAGATTTGCGCAGGCGACGAG > NZ_CP009273/384996‑385161 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |