Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 2,304,328 | 0 | A | C | 52.6% | 14.0 / 20.0 | 19 | V229G (GTA→GGA) | apbE | FAD:protein FMN transferase ApbE |
Reads supporting (aligned to +/- strand): ref base A (5/4); new base C (0/10); total (5/14) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 1.08e-02 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-03 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
TGCCAGAGGTGCTAATCCCATGACCGTTGATATCCACCACAGCCTGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCACCGACACCAGATAGCGGGAAATCCCTTCCTG > NZ_CP009273/2304245‑2304413 | tGCCAGAGGTGCTAATCCCATGCCCTTTAATACCCCCCACACCCTGACCCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTCCCCGcc < 2:203170/90‑1 (MQ=255) ggTGCTAATCCCATGCCCGTTGATACCCCCCACAGCCTGACCCGCGTTTTTTTTATCGGTTGTTTTTTGAATCGCTCCCCGCCACGGCAg < 2:506538/90‑1 (MQ=255) aTCCCATGACCGTTGATATCCACCACAGCCTGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTCCCCGCCACGGCAGGCCTTCa < 1:232202/90‑1 (MQ=255) ccccATGCCCGTTGATCCCCCCCACACCCTGACCCGCTTTTTTTTTATCGGTTGTTTTTTGAATCGCTCCCCGCCACGGCAGGCCTTCAc < 2:121930/89‑1 (MQ=255) ccATGACCGTTGATATCCACCACAGCCTGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCg > 1:108979/1‑90 (MQ=255) tgcccGTTGATCTCCCCCACAGCCGAAACCGGTTTTTTTTTACGGGTTGTTTTTTGAACCGCTCCCCGCCACGGCAGGCCTTCACCGTTc < 2:108979/87‑1 (MQ=255) tgcccGTTGATACCCCCCACAGCCTGACCCGCGTTTTTTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTc < 1:284604/87‑1 (MQ=255) gTTGATATCCACCACAGCCTGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTCATAcc > 1:531530/1‑90 (MQ=255) gTTGATATCCACCACAGCCTGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTCATAcc > 2:18212/1‑90 (MQ=255) tGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGc > 2:340974/1‑51 (MQ=255) tGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGc < 1:340974/51‑1 (MQ=255) gctttttttttATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCAcc < 1:382071/87‑1 (MQ=255) gcgTTTTCTTTATCGGTTGTTTTTTGAATCGCTCCCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCAcc < 1:78309/90‑1 (MQ=255) gcgTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGcgccgc > 2:13851/1‑85 (MQ=255) gcgTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGcgccgc < 1:13851/85‑1 (MQ=255) tttttATCGGTTGGTTTTTGAATCGCTCCCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCACCGACAcc < 1:139604/88‑1 (MQ=255) tggtttggttTTTGAAGCCCTCCCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCACCGACACCAGATAg < 2:173146/86‑1 (MQ=255) tgttttttAAACCCCTCCCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCACCGACACCAGATAGCGGGa < 2:186562/87‑1 (MQ=255) cGCTCCCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCACCGACACCAGATAGCGGGAAATCCCTTCCTg < 1:302698/90‑1 (MQ=255) | TGCCAGAGGTGCTAATCCCATGACCGTTGATATCCACCACAGCCTGAACCGCGTTTTCTTTATCGGTTGGTTTTTGAATCGCTACCCGCCACGGCAGGCCTTCACCGTTCATACCACGGCTGTTCAGCGCGCCGCCCACCGACACCAGATAGCGGGAAATCCCTTCCTG > NZ_CP009273/2304245‑2304413 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |