Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | NZ_CP009273 | 4,527,926 | 0 | A | C | 56.5% | 10.1 / 27.7 | 23 | G219G (GGT→GGG) | nanM | N‑acetylneuraminate epimerase |
Reads supporting (aligned to +/- strand): ref base A (10/0); new base C (0/13); total (10/13) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 8.74e-07 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-04 | |||||||||||
Rejected as consensus: Frequency below/above cutoff threshold. |
AGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATC > NZ_CP009273/4527861‑4528006 | aGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCa > 2:406437/1‑90 (MQ=255) taaaaTCCGGCACCCGTCCGAAATCCTGGTTTGCCTCCGCCA‑TAAAAAAG‑CCGGTTTTTACCCCT‑TTTATTC‑AAAA‑CCGCCGCACCAGcc < 1:320325/88‑1 (MQ=255) taaaaTCCGCCATCCGTCCCCAACCCTGTTTTGGCTTCGCCA‑TTAAAAAC‑CCAGTTTTTATCCCC‑TTTATTC‑ACAA‑CCGCCGCACCAGcc < 1:91/88‑1 (MQ=255) taaaaTACGGCACCCTTTCGCAATCCTGGTTTGGCTTCCCCA‑TTAAAAAG‑CCAGTTTTTATCCCC‑TTTATTC‑ACAA‑CCGCCGCACCAGcc < 1:215097/88‑1 (MQ=255) tCAAATCCGCCCCCCTTTCGCAATCCTGGTTTGGCTCCCCCA‑TTAAAAAG‑CCAGTTTTTATCCCC‑TTTATTC‑ACAA‑CCGCCGCACCAGcc < 1:97572/90‑1 (MQ=255) tCAAAATCGGCCTCCCTTCCCAACTCTTGTTTTGGTTCCCCATTTAATAAC‑CCGGTTTTTACCCCT‑TTTATCC‑ACAA‑CCGCCGCACCAGc < 1:284049/90‑1 (MQ=255) aaaTCCGGCCTCCTTCCCAAACCCTGTTTTGGCTTCGCCA‑TAAAAAAG‑CCAGTTTTTATCCCT‑TTTATTC‑ACAC‑CCGCCGCACCAGCCGt < 2:36860/90‑1 (MQ=255) aaTACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCA‑TTAATAAG‑CCAGGTTTTTTCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGtt > 2:184825/1‑90 (MQ=255) gTTCGCAATCCTGGTTTGGCTTCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGc > 1:447293/1‑90 (MQ=255) tCGCAATCCTGGTTTGGCTTCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGa > 1:382058/1‑90 (MQ=255) acgggTTTGGCTTCCCCC‑TTAAAAAA‑CCAAGTTTTAACCCCTTTTATCC‑ACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGcc < 2:82945/87‑1 (MQ=255) ggTTTGGCTTCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCg > 2:150948/1‑90 (MQ=255) ggCTTCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACt > 2:374699/1‑90 (MQ=255) tCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt > 1:542352/1‑90 (MQ=255) tcccA‑TTAAAAAACCCAGTTTTTACCCCT‑TTTTTCC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAt < 1:41879/88‑1 (MQ=255) cGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTcctgt > 1:82945/1‑87 (MQ=255) cGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCAtt > 2:517752/1‑90 (MQ=255) caataaaaaG‑CCAGGTTTTTACCCC‑TTTTTTCAACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGc < 2:410173/87‑1 (MQ=255) a‑TTAATAAG‑CCAGTTTTTATCCCT‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 1:12971/90‑1 (MQ=255) ataaaaaG‑CCAATTTTTTCCCCC‑TTTTTTC‑CCAACCCGCCGCCCCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCtg < 2:309547/89‑1 (MQ=255) aacccAGGTTTTATCCCC‑TTTATTC‑ACAC‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGa < 1:125803/87‑1 (MQ=255) aG‑CCGGTTTTTATCCCT‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTGCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAg < 1:24851/90‑1 (MQ=255) ggTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATc > 2:140123/1‑90 (MQ=255) | AGTTCAAATACGGCATCCGTTCGCAATCCTGGTTTGGCTTCGCCA‑TTAATAAG‑CCAGGTTTTATCACC‑TTTATTC‑ACAA‑CCGCCGCACCAGCCGTTCCGTACCAGGGCGATTCGCCAGCGTAACTCCATTGCTGTGTTGAGGGATC > NZ_CP009273/4527861‑4528006 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 5 ≤ ATCG/ATCG < 7 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |